2022年门源M_(S)6.9地震的同震滑动分布:联合InSAR、GPS和地表位错的贝叶斯建模

2022年1月8日青海省门源M_(S)6.9地震发生于托莱山断裂东段和冷龙岭断裂西段的交汇部位,其运动学特征关系到青藏高原北缘的动力学行为,同时也是区域地震危险性评估中的基础数据.本研究针对该地震,综合InSAR和GPS观测的同震位移场,以及野外地质调查的同震位错数据,采用贝叶斯反演方法,构建同震滑动分布,断层模型采用三角网格更好地刻画断层的不规则几何形态.结果显示,本次左旋走滑型地震的主要滑动发生在冷龙岭断裂西段,深度范围约2~6 km,计算的地震矩为0.95×10^(19)N·m,对应矩震级M_(W)6.65.利用多数据的联合反演解析出,沿走向存在两个滑动集中区,东侧的最大滑动达约4.8 m.震间亏损能量和库仑应力分析表明,托莱山断裂和冷龙岭断裂仍存在未来发震潜能;该地震造成的应力变化,对长期地震空区金强河—毛毛山段落有加载作用,因此,该段未来地震危险性值得关注.

泸州地区深层页岩气水平井套变成因机理探讨

四川盆地南部泸州地区深层页岩气储量丰富,是继长宁、威远等中浅层页岩气之后重要的天然气资源接替领域。但在泸州地区深层页岩气开发过程中水平井套管变形(以下简称套变)频发,对水平井压裂施工产生明显影响,制约了该区深层页岩气的规模建产。业内普遍认为套管变形是地层中的断裂发生滑移并对水平井产生剪切和挤压作用引起的,但基于该理论的套变模型与诸多实测数据不相吻合。为了进一步明确套管变形机理,提高套变的防控效果,以泸203井区18口套管变形水平井为对象,基于多臂井径成像测井,逐个对比套变形态、套变层位以及平面断裂发育情况,建立了断裂滑移和层理面滑移2种地质模型,并分析了2种地质模型下的不同套管变形特征。研究结果表明:(1)断裂滑移诱发的套变具有以下规律,即平面上套变仅出现在水平井与断裂带的交会处;套管的变形长度理论上应在10 m以内;套管的变形程度应呈近似对称关系,中部变形强,向两侧逐渐减弱。(2)该区多数套变不在预测的断裂带内,且部分套变段变形长度与形态不符合断裂滑移的特征。(3)层理面滑移可诱发多个相似的套变,在该机制下套管变形长度等于井筒与滑移面相交的距离,且层理面滑移的垂向分布与地层的页理发育程度一致。(4)层理面滑移与断裂滑移具有相似的力学机理,通过三维莫尔圆可以分析钻井和压裂过程中诱发断裂/层理面滑移的关键参数变化。结论认为,除断裂滑移外,层理面滑移是诱发页岩气水平井套管变形的另一重要机理,该研究成果可为深层页岩气规模效益开发的套变防控提供理论指导。

穿越多断层破碎带隧道自适应结构走滑错动特性研究

在活动断裂构造带地区,断层运动往往是沿断层带发生,自活动盘到破碎带错动位移表现出逐级传递的特征,因此穿越活动断裂带的隧道结构应力、变形特征也受到断层错动形式的影响。依托西部高烈度区某穿越大型活动断裂带的公路隧道工程,通过数值模拟分析穿越多断层隧道结构在断层多级走滑错动下的位移、应力响应规律,得出隧道结构纵向基于多断层蠕滑错动的影响分区;建立基于柔性初支、刚性二衬和节段铰接的穿越蠕滑断层隧道自适应结构,探究在错动区和缓冲区内设置隧道自适应结构的抗错性能。研究表明:隧道围岩性质的差异是影响穿越断层破碎带隧道结构变形和应力的最重要因素,隧道自适应结构在蠕滑错动过程中,柔性初支能够吸收部分围岩传递至隧道结构的变形,柔性接头的剪切变形和衬砌节段的刚性转动能够承担较大的位错,随着隧道衬砌节段长度的减小、接头数量的增加,各个接头的剪切变形角和衬砌节段偏转角都将减小,大幅提升结构的抗错安全性。

一种多频多系统周跳探测与修复方法

为满足多频多系统精密单点定位(PPP)数据预处理的实际需求,提出了一种周跳探测与修复方法。针对多频多系统PPP数据处理中各卫星的信号数量不同,利用最小二乘模糊度降相关平差(LAMBDA)原理构建模糊度组合。通过仅使用相位观测值估计的倾斜电子总量(STEC)探测不敏感周跳和漏判周跳。对于周跳修复,将伪距和载波相位的观测值在历元间作差求得周跳的浮点解和协方差阵,然后应用LAMBDA算法搜索周跳整数解并进行修复。通过静态、动态以及磁暴环境下周跳探测和修复实验,验证了所提方法的有效性,多频条件下周跳修复的正确率达到100%。

塔里木盆地富满油田东部走滑断裂发育特征与油气成藏

富满油田东部走滑断裂是塔里木盆地奥陶系超深层碳酸盐岩油气藏勘探的重点目标。该走滑断裂活动较弱,识别解译难度较大,需要深入研究其发育、演化及控藏特征。基于富满油田最新采集的三维地震数据,选取并精细解译了3条典型主干断裂(F_(Ⅰ)10,F_(Ⅰ)12和F_(Ⅰ)17),明确了富满油田东部走滑断裂活动特征及主要活动期次。在此基础上,结合油气性质及成藏期次,分析了断裂活动特征与油气成藏的关系,加深了对该区域断裂发育特征与油气分布规律的认识,进一步明确了断裂演化与油气成藏的耦合关系对油气性质的影响。研究表明,富满油田内走滑断裂具有垂向分层差异变形特征,可划分出4个构造变形层,由下往上分别为盐下基底构造层(T∈_(2)之下)、盐构造层(T∈_(2)—T∈_(3))、碳酸盐岩构造层(T∈_(3)—TO∈_(3)t)和碎屑岩构造层(TO∈_(3)t之上)。研究区断裂活动期次整体可划分为加里东早期、加里东中期Ⅲ幕和加里东晚期—海西期3个阶段,其中加里东中期Ⅲ幕是区域内断裂主要活动阶段。结合油气性质差异、油井生产动态、断裂活动特征与油气成藏期的耦合关系认为,断裂活动影响了垂向运移通道的开启,控制了油气垂向运移;断裂持续活动时间长,使油气运移通道保持开启,有利于晚期高成熟油气的持续充注,造成断控油气藏具有油气充注强度大、成熟度高及气/油比高的特征。

多场耦合条件下充填料浆管输壁面滑移特性

为了探究多场耦合作用下高浓度充填料浆管道输送过程中的壁面滑移特性,基于流体动力学理论与非牛顿流体的基本原理对管内料浆流动过程进行分析,建立考虑温度场—化学场—料浆流场的多场耦合数值模型。并通过Comsol进行数值计算,分析初始温度、料浆浓度及管径对壁面滑移速度的影响规律。将模型计算结果与实验结果进行对比,误差小于10%,证明了模型的有效性与可靠性。结果表明:温度升高,加剧尾砂颗粒的布朗运动,促使“絮网”结构向着“液网”结构转换,浆体黏度变小,滑移速度增大;管径增大,浆体所受剪切作用力变小,滑移速度减小,当管径达到180 mm时浆体所受剪切作用力小于发生滑移的临界值,管道输送过程中不存在完整的滑移过程;料浆浓度越大,浆体发生滑移运动所需的临界剪切应力越大,同时料浆浓度较大时,其所形成的滑移层较厚,滑移速度偏大。

基于多传感器融合的拖拉机驱动轮滑转率估算方法

拖拉机驱动轮滑转率的精确估计对提高拖拉机作业效率及安全性,实现拖拉机驱动防滑控制具有重要意义。本文提出了多新息并行扩展卡尔曼滤波算法,融合了包括机器视觉在内的多个传感器信息,并通过在线统计多传感器的新息,引入D-S证据理论进行决策,修正测量噪声矩阵,从而实现拖拉机驱动轮滑转率的精确估计。仿真结果表明,相较于普通的卡尔曼滤波算法,本文提出的融合算法估计滑转率的精度更高,滑转率估计值的均方根误差从2.34%降低到1.45%,且对干扰信号不敏感。试验结果表明,在多组工况下本文所提出的多传感器信息融合算法相较于单一视觉法或雷达法对拖拉机驱动轮滑转率估计的平均绝对误差、均方根误差均有所降低,从而验证了所提出的算法满足拖拉机驱动轮滑转率的精确估计,可对后续实现拖拉机驱动防滑控制提供依据。

基于解耦计算的多步快速成形有限元法

钣金零件是构成飞机机体和汽车车身的重要零部件。为缩短钣金件的生产周期,将有限元法用于零件设计初期评估零件的可制造性。为提高现有板料成形有限元法的计算效率,同时保证模拟精度,在一步快速成形有限元法的理论基础上,研究多步快速成形有限元法,通过引入中间构形的方式来考虑板料成形过程中加载路径与变形历史的影响。中间构形的构造是多步快速成形有限元法的关键,采用解耦思想将中间构形分解为弯曲变形和拉伸变形两个独立的过程进行计算。在弯曲变形阶段,不考虑材料的流动,根据板料与模具之间的位置关系计算获得滑移约束面;在拉伸变形阶段,材料的流动限制在滑移约束面上,通过应力平衡迭代以及节点修正后获得中间构形。以典型的钣金零件为例进行成形模拟,与现有商业有限元软件在计算精度和效率上进行对比,验证了该算法的可行性和有效性。结果表明,所提出的算法能够快速地构造出合理的中间构形,且能够准确地预测零件的成形性和厚度分布。

基于VSG及内模-滑模的MMC控制策略研究

风力发电并网系统的模块化多电平换流器具有零阻尼惯性,负荷突变及网侧故障都可能导致系统频率发生较大偏差,易产生冲击电流及能量堆积。文章提出基于虚拟同步发电机的内模-滑模MMC控制策略。在虚拟同步发电机中增加基于内模控制的滑模控制器,使模块化多电平换流器具有同步发电机一次调频外特性,并提升系统鲁棒性。通过Matlab/Simulink仿真模块搭建模块化多电平换流器并网模型,验证控制策略的正确性与有效性。仿真结果表明,控制策略能在系统发生负荷突变及网侧故障时提供频率支撑,有效抑制系统电流、功率的波动,减小对系统的冲击,提升运行稳定性。

深水脐带缆结构阻尼的数值与试验研究

深水动态脐带缆内部结构复杂,其阻尼特性对脐带缆的舞动以及疲劳损伤有着重要的影响。脐带缆在和流体相互作用过程中,脐带缆内部构件可能会产生滑动,进而影响脐带缆的结构阻尼,这也使得真实海况中脐带缆的舞动和疲劳是一项复杂而值得深入研究的课题。为验证有限元模型能够很好地模拟脐带缆的结构阻尼特性,在考虑内部结构单元间摩擦的基础上,利用ANSYS建立了用于中国南海某油田的脐带缆多层滑移接触摩擦模型。同时,采用自由衰减法对该动态脐带缆开展了原型试验,测得了该缆的结构阻尼比与自振周期。通过对比有限元模型计算结果与原型试验结果发现:数值模型获得的该脐带缆的结构阻尼值与试验值非常接近,这为工程中获得脐带缆结构阻尼比的值提供了好的方向。

多点接触粘滑驱动微动平台复合运动模型构建

为了精确模拟基于多点接触粘滑驱动原理的微动平台动态微位移输出,验证粘滑驱动原理在微动平台应用的可行性,课题组基于4稳态电磁执行器的粘滑驱动微动平台,建立了一种复合动态模型。该动态模型包含4个子动态模型,即磁力计算模型、LuGre摩擦模型、碰撞反弹模型和质心偏移模型;课题组在此基础上搭建了实验平台,研究了不同驱动电流对微动平台动态输出位移影响。结果表明:复合动态模型模拟动态位移输出与实验输出趋势基本一致,验证了所提出和建立的复合动态模型的有效性。

基于局部强度折减法的边坡多滑面分析方法及应用研究

边坡最危险滑动面的搜索方法一直是研究的热点,但边坡内部次级滑动面也可能不满足安全设计要求,因此,考虑边坡多条滑动面的分析方法亦应得到关注。在传统强度折减法中,对边坡整个区域的抗剪强度参数进行折减,此方法仅可得到一个临界滑动面和最小安全系数。提出了一种基于局部强度折减法的多滑面分析方法,即首先定义单元安全系数的概念,并且计算边坡每个单元的安全系数,然后自动搜索出单元安全系数处于不同范围内的单元集合,对各个单元集合的强度参数进行折减计算,即可得到不同安全系数对应的滑动面。通过单台阶和双台阶边坡算例验证了该方法的可行性,结果表明,随着安全系数的增大,潜在滑面的深度和潜在滑动区域亦增大。最后把该方法应用到某隧道进口仰坡的稳定性评价中,通过该方法得到的多级滑动面与现场监测数据吻合较好。

基于岩土材料软化特性的滑坡多级滑动面分析方法研究

在实际工程中,滑坡经常呈现出多级破坏的特征,而一般计算和设计中仅关注最危险的滑动面及对应的最小安全系数,这往往遗留安全隐患。在考虑岩土材料软化特性的基础上,借助FLAC3D和Matlab软件平台构建了无需人工干预便可有效模拟与评价滑坡多级破坏的理论框架。以旬阳县党家坝廉租房小区滑坡为例,通过对塑性剪应变、塑性拉应变、剪切应变增量等特征变量的渐进发展规律分析发现:滑坡多级滑动面的形成,其时间和空间顺序并不一定相对应。各级滑动面在时间上按产生的先后顺序依次为第1级主滑面、第2级主滑面、次级滑面;而在空间上按照从前到后的顺序依次为第1级主滑面、次级滑面、第2级主滑面。从计算所得的滑坡最终破坏形态中发现,坡体内部出现的滑动面条数与现场采集的拉裂缝数基本保持一致;第1级主滑面的入口位置与滑坡前缘拉裂缝位置基本吻合;而第2级主滑面及次级滑面的位置与现场勘查到的拉裂缝位置出现偏差。从材料参数的渐进发展规律中发现,滑面上强度参数的分布及大小均随着滑面的产生、发展而逐渐变化,其变化区间为峰值强度到残余强度,这是此方法可以有效模拟多级滑动面形成过程与各级滑动面之间相互影响的核心所在。借助矢量和法成功地实现了基于滑面上强度参数渐进发展规律的安全系数演化过程的确定。通过对安全系数演化过程的分析发现,在滑坡多级滑动面的形成过程中安全系数大小顺序呈现出3种不同的状态,这很好地揭示了各级滑动面在形成过程中的主、被动关系。

非均质边坡多级稳定性分析方法

为了实现对非均质边坡稳定性的全面表征与系统评价,提出了一种针对非均质边坡的"多级次"评价方法。这种方法以强度折减理论为基础,首先对岩土体进行多次局部强度折减,得到边坡不同深度、不同岩层控制的多级潜在滑动面;然后运用极限平衡法进行各个潜在滑动面的稳定性计算,即可得到边坡不同深度、不同级次的稳定性系数。这种方法在确定多级滑动面的过程中充分利用了强度折减法客观性的优点,避开了极限平衡方法主观性的缺点,折减所得滑动面更能反映边坡在重力作用下的变形破坏过程,优于滑弧法搜索结果。通过规则模型的验证及实际边坡的应用,证明这种方法对于综合评价非均质边坡的稳定状态是有效的,并可为边坡的治理、监测设计及施工提供指导。

多滑面滑坡稳定性分析

大型滑坡往往具有多个滑动面,然而工程中却较少考虑各滑面之间的相互作用。基于塑性力学上限定理,提出了一种能够考虑滑面之间相互作用的滑坡稳定性分析方法。采用竖直条分法将上下滑体分为条块,根据条块间速度间断构造出机动许可的速度场,然后根据外力功率和内能耗散率相等建立平衡条件,通过迭代计算得到最小安全系数。最后,通过将作者文中的安全系数解与他人的计算结果对比证明了该解的可靠性。

基于动态强度折减DDA法的边坡多滑面稳定性分析

为解决岩体结构面损伤的程度不一性问题及边坡多滑面分析中数值模型的选取问题,借鉴DDA计算位移的优势,考虑块体间相对位移,提出动态强度折减DDA法(Dynamic Strength Reduction Discontinuous Deformation Analysis,简称DSR-DDA)的边坡多滑面搜索,通过对不同破坏程度下岩体结构面动态折减不同的强度及动态变化边坡多滑面分析中的数值模型,以解决数值分析中岩体结构面损伤的程度不一性问题及边坡多滑面分析中数值模型的选取问题。利用倾斜平面滑块经典案例测试位移阈值,验证方法的可行性及计算精度,并将其应用到抚顺西露天矿岩质高边坡多滑面稳定性分析中,且与自主研发的GeoSMA-3D系统进行耦合。研究结果表明:在设定位移阈值为1 mm时,DSR-DDA与理论解误差在0.5%以内,满足计算要求。后续得到了抚顺西露天矿边坡最危险滑落面(对应安全系数为1. 136),与分析现场监测位移得出的潜在滑落位置基本一致,并在首次滑坡后的数值边坡模型基础上,得到了次级滑落面位置(对应安全系数为1.189),与常规分析得出的结果大相径庭,佐证了边坡多滑面稳定性分析中数值模型需动态变化的必要性。DSRDDA法与GeoSMA-3D系统识别出的关键块体耦合效果良好,从块体理论角度印证了DSR-DDA法的实用性。抚顺西露天矿边坡破坏模式为牵引式滑动破坏,中下部油母页岩对边坡稳定性起关键作用,故不可继续开采剩余油母页岩,以免引起边坡上部大面积滑坡及矿坑-城市边界滑坡。

重力坝深层抗滑稳定多滑面计算方法探讨

重力坝坝基深层滑裂面出现多滑面或双滑面形式的滑移节点较大偏移坝址垂线,若直接按照规范以双滑面计算式列出各滑移体的算式求解安全度,存在2个问题:建基面上部水平荷载和铅直荷载在滑移体上的分配方式存在多样;计算公式复杂,不利于迭代。为解决上述问题,将建基面上部荷载分别按应力和重力分配,利用等安全度的思想推导,分别得到了简化的分项系数和单一安全系数多滑面计算式。工程实例证明:采用简化的多滑面计算式验算重力坝深层抗滑稳定安全度是可行的;上部荷载按应力分配的计算结果较重力分配法偏于安全;将多滑面简化为双滑面计算得到的安全度偏大。

晚白垩世以来沿阿尔金断裂带的阶段性走滑隆升

年代学和沉积学特征表明,与印度板块约55MaBP以来沿喜马拉雅东、西构造结缝合带碰撞后青藏高原的阶段性隆升相对应,沿阿尔金断裂带阿尔金山地区亦发生阶段性走滑隆升,时间段大致为55～30MaBP、约20MaBP、约10MaBP和1～0.8MaBP。始新世—渐新世期间,阿尔金断裂大型左行走滑活动的发生,对中国西北大陆构造格架产生巨大的影响,不仅改造了古老的盆地和褶皱带,而且控制了新生代的走滑拉分盆地和两侧盆地的沉积速率,并对青藏高原的阶段性隆升和地壳物质东向蠕散起到调节作用。

粘滑过程中的多点错动

在实验室利用分布式多通道瞬态信号采集系统,从近场条件下研究了岩石粘滑过程中沿断层的声发射信号的初动到时和初动方向空间分布,并结合瞬态高频应变记录系统观测沿断层的应变分布特征,对粘滑过程中的震源错动行为进行了研究。实验在双轴伺服加载系统上进行。分析表明,一次突发应力降对应的瞬时震源滑动过程是由多次更短暂的微小错动组成的,每次短小错动都产生自己的声发射事件,并依次对应不同的应变值快速改变。而每一次的微小错动也可能不止一个位错点,而是对应空间上多个启动点,从振动波初动信号的空间分布形式与沿断层应变场不均匀分布中可以找到证据。近场条件下观测到的粘滑过程多点错动现象对于理解震源过程的复杂性,解释现场地震震源不确定性与震源机制解的高矛盾比现象提供了可能的依据。

连接结构接触界面非线性力学建模研究

连接界面上存在的跨尺度、多物理场和非线性行为是引起结构复杂非线性动力学的主要原因.由于连接界面的力学行为的复杂性,以及难以对连接界面进行直接试验观测,连接界面的力学建模一直是非常具有挑战性的科学问题.本文首先从分析结合面的跨尺度物理机理入手,将名义的光滑平面视作凹凸不平的粗糙面,考虑单个微凸体的黏滑摩擦行为,建立接触载荷与变形的非线性关系,然后采用GW(Greenwood and Williamson)模型数理统计方法建立整个粗糙界面的跨尺度力学模型,并与公开文献中试验结果进行对比.考虑连接界面具有典型非线性特征,提出一种改进的Iwan唯象模型,利用精细有限元方法获得非线性特征结果,采用系统辨识理论建立连接结构的降阶力学模型,并利用有限元结果进行模型验证.结果表明,本文提出的粗糙界面跨尺度模型在法向载荷较小时与试验结果吻合较好,改进的Iwan模型能够较好描述连接界面的非线性特征,并与有限元结果吻合较好.