岩芯力学实验的面应变耦合系数及其差值的阈值研究.pdf

第4 3卷 第9期2 0 2 3年9月大 地 测 量 与 地 球 动 力 学J o u r n a l o fG e o d e s ya n dG e o d y n a m i c sV o l.4 3N o.9S e p t.,2 0 2 3收稿日期:2 0 2 2-1 1-1 0项目来源:中国地震局震情跟踪定向工作任务(2 0 2 2 0 1 0 2 1 7);江苏省地震局青年基金(2 0 2 1 0 7)。第一作者简介:何斌,高级工程师,主要从事地球物理及台网形变观测研究,E-m a i l:8 5 3 4 9 7 5 6q q.c o m。D O I:1 0.1 4 0 7 5/j.j g g.2 0 2 3.0 9.0 1 0文章编号:1 6 7 1-5 9 4 2(2 0 2 3)0 9-0 9 3 1-0 5岩芯力学实验的面应变耦合系数及其差值的阈值研究何 斌1 李 锋1 田 韬1 江昊琳1 章 东11 江苏省地震局,南京市卫岗3号,2 1 0 0 1 4摘 要:通过单轴压缩实验获取江苏省体应变仪钻孔岩芯的弹性模量和泊松比,利用面应变理论,计算钻孔应变仪器的理论耦合系数和实测耦合系数,探讨两者的相关性。在此基础上,引入观测精度,研究仪器的观测精度、耦合系数和数据质量的关系,确定耦合系数差值的阈值。结果显示:1)钻孔应变仪的观测精度与理论耦合系数和实测耦合系数呈正相关,钻孔应变仪的理论耦合系数和实测耦合系数一致性越好,仪器的观测精度越高,数据质量越好。2)江苏省各台站面应变的理论耦合系数和实测耦合系数差值的阈值为0.1。差值小于0.1时,观测数据质量较高;差值在0.10.2之间时,井下耦合对观测精度有一定影响,但不显著;差值超过0.2时,井下耦合存在问题,仪器的观测精度较低,数据的观测质量较差。关键词:体应变仪;岩石力学;耦合系数;阈值中图分类号:P 3 1 5 文献标识码:A 钻孔体应变仪传感器的工作环境在地下,仪器与围岩的耦合系数是影响应变观测精度的重要因素。计算耦合系数的常用方法有两类,一类是基于弹性力学理论,建立一个合理的理论模型,诸如双环模型、两层介质模型和双衬套模型等,采用正演的方法计算耦合系数1-2。然而实际计算过程中,钻孔岩芯的弹性模量和泊松比等相关力学参数多是采用岩性的平均弹性模量,并不能真实地反映钻孔应变周围的岩芯特性。另一类方法是基于固体潮汐变形理论,采用钻孔应变的观测值实地标定对耦合系数进行修正3-4。然而,由于实测耦合系数与观测数据的固体潮汐形态呈正相关,很难独立了解钻孔的真实耦合情况。本文对江苏省不同地区的体应变钻孔岩芯进行系统性的单轴压缩实验,获取不同地区钻孔岩芯的弹性模量和泊松系数5-6,再采用不同方法计算钻孔的耦合系数,对仪器观测精度进行相关性分析,以期得到井孔耦合程度的定量化评定指标,进一步探讨耦合系数对观测精度影响的阈值。1 单轴压缩实验1.1 实验样品与制备实验样品为江苏省在运行的T J-2型体应变仪的钻孔岩芯,台站分布见图1。T J-2型体应变仪的传感器为一个充满了硅油的长圆形弹性筒。根据D B/T8.2-2 0 2 0地震台站建设规范 地形变台站第2部分:钻孔地倾斜和地应变台站 要求,体应变仪安装时需要进行钻孔工程的施工:松散层孔段采用正循环无芯钻进,岩石井段采用1 5 9mm金刚石取芯钻头取芯钻进。根据岩芯的致密性进一步确定准确的安装位置,利用水泥将传感器和四周岩石紧密耦合在一起。受到岩石的挤压或拉伸时,筒内的液体压力发生改变,通过液压的增大或缩小,即可获得岩石的应变状态。此外,仪器还可以记录到清晰的固体潮汐变化和地震面波形态。图1 台站分布F i g.1 D i s t r i b u t i o no f t h es t a t i o n s大 地 测 量 与 地 球 动 力 学2 0 2 3年9月取样的钻孔岩芯位于井下传感器安装段,其中徐州台、邳州台、连云港台和溧阳台的平均取芯深度6 0m(表1),占取样总数的6 0%;新沂台和六合台平均取芯深度8 0m;南通台的取芯深度为1 6 0m,取样深度最深。观测层岩性以砂岩为主,占总数5 0%,其次是灰岩。制备实验样品时,首先将7个台站的钻孔岩芯(直径为1 5 0mm)制作成2个圆柱体试样。端面在平面磨床上磨平、磨光,使其长度、外表面光洁度、两端面光滑度以及平行度都达到国际岩石力学学会(I S RM)规定的要求。采用WAW-D 1 0 0 0型电液伺服万能实验系统进行单轴压缩实验,采用荷载速率加载控制,加载速率为0.1k N/s,持续加载直至试样破坏。整个实验在中国地震局地质研究所地震动力学国家重点实验室完成。表1 仪器基础信息T a b.1 T h eb a s i c i n f o r m a t i o no f t h e i n s t r u m e n t s台站井下传感器深度/m观测层岩性试件编号试件尺寸/mm直径高度备注徐州6 0深灰色灰岩13 1.2 3 1 1 1.4 02 0 1 7-1 0故障23 1.2 3 1 0 8.3 52 0 1 9-1 2恢复观测新沂8 7红色砂岩33 1.5 0 1 1 0.8 343 1.5 0 1 1 3.1 5-邳州6 0青灰色石英砂岩53 1.2 35 3.2 763 1.2 3 1 0 1.6 9-连云港6 1混合片麻岩73 1.2 35 9.4 584 1.5 94 6.4 1-南通1 6 0灰色、灰黄色砂岩93 1.2 38 1.7 11 03 1.2 38 0.3 5-溧 阳6 0安山岩1 13 1.5 0 1 1 6.4 61 23 1.5 0 1 1 5.7 7-六合8 4灰岩1 33 1.2 36 1.2 12 0 1 6-0 6故障1 43 1.2 32 8.9 72 0 1 7-1 0恢复观测1.2 实验结果根据单轴压缩条件下岩石试件的应力-应变变化特性,全应力-应变曲线可分为裂纹闭合、线弹性变形、裂纹稳定扩展、裂纹非稳定扩展、峰后破坏5个阶段。本文以分布在茅东断裂带的台站(溧阳台)和郯庐断裂带江苏段的台站(新沂台)为例,分析试件的应力-应变变化特征。因选取的井下传感器安装段岩芯较为致密,在压应力作用下应变量变化微小。尽管2个台站的岩芯均在较短时间内完成了试样的压密阶段,但仍然可以看出新沂台岩芯致密性较溧阳台的差,在压密阶段岩芯的进程出现明显波动(图2、3)。实验进入到弹性阶段,因为岩石在压密阶段完成了微裂隙的闭合,因此该阶段线性表现均较好。在临近峰值荷载的时候,试样完全破坏,岩石强度不断下降,由图2 单轴压缩条件下溧阳台应力-应变曲线F i g.2 S t r e s s-s t r a i nc u r v e so fL i y a n gs t a t i o n i nu n i a x i a lc o m p r e s s i o nt e s t图3 单轴压缩条件下新沂台应力-应变曲线F i g.3 S t r e s s-s t r a i nc u r v e so fX i n y i s t a t i o n i nu n i a x i a lc o m p r e s s i o nt e s t于前期宏观破裂面形成,岩石破裂成几块。根据岩芯的应力-应变曲线的斜率可以得到弹性模量:E=a(1)式中,E为岩石弹性模量(单位G P a),为应力,a为应变。泊松比通过单轴拉伸实验获得:=ba(2)式中,为岩石泊松比,b为应力-应变曲线弹性段横向 应 变,a为 应 力-应 变 曲 线 弹 性 段 纵 向应变。通过单轴压缩实验,获得各试样的力学参数,见表2。为了研究单轴压缩条件下不同台站的岩芯试件物理力学参数特征,对实验过程中的相关数据(表1)进行统计分析可知,同一钻孔内2个岩石样品具有大致相同的矿物颗粒和结构,在宏观上具有一致性,测试曲线具有显著的线弹性特征,但岩石材料本身具有非均质性,局部结构不同,使得测试结果存在差异。由于岩石弹性模量受应力环境的影响是一个动态模量,本文取2个样品的测试平均值。239 第4 3卷第9期何 斌等:岩芯力学实验的面应变耦合系数及其差值的阈值研究表2 各台试件力学参数计算结果T a b.2 C a l c u l a t i o nr e s u l t so fm e c h a n i c a l p a r a m e t e r so f s a m p l e sa t e a c hs t a t i o n台站名称试件编号弹性模量/G P a平均弹性模量/G P a泊松比平均泊松比徐州11 7.322 3.12 0.2 00.2 5 840.3 4 800.3 0 32新沂34 8.643 4.54 1.5 50.2 3 210.2 0 250.2 1 73邳州52 2.163 3.12 7.6 00.2 1 190.1 9 410.2 0 30连云港73 6.784 0.93 8.8 00.2 2 430.1 5 520.1 8 96南通93 7.91 04 0.13 9.0 00.1 6 490.2 0 490.1 8 49溧阳1 12 8.61 24 7.73 8.1 50.2 3 440.2 6 500.2 4 97六合1 33 4.71 45 6.54 5.6 00.2 4 560.1 5 360.1 9 962 面应变耦合系数的计算使用三层复合双衬结构模型1-2计算面应变的理论耦合系数,使用理论固体潮标定法4计算面应变的实测耦合系数。2.1 三层复合双衬结构模型设钢筒、水泥、岩石的弹性模量和泊松比分别为E1和1、E2和2、E3和3。假定采样岩芯介质均为各向同性弹性体,忽略井口及井底对传感器的影响,钻孔的轴向应力3=0。面应变的表达式为:m i=2k(1+2)E3(3)k值的计算参考张凌空等7的双衬套法:k=E3x4E1(4)x4=r22r23(2+2)(3+3)r22(2-3)r21(1+2)+r22(1-2)+r23(3+3)r21(1-2)+r22(2+2)(5)式中,1=(1+1)E1,2=(1+2)/E2,3=(1+3)E3,1=(1-1)/E1,2=(1-2)/E2,3=(1-3)E3,r1、r2分别为钢筒内、外半径,r3为井孔半径。平面应力作用下,面应变观测理论耦合系数计算公式为:m=2k1-3(6)在面应变耦合系数的求解式(6)中,需要知道体应变仪地下传感器安装位置岩芯段的弹性模量和泊松比。根据T J-2体应变仪技术参数可知,E1=2 11 01 0Pa,1=0.3 0,E2=3.01 01 0P a,2=0.2 5,r1=4 2mm,r2=4 4.5mm和r3=6 5mm,E3和3根据单轴压缩实验获取。各台站面应变理论耦合系数m计算结果见表3。表3 各台站的面应变理论耦合系数T a b.3 T h e o r e t i c a l c o u p l i n gc o e f f i c i e n t o fs u r f a c e s t r a i na t e a c hs t a t i o n台站钻孔岩芯力学参数弹性模量/G P a泊松比理论耦合系数m徐州2 0.2 00.3 0 320.8 0 37新沂4 1.5 50.2 1 731.0 6 32邳州2 7.6 00.2 0 300.8 7 22连云港3 8.8 00.1 8 961.0 0 79南通3 9.0 00.1 8 491.0 0 59溧阳3 8.1 50.2 4 971.0 5 86六合4 5.6 00.1 9 961.0 8 452.2 理论固体潮标定法体应变的观测对象是“体积”的相对变化,地下传感器的长度比截面直径大很多,轴向线应变要比径向的小得多,并且观测层一般在6 0m左右,接近地表轴向应力分量3=0,体应变观测值S可以表述为4:S=A(x+y)=A(1+2)(7)式中,A为面应变实测耦合系数,x、y为钻孔体应变仪在水平方向2个分量的观测值。采用球状径向不均匀弹性地球模型3计算理论面应变,考虑太阳和月亮万有引力变化的影响,忽略海潮、地质构造等因素的影响。设和分别为地面上一点的余纬和东经,R为地球半径,g0为地球平均重力加速度,n阶起潮力在地球表面上一点p某一时刻产生的水平方向固体潮理论应变为:=1R g0hnWn+ln2Wn2 (8)=1R g0hnWn+lnc o tWn+lns i n22Wn2 (9)=1R g0s i n2Wn-c o tWn (1 0)式中,hn和ln为勒夫数,hn等于球状径向不均匀弹性地球模型表面在n阶起潮力作用下的垂直位移与平衡潮潮高的比值,ln等于水平位移与平衡潮水平位移的比值,Wn为月亮和地球在地球内部任意一点p的起潮力位,可以近似表示为:Wn=Wm2+Wm3+Ws2(1 1)式中,Wm2、Wm2和Ws2分别为月亮的二阶、三阶起潮力位和太阳的二阶起潮力位。本文使用邱泽华等4的理论应变固体潮程序,通过拟合观测值与理论值的小时差分曲线方法进行实地标定求解A,计算结果见表4。由表4可见,实测耦合系数A随时间变化,339大 地 测 量 与 地 球 动 力 学2 0 2 3年9月总体接近0.4 5;新沂台和邳州台实测耦合系数A明显偏小,邳州台还存在比较明显的波动,这2个台站的井下传感器耦合可能存在一定问题。表4 根据观测值计算的面应变实测耦合系数T a b.4 T h em e a s u r e dc o u p l i n gc o e f f i c i e n t o f s u r f a c es t r a i nc a l c u l a t e df r o mt h eo b s e r v e dv a l u e s时间实测耦合系数A徐州新沂邳州连云港南通溧阳六合2 0 1 5-1 00.3 9 0.2 6 0.1 20.4 70.4 3 0.4 8 0.6 92 0 1 5-1 10.4 6 0.2 2 0.1 30.4 60.4 6 0.4 8 0.6 62 0 1 5-1 20.4 0 0.2 5 0.1 10.4 30.4 5 0.4 9 0.6 62 0 1 6-0 10.4 0 0.2 7 0.1 70.4 70.4 3 0.4 6 0.6 42 0 1 6-0 20.4 1 0.3 1 0.1 20.4 60.4 5 0.4 5 0.6 12 0 1 6-0 30.3 9 0.3 1 0.1 50.4 60.4 2 0.4 8 0.6 9均值0.4 1 0.2 7 0.1 40.4 60.4 4 0.5 1 0.6 63 面应变耦合系数相关性与观测精度 理论耦合系数m和实测耦合系数A存在如下关系7:m=2A(1 2)采用式(1 2),将各个台站体应变仪的理论耦合系数和实测耦合系数代换后,不同方法计算的耦合系数的相关性见图4。可以看出,各个台站实测的面应变耦合系数普遍小于理论耦合系数。实测耦合系数和理论耦合系数符合比较好的有连云港台、南通台、溧阳台和徐州台,2种系数差值小于0.1。其中,徐州台的实测耦合系数和理论耦合系数最为接近,两者仅相差0.0 1 63。实测耦合系数和理论耦合系数符合较差的有新沂台、邳州台和六合台,2种系数差值大于0.1。其中,新沂台和邳州台的理论耦合系数大于实测耦合系数,六合台的实测耦合系数大于理论耦合系数。图4 各台站面应变耦合系数与观测精度F i g.4 S u r f a c es t r a i nc o u p l i n gc o e f f i c i e n t a n do b s e r v a t i o na c c u r a c yo f e a c hs t a t i o n为了进一步定量评价各台站应变仪面应变耦合系数,基于应变仪器可以记录到清晰的固体潮汐变化,采用固体潮汐变化M 2波的潮汐因子相对中误差的比值作为评价体应变观测数据质量的指标,比值越小,说明仪器的观测精度越高、在固体潮频段数据质量越高。采用V e n e d i k o w调和分析方法计算每个台站相对应时间段内观测精度的均值,结果见图4。可以看出,观测精度与面应变理论耦合系数和实测耦合系数呈正相关的关系。面应变理论耦合系数和实测耦合系数一致性越好,观测精度越高(相对误差比值越小)。其中,徐州台、溧阳台、连云港台和南通台面应变耦合系数理论值与实测值的差值较小,相对应台站的观测精度较高。新沂台、邳州台和六合台的观测精度分别为0.0 2 19、0.0 4 38和0.0 1 44,均超过0.0 1,相应的面应变耦合系数理论值和实测值的差异也较大。观测精度可以直观评价观测数据的质量,观测数据质量也反映了传感器与井孔的耦合情况。图5给出观测精度最高的徐州台和观测精度最低的邳州台的原始观测数据。可以看出,徐州台的观测曲线动态变化背景清晰,能清楚地记录到正常的固体潮汐变化,相应的观测精度较高;邳州台日常观测固体潮汐形态较差,面应变耦合系数理论值和实测值相差最大,相应的观测精度较差,可能这个台站的井下传感器存在一定的耦合问题,直接导致观测数据的精度不高。图5 徐州台和邳州台2 0 2 0-0 3-0 13 1原始观测曲线F i g.5 O r i g i n a l o b s e r v a t i o nc u r v ea tX u z h o ua n dP i z h o us t a t i o n s f r o m M a r c h1t o3 1,2 0 2 04 结 语1)井下传感器与围岩的耦合程度,是钻孔体应变观测效能的主要影响因素。体应变钻孔岩芯的弹性参数是非常重要的基础数据,通过单轴压缩实验获取相关的弹性模量和泊松比,结合观测模型、系统各组成部分的几何参数获得较为合理的面应变耦合系数。从计算结果来看,江苏区域内整体的面应变理论耦合系数在0.9左右(表3)。2)徐州台、连云港台、南通台和溧阳台的面应变耦合系数理论值和实测值的差值均小于0.1,说明这4个台站的钻孔安装条件和场地条件较好,井下传感器与围岩耦合理想。新沂台和邳州439 第4 3卷第9期何 斌等:岩芯力学实验的面应变耦合系数及其差值的阈值研究台的面应变耦合系数理论值和实测值的差值均在0.5左右,可能存在井下耦合问题。3)钻孔应变仪的观测精度与理论耦合系数和实测耦合系数呈正相关,钻孔应变仪器的理论耦合系数和实测耦合系数一致性越好,仪器的观测精度越高、数据质量越好。江苏省各台站面应变理论 耦合系数和 实测耦合系 数差值的阈 值为0.1,当理论耦合系数和实测耦合系数差值小于0.1,说明观测数据质量较高;差值在0.10.2之间,说明井下耦合对观测精度有一定影响,但是不显著;当差值超过0.2,说明井下耦合存在问题,仪器的观测精度较低、数据的观测质量较差。参考文献1 张凌空,牛安福.钻孔体应变与面应变观测井孔耦合系数的计算J.地震学报,2 0 1 5,3 7(1):8 0-8 8(Z h a n gL i n g k o n g,N i u A n f u.C a l c u l a t i o n o f B o r e h o l e C o u p l i n g C o e f f i c i e n tB a s e do nB o r e h o l eV o l u m eS t r a i na n dA r e aS t r a i nO b s e r v a-t i o nJ.A c t aS e i s m o l o g i c aS i n i c a,2 0 1 5,3 7(1):8 0-8 8)2 张凌空,牛安福.基于双衬套理论求解地壳岩石应变传递系数的探讨J.地球物理学进展,2 0 1 7,3 2(4):15 2 5-15 3 1(Z h a n gL i n g k o n g,N i uA n f u.D i s c u s s i o no nS o l u t i o no fC r u s t a lR o c kS t r a i nT r a n s f e rC o e f f i c i e n tB a s e do nD o u b l eB u s hT h e o r yJ.P r o g r e s s i nG e o p h y s i c s,2 0 1 7,3 2(4):15 2 5-15 3 1)3 吴庆鹏.重力学与固体潮M.北京:地震出版社,1 9 9 7(W uQ i n g p e n g.G r a v i t ya n dS o l i dT i d eM.B e i j i n g:S e i s m o l o g-i c a lP r e s s,1 9 9 7)4 邱泽华,石耀霖,欧阳祖熙.四分量钻孔应变观测的实地相对标定J.大地测量与地球动力学,2 0 0 5,2 5(1):1 1 8-1 2 2(Q i uZ e h u a,S h iY a o l i n,O u y a n gZ u x i.R e l a t i v eI n-S i t uC a l i b r a t i o no f 4-C o m p o n e n tB o r e h o l eS t r a i n O b s e r v a t i o nJ.J o u r n a lo fG e o d e s ya n dG e o d y n a m i c s,2 0 0 5,2 5(1):1 1 8-1 2 2)5 何斌,江昊琳,章东,等.基于单轴压缩实验的钻孔应变仪岩芯力学性质的研究J.地震学报,2 0 2 0,4 2(6):6 9 7-7 0 6(H eB i n,J i a n gH a o l i n,Z h a n gD o n g,e ta l.S t u d yo nt h eM e-c h a n i c a lP r o p e r t i e so fC o r ea b o u tB o r e h o l eS t r a i n M e t e r sB a s e do nU n i a x i a lC o m p r e s s i o nE x p e r i m e n tJ.A c t aS e i s-m o l o g i c aS i n i c a,2 0 2 0,4 2(6):6 9 7-7 0 6)6 何斌,田韬,王恒知.基于钻孔特性的体应变观测相关性研究J.地震工程学报,2 0 1 8,4 0(1):1 5 3-1 5 8(H eB i n,T i a nT a o,W a n g H e n g z h i.S t u d y o n C o r r e l a t i o n o f V o l u m eS t r a i nO b s e r v a t i o nB a s e do nB o r e h o l eF e a t u r e sJ.C h i n aE a r t h q u a k eE n g i n e e r i n gJ o u r n a l,2 0 1 8,4 0(1):1 5 3-1 5 8)7 张凌空,牛安福.钻孔体应变与面应变观测参数k的计算J.中国地震,2 0 1 3,2 9(3):3 3 5-3 4 6(Z h a n gL i n g k o n g,N i uA n f u.C a l c u l a t i o no fB o r e h o l eV o l u m e t r i cS t r a i na n dA r e aS t r a i nO b s e r v a t i o nP a r a m e t e rkJ.E a r t h q u a k eR e s e a r c hi nC h i n a,2 0 1 3,2 9(3):3 3 5-3 4 6)S t u d yo nt h eB o r e h o l eC o u p l i n gC o e f f i c i e n t sa n dt h eT h r e s h o l do fT h e i rD i f f e r e n c eB a s eo nC o r eM e c h a n i c a lT e s t sHEB i n1 L IF e n g1 T I ANT a o1 J I ANGH a o L i n1 ZHANGD o n g11 J i a n g s uE a r t h q u a k eA g e n c y,3W e i g a n g,N a n j i n g2 1 0 0 1 4,C h i n aA b s t r a c t:T h r o u g hu n i a x i a l c o m p r e s s i o n t e s t,w eo b t a i n t h e e l a s t i cm o d u l u s a n dP o i s s o ns r a t i oo f t h ed r i l l c o r eo f t h eb o r e h o l e s t r a i n m e t e r s i nJ i a n g s up r o v i n c e.W e c a l c u l a t e t h e t h e o r e t i c a l c o u p l i n gc o e f f i-c i e n t a n dt h em e a s u r e dc o u p l i n gc o e f f i c i e n to f t h eb o r e h o l es t r a i n m e t e r su s i n gs u r f a c es t r a i nt h e o r y,a n dd i s c u s s t h e i r c o r r e l a t i o n.O n t h i sb a s i s,w e i n t r o d u c e t h eo b s e r v a t i o na c c u r a c y,s t u d y t h e r e l a t i o n-s h i pb e t w e e n t h eo b s e r v a t i o na c c u r a c yo f t h e i n s t r u m e n t,t h e c o u p l i n gc o e f f i c i e n t a n d t h ed a t aq u a l i t y,a n dd e t e r m i n et h et h r e s h o l do ft h ed i f f e r e n c eo ft h ec o u p l i n gc o e f f i c i e n t.T h er e s u l t ss h o wt h a t:1)T h eo b s e r v a t i o na c c u r a c yo ft h eb o r e h o l es t r a i n m e t e ri sp o s i t i v e l yc o r r e l a t e dw i t ht h et h e o r e t i c a lc o u p l i n gc o e f f i c i e n ta n dt h em e a s u r e dc o u p l i n gc o e f f i c i e n t.T h eb e t t e rt h ec o n s i s t e n c yb e t w e e nt h et h e o r e t i c a l c o u p l i n gc o e f f i c i e n t a n d t h em e a s u r e dc o u p l i n gc o e f f i c i e n t,t h eh i g h e r t h e o b s e r v a t i o na c c u r a-c yo f t h eb o r e h o l es t r a i n m e t e r sa n dt h eb e t t e r t h ed a t aq u a l i t y.2)T h et h r e s h o l do fd i f f e r e n c eb e t w e e nt h et h e o r e t i c a l c o u p l i n gc o e f f i c i e n ta n dt h em e a s u r e dc o u p l i n gc o e f f i c i e n to f t h es u r f a c es t r a i na te a c hs t a t i o ni nJ i a n g s up r o v i n c e i s 0.1.W h e n t h ed i f f e r e n c e i s l e s s t h a n0.1,t h eo b s e r v a t i o nd a t aq u a l i t y i sh i g h e r.W h e nt h ed i f f e r e n c e i s b e t w e e n 0.1a n d 0.2,t h e d o w n h o l e c o u p l i n gh a s a c e r t a i n,b u t n o t s i g n i f i c a n t,e f f e c t o n t h eo b s e r v a t i o na c c u r a c y.Wh e nt h ed i f f e r e n c ee x c e e d s0.2,d o w n h o l ec o u p l i n gp r o b l e m so c c u r,r e s u l t i n gi nl o wi n s t r u m e n to b s e r v a t i o na c c u r a c ya n dp o o rd a t ao b s e r v a t i o nq u a l i t y.K e yw o r d s:b o r e h o l es t r a i n m e t e r;r o c km e c h a n i c s;c o u p l i n gc o e f f i c i e n t;t h r e s h o l dF o u n d a t i o ns u p p o r t:T h eE a r t h q u a k eT r a c k i n gT a s ko fC E A,N o.2 0 2 2 0 1 0 2 1 7;Y o u t hF u n do f J i a n g s uE a r t h q u a k eA g e n c y,N o.2 0 2 1 0 7.A b o u t t h e f i r s t a u t h o r:HEB i n,s e n i o r e n g i n e e r,m a j o r s i ng e o p h y s i c s a n dn e t w o r dd e f o r m a t i o no b s e r v a t i o n s,E-m a i l:8 5 3 4 9 7 5 6q q.c o m.539