程鵬圖,楊 順,王晨曦,馬云舉
(寧夏回族自治區地震局,寧夏 銀川 750001)
前兆觀測是地震監測預報的主要手段之一,具有明確的物理含義以及預測機理,因此對前兆觀測出現的各種異常變化進行實地落實,是對這些異常變化進行分析判斷的主要依據。但在實際觀測中,異常落實具有較強的復雜性,多數異常變化不能做到準確判斷[1]。本文通過對海原小山形變觀測點鉆孔傾斜異常變化進行跟蹤落實,排查了可能影響數據異常的觀測系統、觀測環境等影響因素,結合前期積累的異常應震實例,客觀分析了此次異常變化。
海原地震臺小山形變觀測點位于海原縣海城鎮的小山自然村南西側,位于青藏高原東北弧形構造帶、海原斷裂帶上,地面距離1920年海原81/2級地震主斷裂450 m左右(圖1)。海原斷裂帶展布于米家山北麓,哈思山南麓,黃家洼山南麓,西、南華山北麓和月亮山東麓,南華山及其以西,走向北西西,以東為北西走向,由11條不連續次級剪切斷層組成,各次級斷層均向所在山體傾斜,地貌上是一條醒目的山脈與盆地的界線。

圖1 海原小山形變觀測點地理位置圖Fig.1 The geographic location diagram of Haiyuan Xiaoshan deformation observatory
海原斷裂帶早期為逆沖活動,后期為左旋走滑活動,其活動方式的轉化發生于早更新世末至中更新世初。左旋走滑活動在大多數次級剪切斷層的接續部位形成了拉分型盆地,由斷裂兩側地址體所記錄的水平位移量達幾到幾十公里,穿過斷裂的各級別水系,均發生幅度不等的左旋位錯(圖2)。地貌、水系研究結果顯示,全新世以來斷裂帶的平均走滑速率為6~10mm/a。觀測點附近巖性大多數為灰黑色變質灰巖和泥石混合物,由于1920年海原大地震,地質構造比較復雜,巖層比較破碎。

圖2 海原斷裂帶西華山一帶地質體左旋位移平面圖Fig.2 The left鄄hand displacement plan of the West Huashan area at the Haiyuan fault zone
2017年7月27日,臺站值班人員發現海原鉆孔傾斜資料NS向加速下降,下降幅度102毫角秒,其中7月27日下降幅度65毫角秒,7月28日—8月2日處在低值,且出現新的穩定點;EW向加速上升,上升幅度110毫角秒,其中7月27日上升幅度最大,達74毫角秒,7月28日—8月2日處在上升后穩定期,且出現新的穩定點(圖3)。8月2日,臺站組織技術人員赴海原小山形變觀測點進行鉆孔傾斜資料異常落實。

圖3 海原小山鉆孔傾斜分鐘值曲線(20170720—20170803)Fig.3 The minute value curve of the Haiyuan Xiaoshan borehole tilt anomaly
2017年8月2日,海原地震臺組成異常核實小組,赴現場進行核實調查,并初步制定了異常核實工作方案:(1)對局部觀測環境進行調查;(2)對儀器工作電源進行檢查;(3) 儀器標定記錄、日志記錄等信息的核實;(4)對異常信度進行分析。
異常核實小組對海原小山鉆孔傾斜觀測系統進行了檢查,儀器觀測系統基本正常,主要落實情況如下:
(1) 供電電源檢查
異常核實小組對海原小山鉆孔傾斜供電系統進行了檢查。結果顯示交流供電224.2 V,穩壓電源供電41.17 V,電源供電正常。對觀測電源各接頭進行按壓,保證接頭接觸良好。
(2)數采記錄系統檢查
異常核實小組人員對海原小山鉆孔傾斜儀器數采記錄系統進行檢查,發現NS向有所偏大,粗調后恢復。
(3) 標定格值變化
海原小山鉆孔傾斜格值變化在2013年7月維修之前變化穩定,2013年7月維修之后NS向存在規律的起伏變化,但總體穩定,EW向變化穩定。2016年6月、12月,2017年6月三次標定結果變化穩定。
(4) 日志檢查
查閱了儀器觀測日志,鉆孔傾斜日志記錄完整,記錄系統觀測正常。
3.2.1 儀器所在基巖的巖性情況
海原小山鉆孔傾斜井深60.24 m,儀器探頭安裝在約51 m深處。探頭所處巖性為灰黑色灰巖,巖性完整。
3.2.2 測點周邊的環境變化
海原小山鉆孔傾斜觀測點位于寧夏海原南華山自然保護區內,高程2220 m,附近人煙稀少,對測點周邊5 km范圍進行普查,結果顯示測點周邊無大型的施工作業。
經初步核實,認為2017年7月27日起發生的海原小山鉆孔傾斜NS向加速下降異常、EW向加速上升異常與觀測系統、降雨、氣溫、氣壓等相關性不大;海原小山附近為山區,環境干擾較少,但周邊存在抽水井。異常核實小組赴海原縣自來水公司調查了附近機井抽水情況。海原小山周邊存在抽水機井,主要供縣城用水,其中南部區域機井距離小山觀測點平均4.7 km,中部區域機井距離小山觀測點平均5.5 km,北部區域機井距離小山觀測點6.4 km。小山觀測點鉆孔布置于山腰處,為基巖場區,因此異常組認為抽水對小山形變的干擾較小,排除了機井抽水對形變異常的影響。為了進一步驗證這一推論是否成立,本文對抽水是否影響形變異常進行了以下深入分析:
(1) 從實際觀測可得,2017年7月27日起海原小山鉆孔傾斜EW向加速上升,上升幅度110 ms,其中7月27日上升幅度最大,達74 ms;7月28日—8月2日處在上升后穩定期。且出現新的穩定點;海原小山鉆孔傾斜NS向加速下降,下降幅度102 ms,其中7月27日下降幅度65 ms;7月28日—8月2日處在低值,且出現新的穩定點。根據實地調查,海原小山周圍共有18口機井,距海原小山平均距離為5.5 km,抽水量為800 m3/d·口,18口機井每天抽水總量為18口×800 m3/d·口=1.44×103m3,根據二維載荷均勻分布模型[2],在5km范圍內,每天1.44×103m3抽水量產生的影響僅為0.008 ms,遠小于此次異常的變化幅度,所以可以推定海原小山周圍的機井抽水量變化不會對海原小山鉆孔傾斜產生影響。
(2) 海原小山周邊 18口井抽水量為800 m3/d·口,根據最小“安靜”距離理論的三維集中載荷模型,可以推算抽水量為800m3/d抽水井的最小“安靜”距離小于50 m,而抽水井距海原小山測點的距離最小為4.7 km,遠大于抽水井的最小“安靜”距離,所以可以進一步推斷海原小山周圍的機井抽水量變化不會對海原小山鉆孔傾斜產生影響[3]。
如圖4所示,海原小山鉆孔傾斜異常時段期間,鉆孔氣壓及水位并無明顯的異常變化與之相對應,而在異常時段內,27日至28日出現了連續降雨,但降雨量相對較少,因此同樣可排除降雨對鉆孔傾斜異常的影響。由于海原小山測點距抽水井距離較遠、且抽水量相對較小,因此抽水對海原小山影響很小,可忽略不計,可以推斷海原小山鉆孔傾斜不存在明顯的環境干擾。對海原小山鉆孔傾斜進行潮汐分析和內精度分析,自觀測以來,鉆孔傾斜潮汐因子和內精度變化穩定,觀測數據真實可靠,曾連續在國家局資料評比中獲獎。

圖4 海原小山鉆孔傾斜與鉆孔氣壓、鉆孔水位、降雨的對比曲線Fig.4 A contrast curve between borehole tilt and drilling pressure,drilling water level,and rainfall of Haiyuan Xiaoshan
2017年7月27日起海原小山鉆孔傾斜出現異常,異常形態為NS向加速下降后恢復、EW向加速上升后恢復,期間伴有階躍式變化。類似于此的異常變化,自觀測以來有兩次加速下降、階躍式異常對應地震,具體分析如下:
(1) 2016年3月31日,海原小山鉆孔傾斜NS向出現加速上升,上升幅度約31 ms,EW出現加速下降,下降幅度約43 ms(圖5),2016年4月11日發生寧夏中寧MS3.9級地震,震中距離82 km。

圖5 海原小山鉆孔傾斜分鐘值曲線(中寧MS3.9級地震前)Fig.5 The minute value curve of the Haiyuan Xiaoshan borehole tilt anomaly(Before Zhongning MS3.9)
(2) 2016年11月10日,海原鉆孔傾斜NS向出現加速上升,上升幅度約28 ms,EW向出現加速下降,下降幅度約36 ms(圖6),2016年11月13日發生寧夏海原MS2.3級地震,震中距離28 km。

圖6 海原小山鉆孔傾斜分鐘值曲線(海原MS2.3級地震前)Fig.6 The minute value curve of the Haiyuan Xiaoshan borehole tilt anomaly(Before Haiyuan MS2.3)
由于排除了觀測系統、觀測環境等干擾因素的影響,加之分析了兩次典型的異常應震案例,分析認為此次異常可能為震前異常。
綜合分析認為,海原小山鉆孔傾斜2017年7月27日出現的加速變化真實可靠,無觀測環境、氣象因素干擾,可作為短期異常進行跟蹤,需密切關注海原周邊震情發展,臺站工作人員需做好應急準備。
2016年以來,海原小山鉆孔傾斜出現3月31日、11月10日兩次明顯的加速下降變化,兩次變化之后分別發生了寧夏中寧Ms3.9級、寧夏海原MS2.3級地震,此次異常變化較前兩次略有不同,出現NS向加速下降—恢復、EW向加速上升—恢復的過程,期間伴有階躍式變化,總變化幅度較大,達215 ms,短臨形勢較前兩次嚴峻,需結合地球物理場、地球化學場等其他資料密切跟蹤、分析其變化機理。