李秀英，鄧海巖

(1.大同市地震局，山西 大同 037006;2.同車公司觀測站，山西 大同 037038)

0 引言

地裂縫是在一種由內(nèi)、外引力作用下，巖石和土層發(fā)生變形，當(dāng)作用力大于土層內(nèi)力作用累積時(shí)，巖土層破裂，其連續(xù)性遭到破壞，形成裂縫。我國城市地裂縫及其災(zāi)害研究大致經(jīng)歷以下幾個(gè)階段：早期分散調(diào)查階段(主要有地裂縫場地工程穩(wěn)定性的分散調(diào)查對地裂縫的零散目測調(diào)查)；以落實(shí)震兆為主的調(diào)查研究階段；抗裂與震兆并行的勘察研究高潮階段；多種目的的研究與綜合研究階段。我國專家學(xué)者前后陸續(xù)對大同市地裂縫成因、地下水開采對地裂縫活動(dòng)的影響等展開了一系列的相關(guān)研究，并取得一些成果[1-7]，但大多數(shù)研究都是從定性方面評估，定量和半定量的研究相對較少。

同車公司(原大同機(jī)車廠)地裂縫是大同市唯一監(jiān)測的一條地裂縫，也是大同市最早發(fā)現(xiàn)的地裂縫。文章以同車公司地裂縫觀測數(shù)據(jù)為依據(jù)，收集地裂縫及井水位觀測近20年的資料，運(yùn)用數(shù)學(xué)相關(guān)系數(shù)法探討地下水位變化與地裂縫活動(dòng)之間的相關(guān)性。

1 觀測站基本概況

1.1 構(gòu)造背景

大同盆地屬新生代斷陷盆地，大同市地處大同盆地的北端，NE向的口泉斷裂、NEE走向的六棱山斷裂、恒山北緣斷裂是盆地內(nèi)的3條主要活動(dòng)斷裂，控制著盆地的形成、發(fā)育以及新構(gòu)造格局和地震活動(dòng)。

大同市地裂縫排列有序、展布成帶，具有較強(qiáng)的方向性，與大同市區(qū)斷裂構(gòu)造的主體方向一致。房子村斷裂是夾在口泉斷裂和桑干河斷裂之間的次級斷裂。房子村斷層活動(dòng)是同車公司地裂縫的深部構(gòu)造背景，同車公司地裂縫是房子村斷層在地表的直接反映。據(jù)地質(zhì)勘探表明，房子村斷層與地裂縫的產(chǎn)狀一致，走向NE，傾向SE，傾角約75°。

同車公司地裂縫西起十里河?xùn)|岸，終止南關(guān)南街，全長約5 km，由多條NE走向、右行排列的個(gè)體地裂縫組合而成，優(yōu)勢方位NE 57°，傾向SE，傾角60°～70°，表現(xiàn)為高角度活動(dòng)正斷層,以南東盤下降為主，水平拉張為輔，伴有左旋扭動(dòng),為一張性地裂縫。

1.2 臺(tái)站及儀器安裝概況

地裂縫觀測站位于大同盆地的西北部，地處十里河沖、洪積平原上,臺(tái)站基底為第四系砂、黏土交互層。2012年同車公司地裂縫垂直向、水平向數(shù)字化改造前，使用的是原國家地震研究所研制的DG-1A型靜力水準(zhǔn)儀及SR-1型水平伸縮儀，量程為20 mm，通過調(diào)整可增大到80 mm,儀器格值為0.01 mm,讀數(shù)精度±0.01 mm，兩測頭高差標(biāo)準(zhǔn)偏差小于±0.07 mm。2012年數(shù)字化改造后，使用的是中國地震局地殼應(yīng)力研究所下屬公司(北京震苑迪安防災(zāi)技術(shù)研究中心)研發(fā)的“三維斷層形變觀測系統(tǒng)”(斷層垂直形變分量儀DFG、斷層水平法向分量儀DSG、斷層水平切向分量儀DLG)進(jìn)行觀測,由二維觀測升級為三維觀測。主觀測室為N-S方向，與同車公司地裂縫正交。測線長度為13 m，儀器的分辨力為0.004 mm,量程為40 mm,線性度為0.5%，重復(fù)性誤差為0.5%(見第26頁圖1)。

圖1 觀測站儀器設(shè)備平面布置示意圖Fig.1 Layout of instruments in observation station

2 地裂縫活動(dòng)與地下水位變化相關(guān)分析

據(jù)文獻(xiàn)[8]研究地裂縫活動(dòng)表現(xiàn)出強(qiáng)烈的三維性，以垂直沉降為主，橫向水平拉張次之，左旋水平扭動(dòng)最弱，其比值為垂直沉降量：水平張量：水平扭量為10∶6∶1。因此文章就地裂縫垂向(DFG)觀測值與地下水位變化的相關(guān)性進(jìn)行著重分析。

據(jù)相關(guān)報(bào)告李常保，張曉燕，蘇 丹,等.2012—2015年度大同市及周邊平原區(qū)地下水動(dòng)態(tài)年報(bào).山西省地質(zhì)勘查局二一七地質(zhì)隊(duì)，2012—2015.的研究，大同市城區(qū)及周邊平原區(qū)地下水開采的主要途徑為居民飲水、農(nóng)作物灌溉及占比例較小的工業(yè)用水。

大同市有3個(gè)地下水集中開采區(qū)，形成3個(gè)明顯的地下水降落漏斗，分別為城北、城南、城西水源地降落漏斗區(qū)。觀測站位于大同市城西水源地降落漏斗之中(見圖2)，水源地開采量為7.09萬 m3/d。主要開采第四系淺中層潛水或承壓水以及深部承壓水，含水層巖性以粗砂、中砂、細(xì)砂為主，礫石、粉砂次之。地下水補(bǔ)給來源主要是大氣降雨、十里河水的入滲補(bǔ)給等。

2.1 資料選取

選取地裂縫垂直向(DFG)1998-2016年的資料,D5井也選取此時(shí)段的井水位測值(見第27頁圖3)。

2.2 方法及數(shù)據(jù)分析

相關(guān)關(guān)系是變量之間的非確定關(guān)系，變量之間關(guān)系密切程度通常會(huì)在相關(guān)系數(shù)中顯現(xiàn)出來。文章試圖用地下水位的變化與地裂縫活動(dòng)兩者變量之間的相關(guān)系數(shù)，探討地下水動(dòng)態(tài)變化與地裂縫活動(dòng)間的密切程度。據(jù)文獻(xiàn)[9]，用相關(guān)系數(shù)的變化范圍可判斷兩者之間的相關(guān)程度。R=0時(shí)，完全不相關(guān)；0<∣R∣≤0.3時(shí)，基本不相關(guān)；0.3<∣R∣≤0.5時(shí)，弱相關(guān)；0.5<∣R∣≤0.8時(shí)，較強(qiáng)相關(guān)；0.8<∣R∣<1時(shí)，強(qiáng)相關(guān)；∣R∣=1時(shí)，完全相關(guān)。

圖3 D5井水位與地裂縫垂直向月均值變化對比圖Fig.3 Comparison of water level in well D5 and ground fissures vertical monthly mean value

D5井是山西省地質(zhì)環(huán)境中心的一口長期觀測井，設(shè)立完全符合地下水觀測規(guī)范的要求。D5井位于十里河沖洪積平原，距同車公司地裂縫西端不足1 km，距地裂縫觀測站約3.6 km，與地裂縫觀測站同處城西水源地降落漏斗之中，井水位變化基本反映城西降落漏斗區(qū)的地下水動(dòng)態(tài)。據(jù)相關(guān)報(bào)告①的結(jié)果，地裂縫活動(dòng)并未對地下水動(dòng)態(tài)變化造成影響。

(1) 由D5井水位與地裂縫垂直向(DFG)年動(dòng)態(tài)資料(見圖3、表1)看出，地下水與地裂縫年變幅最大變化時(shí)間并不一致。地下水位在2000年及2016年變化較為顯著。2000年的年變幅最大，達(dá)2.96 m，地裂縫年變幅為3.76 mm；2016年地下水位年變幅為2.37 mm,地裂縫年變幅為0.77 mm；1999年地裂縫年變幅最大，達(dá)7.1 mm,地下水位年變幅為1.46 m。表明地下水位的變化與地裂縫的變化時(shí)間并不同步，兩者無必然的相關(guān)性。

(2) 從地裂縫活動(dòng)與D5井水位月均值相關(guān)程度來看(見表1)，不相關(guān)與弱相關(guān)的年份占53%，分別是1998、2000、2001、2004、2006、2013、2014、2016年，相關(guān)系數(shù)分別是-0.04、-0.27、-0.12、-0.28、0.23、-0.10、0.01、0.10；較強(qiáng)(中等程度)相關(guān)的年份占總年數(shù)的35%，分別是2002、2003、2005、2007、2009、2012年，相關(guān)系數(shù)分別是0.66、0.53、0.69、0.6、-0.7、0.56；強(qiáng)相關(guān)的年份僅占總年數(shù)的12%，分別是2008、2015年，相關(guān)系數(shù)均為0.91，呈正相關(guān)。

表1 地裂縫垂直向與D5井水位兩者年變幅及月均值相關(guān)系數(shù)統(tǒng)計(jì)表Table 1 Correlation coefficient of annual variation amplitude and monthly mean value between vertical variation of ground fissure and water level in well D5

地裂縫垂直向與地下水位動(dòng)態(tài)1998-2016年月均值相關(guān)系數(shù)變化圖(見圖4)看出，兩者的月均值相關(guān)系數(shù)變化極其不穩(wěn)定，波動(dòng)較大，相關(guān)性方向具有隨機(jī)性，只有個(gè)別年份存在強(qiáng)相關(guān)。由此可見，地下水位變化與地裂縫活動(dòng)兩者之間不具有必然的相關(guān)性。

圖4 地裂縫垂直向與D5井水位月均值相關(guān)系數(shù)變化圖Fig.4 Correlation coefficient variation of monthly mean value between vertical variation of ground fissure and water level in well D5

(3) 根據(jù)1998-2009年D5井水位年變幅與地裂縫垂直向年變幅資料統(tǒng)計(jì)計(jì)算，其相關(guān)系數(shù)為0.17(見表1)。說明地下水位變化與地裂縫活動(dòng)基本上不具有相關(guān)性。

(4) 由地裂縫垂直向DFG測值和大同盆地ML1.0以上地震活動(dòng)情況(見第28頁圖5、表2)可看出，地裂縫活動(dòng)經(jīng)歷活躍段(1998-1999年)—平靜段(2000-2005年)—活躍段(2006-2011年)—平靜段(2012-2016年)，表明地裂縫具有間歇活動(dòng)的特征。

圖5 地裂縫垂直向測值與大同盆地ML 1.0以上地震活動(dòng)情況對比圖Fig.5 Comparison of ground fissure vertical variation and activity of earthquakes with ML≥1.0 in Datong Basin

震級地震活動(dòng)次數(shù)統(tǒng)計(jì)1998-1999年2000-2005年2006-2011年2012-2016年1.0～1.98655717266312.0～2.9555120152663.0～3.971722164.0～4.993115.0～5.91000總次數(shù)1 501701901714

由圖5及表2看出：

(1) 地裂縫在第1個(gè)活躍段(1999年)之前，S盤相對于N盤年沉降量在7～9 mm,在此期間,大同盆地小震活動(dòng)密集、頻次高、強(qiáng)度大，呈震級上升趨勢。發(fā)生71次3.0～3.9級地震，9次4.0～4.9級地震，1次5.0～5.9級地震，最大為1999年11月1日大同陽高M(jìn)5.6地震。

(2) 地裂縫在第1個(gè)平靜段(2000-2005年)，S盤相對于N盤年沉降量在2～3 mm,在此期間,大同盆地小震活動(dòng)密集，頻次、強(qiáng)度明顯下降，出現(xiàn)震級下降趨勢。發(fā)生7次3.0～3.9級地震，3次4.0～4.9級地震，最大為2000年11月5日大同縣M4.5地震。

(3) 地裂縫在第2個(gè)活躍段(2006-2011年)，S盤相對于N盤年沉降量最大為5 mm,明顯小于第1個(gè)活躍段。在此期間,大同盆地小震活動(dòng)密集、頻次較高、強(qiáng)度明顯增大，呈震級上升趨勢。發(fā)生22次3.0～3.9級地震，1次4.0～4.9級地震，最大為2010年4月4日大同陽高M(jìn)4.7地震。

(4) 地裂縫在第2個(gè)平靜段(2012-2016年)，S盤相對于N盤年沉降量不足1 mm,在此期間,大同盆地小震活動(dòng)密集，強(qiáng)度明顯下降，出現(xiàn)震級下降趨勢。發(fā)生16次3.0～3.9級地震，1次4.0～4.9級地震，最大為2013年2月22日大同縣、陽高縣與河北陽原交界M4.1地震。

由此可看出,大同盆地小震活動(dòng)變化經(jīng)歷密集(1998-1999年)—平靜(2000-2005年)—密集(2006-2011年)—平靜(2012-2016年)四個(gè)周期。與同車公司地裂縫活動(dòng)在時(shí)間大尺度上具有同步性，表明構(gòu)造活動(dòng)可能是地裂縫活動(dòng)的力源。

3 結(jié)論

綜合上述分析，得出如下結(jié)論：

(1) 從地裂縫垂直向(DFG)和大同市地下水位變化資料分析，兩者變化不同步，無必然的聯(lián)系。

(2) 地下水位與地裂縫活動(dòng)兩者的相關(guān)系數(shù)變化極其不穩(wěn)定，表明兩者不具有必然的相關(guān)性。地裂縫活動(dòng)與地震活動(dòng)在大尺度時(shí)間上的同步性，說明地裂縫的活躍程度取決于地區(qū)的構(gòu)造活動(dòng)，非地下水位變化所致。

(3) 根據(jù)1998-2009年的地下水與地裂縫垂直向年變幅資料，其相關(guān)系數(shù)為0.17，說明地下水位變化與地裂縫活動(dòng)基本上不具有相關(guān)性。

(4) 根據(jù)王景明等研究地下水活動(dòng)，尤其是抽汲地下水水位下降，乃至造成地面沉降，因其發(fā)生在地表之下一定深度，對近地表的區(qū)域微破裂開啟、擴(kuò)寬并成為地裂縫沒有明顯的影響，但對地下斷層蠕滑及其兩側(cè)土體差異沉降有一定的控制作用[10]。

(5) 根據(jù)大同市內(nèi)地裂縫的空間分布，地裂縫在地下水集中開采區(qū)有分布，在非集中開采區(qū)也有發(fā)育。

(6) 通過地震物探勘測和鉆探驗(yàn)證發(fā)現(xiàn)，同車公司地裂縫在深部發(fā)育有房子村活斷層，形成地裂縫的應(yīng)力來源是由深部的房子村活斷層的活動(dòng)向上傳遞，不是由地下水位變化引起。

(7) 通過同車公司地裂縫活動(dòng)與大同盆地地震活動(dòng)情況對比發(fā)現(xiàn),地裂縫活動(dòng)與大同盆地地震活動(dòng)在時(shí)間大尺度上具有同步性。

由上分析認(rèn)為，同車公司地裂縫活動(dòng)的主要原因不是地下水位變化，深部基底斷層的活動(dòng)可能是地裂縫活動(dòng)的力源。