程鵬圖,楊 順,王晨曦,馬云舉
(寧夏回族自治區(qū)地震局,寧夏 銀川 750001)
前兆觀測(cè)是地震監(jiān)測(cè)預(yù)報(bào)的主要手段之一,具有明確的物理含義以及預(yù)測(cè)機(jī)理,因此對(duì)前兆觀測(cè)出現(xiàn)的各種異常變化進(jìn)行實(shí)地落實(shí),是對(duì)這些異常變化進(jìn)行分析判斷的主要依據(jù)。但在實(shí)際觀測(cè)中,異常落實(shí)具有較強(qiáng)的復(fù)雜性,多數(shù)異常變化不能做到準(zhǔn)確判斷[1]。本文通過(guò)對(duì)海原小山形變觀測(cè)點(diǎn)鉆孔傾斜異常變化進(jìn)行跟蹤落實(shí),排查了可能影響數(shù)據(jù)異常的觀測(cè)系統(tǒng)、觀測(cè)環(huán)境等影響因素,結(jié)合前期積累的異常應(yīng)震實(shí)例,客觀分析了此次異常變化。
海原地震臺(tái)小山形變觀測(cè)點(diǎn)位于海原縣海城鎮(zhèn)的小山自然村南西側(cè),位于青藏高原東北弧形構(gòu)造帶、海原斷裂帶上,地面距離1920年海原81/2級(jí)地震主斷裂450 m左右(圖1)。海原斷裂帶展布于米家山北麓,哈思山南麓,黃家洼山南麓,西、南華山北麓和月亮山東麓,南華山及其以西,走向北西西,以東為北西走向,由11條不連續(xù)次級(jí)剪切斷層組成,各次級(jí)斷層均向所在山體傾斜,地貌上是一條醒目的山脈與盆地的界線。
圖1 海原小山形變觀測(cè)點(diǎn)地理位置圖Fig.1 The geographic location diagram of Haiyuan Xiaoshan deformation observatory
海原斷裂帶早期為逆沖活動(dòng),后期為左旋走滑活動(dòng),其活動(dòng)方式的轉(zhuǎn)化發(fā)生于早更新世末至中更新世初。左旋走滑活動(dòng)在大多數(shù)次級(jí)剪切斷層的接續(xù)部位形成了拉分型盆地,由斷裂兩側(cè)地址體所記錄的水平位移量達(dá)幾到幾十公里,穿過(guò)斷裂的各級(jí)別水系,均發(fā)生幅度不等的左旋位錯(cuò)(圖2)。地貌、水系研究結(jié)果顯示,全新世以來(lái)斷裂帶的平均走滑速率為6~10mm/a。觀測(cè)點(diǎn)附近巖性大多數(shù)為灰黑色變質(zhì)灰?guī)r和泥石混合物,由于1920年海原大地震,地質(zhì)構(gòu)造比較復(fù)雜,巖層比較破碎。
圖2 海原斷裂帶西華山一帶地質(zhì)體左旋位移平面圖Fig.2 The left鄄hand displacement plan of the West Huashan area at the Haiyuan fault zone
2017年7月27日,臺(tái)站值班人員發(fā)現(xiàn)海原鉆孔傾斜資料NS向加速下降,下降幅度102毫角秒,其中7月27日下降幅度65毫角秒,7月28日—8月2日處在低值,且出現(xiàn)新的穩(wěn)定點(diǎn);EW向加速上升,上升幅度110毫角秒,其中7月27日上升幅度最大,達(dá)74毫角秒,7月28日—8月2日處在上升后穩(wěn)定期,且出現(xiàn)新的穩(wěn)定點(diǎn)(圖3)。8月2日,臺(tái)站組織技術(shù)人員赴海原小山形變觀測(cè)點(diǎn)進(jìn)行鉆孔傾斜資料異常落實(shí)。
圖3 海原小山鉆孔傾斜分鐘值曲線(20170720—20170803)Fig.3 The minute value curve of the Haiyuan Xiaoshan borehole tilt anomaly
2017年8月2日,海原地震臺(tái)組成異常核實(shí)小組,赴現(xiàn)場(chǎng)進(jìn)行核實(shí)調(diào)查,并初步制定了異常核實(shí)工作方案:(1)對(duì)局部觀測(cè)環(huán)境進(jìn)行調(diào)查;(2)對(duì)儀器工作電源進(jìn)行檢查;(3) 儀器標(biāo)定記錄、日志記錄等信息的核實(shí);(4)對(duì)異常信度進(jìn)行分析。
異常核實(shí)小組對(duì)海原小山鉆孔傾斜觀測(cè)系統(tǒng)進(jìn)行了檢查,儀器觀測(cè)系統(tǒng)基本正常,主要落實(shí)情況如下:
(1) 供電電源檢查
異常核實(shí)小組對(duì)海原小山鉆孔傾斜供電系統(tǒng)進(jìn)行了檢查。結(jié)果顯示交流供電224.2 V,穩(wěn)壓電源供電41.17 V,電源供電正常。對(duì)觀測(cè)電源各接頭進(jìn)行按壓,保證接頭接觸良好。
(2)數(shù)采記錄系統(tǒng)檢查
異常核實(shí)小組人員對(duì)海原小山鉆孔傾斜儀器數(shù)采記錄系統(tǒng)進(jìn)行檢查,發(fā)現(xiàn)NS向有所偏大,粗調(diào)后恢復(fù)。
(3) 標(biāo)定格值變化
海原小山鉆孔傾斜格值變化在2013年7月維修之前變化穩(wěn)定,2013年7月維修之后NS向存在規(guī)律的起伏變化,但總體穩(wěn)定,EW向變化穩(wěn)定。2016年6月、12月,2017年6月三次標(biāo)定結(jié)果變化穩(wěn)定。
(4) 日志檢查
查閱了儀器觀測(cè)日志,鉆孔傾斜日志記錄完整,記錄系統(tǒng)觀測(cè)正常。
3.2.1 儀器所在基巖的巖性情況
海原小山鉆孔傾斜井深60.24 m,儀器探頭安裝在約51 m深處。探頭所處巖性為灰黑色灰?guī)r,巖性完整。
3.2.2 測(cè)點(diǎn)周邊的環(huán)境變化
海原小山鉆孔傾斜觀測(cè)點(diǎn)位于寧夏海原南華山自然保護(hù)區(qū)內(nèi),高程2220 m,附近人煙稀少,對(duì)測(cè)點(diǎn)周邊5 km范圍進(jìn)行普查,結(jié)果顯示測(cè)點(diǎn)周邊無(wú)大型的施工作業(yè)。
經(jīng)初步核實(shí),認(rèn)為2017年7月27日起發(fā)生的海原小山鉆孔傾斜NS向加速下降異常、EW向加速上升異常與觀測(cè)系統(tǒng)、降雨、氣溫、氣壓等相關(guān)性不大;海原小山附近為山區(qū),環(huán)境干擾較少,但周邊存在抽水井。異常核實(shí)小組赴海原縣自來(lái)水公司調(diào)查了附近機(jī)井抽水情況。海原小山周邊存在抽水機(jī)井,主要供縣城用水,其中南部區(qū)域機(jī)井距離小山觀測(cè)點(diǎn)平均4.7 km,中部區(qū)域機(jī)井距離小山觀測(cè)點(diǎn)平均5.5 km,北部區(qū)域機(jī)井距離小山觀測(cè)點(diǎn)6.4 km。小山觀測(cè)點(diǎn)鉆孔布置于山腰處,為基巖場(chǎng)區(qū),因此異常組認(rèn)為抽水對(duì)小山形變的干擾較小,排除了機(jī)井抽水對(duì)形變異常的影響。為了進(jìn)一步驗(yàn)證這一推論是否成立,本文對(duì)抽水是否影響形變異常進(jìn)行了以下深入分析:
(1) 從實(shí)際觀測(cè)可得,2017年7月27日起海原小山鉆孔傾斜EW向加速上升,上升幅度110 ms,其中7月27日上升幅度最大,達(dá)74 ms;7月28日—8月2日處在上升后穩(wěn)定期。且出現(xiàn)新的穩(wěn)定點(diǎn);海原小山鉆孔傾斜NS向加速下降,下降幅度102 ms,其中7月27日下降幅度65 ms;7月28日—8月2日處在低值,且出現(xiàn)新的穩(wěn)定點(diǎn)。根據(jù)實(shí)地調(diào)查,海原小山周?chē)灿?8口機(jī)井,距海原小山平均距離為5.5 km,抽水量為800 m3/d·口,18口機(jī)井每天抽水總量為18口×800 m3/d·口=1.44×103m3,根據(jù)二維載荷均勻分布模型[2],在5km范圍內(nèi),每天1.44×103m3抽水量產(chǎn)生的影響僅為0.008 ms,遠(yuǎn)小于此次異常的變化幅度,所以可以推定海原小山周?chē)臋C(jī)井抽水量變化不會(huì)對(duì)海原小山鉆孔傾斜產(chǎn)生影響。
(2) 海原小山周邊 18口井抽水量為800 m3/d·口,根據(jù)最小“安靜”距離理論的三維集中載荷模型,可以推算抽水量為800m3/d抽水井的最小“安靜”距離小于50 m,而抽水井距海原小山測(cè)點(diǎn)的距離最小為4.7 km,遠(yuǎn)大于抽水井的最小“安靜”距離,所以可以進(jìn)一步推斷海原小山周?chē)臋C(jī)井抽水量變化不會(huì)對(duì)海原小山鉆孔傾斜產(chǎn)生影響[3]。
如圖4所示,海原小山鉆孔傾斜異常時(shí)段期間,鉆孔氣壓及水位并無(wú)明顯的異常變化與之相對(duì)應(yīng),而在異常時(shí)段內(nèi),27日至28日出現(xiàn)了連續(xù)降雨,但降雨量相對(duì)較少,因此同樣可排除降雨對(duì)鉆孔傾斜異常的影響。由于海原小山測(cè)點(diǎn)距抽水井距離較遠(yuǎn)、且抽水量相對(duì)較小,因此抽水對(duì)海原小山影響很小,可忽略不計(jì),可以推斷海原小山鉆孔傾斜不存在明顯的環(huán)境干擾。對(duì)海原小山鉆孔傾斜進(jìn)行潮汐分析和內(nèi)精度分析,自觀測(cè)以來(lái),鉆孔傾斜潮汐因子和內(nèi)精度變化穩(wěn)定,觀測(cè)數(shù)據(jù)真實(shí)可靠,曾連續(xù)在國(guó)家局資料評(píng)比中獲獎(jiǎng)。
圖4 海原小山鉆孔傾斜與鉆孔氣壓、鉆孔水位、降雨的對(duì)比曲線Fig.4 A contrast curve between borehole tilt and drilling pressure,drilling water level,and rainfall of Haiyuan Xiaoshan
2017年7月27日起海原小山鉆孔傾斜出現(xiàn)異常,異常形態(tài)為NS向加速下降后恢復(fù)、EW向加速上升后恢復(fù),期間伴有階躍式變化。類(lèi)似于此的異常變化,自觀測(cè)以來(lái)有兩次加速下降、階躍式異常對(duì)應(yīng)地震,具體分析如下:
(1) 2016年3月31日,海原小山鉆孔傾斜NS向出現(xiàn)加速上升,上升幅度約31 ms,EW出現(xiàn)加速下降,下降幅度約43 ms(圖5),2016年4月11日發(fā)生寧夏中寧MS3.9級(jí)地震,震中距離82 km。
圖5 海原小山鉆孔傾斜分鐘值曲線(中寧MS3.9級(jí)地震前)Fig.5 The minute value curve of the Haiyuan Xiaoshan borehole tilt anomaly(Before Zhongning MS3.9)
(2) 2016年11月10日,海原鉆孔傾斜NS向出現(xiàn)加速上升,上升幅度約28 ms,EW向出現(xiàn)加速下降,下降幅度約36 ms(圖6),2016年11月13日發(fā)生寧夏海原MS2.3級(jí)地震,震中距離28 km。
圖6 海原小山鉆孔傾斜分鐘值曲線(海原MS2.3級(jí)地震前)Fig.6 The minute value curve of the Haiyuan Xiaoshan borehole tilt anomaly(Before Haiyuan MS2.3)
由于排除了觀測(cè)系統(tǒng)、觀測(cè)環(huán)境等干擾因素的影響,加之分析了兩次典型的異常應(yīng)震案例,分析認(rèn)為此次異常可能為震前異常。
綜合分析認(rèn)為,海原小山鉆孔傾斜2017年7月27日出現(xiàn)的加速變化真實(shí)可靠,無(wú)觀測(cè)環(huán)境、氣象因素干擾,可作為短期異常進(jìn)行跟蹤,需密切關(guān)注海原周邊震情發(fā)展,臺(tái)站工作人員需做好應(yīng)急準(zhǔn)備。
2016年以來(lái),海原小山鉆孔傾斜出現(xiàn)3月31日、11月10日兩次明顯的加速下降變化,兩次變化之后分別發(fā)生了寧夏中寧Ms3.9級(jí)、寧夏海原MS2.3級(jí)地震,此次異常變化較前兩次略有不同,出現(xiàn)NS向加速下降—恢復(fù)、EW向加速上升—恢復(fù)的過(guò)程,期間伴有階躍式變化,總變化幅度較大,達(dá)215 ms,短臨形勢(shì)較前兩次嚴(yán)峻,需結(jié)合地球物理場(chǎng)、地球化學(xué)場(chǎng)等其他資料密切跟蹤、分析其變化機(jī)理。