中国地震  2022, Vol. 38 Issue (4): 772-784
工业开采诱发地震的应力降研究现状和主要科学议题
姜丛1, 尹欣欣1,2, 蒋长胜1, 翟鸿宇1, 王骞1, 张延保1, 来贵娟1, 尹凤玲1     
1. 中国地震局地球物理研究所,北京 100081;
2. 甘肃省地震局,兰州 730000
摘要:诱发地震及其灾害风险已成为非常规油气资源开发和废水回注等新型工业活动顺利实施的重要威胁。应力降是反映地震破裂前后断层上平均应力差的物理量,对诱发地震中应力降的系统研究有助于揭示诱发地震的震源特征、高频地震动特征、构造背景和地震危险性。本文系统地总结了近20年来工业开采诱发地震中应力降的相关研究进展,分别从计算所用理论方法、获得的平均应力降数值、影响应力降数值的主客观因素、围绕应力降的重要科学讨论等方面做了归纳分析。介绍了直接进行谱拟合、经验格林函数(EGF)、谱叠加和广义反演3种适用性的计算方法,通过案例考察了工业开采诱发地震的应力降平均值,分析影响应力降测定数值的主客观因素,并总结应力降空间分布依赖性、与震源机制类型相依性、自相似性、深度依赖性等4个主要科学问题。同时,阐述了诱发地震应力降的测定可靠性仍未突破、主要科学问题存在广泛争议、在新型工业开采活动中具有重要应用潜力等现状。本文的总结分析可为从事新型能源开发、地震安全监管和科学研究领域的企业、管理人员和科学研究人员提供参考。
关键词水力压裂    诱发地震    震源参数    应力降    减轻地震灾害风险    
Research Status and Major Scientific Issues in Stress Drop in Mining-induced Earthquake
Jiang Cong1, Yin Xinxin1,2, Jiang Changsheng1, Zhai Hongyu1, Wang Qian1, Zhang Yanbao1, Lai Guijuan1, Yin Fengling1     
1. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
2. Gansu Earthquake Agency, Lanzhou 730000, China
Abstract: Induced earthquake and the disaster risk have become serious threats to the smooth implementations of newly industrial activities such as unconventional oil and gas extraction and wastewater reinjection. The stress drop, which reflects the average stress difference on the fault before and after the earthquake rupture, is helpful to reveal the source characteristics, high-frequency ground motion characteristics, tectonic setting and seismic risk of induced earthquake. In this study, we systematically summarize the research progress of stress drop that is related to mining-induced earthquake during the past 20 years, referring to the theory and method used in calculation, the estimated average stress drop values, the subjective and objective factors that affect the stress drop, and important scientific issues. We introduce three applicable calculation methods, including direct spectral fitting, the Empirical Green's function(EGF), and spectrum superposition and generalized inversion. We study the average stress drop that is caused by mining-induced earthquake through case investigation, analyze the subjective and objective factors affecting the estimation of stress drop, and summarize four major scientific issues, including the space dependency, focal mechanism type dependency, self-similarity, and depth dependency. We also mention that there is no breakthrough in the reliability of measured stress drop; the main scientific issues are widely disputed; and it has important application potential in newly mining activities. We hope the summary and analysis in this study will provide some scientific references for enterprises, managers and researchers who are engaged in the fields of new energy development, seismic safety supervision and scientific research.
Key words: Hydraulic fracturing     Induced earthquake     Source parameters     Stress drop     Earthquake disaster risk reduction    
0 引言

在地壳浅部采用水力压裂等方式进行的页岩气和干热岩等非常规能源开采、废水回注以及抽取方式进行的矿山开采,均可造成地壳应力状态的扰动并产生诱发地震,给这些工业活动的顺利实施以及地表社区与重大基础设施的地震安全带来较大挑战。应力降Δσ是描述地震震源特征的重要参量,被定义为地震破裂前后断层上的平均应力差(Hanks,1979)。由于应力降可被用来考察高频地震动特征(Boore,1983)、区分所在地区的构造背景(Allmann et al,2009; Boyd et al,2017)、根据其自相似性给出关于地震危险性的认识(Abercrombie,1995; Oth et al,2010),因此其在天然地震的地球动力学研究和减轻地震灾害风险研究中广受青睐。在工业开采诱发地震活动中,应力降还可用来了解水力压裂等工业活动与诱发地震震源性质之间因果关系(Yu et al,2020),以及作为重要参量来揭示断层活化现象(Chen et al,2020),这使其成为目前诱发地震研究中的重点内容之一。尤其是开采活动中布设大量的密集地震观测网络,为更深层次利用应力降提供了新机遇(Shearer et al,2019)。

最初的Δσ是根据断层力学的静态理论定义的,这是由于断层可以用弹性介质中的剪切位错来表示,地震释放的断层附近累积的应变变化量约为$,则断层上的平均应力降(Stein et al,2009)为

$ \Delta \sigma \approx \mu \bar{U} / L $ (1)

其中,L为断层特征尺度,μ为剪切模量,$ \bar{U}$为断层的平均位错,表明应力降仅与断层的物理参数有关(Holmgren et al,2019)。在弹性位错理论发展起来后,地震矩的概念被提出并被定义为断层面上的平均位错$ \bar{U}$与断层面面积S的乘积,即$ M_0=\mu \bar{U} S$。由此,应力降与观测获得的标量地震矩建立起联系,即

$\Delta \sigma=c M_0 / S / L $ (2)

从而实现了直接由地震波测量应力降的过程。由上述定义可见: ①受到震源模型假设的影响,Δσ的具体表达形式依赖于震源模型,例如半径为R的圆形断层上的应力降为Δσ=7M0/16R3,长为L宽为w的矩形断层上的应力降为Δσ=2M0/πw2L; ②由于断层面积SL2,应力降Δσ与1/L3相关,其测定结果的不确定性很大(Allmann et al,2007); ③由于应力降的概念基于静态位错理论,仅考虑了断层的相对移动,没有考虑断层的破裂,因此通常所说的静态应力降要比实际的应力降小。

鉴于上述关于应力降Δσ的定义和物理含义,以及其在新型工业活动和诱发地震灾害风险管控中的重要价值,有必要系统全面地总结目前关于诱发地震应力降研究的新进展、技术发展趋势和已获得的科学认识。而在人为工业活动中的这些进展和认识,对天然地震的研究也同样有重要启发意义。本文尝试对近20年的相关研究进行总结,以期为新型工业活动伴随的诱发地震的减灾研究以及天然地震的研究提供参考。

1 工业开采诱发地震应力降的计算

包括应力降在内的震源参数的测定,需要分离或同时测定几何扩散、介质衰减、场地响应等波形信号的耦合信息。新型工业活动诱发地震多为震级较小的微震事件,波形信号的低信噪比给测定工作带来较大难度。目前已发展的测定方法有如下三类:

(1) 谱拟合方法是计算诱发地震的应力降的最为直接的方法。一般首先利用多锥形窗方法(Percival et al,1993)分别对三分量的地动速度谱进行平滑处理,并利用Boatwright点源模型拟合地动速度谱(Boatwright,1980),即

$\dot{u}^{t h}\left(f ; M_0, f_c, Q_c\right)=\frac{R_c}{\rho V_c^3 R} \frac{f M_0}{\left(1+\left(f / f_c\right)^4\right)^{0.5}} \exp \left(-\frac{\pi R f}{Q_c V_c}\right) $ (3)

其中,f为频率,M0为地震矩,fc为拐角频率,Qc为品质因子。式(3)中的震源参数[M0fcQc]的初值可采用网格搜索法,以及求解理论地动速度谱与实际观测谱差值的L1范数最小值来获得,一些研究还采用郊狼算法(Coyote Optimization Algorithm,COA,Pierezan et al,2018)等智能全局搜索方法进一步获得更为准确的震源参数估计值(Kwiatek et al,2011),然后利用Brune(1970)圆盘破裂模型Δσ=7M0/(16r03)以及震源破裂半径r0=kVc/fc的定义进一步计算应力降(Eshelby,1957; Madariaga,1979)。谱拟合方法可分别对P波和S波求解震源参数,其中对波速k的取值合理性已有大量的研究(Abercrombie,1995; Kwiatek et al,2011; Chaves et al,2020)。由于谱拟合方法并未分离震源项、衰减项和场地响应,若采用单事件单台站的直接拟合,计算结果有可能具有较大的不确定性。

(2) 单事件的经验格林函数方法(EGF)是诱发地震应力降计算的重要方法,其优点是可一定程度提高拐角频率fc和标量地震矩M0的测量精度(Kwiatek et al,2015)。EGF方法的原理是,对拟求解震源参数的地震,找到位置接近的小震级事件(震级差ΔM>1.0)作为经验格林函数,利用谱阈的反卷积方法,消除待求解事件与EGF事件相同的场地响应、传播路径衰减项,从而获得震源谱。其中拟求解事件i与EGF事件j的地动速度谱的比值为

$ \varPsi^{i j}\left(f ; M_0^i, M_0^j, f_c^i, f_c^j\right)=\frac{\dot{u}^{t h}\left(f ; M_0^i, f_c^i\right)}{\dot{u}^{t h}\left(f ; M_0^j, f_c^j\right)}=\frac{M_0^i}{M_0^j}\left(\frac{1+\left(f / f_c^i\right)^4}{1+\left(f / f_c^j\right)^4}\right)^{0.5} $ (4)

对式(4)一般可采用最小二乘法拟合理论谱与实际的观测比值谱的L2范数最小值来获得参数结果(Holmgren et al,2019; Wang et al,2020)。在实际应用中,能否找到严格吻合EGF定义的小地震是制约经验格林函数方法的重要因素(Huang et al,2019),此外对EGF方法计算结果分辨率的限度也有较多讨论(Abercrombie,2015; Shearer et al,2019)。

(3) 利用谱叠加和广义反演方法(Castro et al,1990; Shearer et al,2006; Oth et al,2011)获得更具普遍意义的经验格林函数(EGF)是计算应力降的重要方法。例如,Warren等(2002)以及Shearer等(2006)发展了一种P波谱叠加方法,可以将震源谱与几何扩散、传播路径的衰减、场地响应分离出来。该方法通过将同一台站大量的地震记录进行叠加,迭代后消除几何扩散、路径衰减等差异项,再利用震级分组的震源项与理论震源谱的拟合获得修正后的经验格林函数,从而更为完整地去除几何扩散、路径衰减和场地响应,获得每个地震的震源谱,并利用Brune圆盘震源模型计算应力降等震源参数。Trugman等(2017)在该方法的基础上,发展了非参数重采样方法来估计震源参数的不确定性。为提高谱叠加和广义反演方法中的震源谱参数反演的可靠性,Picozzi等(2017)还发展了一种基于数据驱动策略的无参数计算方法,其核心是设计了震源谱lg10(S(f))对各个震源参数的无量纲的偏导数作为频域相依的敏感度函数,作为反演的权重函数。

2 工业开采诱发地震应力降的平均值

目前在应力降的数值上已取得一定共识。工业开采诱发地震的应力降通常较震源深度更深的天然地震的应力降值偏小(Fehler et al,1991; Abercrombie et al,1993; Hua et al,2013; Hough,2014)。已有的对全球天然地震的研究表明,应力降至少在三个数量级上变化(例如,Shearer et al,2006; Allmann et al,20072009),数值通常在1~100MPa之间(Abercrombie,1995),且符合正态分布(Allmann et al,2009)。而对工业开采诱发地震,应力降数值则通常在0.01~10MPa之间。本文罗列了一些工业开采区诱发地震应力降的测定结果,以便进行比较研究。

在北美有大量的诱发地震应力降测定案例。例如,Wu等(2019)对2015年美国俄克拉荷马州加斯里(Guthrie)MW4.0诱发地震应力降的测定结果为3.4MPa,这一数值接近Wu等(2018)测定的俄克拉荷马州其他诱发地震应力降的中值(约2MPa)。Kwiatek等(2015)采用Madariaga震源模型观测到美国加利福尼亚州盖尔瑟斯地热田Prati-9井附近的一组诱发地震,震级分布在MW2.0~3.8之间,这些地震的平均应力降约为7MPa,但明显低于Viegas等(2011)采用经验格林函数方法对同一水力压裂注入井测定的应力降中值(28MPa)。Jeong等(2020)对美国得克萨斯州沃思堡盆地的90次诱发地震(ML2.0~3.9)的应力降测定结果显示,平均应力降约为4.46MPa。在加拿大不列颠哥伦比亚省东北部的蒙特尼区块,通过P波估算的应力降中值为0.08MPa,对数标准差为0.8(以lg为单位),同时利用S波估算的应力降中值为0.03MPa,对数标准差为0.7(Yu et al,2020)。在加拿大西部沉积盆地(WCSB),Holmgren等(2019)利用S波测定了87次震级为M2.3~4.2的诱发地震的应力降,发现中小震级事件的平均应力降为7.5±0.5MPa,但测定结果在0.2~98MPa之间变化,且有一个极端事件的应力降测定值为370MPa。

在欧洲和亚洲,部分页岩气和干热岩开采诱发地震的应力降得到了测定和研究。其中,Kettlety等(2020)测定了英国兰开夏郡普雷斯顿新路PNR-1z井2018年发生的震级为ML1.5的诱发地震,平均应力降值为1.0MPa。对于瑞士巴塞尔地热田附近2006年发生的MW 3.2±0.1诱发地震,以及圣加仑地区2013年发生的MW 3.4±0.1诱发地震事件,Edwards等(2015)测得的平均应力降分别为3.5MPa和3.0MPa。在芬兰赫尔辛基,Kwiatek等(2019)对OTN-3井附近震级为MW2.0的诱发地震测得的平均应力降为1.6MPa。在亚洲的韩国浦项地区,Song等(2019)对PH02井附近的诱发地震计算获得的平均应力降约为2MPa,与Chai等(2020)得到的1.92MPa结果一致。

表 1给出了这些诱发地震应力降测定结果。由上述案例可见,对于多数震级在1.5~4.0的诱发地震,平均应力降分布在0.03~7MPa之间,且多数主要分布在1.0~7MPa之间。目前认为诱发地震应力降的数值相比天然地震更小的现象,主要是由诱发地震发生在地壳更浅部造成的。

表 1 全球部分地区工业开采诱发地震的平均应力降测定结果
3 影响应力降测定数值的因素

诱发地震的发生机制与天然地震明显不同,至少包括注入流体影响岩石孔隙压力、注入的较冷流体与高温岩石相互作用引起热弹性应变和裂缝滑动(Rawal et al,2014; 尹欣欣等,2021)等较为独特的机制,因此影响应力降数值的因素可能也有所不同。下面分别从主观因素造成测定结果的不可靠性以及客观因素造成测定结果的复杂性做出简要总结。

震源模型和常数参数的设置,是应力降计算的重要假设,可直接影响应力降的计算结果。例如,Yu等(2020)使用圆盘震源模型(Brune,1970)和取值S波的k =0.37,相比于使用Madariaga模型(Madariaga,1976)以及取值k =0.21,其计算得到的应力降较后者小约5.5MPa。Yamada等(2007)发现使用圆盘震源模型(Brune,1970)要比使用Madariaga模型(Madariaga,1976)获得的应力降小,并分析认为使用Madariaga模型更为合理。破裂速度的取值也是应力降计算中的重要影响因素,较慢的破裂速度将导致震源半径的估计值较小、应力降的估计值较大。而破裂速度的取值通常带有明显的主观性,例如在Madariaga模型中常取剪切波速3.9km/s的90%,这在采用中值计算平均应力降的研究中影响更为明显(Allmann et al,2009)。地震波传播介质的品质因子Q值的计算和取值,也会显著影响应力降测定的可靠性(Ko et al,2012)。Q值对岩石的力学性质非常敏感,其空间变化可以作为水力压裂影响局部应力状态的指标(Yu et al,2020),但当对衰减校正不充分时常常可得到较低的应力降测定结果(Tomic et al,2009)。由于上述震源模型和参数设置的人为主观因素影响,Holmgren等(2019)Yu等(2020)Wang等(2020)建议,应力降的绝对数值只有在使用相同的震源模型和参数设置的条件下才可以进行比较。

一些客观因素也显著影响了应力降的分布差异,包括构造环境的差异、存在的水量、破裂速度的变化、断层面正应力的变化等(Tomic et al,2009)。这些因素导致应力降数值时空分布出现差异,另一方面也使得根据应力降数值反推地下施工区或储层应力状态变得更为复杂。这些影响因素至少包括: ①水热型地热田和增强型地热田(干热岩)的高温特征,对地壳的摩擦强度等力学行为有重要作用,例如Oth(2013)对日本地区的应力降的研究表明,低应力降的数值在空间上与高热流区有关; ②诱发地震在深度分布上的差别以及流体介入地震破裂过程的程度也能明显影响应力降的数值,这是由于震源深度较浅和流体注入时,可能会降低有效应力,从而导致剪切应力更容易克服莫尔-库仑破坏准则进而产生诱发地震(Goertz-Allmann et al,2013); ③诱发地震的应力降在空间上具有方向性差异,且对中等大小的地震(Lui et al,2019)和微震(Folesky et al,2016)均可观测到。而当地震记录的台站覆盖较为稀疏时,应力降测量方式不当则可能带来结果的较大不确定性,例如仅从断层的前向或后向测量,可能会将拐角频率fc的测定结果带来2倍误差,并导致应力降结果达到8倍偏差(Holmgren et al,2019)。

4 与应力降有关的重要科学议题 4.1 诱发地震应力降的空间相关性

探讨诱发地震应力降的空间相关性问题具有重要的科学价值。这是由于断层带和裂隙上可能存在应力的不均匀性(Chen et al,2020; 姜丛等,2022),应力降的空间变化可以揭示断层强度的非均质性,以及反映断层带的水文地质性质,这有助于理解诱发地震的发生机制。

应力降是否具有“近低、远高”分布现象近年来在诱发地震研究中成为热点问题。Yu等(2020)在对加拿大不列颠哥伦比亚省东北部的蒙特尼区块水力压裂井周边3km范围内的484次地震(M-1.0~3.0)的研究中发现,近井端的应力降(0.1~1MPa)要低于远井端的应力降(1~10MPa),且约小一个数量级。类似的现象也在加拿大阿尔伯塔省克鲁克湖(Crooked Lake)附近的注水诱发地震(Clerc et al,2016),以及瑞士巴塞尔和萨尔瓦多柏林地区的增强型地热系统(EGS)开采区诱发地震研究(Goertz-Allmann et al,2011; Kwiatek et al,2014)中被证实。上述例子中应力降数值与井的距离相关的现象,在统计学意义上也得到了证明(Yu et al,2020)。Yu等(2020)认为,应力降的这种空间分布特点可能是由于井筒附近水力压裂受损岩石的破裂速度较慢,裂隙密度高或孔隙压力较高,阻止了井筒附近地壳岩石存储和释放较大应力,而远端地震事件的应力扰动受水力压裂的影响则相对有限或几乎可以忽略。

应力降是否具有空间分布的恒定性也是诱发地震研究的关注点之一。在对加拿大不列颠哥伦比亚省东北部的蒙特尼区块水力压裂诱发地震的研究中,Yu等(2020)还发现无论是发生在距离井筒近处还是远处的地震,应力降的数值随着时间变化而大致不变。Kwiatek等(2015)Picozzi等(2017)分别采用不同的应力降计算方法,在对美国Geysers地热田Prati-9井诱发地震的分析中,也表示并未观察到静态应力降随时间的明显变化。这种空间分布的恒定性,可能表明具有相似流变特性的岩石中的动态破裂通常具有相对恒定的应力降(Yu et al,2020)。还有一些研究发现,地震序列的应力降具有随时间变化的特点,主要表现在主震后应力降随时间增加,在余震序列的后期观察到的应力降总体较高(Sumy et al,2017; Chen et al,2020; Goertz-Allmann et al,2011),Goertz-Allmann等(2011)认为这是孔隙压力过高导致的。但这些序列的应力降时序变化是否在同一位置、余震序列活动的短期调整变化是否能影响空间分布的恒定性,还缺乏更深入的研究。

4.2 诱发地震应力降与震源机制类型的相依性

即使在天然构造地震研究中,应力降与震源机制类型的依赖性也是重要的研究课题。根据安德森断层理论,逆断层作用的剪应力最高,正断层作用的剪应力最低,走滑型断层的剪应力介于前两者之间,尽管应力降和绝对应力之间的关系尚不清楚,但传统观点认为应力降和绝对应力在与震源机制类型的关系上相同(McGarr,1984; McGarr et al,2002)。但也正因为应力降和绝对应力之间的关系不明确,一些新的观测现象在逐渐揭示其复杂性。例如Allmann等(2009)在对地震应力降的全球变化的研究中观察到了应力降中值对震源机制类型有依赖性,其中走滑事件的应力降最高,中值约为10MPa,而正断层和逆断层事件的应力降较低,中值为2~3MPa,并未显示逆冲型地震具有更高的应力降。

对于诱发地震,Kwiatek等(2015)对美国盖尔瑟斯地热田Prati-9井水力压裂诱发地震的研究表明,尽管正断层型地震的应力降略低于走滑型地震的应力降,但由于走滑型地震在储层深部分布的数量更多,因此这种应力降与震源机制类型的相关性难于验证。研究这一问题的难点在于影响应力降数值大小的因素,例如在一些应力降的计算中假设断层破裂速度是恒定的,而Allmann等(2009)采用可变的断层破裂速度来计算应力降时,可以解释走滑型地震具有较高的应力降可能是由更快的断层破裂速度造成的。

4.3 诱发地震应力降的非自相似性

应力降是否具有自相似性在物理上具有特殊意义,其实质是大地震和小地震是否具有相似的动态破裂过程。在实际应用中,同一地区小地震的应力降空间分布是否与中等以上地震存在相关性、能否用平均的应力降代表这一地区的总体应力降水平,直接关系到是否可用数量更多的小震级的诱发地震来预测未来破坏性诱发地震的应力降,并减少概率性地震危险性评估和地震预警中的地面运动不确定性(Hardebeck,2020)。

对应力降的自相似性/非自相似性研究,最重要的研究方式是考察应力降是否与震级或标量地震矩存在相关性。Holmgren等(2019)在对加拿大西部沉积盆地87次M 2.3~4.2级诱发地震应力降的测定中,发现应力降随震级增大呈上升趋势。Novakovic等(2018)Chen等(2020)在对美国俄克拉荷马州以及Picozzi等(2017)对加利福尼亚州盖尔瑟斯地热田进行诱发地震研究中,也均发现较大震级的诱发地震具有相对更高的应力降。Chai等(2020)对韩国浦项诱发地震序列的应力降研究中也获得了同样结论。其他案例还包括Mayeda等(2005)Venkataraman等(2006)等的研究。概括起来,目前支持工业开采诱发地震应力降具有自相似性的案例总体上较少,仅有Yamada等(2007)对南非金矿开采诱发的20次微震的静态应力降的研究等少数案例,此外,尽管Picozzi等(2017)在美国加州盖尔瑟斯地热田诱发地震的研究中观察到了大地震群体具有自相似趋势,但对单个地震序列的研究往往发现强烈的非自相似行为。

4.4 诱发地震应力降的深度依赖性

在天然地震的应力降研究中,许多案例显示出应力降与震源深度存在依赖性(Hardebeck et al,2009; Allmann et al,2007; Baltay et al,2017),并认为这可能与断层的摩擦强度沿深度的变化有关。工业开采诱发地震多发生在地下3~10km深度的地壳浅部,表 1给出的大量案例也表明诱发地震的应力降总体相比天然地震偏小,因此目前普遍持有的观点是造成诱发地震应力降比天然地震偏小的原因也在于震源深度的影响(Satoh,2006; Hough,2014; Yenier et al,2015; Zhang et al,2016; Atkinson et al,2017)。

但对无论是诱发地震还是天然地震是否存在应力降大小的深度依赖性目前还有较大争议,例如Wang等(2020)认为在计算应力降时假设剪切波速恒定可能导致人为的深度依赖性,特别是发生在地壳浅部5km以上的地震,Allmann等(2007)也持有类似观点。Holmgren等(2019)进一步解释认为,由于应力降与速度的立方呈正相关,也认为由于剪切波速随深度增加,应力降的估计对假定的震源深度很敏感。在具体的应力降测定研究中,Allmann等(2009)对全球范围的天然地震研究表明,尽管存在较大的分散性,但平均应力降估计值随深度的变化很小。进一步从不同震源机制类型的角度分析,Allmann等(2009)还发现正断层事件和走滑事件的平均应力降均没有深度依赖性,而逆断层作用下平均应力降随深度的增加与Bilek等(19981999)在俯冲带地震中发现的结果一致,可以用深度相关的刚度变化来解释。左可桢等(2021)使用谱比法分析四川长宁地区微震震源参数之间的相互关系以及应力降的时空分布特征,也得到应力降总体与震源深度关系并不明显的结果。

尽管正反两种观点的争议较大,表 1给出的大量案例至少给出了诱发地震的应力降数值分布范围更窄、总体比天然地震偏小的初步认识。由于应力降测定结果的主观可靠性、客观影响因素较为复杂,这也使得对这一问题的争论和研究在未来仍具有开放性。

5 结论和讨论

鉴于应力降在研究工业开采诱发地震的震源特征、高频地震动特征、区分构造背景和地震危险性分析等诸多领域的重要价值,本文系统地总结了近20年来相关方向的研究进展。分别从计算所用理论方法、获得的平均应力降数值、影响应力降数值的主客观因素、围绕应力降的重要科学议题4个方面做了总结。获得的主要认识如下:

(1) 在工业开采诱发地震的理论计算方法上,针对诱发地震微震较多、信噪比较低、实时快速计算需求明显等特点,目前主要发展了三类方法,包括直接采用谱拟合方法求取震源参数、利用小震级事件作为经验格林函数(EGF)提高拐角频率fc和标量地震矩M0测量精度以及利用谱叠加和广义反演方法对多台站多事件同时计算震源参数。不同方法尚存在明显的不足之处,这也限制了对应力降数值本身和分布规律探索的过高期待。

(2) 对于工业开采诱发地震的应力降平均值,本文收集了相关案例,结果显示多数震级在1.5~4.0的诱发地震平均应力降分布在0.03~7MPa之间,且多数主要分布在1.0~7MPa之间,应力降数值相比天然地震分布的范围更窄、总体偏小。对这一现象的主要解释集中在诱发地震往往发生在地壳更浅部,以及流体更多地参与诱发地震发生过程等因素上。

(3) 影响工业开采诱发地震测定数值的因素,在测定方法的可靠性上,主要受到采用相对简单的圆盘模型等震源模型假定,以及对波速、断层破裂速度、地壳介质品质因子Q值的具体取值差异等因素影响; 此外,构造环境的差异、存在的水量、破裂速度的变化、断层面正应力的变化等客观因素也影响了应力降的数值本身。

(4) 围绕工业开采诱发地震应力降的科学问题,总结了应力降是否具有“近低、远高”分布特征以及是否存在空间分布恒定性等空间分布依赖性问题,以及诱发地震的应力降是否与震源机制类型存在相依性、是否存在相似性或非自相似性、是否存在深度依赖性等4个问题。其中对应力降存在距离依赖性、存在非自相似性有一定程度的共识,但鉴于应力降测定可靠性和多种客观因素的影响,这些科学问题未来均有较大的开放性、需要持续探索。

工业开采区通常布设更为密集、覆盖程度更好的现代化地震台网,这为诱发地震获得更为可靠的包括应力降在内的震源参数提供了便利。但通过本文的总结可见,受观测技术和测定理论方法的局限,诱发地震的应力降测定的可靠性问题并未获得实质性突破。对诱发地震的应力降很难精确测量似乎成了为数不多的共识性问题(Shearer et al,2019),Stein等(2009)甚至认为当前的静态应力降研究仅反映了观测谱或震源谱的一般特征、无法直接反映震源的物理特性,甚至会误导对真实的断层活动的认识。由于应力降是由拐角频率通过假设的震源模型和参数间接计算得到的,对此Holmgren等(2019)Shearer等(2019)建议应该转而关注拐角频率fc,这是相比应力降更基本、有用且更接近实际数据的震源参数。

但与此同时,应力降在工业开采活动施工和诱发地震风险管控中的独特价值,使得我们对应力降的研究“无法拒绝”。例如在实际应用中,应力降是自发性停滞破裂(self-arrested ruptures)最大震级Mwmax-arr预测中的重要参量(Orlecka-Sikora et al,2020)。而那些相比多数诱发地震应力降数值更大或更小的“异常”地震,可能揭示了更多重要信息,例如局部的低应力降的分布,可能是沿弱化断层段或在断层-流体相互作用明显的地方发生破裂的指标(Chen et al,2013)。Picozzi等(2017)对美国盖尔瑟斯诱发地震的研究也表明,一些具有更大的应力降数值的诱发地震,主要分布在储层深部,且可能与深部断层的重新活化以及区域构造活动有关。陡倾角的断层上分布的低应力降和高地震效率ηSW的事件,则帮助揭示了这些断层在流体压力扩散中扮演的重要作用。这种伴随着不完善和争议的总体现状,也决定了工业开采诱发地震应力降的测定和应用研究在未来仍是重要的前沿科研领域。

致谢: 中国地震局地球物理研究所博士研究生张琰在文献检索上予以帮助,郭祥云高级工程师等专家提供了诸多建议,在此一并表示感谢。
参考文献
姜丛、蒋长胜、尹欣欣等, 2022, 水力压裂诱发地震活动中的b值时空异质性及其应用, 地震地质, 44(5): 1333-1349. DOI:10.3969/j.issn.0253-4967.2022.05.015
尹欣欣、蒋长胜、翟鸿宇等, 2021, 全球干热岩资源开发诱发地震活动和灾害风险管控, 地球物理学报, 64(11): 3817-3836. DOI:10.6038/cjg2021O0448
左可桢、赵翠萍, 2021, 四川长宁地区地震震源参数的时空分布特征, 中国地震, 37(2): 472-482. DOI:10.3969/j.issn.1001-4683.2021.02.019
Abercrombie R, Leary P, 1993, Source parameters of small earthquakes recorded at 2.5km depth, Cajon Pass, southern California: Implications for earthquake scaling, Geophys Res Lett, 20(14): 1511-1514. DOI:10.1029/93GL00367
Abercrombie R E, 1995, Earthquake source scaling relationships from -1 to 5 ML using seismograms recorded at 2.5-km depth, J Geophys Res: Solid Earth, 100(B12): 24015-24036. DOI:10.1029/95JB02397
Abercrombie R E, 2015, Investigating uncertainties in empirical Green's function analysis of earthquake source parameters, J Geophys Res: Solid Earth, 120(6): 4263-4277. DOI:10.1002/2015JB011984
Allmann B P, Shearer P M, 2007, Spatial and temporal stress drop variations in small earthquakes near Parkfield, California, J Geophys Res: Solid Earth, 112(B4): B04305.
Allmann B P, Shearer P M, 2009, Global variations of stress drop for moderate to large earthquakes, J Geophys Res: Solid Earth, 114(B1): B01310.
Atkinson G M, Assatourians K, 2017, Are ground-motion models derived from natural events applicable to the estimation of expected motions for induced earthquakes?, Seismol Res Lett, 88(2A): 430-441. DOI:10.1785/0220160153
Baltay A S, Hanks T C, Abrahamson N A, 2017, Uncertainty, variability, and earthquake physics in ground-motion prediction equations, Bull Seismol Soc Am, 107(4): 1754-1772.
Bilek S L, Lay T, 1998, Variation of interplate fault zone properties with depth in the Japan subduction zone, Science, 281(5380): 1175-1178. DOI:10.1126/science.281.5380.1175
Bilek S L, Lay T, 1999, Rigidity variations with depth along interplate megathrust faults in subduction zones, Nature, 400(6743): 443-446. DOI:10.1038/22739
Boatwright J, 1980, A spectral theory for circular seismic sources; simple estimates of source dimension, dynamic stress drop, and radiated seismic energy, Bull Seismol Soc Am, 70(1): 1-27.
Boore D M, 1983, Stochastic simulation of high-frequency ground motions based on seismological models of the radiated spectra, Bull Seismol Soc Am, 73(6A): 1865-1894.
Boyd O S, McNamara D E, Hartzell S, et al, 2017, Influence of lithostatic stress on earthquake stress drops in North America, Bull Seismol Soc Am, 107(2): 856-868. DOI:10.1785/0120160219
Brune J N, 1970, Tectonic stress and the spectra of seismic shear waves from earthquakes, J Geophys Res, 75(26): 4997-5009. DOI:10.1029/JB075i026p04997
Castro R R, Anderson J G, Singh S K, 1990, Site response, attenuation and source spectra of S waves along the Guerrero, Mexico, subduction zone, Bull Seismol Soc Am, 80(6A): 1481-1503.
Chai G, Yoo S H, Rhie J, et al, 2020, Stress-drop scaling of the 2016 Gyeongju and 2017 Pohang earthquake sequences using coda-based methods, Bull Seismol Soc Am, 110(5): 2047-2057. DOI:10.1785/0120200132
Chaves E J, Schwartz S Y, Abercrombie R E, 2020, Repeating earthquakes record fault weakening and healing in areas of megathrust postseismic slip, Sci Adv, 6(32): eaaz9317. DOI:10.1126/sciadv.aaz9317
Chen X, Abercrombie R E, 2020, Improved approach for stress drop estimation and its application to an induced earthquake sequence in Oklahoma, Geophys J Int, 223(1): 233-253. DOI:10.1093/gji/ggaa316
Chen X W, Shearer P M, 2013, California foreshock sequences suggest aseismic triggering process, Geophys Res Lett, 40(11): 2602-2607. DOI:10.1002/grl.50444
Clerc F, Harrington R M, Liu Y J, et al, 2016, Stress drop estimates and hypocenter relocations of induced seismicity near Crooked Lake, Alberta, Geophys Res Lett, 43(13): 6942-6951. DOI:10.1002/2016GL069800
Edwards B, Kraft T, Cauzzi C, et al, 2015, Seismic monitoring and analysis of deep geothermal projects in St Gallen and Basel, Switzerland, Geophys J Int, 201(2): 1022-1039. DOI:10.1093/gji/ggv059
Eshelby J D, 1957, The determination of the elastic field of an ellipsoidal inclusion, and related problems, Proc Roy Soc Lond. Ser A, Math Phys Sci, 241(1226): 376-396.
Fehler M, Phillips W S, 1991, Simultaneous inversion for Q and source parameters of microearthquakes accompanying hydraulic fracturing in granitic rock, Bull Seismol Soc Am, 81(2): 553-575. DOI:10.1785/BSSA0810020553
Folesky J, Kummerow J, Shapiro S A, et al, 2016, Rupture directivity of fluid-induced microseismic events: Observations from an enhanced geothermal system, J Geophys Res: Solid Earth, 121(11): 8034-8047. DOI:10.1002/2016JB013078
Goertz-Allmann B P, Goertz A, Wiemer S, 2011, Stress drop variations of induced earthquakes at the Basel geothermal site, Geophys Res Lett, 38(9): L09308.
Goertz-Allmann B P, Wiemer S, 2013, Geomechanical modeling of induced seismicity source parameters and implications for seismic hazard assessment, Geophysics, 78(1): KS25-KS39. DOI:10.1190/geo2012-0102.1
Hanks T C, 1979, b values and ω seismic source models: Implications for tectonic stress variations along active crustal fault zones and the estimation of high-frequency strong ground motion, J Geophys Res: Solid Earth, 84(B5): 2235-2242. DOI:10.1029/JB084iB05p02235
Hardebeck J L, 2020, Are the stress drops of small earthquakes good predictors of the stress drops of moderate-to-large earthquakes?, J Geophys Res: Solid Earth, 125(3): e2019JB018831.
Hardebeck J L, Aron A, 2009, Earthquake stress drops and inferred fault strength on the Hayward fault, east San Francisco Bay, California, Bull Seismol Soc Am, 99(3): 1801-1814. DOI:10.1785/0120080242
Holmgren J M, Atkinson G M, Ghofrani H, 2019, Stress drops and directivity of induced earthquakes in the Western Canada sedimentary basin, Bull Seismol Soc Am, 109(5): 1635-1652. DOI:10.1785/0120190035
Hough S E, 2014, Shaking from injection-induced earthquakes in the central and eastern United States, Bull Seismol Soc Am, 104(5): 2619-2626. DOI:10.1785/0120140099
Hua W, Chen Z L, Zheng S H, 2013, Source parameters and scaling relations for reservoir induced seismicity in the Longtan reservoir area, Pure Appl Geophys, 170(5): 767-783. DOI:10.1007/s00024-012-0459-7
Huang Y H, Beroza G C, Ellsworth W L, 2016, Stress drop estimates of potentially induced earthquakes in the Guy-Greenbrier sequence, J Geophys Res: Solid Earth, 121(9): 6597-6607. DOI:10.1002/2016JB013067
Huang Y H, De Barros L, Cappa F, 2019, Illuminating the rupturing of microseismic sources in an injection-induced earthquake experiment, Geophys Res Lett, 46(16): 9563-9572. DOI:10.1029/2019GL083856
Huang Y H, Ellsworth W L, Beroza G C, 2017, Stress drops of induced and tectonic earthquakes in the central United States are indistinguishable, Sci Adv, 3(8): e1700772. DOI:10.1126/sciadv.1700772
Jeong S, Stump B W, DeShon H R, 2020, Spectral characteristics of ground motion from induced earthquakes in the Fort Worth Basin, Texas, using the generalized inversion technique, Bull Seismol Soc Am, 110(5): 2058-2076. DOI:10.1785/0120200097
Kettlety T, Verdon J P, Werner M J, et al, 2020, Stress transfer from opening hydraulic Fractures controls the distribution of induced seismicity, J Geophys Res: Solid Earth, 125(1): e2019JB018794.
Ko Y T, Kuo B Y, Hung S H, 2012, Robust determination of earthquake source parameters and mantle attenuation, J Geophys Res: Solid Earth, 117(4): B04304.
Kwiatek G, Bulut F, Bohnhoff M, et al, 2014, High-resolution analysis of seismicity induced at Berlín geothermal field, El Salvador, Geothermics, 52: 98-111. DOI:10.1016/j.geothermics.2013.09.008
Kwiatek G, Martínez-Garzón P, Dresen G, et al, 2015, Effects of long-term fluid injection on induced seismicity parameters and maximum magnitude in northwestern part of The Geysers geothermal field, J Geophys Res: Solid Earth, 120(10): 7085-7101. DOI:10.1002/2015JB012362
Kwiatek G, Plenkers K, Dresen G, et al, 2011, Source parameters of picoseismicity recorded at Mponeng Deep Gold Mine, South Africa: Implications for scaling relations, Bull Seismol Soc Am, 101(6): 2592-2608. DOI:10.1785/0120110094
Kwiatek G, Saarno T, Ader T, et al, 2019, Controlling fluid-induced seismicity during a 6.1-km-deep geothermal stimulation in Finland, Sci Adv, 5(5): eaav7224.
Lui S K Y, Huang Y H, 2019, Do injection-induced earthquakes rupture away from injection wells due to fluid pressure change?, Bull Seismol Soc Am, 109(1): 358-371. DOI:10.1785/0120180233
Madariaga R, 1976, Dynamics of an expanding circular fault, Bull Seismol Soc Am, 66(3): 639-666. DOI:10.1785/BSSA0660030639
Madariaga R, 1979, On the relation between seismic moment and stress drop in the presence of stress and strength heterogeneity, J Geophys Res: Solid Earth, 84(B5): 2243-2250. DOI:10.1029/JB084iB05p02243
Mayeda K, Gök R, Walter W R, et al, 2005, Evidence for non-constant energy/moment scaling from coda-derived source spectra, Geophys Res Lett, 32(10): L10306. DOI:10.1029/2005GL022405
McGarr A, 1984, Scaling of ground motion parameters, state of stress, and focal depth, J Geophys Res: Solid Earth, 89(B8): 6969-6979. DOI:10.1029/JB089iB08p06969
McGarr A, Fletcher J B, 2002, Mapping apparent stress and energy radiation over fault zones of major earthquakes, Bull Seismol Soc Am, 92(5): 1633-1646. DOI:10.1785/0120010129
Novakovic M, Atkinson G M, Assatourians K, 2018, Empirically calibrated ground-motion prediction equation for Oklahoma, Bull Seismol Soc Am, 108(5A): 2444-2461. DOI:10.1785/0120170331
Orlecka-Sikora B, Cielesta S, 2020, Evidence for subcritical rupture of injection-induced earthquakes, Sci Rep, 10: 4016. DOI:10.1038/s41598-020-60928-0
Oth A, 2013, On the characteristics of earthquake stress release variations in Japan, Earth Planet Sci Lett, 377-378: 132-141. DOI:10.1016/j.epsl.2013.06.037
Oth A, Bindi D, Parolai S, et al, 2010, Earthquake scaling characteristics and the scale-(in)dependence of seismic energy-to-moment ratio: Insights from KiK-net data in Japan, Geophys Res Lett, 37(19): L19304.
Oth A, Bindi D, Parolai S, et al, 2011, Spectral analysis of K-net and KiK-net data in Japan, Part Ⅱ: On attenuation characteristics, source Spectra, and site response of borehole and surface stations, Bull Seismol Soc Am, 101(2): 667-687. DOI:10.1785/0120100135
Percival D B, Walden A T, 1993, Spectral Analysis for Physical Applications, Cambridge: Cambridge University Press.
Picozzi M, Oth A, Parolai S, et al, 2017, Accurate estimation of seismic source parameters of induced seismicity by a combined approach of generalized inversion and genetic algorithm: Application to The Geysers geothermal area, California, J Geophys Res: Solid Earth, 122(5): 3916-3933. DOI:10.1002/2016JB013690
Pierezan J, Coelho L D S, 2018. Coyote optimization algorithm: A new metaheuristic for global optimization problems. In: 2018 IEEE Congress on Evolutionary Computation(CEC). Rio de Janeiro, Brazil: IEEE, 1~8.
Rawal C, Ghassemi A, 2014, A reactive thermo-poroelastic analysis of water injection into an enhanced geothermal reservoir, Geothermics, 50: 10-23. DOI:10.1016/j.geothermics.2013.05.007
Satoh T, 2006, Influence of fault mechanism, depth, and region on stress drops of small and moderate earthquakes in Japan, Struct Eng/Earthq Eng, 23(1): 125s-134s.
Shearer P M, Abercrombie R E, Trugman D T, et al, 2019, Comparing EGF methods for estimating corner frequency and stress drop from P-wave spectra, J Geophys Res: Solid Earth, 124(4): 3966-3986. DOI:10.1029/2018JB016957
Shearer P M, Prieto G A, Hauksson E, 2006, Comprehensive analysis of earthquake source spectra in southern California, J Geophys Res: Solid Earth, 111(B6): B06303.
Song S G, Lee H, 2019, Static slip model of the 2017 MW5.4 Pohang, South Korea, earthquake constrained by the InSAR data, Seismol Res Lett, 90(1): 140-148. DOI:10.1785/0220180156
Stein S, Wysession M, 2009, An Introduction to Seismology, Earthquakes, and Earth Structure, New York: John Wiley & Sons.
Sumy D F, Neighbors C J, Cochran E S, et al, 2017, Low stress drops observed for aftershocks of the 2011 MW5.7 Prague, Oklahoma, earthquake, J Geophys Res: Solid Earth, 122(5): 3813-3834. DOI:10.1002/2016JB013153
Tomic J, Abercrombie R E, Do Nascimento A F, 2009, Source parameters and rupture velocity of small M≤2.1 reservoir induced earthquakes, Geophys J Int, 179(2): 1013-1023. DOI:10.1111/j.1365-246X.2009.04233.x
Trugman D T, Shearer P M, 2017, Application of an improved spectral decomposition method to examine earthquake source scaling in Southern California, J Geophys Res: Solid Earth, 122(4): 2890-2910. DOI:10.1002/2017JB013971
Venkataraman A, Beroza G C, Ide S, et al, 2006, Measurements of spectral similarity for microearthquakes in western Nagano, Japan, J Geophys Res, 111(B3): B03303.
Viegas G, Hutchings L, 2011, Characterization of induced seismicity near an injection well at the northwest Geysers geothermal field, California, Trans Geotherm Resour Council, 35(2): 1773-1780.
Wang B, Harrington R M, Liu Y J, et al, 2020, A study on the largest hydraulic-fracturing-induced earthquake in Canada: Observations and static stress-drop estimation, Bull Seismol Soc Am, 110(5): 2283-2294. DOI:10.1785/0120190261
Wang R J, Gu Y J, Schultz R, et al, 2016, Source analysis of a potential hydraulic-fracturing-induced earthquake near Fox Creek, Alberta, Geophys Res Lett, 43(2): 564-573. DOI:10.1002/2015GL066917
Warren L M, Shearer P M, 2002, Mapping lateral variations in upper mantle attenuation by stacking P and PP spectra, J Geophys Res: Solid Earth, 107(B12): 2342.
Wu Q M, Chapman M, Chen X W, 2018, Stress-drop variations of induced earthquakes in Oklahoma, Bull Seismol Soc Am, 108(3A): 1107-1123. DOI:10.1785/0120170335
Wu Q M, Chen X W, Abercrombie R E, 2019, Source complexity of the 2015 MW4.0 Guthrie, Oklahoma earthquake, Geophys Res Lett, 46(9): 4674-4684. DOI:10.1029/2019GL082690
Yamada T, Mori J J, Ide S, et al, 2007, Stress drops and radiated seismic energies of microearthquakes in a South African gold mine, J Geophys Res: Solid Earth, 112(B3): B03305.
Yenier E, Atkinson G M, 2015, An equivalent point-source model for stochastic simulation of earthquake ground motions in California, Bull Seismol Soc Am, 105(3): 1435-1455. DOI:10.1785/0120140254
Yoshimitsu N, Ellsworth W L, Beroza G C, 2019, Robust stress drop estimates of potentially induced earthquakes in Oklahoma: Evaluation of empirical Green's function, J Geophys Res: Solid Earth, 124(6): 5854-5866. DOI:10.1029/2019JB017483
Yu H Y, Harrington R M, Kao H, et al, 2020, Well proximity governing stress drop variation and seismic attenuation associated with hydraulic fracturing induced earthquakes, J Geophys Res: Solid Earth, 125(9): e2020JB020103.
Zhang H L, Eaton D W, Li G, et al, 2016, Discriminating induced seismicity from natural earthquakes using moment tensors and source spectra, J Geophys Res: Solid Earth, 121(2): 972-993. DOI:10.1002/2015JB012603