Anti-seismic Analysis of Gas Pipeline Crossing through the Kezil Thrustfault Fault
-
摘要: 输气管道作为1种薄壁壳体结构,逆冲断层引起的管道压缩变形容易使其破坏。本文以大北南疆输气管道工程为例,探讨了穿越克孜尔逆冲断层的输气管道地震安全问题。在确定管道穿越处的断层倾角、设防断层位错量、表征管土相互作用的土弹簧参数以及钢管容许应变等参数后,采用壳有限元方法,分析了穿越克孜尔逆冲断层的输气管道变形反应。分析结果显示,管道在逆冲断层作用下以压缩应变为主,管道内的最大轴向压缩应变的幅值随着交角的减小而减小。在通过探槽等方法确定断层活动位置后,该管道若以小于或等于11°的交角通过克孜尔断裂,断层引起的最大轴向压缩应变和拉伸应变均在管道相应的容许应变范围内,满足相关规范的抗震要求。Abstract: As a thin-walled shell structure, the gas supply pipeline is easy to be damaged with large compression deformation caused by thrust fault movement. Taking the Dabei-South xinjiang pipeline project as an example, this paper discusses the seismic safety of gas pipeline crossing the Kezil thrust fault. The shell finite element modal was used to analyze the large deformation reaction of the gas pipeline under this thrust fault after determining the fault plane inclination angle, the fortified fault displacement, the soil spring parameters representing the pipe-soil interaction and the allowable strain of the steel pipe. The FEM results show that the main strain response of the pipe under thrust fault movement is the compression strain, and the amplitude of the maximum axial compression strain in the pipeline decreases with the decrease of the intersection angle.We found that in order to meet the requirements by the national seismic code of oil/gas pipeline, the maximum axial compression strain and tensile strain caused by the fault have to be within the permissble strain range of this pipe when the crossing angle is less than or equal to 11°.
-
Key words:
- Gas pipeline /
- Thrust fault /
- Shell finite element model /
- Seismic analysis
-
引言
随着中国城镇化建设的发展, 地震观测技术由地表观测向井下观测发展, 井下观测技术将传感器密封在筒内, 在几十米至几百米的钻孔中耦合安装, 可较好地屏蔽环境(降雨、雷电等)及人类活动的干扰, 有利于获得高精度的观测数据, 可用于记录地倾斜固体潮汐、地形变积累过程等丰富的地壳形变信息(董云开等, 2014)。CBT型钻孔倾斜传感器为井下安装使用的测地倾斜仪器, 此外还可与RZB系列的钻孔应变仪多形式组合为综合钻孔观测系统, 获得多测项更丰富的信息, 为地震监测和地球物理研究提供基础数据。
随着CBT型钻孔倾斜传感器的不断推广应用, 台站连续观测数据给仪器研发人员提供了进一步了解和提高仪器性能的机会。CBT型钻孔倾斜传感器观测数据质量受仪器自身性能、台站地质条件、钻孔条件、安装耦合情况等诸多因素的综合影响, 这些影响因素均为控制钻孔倾斜传感器观测质量的必要条件, 缺一不可。
CBT型钻孔倾斜传感器灵敏度达0.0002角秒(欧阳祖熙等, 2009), 为使观测系统能以万分之几角秒量级稳定工作, 倾斜传感器装配后需较长时间的预老化放置, 以消除加工和安装过程中产生的残余应力, 保证传感器的自身稳定性, 保证CBT型钻孔倾斜传感器顺利安装并投入使用。由于摆式倾斜仪器的特殊悬挂结构, 观测数据在较长一段时间内会存在突跳台阶及尖峰(李明等, 2011)、漂移量大(王梅等, 2003)等问题。目前为获得高灵敏度和高稳定性, 依靠热处理和自然时效处理摆式传感器弹性部件的方法非常普遍, 这种说法存在一定局限性, 一方面在工业环保的高要求下, 能进行热处理的企业越来越少;另一方面, 自然时效的周期太长, 部分仪器弹性部件需自然放置几个月, 甚至长达几年。倾斜传感器加工稳定的长周期与快速发展的地震监测网络台站建设之间的矛盾日益凸显, 如何快速消除倾斜传感器, 特别是摆体悬挂弹性部件的残余应力, 使其能在加工装配后快速稳定投入使用是解决问题的关键。
针对CBT型钻孔倾斜传感器存在的上述实际情况, 在实验室开展了倾斜传感器整机预老化工艺研究和实践, 取得良好的应用效果。
1. CBT型钻孔倾斜传感器简介
1.1 工作原理
CBT型钻孔倾斜传感器采用电容式倾斜传感原理(吴立恒等, 2010), 传感器为刚性支架上竖直悬挂的重摆, 如图 1(a)所示, 固定极板A、固定极板B和重摆M构成了三端电容差动式倾斜传感器。重摆M受到地球重力作用, 始终保持铅垂方向。当地面向某一方向倾斜角度ΔΨ时, 摆支架随之产生ΔΨ的倾斜量, 而重摆M在重力作用下仍保持铅垂方向, 重摆与摆支架必然产生相对位移, 如图 1(b)所示, 地面倾斜ΔΨ使重摆与固定极板产生相对位移Δδ, 关系式为:
$$ \Delta \mathit{\Psi }=\Delta \delta / L \text { (弧度 }) $$ (1) 物体间的电容量C与构成电容元件的2个极板面积S、相互距离δ、极板间介电常数ε有关, 关系式为:
$$ C=\varepsilon S / \delta $$ (2) 结合式(1)、式(2)可知, 重摆位移量仅与地面倾斜角成正比关系, 因此将测角问题转化为测位移问题。电容量仅与重摆位移量成正比关系, 进而将测位移问题转化为测电容问题。因此建立被测物理量ΔΨ和测量电容的关系, 即为CBT型电容式倾斜传感器的工作原理。
1.2 残余应力分析
根据工作原理设计结构尺寸合理的倾斜传感器, 其中悬挂部件为倾斜传感器成功与否的关键, CBT型钻孔倾斜传感器选用0.1mm厚的铍青铜弹簧片悬挂重摆M, 弹簧片性能尤其重要。传感器各部件通过加工、装配成有机整体, 实现倾斜测量。传感器在加工、装配过程中不可避免地产生残余应力, 主要包括加工材料自身的残余应力;加工过程中产生的残余应力;装配过程中产生的残余应力。对于普通观测系统, 残余应力量级很小, 不影响观测系统的性能, 可忽略不计。但对于具有0.0002角秒灵敏度的测地倾斜仪器, 将在很长一段时间内存在残余应力, 传感器会不断释放和调整。而调整为几角秒量级, 在一般的实验室环境下被背景噪声淹没, 因而无法观测。将新组装的倾斜传感器组(2个)密封后置于河北怀来地震台山洞中进行观测, 观测数据如图 2所示, 由图 2可知, 固体潮汐形态不平滑, 突跳多, 稳定周期长, 甚至长时间无法稳定。
1.3 以往残余应力的解决方法
通过多年的实践, 初步认为CBT型钻孔倾斜传感器摆体悬挂部件即弹簧片是关系传感器快速稳定、质量优良的核心部件, 以往对其采用热处理的方式消除材料自身及加工过程中的残余应力。此外, 当传感器组装完成后, 通过数周到数个月的悬挂放置消除装配残余应力, 此方法在实践中得到检验, 有利于钻孔倾斜传感器现场安装应用的快速稳定, 如漳州地震台深井综合观测系统于2008年12月安装, 倾斜单元在一星期内测到了平滑的固体潮汐, 图 3所示为2009年2月的数据, 由图 3可知数据质量良好, 具备入网工作条件。此套传感器悬挂放置时间长达8个月。
2. CBT型钻孔倾斜传感器振动时效工艺及试验
2.1 振动时效定义及原理
《振动时效工艺参数选择及效果评定方法》(GB/T 25712—2010)(中华人民共和国国家质量监督检验检疫总局等, 2011)中定义振动时效为物理过程, 即利用严格受控的振动能量, 对金属工件进行处理, 以解决工件加工过程中和加工后出现的内部残余应力导致的尺寸变化及抗荷载能力变化问题, 对消除、减少或均化金属工件内的残余应力, 提高工件抗动静荷载、抗变形能力, 稳定尺寸精度有超卓功效。从微观方面分析, 振动时效可视为一种以循环荷载的形式施加于零件上的附加应力。当工件受到振动, 施加于零件上的交变应力与零件中的残余应力叠加, 当应力叠加结果达到一定数值后, 在应力集中最严重的部位因超过材料的屈服极限而发生塑性变形。塑性变形降低了该处残余应力峰值, 并强化了金属基体, 而后振动又在其他应力集中较严重的部位产生同样作用, 直至振动附加应力与残余应力叠加的代数和不能引起任何部位发生塑性变形为止, 此时, 振动消除、均化残余应力及强化金属的过程结束。
2.2 传感器振动时效系统构建
根据振动时效原理, 结合倾斜传感器工作原理, 构建实验室倾斜传感器整机振动时效试验系统, 如图 4所示。该系统主要包括倾斜传感器、二维倾斜调平台、水平振动台、位移传感器、振动台采集控制单元及PC控制终端。
将倾斜传感器安装在振动平台上, 通过倾斜调平台调平传感器, 使传感器悬挂弹簧片在受力均衡的条件下开展试验。与此同时, PC控制终端能实时采集传感器输出, 保证振动时效试验在传感器量程内受控进行, 进而保证试验过程中传感器的安全。通过PC机控制界面, 将振型、幅值、频率等参数输出给振动台采集控制单元, 使水平振动台按照设定的模式带动传感器一维振动, 在一定时间内达到释放传感器残余应力、实现预老化的目的。此外, 位移传感器全程监测水平振动台振动幅度, 并作为反馈信号, 可更好地保证试验参数的准确性及试验的安全性。
2.3 传感器振动时效工艺
振动时效工艺目前主要应用于黑色金属领域, 对铸铁件和焊接件中应力集中的部位进行应力消减与均化, 以达到工程要求(廖凯等, 2019), 具有较成熟的工艺可供借鉴。然而对于不同材质装配而成的复杂工件组(如倾斜传感器)无法找到相关的工艺参数作为指导。合理的振动时效参数设计尤为重要, 由于很多台站应用已证实井下安装后的倾斜传感器逐渐趋于稳定, 因此可通过模拟传感器在井下的工作状态, 使其在实验室条件下提前老化稳定, 作为选择振动时效参数的重要思路。
CBT型钻孔倾斜传感器已在实验室和现场应用中积累了丰富经验, 无论是传感器在实验室的自然时效, 还是在井下工作初期的自然时效, 均受地球固体潮汐的作用, 产生24小时2波峰2波谷的运动(视为2个周期的正弦波), 因此将振动运动形式选为正弦波型。
为获得高灵敏度, 倾斜传感器摆体与左右极板的间隙仅为0.2mm, 因此振幅可选0.15mm、0.1mm等, 根据传感器自振周期, 选择振动频率为10Hz、8Hz等。此外, 结合多台站自稳定周期, 按3个月计算需要振动时效确定老化试验时间, 在实验室开展参数组合, 优化组合得到最佳方案。
2.4 传感器振动时效预老化应用实例
选用16-4、16-5号钻孔倾斜传感器作为应用实例样本, 传感器装配时间为2016年6月29日, 根据实验室优化参数, 进行实验室振动时效预老化工作, 工艺参数如表 1所示。
表 1 倾斜传感器振动时效预老化工艺参数Table 1. Vibrating stress relief system process parameters of the tilt sensor传感器号 振动类型 振动幅值/mm 振动频率/Hz 振动时间/h 16-4号传感器 正弦波 0.1 10 4 16-5号传感器 正弦波 0.1 10 4 2016年7月31日, 16-4号钻孔倾斜传感器在新疆阿图什下井安装, 安装后的早期数据如图 5所示。
2016年8月3日, 16-5号钻孔倾斜传感器在新疆新源地震观测点下井安装, 安装后的早期数据如图 6所示。
由图 5、6可知, 仪器安装后传感器快速稳定, 很快便记录到固体潮汐, 形态光滑, 漂移量较小, 说明振动时效预老化起到了积极作用。
需指出的是, 钻孔倾斜传感器的稳定工作与传感器自身稳定、钻孔条件、台站地质条件、耦合可靠性等诸多因素密不可分, 多因素共同作用才能保证钻孔倾斜传感器稳定工作, 产出高质量的观测数据。而传感器自身稳定可靠性仅为其中一个影响因素, 是产出高质量观测数据的必要条件, 而非充要条件。
由于残余应力的检测只能针对某个零部件, 一般还会产生损伤, 因此传感器振动时效预老化前后难以开展整体检测工作, 在实践中的应用效果检验无疑是一种有效手段。CBT型钻孔倾斜传感器在阿图什地震台快速稳定并投入观测的实例, 充分说明实验室振动时效预老化工艺对倾斜传感器快速稳定起到了积极作用, 此预老化系统及工艺对地震监测领域的摆式传感器具有较大借鉴意义。
3. 结语
CBT型钻孔倾斜传感器采用振动时效预老化工艺, 消除传感器加工、装配过程中的残余应力, 满足传感器快速稳定应用的要求, 取得良好的实践应用效果, 为地震监测领域的摆式传感器提供了新的预老化思路。与此同此, 应认识到影响传感器稳定工作的因素较多, 需诸多因素的共同配合, 且预老化工艺的积极作用需在更多的台站应用中加以检验, 该工艺也需进一步探索和完善, 才能更好地服务于地震监测事业。
-
表 1 3个方向土弹簧参数
Table 1. Three-direction soil spring parameters
土弹簧参数 管轴方向 水平横向 垂直方向(向上) 垂直方向(向下) 最大作用力/N·m-1 fs=1.1×104 Pu=8.8×105 qu=4.1×104 qul=2.5×105 屈服位移/m Zu=0.004 Xu=0.058 Yu=0.018 Yul=0.051 表 2 管道以不同的交角穿越逆冲断层的分析结果
Table 2. Analysis result of pipeline crossing thrust fault with different crossing angles
工况 交角/° 最大轴向拉伸应变/% 容许拉伸应变/% 最大轴向压缩应变/% 容许压缩应变/% 1 0.02 0.432 1.29 -0.342 -0.75 2 2 0.421 1.29 -0.381 -0.75 3 6 0.399 1.29 -0.518 -0.75 4 10 0.388 1.29 -0.661 -0.75 5 11 0.386 1.29 -0.703 -0.75 6 12 0.384 1.29 -0.7502 -0.75 7 13 0.383 1.29 -0.7915 -0.75 8 16 0.382 1.29 -0.887 -0.75 9 20 0.383 1.29 -1.01 -0.75 10 25 0.384 1.29 -1.18 -0.75 11 30 0.383 1.29 -1.36 -0.75 12 60 0.357 1.29 -2.09 -0.75 13 90 0.340 1.29 -2.31 -0.75 -
常想德, 胡伟华, 2016.大北至南疆利民3号阀室输气管道工程场地地震安全性评价报告.乌鲁木齐: 新疆防御自然灾害研究所. 冯启民, 赵林, 2001.跨越断层埋地管道屈曲分析.地震工程与工程振动, 21(4):81-87. http://cdmd.cnki.com.cn/Article/CDMD-10615-1016098645.htm 郭恩栋, 冯启民, 1999.跨断层埋地钢管道抗震计算方法研究.地震工程与工程振动, 19(4):43-47. doi: 10.3969/j.issn.1000-1301.1999.04.007 国家发展和改革委员会, 2004.SY/T 0450-2004输油(气)钢质管道抗震设计规范.北京:石油工业出版社. 侯忠良, 1990.地下管线抗震.北京:学术书刊出版社. 金浏, 李鸿晶, 2010.穿越逆冲断层的埋地管道非线性反应分析.防灾减灾工程学报, 30(2):130-134. http://d.old.wanfangdata.com.cn/Periodical/dzxk201002003 刘爱文, 张素灵, 胡聿贤等, 2002.地震断层作用下埋地管线的反应分析.地震工程与工程振动, 22(2):22-27. doi: 10.3969/j.issn.1000-1301.2002.02.004 刘爱文, 胡聿贤, 赵凤新等, 2004.地震断层作用下埋地管线壳有限元分析的等效边界方法.地震学报, 26(增刊):141-147. 刘爱文, 胡聿贤, 李小军等, 2005.大口径埋地钢管在地震断层作用下破坏模式的研究.工程力学, 22(3):82-87. doi: 10.3969/j.issn.1000-4750.2005.03.016 刘学杰, 孙绍平, 2005.地下管道穿越断层的应变设计方法.特种结构, 22(2):81-85. doi: 10.3969/j.issn.1001-3598.2005.02.028 史航, 王丽, 栾鲁滨, 2009.兰成渝管道的抗震设计.油气储运, 28(10):57-59. http://d.old.wanfangdata.com.cn/Periodical/yqcy200910016 田勤俭, 丁国瑜, 郝平, 2006.南天山及塔里木北缘构造带西段地震构造研究.地震地质, 28(2):213-223. doi: 10.3969/j.issn.0253-4967.2006.02.005 徐锡伟, 闻学泽, 叶建青等, 2008.汶川MS 8.0地震地表破裂带及其发震构造.地震地质, 30(3):597-629. http://d.old.wanfangdata.com.cn/Periodical/dzdz200803003 赵雷, 唐晖, 彭小波等, 2010.埋地钢管在逆断层作用下失效模式研究.应用基础与工程科学学报, 18(增刊):111-118. http://www.cnki.com.cn/Article/CJFDTOTAL-YJGX2010S1018.htm 中华人民共和国住房和城乡建设部, 2009.GB 50470-2008油气输送管道线路工程抗震技术规范.北京:中国计划出版社. Kennedy R. P., Chow A. W., Williamson R. A., 1977. Fault movement effects on buried oil pipeline. Transportation Engineering Journal of ASCE, 103(5):617-633. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=Doaj000004190648 Newmark N. M., Hall W. J., 1975. Pipeline design to resist large fault displacement. Proceedings of U.S. National Conference on Earthquake Engineering. Ann Arbor MI, 416-425. http://cn.bing.com/academic/profile?id=cd01ed1004275d92058805477f7c872c&encoded=0&v=paper_preview&mkt=zh-cn Takada S., Liang J. W., Li T. Y., 1998. Shell model response of buried pipelines to large fault movements. Journal of Structural Engineering, JSCE, 44A:1637-1646. http://cn.bing.com/academic/profile?id=4968f02614fabf0188ca337191357d9f&encoded=0&v=paper_preview&mkt=zh-cn Wang L. R. L., Yeh Y. H., 1985. A refined seismic analysis and design of buried pipeline for fault movement. Earthquake Engineering & Structural Dynamics, 13(1):75-96. http://cn.bing.com/academic/profile?id=f8aa415db8fb566e712da4dee2b8922e&encoded=0&v=paper_preview&mkt=zh-cn 期刊类型引用(1)
1. 张涛,邹进贵. 基于万有引力的普适定位. 全球定位系统. 2023(06): 91-97 . 百度学术
其他类型引用(0)
-