碳酸盐岩洞穴型储层双侧向测井数值模拟和响应特征(英文)

洞穴型碳酸盐岩地层具有显著的非均质性和测井预测难度大的特点,如何利用双侧向测井对识洞穴进行准确识别与定量评价,是当前的一个研究难点与热点。利用数值模拟方法计算电测井响应,可为洞穴识别与评价提供理论依据。本文根据双侧向测井原理,采用有限元素法,研究了球状洞穴模型的测井响应特征,并考察了不同尺寸的洞穴对电测井响应的影响;然后,系统研究了不同充填程度的洞穴模型的测井响应,得出了电阻率与充填程度之间的关系;最后,研究了洞穴中充填泥质、灰质、砾岩和砂泥岩薄互层等多种充填物的测井响应特征。研究结果对碳酸盐岩洞穴型储层的识别与评价具有较高的应用价值和研究意义。     

神经网络反演双侧向电阻率测井曲线的物理约束

以一种新型高分辨率双侧向测井仪器在二维轴对称地层模型中的模拟响应为训练集训练BP神经网络，得到了针对该双侧向测井仪的反演网络模型。在训练中，神经网络结构的确定一般采用交叉验证法，但这种验证法经验性的成分偏重，不能完全解决网络结构的范化问题。为此，本文在模型的训练中，在交叉验证的基础上，根据双侧向测井的原理和仪器响应特性，提出了一种新的物理约束方法；反演地层电阻率的误差应随着侵入半径的增加而加大，违反此规律的模型不予采纳。理论研究结果表明，由此得到的神经网络模型具有很好的范化能力。     

双侧向电阻率测井反演算法研究

双侧向测井仪因其激励源的复杂性给其测量数据的反演带来了极大的困难．本文利用变形玻恩（Born）迭代法对轴对称二维非均匀介质中双侧向测井仪的位场数据进行了严格反演．在每一次迭代过程中采用了高效数值模式匹配的正演计算方法,它可半解析地求出反演中所需计算的格林（Green）函数及其偏导数,同时利用所推导出的位场非线性积分方程中积分运算的半解析形式,大大提高了计算效率和反演质量．通过反演实例表明,反演结果与电导率原状剖面吻合较好,并且受噪声影响较小．     

双侧向电阻率测井反演算法研究

双侧向测井仪因其激励源的复杂性给其测量数据的反演带来了极大的困难．本文利用变形玻恩（Born）迭代法对轴对称二维非均匀介质中双侧向测井仪的位场数据进行了严格反演．在每一次迭代过程中采用了高效数值模式匹配的正演计算方法,它可半解析地求出反演中所需计算的格林（Green）函数及其偏导数,同时利用所推导出的位场非线性积分方程中积分运算的半解析形式,大大提高了计算效率和反演质量．通过反演实例表明,反演结果与电导率原状剖面吻合较好,并且受噪声影响较小．     

大岩心电阻率测量响应的数值模拟

根据电极法岩心电阻率测量原理及准静态电场理论,用三维有限元法编制了大直径岩心电阻率测量数值模拟程序,并运用该程序计算考察了裂缝对大岩心电阻率的影响。得出大岩心的电阻率随着裂缝与电场方向夹角的增大而增大,随着裂缝宽度的增加而降低。该结果可供大岩心实验系统和大岩心实验的研发参考。     

双侧向测井仪器响应的数值分析

在简化的二维全非均匀介质模型情况下，利用有限差分方法，编制了计算双侧井仪器响应的数值模拟程序。利用模拟程序，考查并分析了模型中各参数（侵入带半径、侵入带电阻率围岩电阻率、目的层厚度等）对仪器响应的数值分析，进行实际资料反演处理的初始参数选取，以提高反演的收敛速度。     

用模式匹配法计算电阻率测井响应研究进展

数值模式匹配法是地球物理电磁场计算中的一种半数值、半解析的高效算法,相对于别的数值算法具有明显的优越性,近年来在电法测井数值模拟中得到了成功应用.本文详细地介绍了数值模式匹配法的原理和模式匹配法在电法测井的发展历程和最新应用情况,系统总结了用模式匹配法计算电阻率测井响应的思路和基本步骤,并展望了其在电法测井三维数值模拟中的应用前景.     

过套管电阻率测井响应数值模拟研究进展

套管井电阻率测井仪的研制成功是套管井测井技术的重要进展,在老井地层评价中起到了越来越重要的作用.对过套管电阻率测井响应数值模拟,本文综述了近十几年来的研究现状和进展,分析了过套管电阻率测井的测量方法,比较了过套管电阻率测井响应的计算方法,讨论了多层状地层、有限长套管测井响应计算的结果,对我们最近计算的套管接箍测井响进行了分析.针对过套管电阻率测井响应计算问题,指出了在实际应用中尚需解决的实际问题和解决的基本思路,并力图展望今后的研究发展趋势和解决问题的方法.     

自适应hp-FEM在随钻电阻率测井仪器响应数值模拟中的应用

随钻电阻率测井仪器响应数值模拟在电法测井仪器研制和资料处理与解释中具有重要的意义.本文利用自适应hp有限元算法(hp-FEM)模拟和分析了地层各向异性以及相对倾斜对仪器接收信号幅值比和相位差的影响,并对仪器响应特性进行了评价;同时,利用HERMES有限元软件包与VC++开发平台相结合,可以对随钻电阻率测井仪器响应数值模拟过程中,仪器产生的感应电场在地层中和地层分界面处的分布进行显示,从而为钻井地质导向提供有效的信息.计算结果表明,自适应hp-FEM算法收敛速度快且计算精度高,非常适合于对随钻电阻率测井仪器响应进行数值模拟计算.     

长偏移距瞬变电磁测深甚晚期响应及视电阻率的数值计算

文中讨论了提高长偏移距瞬变电磁测深甚晚期响应计算精度及相应的全区视电阻率的计算方法。在频率域 ,利用直接数值积分结合连分式展开提高低频段响应的计算精度 ,从而提高时间域甚晚期响应的计算精度。同时 ,利用连分式渐进展开 ,将均匀半空间电阻率表示为该模型长偏移距瞬变电磁测深响应的反函数 ,利用迭代方法求出适合全时间段的全区视电阻率。数值结果表明 ,文中的方法可有效地提高长偏移距瞬变电磁测深 (甚 )晚期数据的模拟与解释精度     

有限元数值模拟方法在华北地区地震地质研究中的应用进展

文中综述了近30年来有限元数值模拟方法在华北地区地震地质研究中的应用和发展情况,从华北地区地震地质数值模拟的有限元模型、有限元数值模拟所采用的实测资料以及模拟结果3方面进行了讨论,综合分析了华北地区有限元数值模拟从弹性到弹塑性、黏弹性和蠕变性,从二维线性到三维非线性,从连续变形到非连续变形,从个别资料的应用到GPS网络观测资料拟合的发展过程。简单阐述了华北地区有限元数值模拟工作的意义、存在问题和发展趋势     

云南鲁甸M_S6.5地震红石岩滑坡稳定性的数值模拟

2014年8月3日鲁甸MS6.5地震触发了大量的滑坡崩塌,其中,位于鲁甸县李家山村和巧家县红石岩村交界处的牛栏江干流北岸的红石岩滑坡规模巨大,与此处位于左岸的红石岩古滑坡体的前缘部分一起堵塞了牛栏江而形成高达120m、体积达1 200×104m3的大型堰塞体。通过震后开展的野外实地调查,获得了红石岩滑坡发生处的地形地貌、地质构造、岩体结构及物质组成等资料。以这些第一手资料为基础,构建了红石岩滑坡的边坡模型,并应用边坡稳定性分析软件Geo Studio中Slope/W模块分别计算了红石岩滑坡体震前坡体安全系数和地震作用下的坡体安全系数。结果表明,红石岩滑坡体发生处的坡体安全系数在地震前为1.450,处于相对稳定状态,而鲁甸地震的地震动作用则使坡体的安全系数降低至0.962,直接导致红石岩坡体的失稳。文中进一步讨论了坡体滑动面的存在与否对坡体稳定性的影响:安全系数计算的结果表明,在中强地震作用下,先存滑动面的存在是导致大型滑坡形成的重要条件;对于高陡岩质边坡,如果没有先存滑动面,只可能形成浅表性滑坡。     

马17井映震灵敏的应力条件数学模拟研究

本文利用线弹性理论的二维有限元方法,计算了马17井的区域应力及井区平面与剖面的应力状态,研究该井映震灵敏的应力条件。计算结果表明井区在华北地区中强以上地震活动时具有一定的应力集中     

1966年以来华北地区一系列七级大震破裂过程的数值模拟

用非连续变形分析方法(DDA+FEM)数值模拟，在华北地区各地块相互制约的块体系统环境中，地块边界断层上发生1966年邢台地震、 1969年渤海地震、 1975年海城地震、 1976年唐山地震等不同类型破裂模式大震的破裂过程. 数值模拟结果给出每个大震释放的主应力场，最大剪应力变化等值线图，地震前后位移变化矢量图，发震断层滑移随时间变化以及走滑错距和应力降等震源参数. 这些结果与地震的震源机制，用地震波资料研究得到的震源参数,宏观等震线,地表观测的水平位移矢量图基本一致. 其中1969年渤海地震正交破裂模式的结果与宏观等震线及小震分布图像更接近. 1976年唐山地震复杂震源模式与该震早期余震分布图像更相符. 表明用DDA+FEM方法进行数值模拟研究，能较好地模拟地震破裂过程.     

中强震前新沂台电阻率各向异性值（S）变化特征

对新沂地电阻率观测资料进行处理的基础上,研究了地震地电阻率各向异性的前兆异常特征,结果表明：1．在中强震前,震区一定范围内的地电阻率各向异性值（S）会出现上升型异常变化,这种变化表现为在原有平稳变化的背景上的上升变化或者是在原有上升趋势背景上的加速上升变化、亦或是在原有下降趋势背景上的上升变化;2．台站周围400km内5级以上地震发生前均存在明显的地电阻率各向异性值（S）前兆异常。     

NUMERICAL SIMULATION OF DEFORMATION MOVEMENT AND STRESS ACCUMULATION IN LONGMENSHAN AND ITS ADJACENT FAULTS BEFORE WENCHUAN EARTHQUAKE

In order to reveal the deformation and cumulative stress state in Longmenshan and its adjacent faults before Wenchuan earthquake,a 3D viscoelastic finite element model,which includes Longmenshan,Longriba,Minjiang and Huya faults is built in this paper.Using the GPS measurement results of 1999-2004 as the boundary constraints,the deformation and movement of Longmenshan fault zone and its adjacent zones before Wenchuan earthquake are simulated.The conclusions are drawn in this paper as follows:First,velocity component parallel to Longmenshan Fault is mainly absorbed by Longriba Fault and velocity component perpendicular to the Longmenshan Fault is mainly absorbed by itself.Because of the barrier effect of Minjiang and Huya faults on the north section of Longmenshan Fault,the compression rate in the northern part of Longmenshan Fault is lower than that in the southern part.Second,extending from SW to NE direction along Longmenshan Fault,the angle between the main compressive stress and the direction of the fault changes gradually from the nearly vertical to 45 degrees. Compressive stress and shear stress accumulation rate is high in southwest segment of Longmenshan Fault and compressive stress is greater;the stress accumulation rate is low and the compressive stress is close to shear stress in the northeast segment of the fault.This is coincident with the fact that small and medium-sized earthquakes occurred frequently and seismic activity is strong in the southwest of the fault,and that there are only occasional small earthquakes and the seismic activity is weak in the northeast of the fault.It is also coincident with the rupture type of thrust and right-lateral strike-slip of the Wenchuan earthquake and thrust of the Lushan earthquake.Third,assuming that the same type and magnitude of earthquake requires the same amount of stress accumulation,the rupture of Minjiang Fault,the southern segment of Longmenshan Fault and the Huya Fault are mainly of thrust movement and the earthquake recurrence period of the three faults increases gradually.In the northern segment of Longriba Fault and Longmenshan Fault,earthquake rupture is of thrusting and right-lateral strike-slip. The earthquake recurrence period of former is shorter than the latter.In the southern segment of Longriba Fault,earthquake rupture is purely of right-lateral strike-slip,it is possible that the earthquake recurrence period on the fault is the shortest in the study region.     

中卫活动断裂带地震破裂危险区演化特征模拟

通过建立的岩石破裂危险度KR及断层滑动危险系数KF2个判别指标,对中卫一同心活动断裂带进行粘弹性有限元数值模拟,得到该断裂带在现代构造应力场作用下大震后的应力演变及地震破裂危险区长度、范围、峰值、衰减规律等的演化特征。结果表明,该断裂带大震后应力场的调整主要产生于震后300～400年间,震后地震破裂危险区转移到断裂带的中西段,并逐渐向东迁移,范围越来越集中,最终转移到断裂带的弧形顶点部位(红谷梁附近)。     

双相各向同性介质弹性波高精度波场分离数值模拟方法

为了便于研究双相介质固流相混合弹性波场中纵横波波场的传播规律,提出了基于交错网格的Biot双相各向同性介质弹性波动方程高精度波场分离正演数值模拟方法.采用高阶交错网格有限差分法来构建一阶双曲型双相各向同性介质弹性波动方程正演算子实现波场正演,并在每一步递推过程中,分别计算出同相和流相分量相应的散度场(纯纵波场)和旋度场...     

孕震过程中波速比及S波分裂时空变化的数值模拟

本文利用有限元方法模拟计算了强震发生前孕震区的波速比及Ｓ波分裂的时空变化，结果表明，孕震过程中，随着应力的局部集中，岩体发生扩容，引起裂隙的增生及定向排列，从而使孕震区产生各向异性，引起波速比的变化及Ｓ波分裂现象，在不同阶段，穿过同一路径的地震射线其波速比和Ｓ波分裂特征不同；在同一时刻，沿不同路径传播的地震射线其波速比及Ｓ波分裂也不同。总体看来，平行于孕震断层传播的地震波速比变化比垂直于其波速比及