• ISSN 1673-5722
  • CN 11-5429/P

基于数值模拟与含水层垂向应力反演的静乐井水位异常分析

张淑亮 李宏伟 吕芳 陈慧

崔光耀, 伍修刚, 王明年, 荆鸿飞. 黏滑断层隧道减错措施参数对减错效果的影响分析[J]. 震灾防御技术, 2018, 13(3): 502-511. doi: 10.11899/zzfy20180302
引用本文: 张淑亮, 李宏伟, 吕芳, 陈慧. 基于数值模拟与含水层垂向应力反演的静乐井水位异常分析[J]. 震灾防御技术, 2019, 14(4): 854-868. doi: 10.11899/zzfy20190417
Cui Guangyao, Wu Xiugang, Wang Mingnian, Jing Hongfei. Analysis of Influence of Parameters of Reducing Measures on the Effect of Reducing Dislocation on Tunnels Crossing Stick-slip Faults[J]. Technology for Earthquake Disaster Prevention, 2018, 13(3): 502-511. doi: 10.11899/zzfy20180302
Citation: Zhang Shuliang, Li Hongwei, Lv Fang, Chen Hui. Analysis of the Jingle Well Water-level Anomaly Based on Numerical Simulation and Aquifer Stress Inversion[J]. Technology for Earthquake Disaster Prevention, 2019, 14(4): 854-868. doi: 10.11899/zzfy20190417

基于数值模拟与含水层垂向应力反演的静乐井水位异常分析

doi: 10.11899/zzfy20190417
基金项目: 

山西省地震局科技攻关项目 201903D321013

山西省地震局科技攻关项目 SBK-1936;SBK-1936-3

中国地震局震情跟踪项目 2018020304

详细信息
    作者简介:

    张淑亮, 女, 生于1963年。正研级高级工程师。主要从事地震综合预测与研究工作。E-mail:zzsl_009@163.com

Analysis of the Jingle Well Water-level Anomaly Based on Numerical Simulation and Aquifer Stress Inversion

  • 摘要: 利用数值模拟与含水层垂向应力反演的方法,对降雨与河流荷载作用下静乐井区应力进行定量计算,并从区域应力场与同一构造区其他前兆测项变化特征对2种方法计算结果进行分析。数值模拟结果显示,降雨与河流荷载的共同作用使静乐井所在区域产生的最大垂向位移约为0.7mm,最大垂向应力约为30hPa;同时段静乐井-含水层垂向应力反演结果显示静乐井-含水层最大垂向应力约为79hPa;数值模拟与静乐井-含水层垂向应力反演结果均表明,静乐井水位高值异常与同时段降雨量增多、河流荷载效应增强关系密切,但2种结果存在一定差异。通过可反映井区应力场变化特征的山西地震带3级地震缺震和b值、穿过静乐井的GPS基线、地震矩释放、同一构造单元同时段其他前兆测项与静乐井水位高值异常对比等分析结果可知,2017年静乐井水位高值异常期间井区应力场存在显著异常,推测静乐井水位高值异常除受降雨与河流荷载作用的影响外,也可能受构造活动增强的影响。
  • 随着我国交通基础建设的持续深入发展,穿越黏滑断层的隧道不断涌现,如在建的雅康高速二郎山隧道穿越保凰活动断裂,成兰铁路诸隧道穿越龙门山活动断裂等。在2008年汶川地震中,都映高速部分隧道,其所穿越的黏滑断层在强震诱发下发生黏滑错动,造成错动段隧道结构和围岩严重破坏,甚至出现了隧道整体垮塌的严重震害(图 1),这对灾后救援和灾后重建工作造成了严重影响(四川省交通厅公路规划勘察设计研究院,2008于海英等,2008)。如何提高断层黏滑错动时隧道结构的安全性和稳定性是亟待研究和解决的关键技术问题之一。

    图 1  龙溪隧道黏滑断层段隧道整体垮塌
    Figure 1.  Collapse of stick-slip fracture section in Longxi tunnel

    目前,黏滑断层隧道的减错措施主要有二衬设置减错缝和初支与二衬之间设置减错层2种。设置减错缝可减小隧道结构的纵向刚度,提高了隧道结构适应断层黏滑错动所产生相对变位的性能;设置减错层可减小隧道结构的横向刚度,消减了传递至二衬结构的强制位移。

    国内外专家、学者对穿越黏滑断层隧道的减错技术做了部分研究,主要有:穿越活动断层隧道震害特征及震害机理(高孟潭等,2008崔光耀等,2013a何川等,2014);断层错动下围岩与隧道结构的位移和应力变化(赵伯明等,2009李玉江等,2013);不同倾角、错距下断层黏滑错动对山岭隧道洞身结构的受力特性(熊炜等,2010蒋建平等,2011刘学增等,2014张理平等,2017);设置不同减错缝的减错效果(崔光耀等,2013b李学峰等,2014信春雷等,2015);设置不同减错层的减错效果(崔光耀等,2013c)等。以上研究在穿越黏滑断层隧道断层错动机制,设置减错缝不同间距、减错层不同厚度减错效果方面做了部分工作,但均未对穿越黏滑断层隧道设置减错缝不同缝宽、减错层不同刚度的减错效果进行研究。本文依托都汶高速友谊隧道F1黏滑断层段,对黏滑断层隧道设置不同缝宽减错缝、不同刚度减错层的减错效果进行研究,以期能够促进穿越黏滑断层隧道减错技术的发展。

    1.1.1   工程地质

    隧址区位于纸厂沟与小桃沟之间的山体内,受断层和岩性影响,节理裂隙较发育,岩体较破碎。

    F1黏滑断层为挤压性逆断层,隧道洞身穿F1黏滑断层而过,走向北东,倾向北西,倾角60°,破碎带宽度0.5—2.0m,断层破碎带主要由糜棱岩、碎块岩组成,Ⅴ级围岩。上下盘岩体均为Ⅴ级围岩。

    1.1.2   支护结构设计

    隧道支护结构为复合式衬砌,跨度×高度=9.4m×8.0m;隧道初支采用C25喷射混凝土,厚度25cm,二衬采用C25模筑混凝土,厚度40cm。

    1.2.1   计算模型

    以友谊隧道F1黏滑断层段为研究背景,利用FLAC3D有限差分数值模拟软件建立计算模型,隧道埋深为40m,隧底围岩厚40m,左右两侧宽度取5倍隧道跨度,约60m,纵向开挖长度100m。本次计算采用brick单元,其中围岩单元47550个,衬砌单元6820个。计算时围岩按弹塑性材料考虑,衬砌则按弹性材料考虑。减错缝采用12m设缝间距(崔光耀等,2013b),减错层厚度采用10cm(崔光耀等,2013c)。通过在上盘围岩施加位移荷载来模拟断层错动,位移荷载通过施加速度来实现,错距10cm。

    1.2.2   计算参数

    施工中减错缝主要采用止水带进行连接,减错层一般采用橡胶板(王帅帅等,2016),围岩根据勘测资料选取Ⅴ级围岩。计算模型物理力学参数如表 1所示。

    表 1  计算模型参数
    Table 1.  Parameters of calculation model
    参数 重度/kN·m-3 弹性模量/MPa 泊松比 黏聚力/MPa 内摩擦角/°
    围岩 20 15000 0.40 0.10 24
    C25喷射混凝土 22 23000 0.20
    C25模注混凝土 25 28000 0.20
    减错缝 6 10 0.45
    减错层 10 0—2000 0.45
    下载: 导出CSV 
    | 显示表格
    1.2.3   计算工况

    本文主要研究减错缝不同缝宽、减错层不同刚度对减错效果的影响规律,计算工况如表 2所示。

    表 2  计算工况
    Table 2.  Calculation condition
    工况 计算内容 备注
    1 无减错措施
    2 施设减错缝 缝宽分别为:5cm、10cm、15cm、20cm
    3 施设减错层 弹模分别为:2.5MPa、5MPa、10MPa、20MPa、50MPa、100MPa、200MPa、500MPa、1000MPa、2000MPa
    下载: 导出CSV 
    | 显示表格
    1.2.4   监测布置

    计算模型上盘设4个监测断面、下盘设3个监测断面,每个监测断面设置8个监测点,监测布置如图 2所示。

    图 2  监测布置
    Figure 2.  Arrangement of monitoring

    提取各工况错动完成时的各监测断面计算数据,以缝宽10cm为例进行内力分析。轴力、弯矩的计算按照公式(1)、(2)计算(中华人民共和国交通部,2004)。监测断面最大轴力值见表 3,监测断面最大弯矩值见表 4

    表 3  监测断面最大轴力(单位:kN)
    Table 3.  Maximum axial force of monitoring section (unit: kN)
    工况 距断层距离/m
    -22 -10 -2 2 10 22 34
    无减错措施 -6468 -9754 -14884 -13376 -13341 -9551 -3702
    缝宽10cm -5917 -9133 -14857 -13395 -12558 -8646 -3153
    注:距断层距离正值为上盘,负值为下盘,其它同理。
    下载: 导出CSV 
    | 显示表格
    表 4  监测断面最大弯矩(单位:kN·m)
    Table 4.  Maximum bending moment of monitoring section (unit: kN·m)
    工况 距断层距离/m
    -22 -10 -2 2 10 22 34
    无减错措施 -60.2 -142.1 -217.6 -185.0 -121.9 -45.1 -126.7
    缝宽10cm -61.3 -140.2 276.1 -200.5 -127.7 48.5 -116.6
    下载: 导出CSV 
    | 显示表格

    轴力、弯矩计算公式为:

    $$ N=\frac{1}{2}E({\varepsilon _内}+{\varepsilon _外})bh $$ (1)
    $$ M=\frac{1}{{12}}E({\varepsilon _内}-{\varepsilon _外})b{h^2} $$ (2)

    式中,b为截面宽度(m),取1m;h为截面厚度(m);E为弹性模量(MPa);${\varepsilon _内}$、${\varepsilon _外}$为结构内、外侧应变;N为结构轴力(N);M为弯矩(N·m);后同。

    表 3可知,施设减错缝与无减错措施工况监测断面最大轴力值均呈现压应力状态,且沿隧道纵向分布趋势相似,最大值均出现在下盘靠近断层的断面,其中无减错措施工况值为-14884kN,施设减错缝工况值为-14857kN;随着监测断面与断层的距离增大,监测断面最大轴力呈下降趋势。由各监测断面最大轴力整体分析可知,施设减错缝工况监测断面最大轴力略小于无减错措施工况。

    表 4可知,无减错措施工况监测断面最大弯矩均表现为外侧受拉,施设减错缝改变了隧道结构受力,出现内侧受拉部位;监测断面最大弯矩值均出现在下盘靠近断层的断面,其中不设减错缝工况值为-217.6kN·m,施设减错缝工况值为276.1kN·m。由各监测断面最大弯矩整体分析可知,施设减错缝工况监测断面最大弯矩略大于无减错措施工况。

    轴力、弯矩指标主要用于结构设计,单一分析轴力、弯矩不能综合体现结构的力学性能。隧道支护的结构安全性可由安全系数(综合考虑了轴力、弯矩等指标)体现。

    提取各工况错动完成时的内力计算数据进行安全系数分析,安全系数K按照公式(3)、(4)计算(中华人民共和国交通部,2004)。以缝宽10cm进行安全系数分析,监测断面最小安全系数见表 5

    $$ KN \le \varphi \alpha{R_{\rm{a}}}bh $$ (3)
    $$ KN \le \varphi \frac{{1.75{R_1}bh}}{{6{e_0}/h - 1}} $$ (4)
    表 5  监测断面最小安全系数
    Table 5.  Minimum safety factor of monitoring section
    工况 距断层距离/m
    -22 -10 -2 2 10 22 34
    无减错措施 1.933 1.282 0.840 0.935 0.937 1.309 3.377
    缝宽10cm 2.112 1.326 0.841 0.933 0.995 1.446 3.964
    下载: 导出CSV 
    | 显示表格

    式中,Ra为混凝土抗压极限强度(MPa);Rl为混凝土抗拉极限强度(MPa);φ为构件纵向弯曲系数;α为轴向力偏心影响系数;e0为轴向力对截面的偏心距(m)。

    表 5可知,施设减错缝与无减错措施工况监测断面最小安全系数沿隧道纵向分布趋势相似,最小值均出现在下盘靠近断层监测断面,其中无减错措施工况安全系数最小值为0.840,施设减错缝工况安全系数最小值为0.841;随着监测断面与断层距离的增大,监测断面最小安全系数呈增大趋势,上盘安全系数整体较下盘小。由各监测断面最小安全系数整体分析可知,施设减错缝工况监测断面最小安全系数略大于无减错措施工况。

    根据安全系数分析各工况减错效果,减错效果可按式(5)进行计算,各工况减错效果如图 3所示。

    $$ \lambda= \frac{{{K_1}-{K_0}}}{{{K_0}}} \times 100\% $$ (5)
    图 3  减错效果
    Figure 3.  Effect of reducing dislocation

    式中:λ为减错效果;K0为无减错措施工况监测断面最小安全系数;K1为施设减错缝工况监测断面最小安全系数。

    图 3可知,各工况减错效果沿隧道纵向分布趋势相似,随着监测断面与断层的距离增大,上下盘减错效果整体呈上升趋势,且上盘减错效果较下盘整体更好,其中上盘减错效果最大为24.50%,下盘减错效果最大为9.26%。随着减错缝宽度的增加,隧道上盘减错效果变好,下盘缝宽10—15cm时减错效果最好,主要由于断层黏滑错动对上盘的影响远大于下盘。减错缝宽度的增加适应了断层黏滑错动造成的上盘强制位移,断层黏滑错动的能量得到一定程度的释放,从而提高了上盘隧道的结构安全性,可以设想,对于隧道上盘部分存在最优减错缝宽度;减错缝宽度的增加虽然适应了断层黏滑错动造成的下盘强制位移,但由于抵御断层黏滑错动对下盘影响所需减错缝缝宽为10—15cm,当缝宽增加到20cm时,下盘隧道整体刚度的下降对下盘隧道的结构安全性造成了不利影响,致使减错效果下降。

    提取各工况错动完成时各监测断面的计算数据,以减错层弹性模量为2.5MPa为例进行内力分析。监测断面最大轴力值见表 6,监测断面最大弯矩值见表 7

    表 6  监测断面最大轴力(单位:kN)
    Table 6.  Maximum axial force of monitoring section (unit: kN)
    工况 距断层距离/m
    -21 -11 -2 2 11 21 32
    无减错措施 -6515 -10663 -17049 -14062 -12114 -7572 -3777
    减错层弹模(2.5MPa) -5557 -7014 -8375 -8981 -8490 -6154 -3788
    下载: 导出CSV 
    | 显示表格
    表 7  监测断面最大弯矩(单位:kN·m)
    Table 7.  Maximum bending moment of monitoring section (unit: kN·m)
    工况 距断层距离/m
    -21 -11 -2 2 11 21 32
    无减错措施 -78.5 -155.9 -218 -187.2 -100.3 -46.8 -35.9
    减错层弹模(2.5MPa) -173.4 -269 -254.1 -217.1 -108.0 -59.3 -42.7
    下载: 导出CSV 
    | 显示表格

    表 6可知,施设减错层与无减错措施工况监测断面最大轴力值均呈现压应力状态,且沿隧道纵向分布趋势相似,最大值均出现在上下盘靠近断层的断面,其中无减错措施工况值为-17049kN,施设减错层工况值为-8981kN;随着监测断面与断层距离的增大,监测断面最大轴力呈下降趋势。由各监测断面最大轴力整体分析可知,施设减错层工况监测断面最大轴力小于无减错措施工况。

    表 7可知,施设减错层与无减错措施工况监测断面最大弯矩值均表现为外侧受拉;无减错措施工况监测断面最大弯矩值出现在下盘靠近断层的断面,其值为-218kN·m,施设减错层监测断面最大弯矩值出现在下盘距离断层11m处,其值为-269kN·m。由各监测断面最大弯矩整体分析可知,施设减错层工况监测断面最大弯矩大于无减错措施工况。

    提取各工况错动完成时各监测断面的计算数据,以减错层弹性模量为2.5MPa为例进行安全系数分析,监测断面最小安全系数见表 8

    表 8  监测断面最小安全系数
    Table 8.  Minimum safety factor of monitoring section
    工况 距断层距离/m
    -21 -11 -2 2 11 21 32
    无减错措施 1.919 1.172 0.733 0.889 1.032 1.651 3.310
    减错层弹模(2.5MPa) 2.455 1.798 1.505 1.583 1.592 2.031 3.399
    下载: 导出CSV 
    | 显示表格

    表 8可知,施设减错层与无减错措施工况监测断面最小安全系数沿隧道纵向分布趋势相似,最小值均出现在下盘靠近断层监测的断面,其中无减错措施工况最小安全系数最小值为0.733,施设减错层工况最小安全系数最小值为1.505;随着监测断面与断层距离的增大,监测断面最小安全系数呈增大趋势,上盘安全系数整体较下盘小。由各监测断面最小安全系数整体分析可知,施设减错层工况监测断面最小安全系数大于无减错措施工况。

    根据安全系数分析各工况减错效果,减错效果计算方法同式(5),各工况减错效果如图 4所示。

    图 4  减错效果
    Figure 4.  Effect of reducing dislocation

    图 4可知,各工况减错效果沿隧道纵向分布趋势相似,随着监测断面与断层距离的增大,上下盘减错效果整体呈下降趋势,且下盘减错效果较上盘整体略好,其中下盘减错效果最大为105.32%,上盘减错效果最大为78.07%;当减错层弹性模量增加到一定程度(约100MPa),减错效果趋于稳定。

    (1)断层黏滑错动对隧道上盘的影响远大于下盘,上盘部分是黏滑断层隧道减错设防设计的重点段落。

    (2)减错缝对上盘部分隧道结构的减错效果优于下盘,其中上盘减错效果最大为24.50%,下盘减错效果最大为9.26%;减错层对下盘部分隧道结构的减错效果略优于上盘,其中下盘减错效果最大为105.32%,上盘减错效果最大为78.07%。

    (3)随着减错缝宽度的增加,隧道上盘减错效果变好,下盘缝宽10—15cm减错效果最好,主要是由于断层黏滑错动对上盘的影响远大于下盘。减错缝宽度的增加适应了断层黏滑错动造成的上盘强制位移,断层黏滑错动的能量得到一定程度的释放,从而提高了上盘隧道的结构安全性,可以设想,对于隧道上盘部分存在最优减错缝宽度;减错缝宽度的增加虽然适应了断层黏滑错动造成的下盘强制位移,由于抵御断层黏滑错动对下盘影响所需减错缝缝宽为10—15cm,当缝宽增加到20cm时,下盘隧道整体刚度的下降对下盘隧道的结构安全性造成了不利影响,致使减错效果下降。

    (4)随着减错层弹性模量的增加,隧道上下盘减错效果降低;当减错层弹性模量增加到一定程度(约100MPa),减错效果趋于稳定。

  • 图  1  井区地质构造

    Figure  1.  Well area tectonic map

    图  2  静乐井水位与降雨量和东碾河流量的对比

    Figure  2.  Comparison of Jingle well level with rainfall and Dongnian river discharge

    图  3  单降雨荷载作用下静乐井垂向位移和垂向应力

    Figure  3.  Vertical displacement and vertical stress of Jingle well under single rainfall load

    图  4  单河流荷载作用下静乐井垂向位移和垂向应力

    Figure  4.  Vertical displacement and vertical stress of Jingle well under single river load

    图  5  河流与降雨荷载共同作用下静乐井垂向位移和垂向应力

    Figure  5.  Vertical displacement and vertical stress of Jingle well under combined action of river and rainfall loads

    图  6  东碾河流量、静乐井-含水层垂向应力与降雨量的对比

    Figure  6.  Flow rate of Dongnian river, vertical stress of Jingle well-aquifer and rainfall comparison chart

    图  7  静乐井水位月均值与山西地震带3级地震缺震和b值的对比

    Figure  7.  Comparison chart of monthly average water level of jingle well with magnitude 3 earthquake deficiency and b value in Shanxi seismic belt

    图  8  山西岢岚-山西太原GPS基线时序曲线

    Figure  8.  GPS baseline timing curve of Kelan, Shanxi-Taiyuan, Shanxi

    图  9  静乐井所在区域地震矩加速释放时空图

    Figure  9.  Accelerated release diagram of seismic moment in jingle well area

    图  10  静乐井水位与同一构造区其它前兆测项对比图

    Figure  10.  Comparison between Jingle well water level and other precursor measurements in the same structure area

    表  1  静乐井-含水层多种参数计算结果

    Table  1.   Calculation results of various parameters of Jingle well-aquifer

    时间 气压系数/(mm·hPa−1 潮汐效率/(mm·10−9 体积压缩系数/pa−1 孔隙度/% 固体体积压缩系数/pa−1 水位变化量/m 含水层垂向应力/hPa
    2016-01 4.618 0.191 69.285 1.806 5.117 0.013 1.4019
    2016-02 2.522 0.330 9.766 2.442 0.751 0.007 0.5703
    2016-03 3.267 0.308 167.63 1.008 6.112 −0.004 −0.8453
    2016-04 5.997 0.180 960.22 0.434 18.121 −0.014 −5.9814
    2016-05 4.328 0.342 530.25 0.516 10.970 0.005 1.8874
    2016-06 5.032 0.264 4926.8 0.053 10.884 −0.0009 −3.2224
    2016-07 5.742 0.316 221.25 0.657 6.260 −0.003 −0.8655
    2016-08
    9.897 0.238 592.66 0.203 5.651 −0.017 −15.5497
    2016-09 3.960 0.167 170.08 1.215 8.059 −0.006 −0.9238
    2016-10 4.642 0.148 90.599 1.785 6.621 −0.007 −0.8333
    2016-11 2.459 0.269 162.6 1.366 6.875 0.002 0.2584
    2016-12 3.767 0.147 1709.5 0.106 6.9188 0.008 14.8921
    2017-01 4.007 0.154 123.43 1.52 7.348 0.007 0.8916
    2017-02 2.861 0.214 254.43 1.089 9.405 0.001 0.1522
    2017-03 2.784 0.227 1454.1 0.071 3.449 0.001 1.3134
    2017-04 3.808 0.292 148.14 0.997 5.684 −0.007 −1.3365
    2017-05 6.093 0.206 600.47 0.494 12.952 0.011 4.2036
    2017-06 4.756 0.285 84.983 1.233 4.318 0.005 0.8134
    2017-07 6.225 0.249 114.4 1.032 5.173 −0.004 −0.7851
    2017-08 18.400 0.191 4233.5 0.018 36.619 −0.03 −323.3068
    2017-09 15.148 0.164 751.41 0.388 14.207 −0.015 −7.5280
    2017-10 6.703 0.269 3324.2 0.055 8.1277 −0.012 −42.1179
    2017-11 3.473 0.161 813.31 0.676 20.424 0.010 2.9897
    2017-12 0.246
    下载: 导出CSV
  • 车用太, 鱼金子, 刘五洲, 1999.华北北部地区3次强震前地下流体异常场及其形成与演化机理.中国地震, 15(2):139-150. http://www.cqvip.com/Main/Detail.aspx?id=3579726
    车用太, 鱼金子, 2004.地下流体典型异常的调查与研究.北京:气象出版社.
    丁风和, 哈媛媛, 王勇等, 2015.基于数字化水位的张渤带地区构造应力场时序特征分析.地震, 35(2):133-138. doi: 10.3969/j.issn.1000-3274.2015.02.014
    黄辅琼, 晏锐, 陈颙等, 2004.利用深井地下水位动态研究大华北地区现今构造应力场状态.地震, 24(1):112-118. http://d.old.wanfangdata.com.cn/Periodical/diz200401017
    李春洪, 陈益惠, 田竹君, 1990.井-含水层系统对固体潮的动态响应及其影响因素.中国地震, 6(2):37-45. http://www.cqvip.com/Main/Detail.aspx?id=410642
    刘雁冰, 裴顺平, 2017.汶川地震前后b值的时空变化及构造意义.地球物理学报, 60(6):2104-2112. http://www.cqvip.com/QK/94718X/20176/672463573.html
    刘峡, 马瑾, 占伟等, 2013.汶川地震前后山西断陷带的地壳运动.大地测量与地球动力学, 33(3):5-10. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dkxbydz201303002
    刘瑞春, 李自红, 赵文星等, 2014.汶川M8.0地震前后山西地震带水平形变场变化特征研究.地震工程学报, 36(3):634-638. http://d.old.wanfangdata.com.cn/Periodical/xbdzxb201403034
    孙小龙, 刘耀炜, 晏锐, 2011.利用水位资料反演华北地区构造应力场变化.地震, 31(2):42-49. http://d.old.wanfangdata.com.cn/Periodical/diz201102005
    王爱英, 张淑亮, 杨占山, 1998.山西静乐井水位高值异常及其与华北地震活动的关系.地震, 18(2):201-205. doi: 10.3321/j.issn:0253-3782.1998.02.012
    吴小平, 1990.b值物理机制的再探讨.西北地震学报, 12(3):1-13. http://d.old.wanfangdata.com.cn/Thesis/Y065065
    鱼金子, 车用太, 张大维等, 1994.降雨对汤坑群井水位动态影响的剖析.地震研究, 17(3):273-280. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400058071
    张淑亮, 刘巍, 王爱英, 1997.从地下流体附加力源看晋5-1井的水位变化.山西地震, (4):27-31. https://www.zhangqiaokeyan.com/academic-journal-cn_earthquake-research-shanxi_thesis/0201254234742.html
    张淑亮, 李冬梅, 范雪芳, 2005.井水位前驱波与气压、风、降雨及强震关系的分析.地震, 25(3):69-77. http://d.old.wanfangdata.com.cn/Periodical/diz200503010
    张昭栋, 张广诚, 1987.利用水位阶变资料反演震时应力的调整变化.地震研究, 10(6):693-702. http://www.cnki.com.cn/Article/CJFDTotal-DZYJ198706004.htm
    张昭栋, 郑金涵, 冯初刚, 1989.水井水位的气压效率和降水荷载效率之间的定量关系.地震, (6):38-44. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000000213971
    张昭栋, 郑金涵, 冯初刚等, 1990.深井水位的地表荷载效应及其与含水层物理力学参数关系初探.地震学刊, (1):27-32. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000000378638
    张昭栋, 郑金涵, 冯初刚等, 1993.根据井水位阶变反演含水层应力变化的三种方法的对比.华北地震科学, 11(1):39-44. https://www.ixueshu.com/document/946e317b94c487682ce096f8b8e9acd0318947a18e7f9386.html
    朱艾斓, 徐锡伟, 胡平等, 2005.首都圈地区b值随震源深度的变化:对地震成核的意义.科学通报, 50(8):788-792. doi: 10.3321/j.issn:0023-074X.2005.08.012
    朱艾斓, 解朝娣, 徐锡伟等, 2010.鄂尔多斯块体周缘地区近期地震活动性与汶川地震应力触发作用的关系.地学前缘, 17(5):206-214. http://d.old.wanfangdata.com.cn/Periodical/dxqy201005018
    Bredehoeft J. D., 1967.Response of well-aquifer systems to Earth tides. Journal of Geophysical Research(72):3075-3087. doi: 10.1029-JZ072i012p03075/
  • 加载中
图(10) / 表(1)
计量
  • 文章访问数:  99
  • HTML全文浏览量:  52
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-07-05
  • 刊出日期:  2019-12-01

目录

/

返回文章
返回