基于电张量矩阵分析裂缝形态对页岩模型等效电特性的影响     

范珍珍  郭晨  刘婉萍  杨志芳 《地球物理学进展》2021,36(5):2168-2177

对具有裂缝结构的各向异性页岩模型等效电特性的研究是非常规油气勘探、页岩储层评价和测井解释的重要基础.针对具有强各向异性的裂缝页岩模型,本文提出使用二阶电学张量矩阵来对裂缝模型的电各向异性进行分析.论文工作采用有限元方法对含针状裂缝和圆盘状裂缝的理想页岩模型进行数值模拟,系统地分析了裂缝的形状、体积以及在岩体中的不同位置形态对页岩等效电参数张量的影响.结果 表明:裂缝的存在是导致岩石模型表现出电各向异性的原因之一;当裂缝的体积相同时,裂缝的形状以及裂缝的位置不同,都会使页岩模型的介电张量出现较大差异,并随裂缝相对位置及角度改变而呈现一定规律.这一研究结果有望为电磁测井和页岩油页岩气参数的反演提供关键的理论参考依据.  相似文献   

某电厂铁路边坡稳定性分析     

侯连成  王笑二  徐立军 《中国地质灾害与防治学报》2010,21(3):33-37

本文以某电厂铁路沿线边坡为研究对象,采用有限元法结合极限平衡法分析了边坡稳定性。通过有限元分析,给出了边坡岩体应力及应变分布,分析表明,边坡岩体加荷后未发现大剪应变的贯通区。刚体极限平衡法评价了沿节理滑动的边坡稳定性,三个剖面边坡稳定性系数均大于1.4。安全系数均符合相关规范要求,计算结果为边坡加固设计提供了重要依据。  相似文献   

三维井地电阻率有限元数值模拟及反演     

刘申芬  李天成  慕洪涛 《物探与化探》2009,33(1):88-90

油气资源和矿产资源勘探中,常需要较大的勘探深度和高的深部分辨能力,充分利用现有钻孔的三维井地电阻率成像技术能够适应于这种要求。通过设计一个高阻水平板状体模型,利用Ansys有限元软件对模型进行三维井地电阻率法有限元正演计算,并将模拟的响应视电阻率作为反演的输入,结合平滑约束最小二乘反演,分析了异常特征和分布规律,其结论和认识具有一定实用价值。  相似文献   

缓坡方程的有限元解法及应用          下载免费PDF全文

郑俊  李瑞杰  江森汇  罗锋 《水科学进展》2009,20(2):275-280

利用有限元方法离散椭圆型缓坡方程,能适用于复杂区域,并很好地拟合不规则边界;采用改进共轭梯度法求解离散方程组,可以大大降低计算内存要求,提高计算效率。利用结合上述两种方法的模式对规划的日照港区水域进行了波浪数值计算,并将计算结果与物理模型试验值进行比较,结果表明:该模式能适用于较大区域的波浪场计算,并可以得到较好的计算结果。  相似文献   

基于二次场电导率分块连续变化的三维可控源电磁有限元数值模拟   总被引：5，自引：5，他引：0       下载免费PDF全文

李勇  吴小平  林品荣 《地球物理学报》2015,58(3):1072-1087

从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态.  相似文献   

基于紧凑存储的任意四边形网格有限元法地震波场数值模拟     

刘瑞合  周建科  周学锋  印兴耀  曹丹平 《CT理论与应用研究》2015,24(5):667-680

采用有限元法精确模拟复杂介质条件下的地震波场,一般采用三角单元,但在节点数相同的情况下,三角单元的计算精度不如矩形单元高,采用三角单元模拟复杂界面时,编制确定结构刚度矩阵非零元素位置的程序也较为麻烦。采用矩形单元离散含有倾斜或起伏界面的地质模型时,无法避免绕射噪声,若加密网格又会增加计算量。为此,本文基于任意四边形单元模拟声波的传播,在倾斜或起伏界面条件下,可以有效避免因离散引起的“阶梯状”界面,在不增加计算量以及内存占用的前提下,有效地消除离散绕射噪声。采用对角的集中质量矩阵代替一致质量矩阵,避免矩阵的求逆运算,从而提高显式有限元法的计算效率;对结构刚度矩阵采用紧凑存储格式,每一行需要存储的元素最多为5个,同时零元素不参与运算,既减少内存的占用,又极大地提高计算效率。   相似文献   

泉州盆地地震效应特征的一维等效线性和二维非线性分析的比较     

董菲蕃  陈国兴  金丹丹 《地震学刊》2014,(2):154-160

采用等效线性动粘弹性模型描述土的动力非线性特性,基于一维等效线性波传法,对泉州盆地地震效应进行了分析;同时,采用修正Martin-Seed-Davidenkov动粘弹塑性模型描述土的动力非线性特性,对泉州盆地非线性地震效应进行了大尺度二维精细化有限元分析,研究了地形地貌和土层横向不均匀性对地震效应的影响。将两种分析结果进行对比,结果表明：①随着基岩输入地震动强度增大,地表峰值加速度PGA放大效应总体呈现减小趋势,中震与小震、大震与小震的地表PGA放大系数之比依次为0.83~0.99、0.72~0.97;②该盆地Ⅲ类场地处,基岩、地表起伏不大,且土层横向分布较均匀,两种方法计算得到的地震效应特征类似;基岩或地表起伏剧烈、土层横向分布明显不均匀的Ⅱ类场地上,二维非线性分析给出的地表PGA放大系数明显大于一维等效线性结果,两种方法得到的地表加速度反应谱及PGA随土层深度的变化特征存在显著差异,二维非线性分析给出的地表加速度反应谱大多呈现双峰甚至多峰现象,且PGA在土层特定深度处存在聚集效应,使PGA随土层深度的变化呈现非单调性。  相似文献   

三角网格谱元法地震波场数值模拟   总被引：5，自引：0，他引：5  

刘有山  滕吉文  徐涛  刘少林  司芗  马学英 《地球物理学进展》2014,(4)

谱元法结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度,是近年来发展的一种重要的地震波场数值模拟方法.经典的谱元法采用四边形(六面体)网格,利用一维Gauss-Legendre-Lobatto(GLL)积分的张量积得到对角的质量矩阵,以大大提高计算效率,但是四边形(六面体)网格不能够灵活地刻画复杂的几何模型的弯曲界面.为此,在谱元法中引入三角形(四面体)网格到二维(三维)是十分必要的.不同于经典的谱元法,在非结构化网格中不能使用GLL积分的张量积,使得非结构化网格的谱元法的实现存在着诸多的困难.目前,比较流行的三角网格谱元法,通过使用KoornwinderDubiner(KD)正交多项式,并正交化这些KD多项式构建基函数,同时利用重合的插值节点和积分节点以获取对角的质量矩阵;它所使用的积分点为优化的点集——Fekete点,且这些积分点能与四边形网格完全耦合.相比于四边形,三角网格谱元法能显著提高复杂模型的描述能力,对起伏地表模型有很大优势.本文引入高效的最佳匹配层(PML)吸收边界条件,并通过数值试验将三角网格谱元法与经典的谱元法进行对比研究.相比于经典的谱元法,三角网格谱元法显著缺点为较低的计算精度.对于7阶谱元,为了能够精确地模拟面波,三角网格谱元法需要在每个最短的面波波长内至少有11个采样点,然而经典的谱元法仅需4个采样点,并且前者所需的内存量约为后者的5.5倍.  相似文献   

起伏地表弹性波传播的间断Galerkin有限元数值模拟方法   总被引：2，自引：1，他引：1       下载免费PDF全文

薛昭  董良国  李晓波  刘玉柱 《地球物理学报》2014,57(4):1209-1223

间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.  相似文献   

水力劈裂对岩体中自然裂纹的影响研究     

石路杨  李建  许晓瑞  余天堂 《岩土力学》2016,37(10):3003-3010

建立了求解自然裂纹和水力裂纹扩展的扩展有限元法,对裂纹附近区域的节点采用广义形函数,并采用线增函数消除混合单元,以提高裂纹附近的精度。引入水力劈裂的非耦合模型,即假设裂纹中的水压力为均布力;用砂浆法（线段-线段接触法）结合增广型拉格朗日乘子法处理受压裂纹段的接触条件。并通过算例分析了以下内容：计算了受压裂纹和裂纹面分布均布水压力的水力裂纹的应力强度因子,并与解析解进行了比较,结果表明,提出的方法具有很高的精度;模拟了水力裂纹对自然裂纹面的影响,并分析了自然裂纹面上的接触力和接触状态。  相似文献