岩芯力学实验的面应变耦合系数及其差值的阈值研究.pdf
《岩芯力学实验的面应变耦合系数及其差值的阈值研究.pdf》由会员分享,可在线阅读,更多相关《岩芯力学实验的面应变耦合系数及其差值的阈值研究.pdf(5页珍藏版)》请在咨信网上搜索。
1、第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
2、.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)钻孔应变仪的观测精度与理论耦合系数和实测耦合系数呈正相
3、关,钻孔应变仪的理论耦合系数和实测耦合系数一致性越好,仪器的观测精度越高,数据质量越好。2)江苏省各台站面应变的理论耦合系数和实测耦合系数差值的阈值为0.1。差值小于0.1时,观测数据质量较高;差值在0.10.2之间时,井下耦合对观测精度有一定影响,但不显著;差值超过0.2时,井下耦合存在问题,仪器的观测精度较低,数据的观测质量较差。关键词:体应变仪;岩石力学;耦合系数;阈值中图分类号:P 3 1 5 文献标识码:A 钻孔体应变仪传感器的工作环境在地下,仪器与围岩的耦合系数是影响应变观测精度的重要因素。计算耦合系数的常用方法有两类,一类是基于弹性力学理论,建立一个合理的理论模型,诸如双环模型、
4、两层介质模型和双衬套模型等,采用正演的方法计算耦合系数1-2。然而实际计算过程中,钻孔岩芯的弹性模量和泊松比等相关力学参数多是采用岩性的平均弹性模量,并不能真实地反映钻孔应变周围的岩芯特性。另一类方法是基于固体潮汐变形理论,采用钻孔应变的观测值实地标定对耦合系数进行修正3-4。然而,由于实测耦合系数与观测数据的固体潮汐形态呈正相关,很难独立了解钻孔的真实耦合情况。本文对江苏省不同地区的体应变钻孔岩芯进行系统性的单轴压缩实验,获取不同地区钻孔岩芯的弹性模量和泊松系数5-6,再采用不同方法计算钻孔的耦合系数,对仪器观测精度进行相关性分析,以期得到井孔耦合程度的定量化评定指标,进一步探讨耦合系数对观
5、测精度影响的阈值。1 单轴压缩实验1.1 实验样品与制备实验样品为江苏省在运行的T J-2型体应变仪的钻孔岩芯,台站分布见图1。T J-2型体应变仪的传感器为一个充满了硅油的长圆形弹性筒。根据D B/T8.2-2 0 2 0地震台站建设规范 地形变台站第2部分:钻孔地倾斜和地应变台站 要求,体应变仪安装时需要进行钻孔工程的施工:松散层孔段采用正循环无芯钻进,岩石井段采用1 5 9mm金刚石取芯钻头取芯钻进。根据岩芯的致密性进一步确定准确的安装位置,利用水泥将传感器和四周岩石紧密耦合在一起。受到岩石的挤压或拉伸时,筒内的液体压力发生改变,通过液压的增大或缩小,即可获得岩石的应变状态。此外,仪器还
6、可以记录到清晰的固体潮汐变化和地震面波形态。图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个圆柱体试样。端面在平面磨床上磨平、磨光,使
7、其长度、外表面光洁度、两端面光滑度以及平行度都达到国际岩石力学学会(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
8、-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
9、 6-0 6故障1 43 1.2 32 8.9 72 0 1 7-1 0恢复观测1.2 实验结果根据单轴压缩条件下岩石试件的应力-应变变化特性,全应力-应变曲线可分为裂纹闭合、线弹性变形、裂纹稳定扩展、裂纹非稳定扩展、峰后破坏5个阶段。本文以分布在茅东断裂带的台站(溧阳台)和郯庐断裂带江苏段的台站(新沂台)为例,分析试件的应力-应变变化特征。因选取的井下传感器安装段岩芯较为致密,在压应力作用下应变量变化微小。尽管2个台站的岩芯均在较短时间内完成了试样的压密阶段,但仍然可以看出新沂台岩芯致密性较溧阳台的差,在压密阶段岩芯的进程出现明显波动(图2、3)。实验进入到弹性阶段,因为岩石在压密阶段完成了
10、微裂隙的闭合,因此该阶段线性表现均较好。在临近峰值荷载的时候,试样完全破坏,岩石强度不断下降,由图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
11、 e s s i o nt e s t于前期宏观破裂面形成,岩石破裂成几块。根据岩芯的应力-应变曲线的斜率可以得到弹性模量:E=a(1)式中,E为岩石弹性模量(单位G P a),为应力,a为应变。泊松比通过单轴拉伸实验获得:=ba(2)式中,为岩石泊松比,b为应力-应变曲线弹性段横向 应 变,a为 应 力-应 变 曲 线 弹 性 段 纵 向应变。通过单轴压缩实验,获得各试样的力学参数,见表2。为了研究单轴压缩条件下不同台站的岩芯试件物理力学参数特征,对实验过程中的相关数据(表1)进行统计分析可知,同一钻孔内2个岩石样品具有大致相同的矿物颗粒和结构,在宏观上具有一致性,测试曲线具有显著的线弹性特
12、征,但岩石材料本身具有非均质性,局部结构不同,使得测试结果存在差异。由于岩石弹性模量受应力环境的影响是一个动态模量,本文取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
13、 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 56
14、0.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
15、)式中,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计算结果见
16、表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.
17、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为地球平均重力加速度,
18、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的理论应
19、变固体潮程序,通过拟合观测值与理论值的小时差分曲线方法进行实地标定求解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 m
20、t 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
21、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。其中,徐州台的实测耦合系数和理论耦合系数
22、最为接近,两者仅相差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为了进一步定量评价各台站应变仪面应变耦合系数,基于应变仪器可以记录到清晰的固体潮汐
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 力学 实验 应变 耦合 系数 及其 差值 阈值 研究
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。