甲烷水合物形成对蒙脱石表面区域电学特性的影响:分子动力学和有限元数值模拟

含天然气水合物泥质沉积物的电学特性是反演水合物饱和度、渗透率等储层参数的重要物理性质。甲烷水合物形成对黏土矿物颗粒表面电学性质的影响机制尚不明晰,制约了对骨架颗粒表面电导及沉积物电性参数的准确建模和计算。提出分子动力学(MD)与有限元(FE)数值模拟相结合的新方法,基于Gouy-Chapman-Stern(GCS)理论建立无/含水合物-蒙脱石-Na Cl溶液体系的MD和FE数值模型,探讨甲烷水合物形成对蒙脱石表面双电层性质和表面区域电势分布的影响规律及机理。研究结果表明:(1)利用MD模型可获得含水合物条件下蒙脱石表面双电层中Stern层和扩散层的厚度以及溶液中离子的扩散系数,为FE模型提供双电层结构参数与离子运动参数,基于FE模型可分析蒙脱石表面区域电势分布特性;(2)蒙脱石表面形成稳定致密的钠离子层与界面水层,甲烷分子无法在蒙脱石表面生成水合物,因此水合物不影响Stern层而仅侵占部分扩散层区域;水合物形成过程中“排盐效应”使溶液的盐离子浓度增大,对扩散层中Na+产生排斥作用,使扩散层中的部分Na+进入Stern层,导致扩散层变薄;水合物的生长降低了溶液中离子的扩散系数;(3)水合物对扩散层的侵占使得扩散层部分区域的相对介电常数发生变化,引起该区域中电势的非均匀分布;水合物的存在不影响Stern层内的电势分布,但被水合物侵占的扩散层部分相对介电常数变小,引起扩散层内平均电势的升高;随着水合物与黏土颗粒表面距离的增加,平均电势逐渐降低。研究结果可为认识含水合物泥质沉积物中骨架颗粒的电学特性、量化颗粒表面电性参数提供理论和模型基础。

基于有限元数值模拟的直流电阻率激电找矿探测技术

为提高直流电阻率激电探测成果地质解释准确性,促进直流电阻率激电在找矿勘查中的应用,采用有限元数值模拟对直流电阻率激电探测中低阻高极化体、直立板状体、低阻盖层、高阻盖层、地表不均匀体、地形起伏及实测数据反演解释几个主要方面进行了对比分析。研究结果表明,对于实际勘查中的反演解释结果,根据异常特征构建地质体模型并进行有限元数值模拟,能快速判别探测目标体异常属性、区分虚假异常,提高定量解释水平及找矿勘查效率。研究结果为直流电阻率激电找矿探测定量解释提供了参考。

NaF-AlF_(3)体系热导率的有限元数值模拟

固态/液态铝电解质热导率是开展铝电解槽热平衡计算的重要参数,对预测电解槽侧壁上的电解质炉帮厚度变化具有重要意义。本文采取传感器测量与有限元数值模拟相结合的方法,开发研究了一种确定熔盐热导率的新方法。使用ANSYS/Fluent软件中的流体体积模块、离散坐标辐射模块和熔化/凝固模块对传感器冷却过程的热行为进行数值模拟,以测试值和模拟值的残差分析为标准来确定铝电解质的热导率。研究结果表明,该方法对单组元熔盐的热导率和双组元熔盐的热导率的计算都有效。利用该方法,建立了NaF-AlF_(3)体系的固态热导率和熔融态热导率的温度分段函数方程,用于计算该体系随着分子比变化和温度变化的热导率,计算结果与文献数据良好吻合。

有限元数值模拟技术在潜山裂缝定量预测中的应用——以珠江口盆地惠州凹陷惠州26构造为例

潜山储层的裂缝发育预测是勘探开发领域一个技术难题。目前多采用地震刻画、地震属性提取、周边钻井资料揭示以及岩芯观察等方法对地下裂缝的分布情况进行预测,但存在不少问题。笔者等通过有限元数值模拟方法对珠江口盆地惠州26构造裂缝发育情况进行预测,并提出预测流程。该方法首先采用断层生长连接分析结合演化史分析建立合适的几何模型和数理模型,通过给定力学参数和边界条件,模拟裂缝主要形成期文昌期、恩平期的古构造应力场,最终通过对张应力和剪应力大小及分布的分析来预测储层构造裂缝的发育程度。模拟结果显示,文昌期裂缝强烈发育区位于F2断层西段和F1断层西段,在惠州26构造处,构造的高部位裂缝发育强度较强,在4井位置裂缝发育较弱,在7井位置发育较好的裂缝。恩平期受应力方向和大小的改变影响,主要的裂缝发育区集中在F2断层下盘惠州26构造高部位。将预测结果与钻井结果进行比较,预测结果和钻探结果较一致,印证了该方法在古潜山裂缝预测中的实用性。有限元数值模拟方法不仅能够恢复构造裂缝形成期构造应力场的分布情况,也可以对储层构造裂缝的分布及发育程度进行预测,为裂缝型潜山油气藏勘探开发提供技术支持。

连续电性介质线源井-地电位三维有限元数值模拟

为了更好的研究油藏开发过程中剩余油的分布情况以及注水前沿的推进位置,给出了三维空间任意线电流源且电导率连续变化的井-地电位三维有限元数值模拟方法。在模拟过程中,采用异常电位法消除线电流源奇异性的影响,并且将网格节点间的电导率设计为连续变化,使其能够模拟实际中更为复杂的地质情况;考虑到采油井或注水井可能存在垂直、倾斜及水平等任意情况,推导了任意线电流源正常电位的计算公式;采用对称超松弛预条件共轭梯度法解大型线性方程组,解决了井-地电位三维模拟速度慢的问题,并编制了实用、高效的井-地电位三维模拟软件。通过对3例地电模型进行模拟计算,分析了水泥环胶结质量、回流井到观测井的距离以及目标层的深度对井-地电位观测结果的影响。

各向异性弹性波场的有限元数值模拟

研究了各向异性弹性波有限元正演系统的精度和效率问题 ,提出了一种透射加衰减的组合人工边界方案 (吸收边界条件 ) ,它对各向异性波具有较好的吸收效果 ,并且有较好的稳定性能 .均匀TI介质中的模拟获得了非常清晰的波场快照 ,其波场特征与理论分析能够准确吻合 .各向异性介质模型的地表地震记录表明 ,各向异性波炮集记录在波的类型、同相轴形态、能量分布和相位等方面与各向同性波都有很大差别 .

辽宁及邻区地壳构造应力场及其与地震活动关系的三维有限元数值模拟研究

依据最新的辽宁及其邻近地区地壳上地幔三维地震波速结构 ,并考虑地质构造、断裂以及地壳深部的构造环境等因素 ,建立研究区的三维地质模型 ;结合震源机制解与原地应力测量资料 ,分析区域现代地应力场特征 ,确定了研究区的边界约束和动力边界条件 .在以上基础上 ,采用线弹性介质模型进行研究区三维构造应力场的有限元数值模拟 ,探讨研究区地壳应力场的基本特征 .对比分析表明 ,本区地壳内应力场的三维有限元数值模拟结果与已发生的地震有着较好的对应关系 。

电子封装件受热载荷作用有限元数值模拟分析

针对典型塑料方形扁平封装体PQFP在工作过程中的受热分析问题,基于热弹性力学理论建立了二维和三维有限元数值模拟分析模型.研究了封装体在功率耗散情况下受均匀和非均匀热载作用时其材料的热膨胀和导热性质.有限元数值模拟取得了与实验一致的结果.数值结果表明,采用较小弹性模量和热膨胀率的材料可以有效地减小热应力,基板和芯片的厚度是影响封装体变形的主要参数.数值分析结果为提高封装件的可靠性和优化设计提供了理论依据.

高速冷滚打成形过程的有限元数值模拟

通过对高速冷滚打成形原理进行分析,利用有限元软件Abaqus对LY12铝合金高速冷滚打成形过程进行数值模拟;计算了成形过程中的变形击打力;并利用自行设计的试验装置进行了连续成形过程中变形击打力的试验验证。结果表明:高速冷打滚成形的数值模拟结果与实际试验结果基本一致,验证了数值模拟的有效性,为进一步深入研究高速冷滚打成形技术奠定了理论及应用基础。

起伏地形电导率连续变化的三维激电数据有限元数值模拟

为了模拟起伏地形条件下的地电体,这里给出了电导率连续变化的三维激电数据有限元数值模拟方法。在模拟过程中,采用超松弛预条件共轭梯度法解线性方程组,并探讨了模拟的计算精度和效率问题。最后,对起伏地形下的地电体进行模拟,验证了本文算法的正确性和可行性。

喜马拉雅西构造结及邻区岩石圈演化三维有限元数值模拟

用数值模拟的方法探讨了自10 Ma以来喜马拉雅弧形造山带的西构造结及邻区岩石圈的形变和位移变化特征，并采用不同的流变参数和流变结构模拟了前2 Ma内的形变和应力场的演化过程，通过对不同黏性系数、本构关系和初始地壳厚度的改变在演化过程中所起的作用研究表明：（1）采用黏弹性模型时，西构造结区黏性参数的选取对西构造结及邻区隆升高度和隆升范围有重要影响，计算结果给出西构造结区地壳黏性系数应小于1023Pa·s；（2）相对刚性的塔里木盆地采用弹性结构时，对整个西构造结区域的应力和应变场均能引致明显变化；（3）力学性质稳定块体的存在并没有使应力场的传播出现明显的解耦现象.

低渗透气藏介质变形-渗流耦合有限元数值模拟

低渗透气藏开发过程中介质变形对储层孔隙度、渗透率及产能的影响不可忽视。综合考虑储层变形、地应力、流体渗流等因素，建立了低渗透气藏介质变形一渗流耦合模型；采用有限元数值模拟技术并开发相应程序对耦合模型进行求解，定量研究开发过程中储层变形及其对储层物性和气藏产能的影响。数值模拟结果表明：非均匀水平地应力造成近井壁区域储层变形及渗透率变化各向异性，产能评价中应综合考虑介质变形的影响。

长周期地震动三维有限元数值模拟方法

采用数值方法对长周期地震动进行模拟已成为当前强地面运动预测的研究方向之一。本文提出了一种有限元并行计算的高效地震动数值模拟方法,可以同时考虑震源、传播途径和场地条件对地震动的影响,并根据昆明地区未来可能发生的设定地震对昆明市的影响地震动进行了预测。这一方法可以应用于国内已探测的发震活断层进行地震危害性评价工作。

有限元数值模拟技术应用于基坑降水设计和结构防浮分析

利用有限元数值模拟技术进行建筑基坑开挖工程降水设计,在首都机场新航站楼工程中,用此项技术进行基坑降水和结构防浮分析,取得良好了效果。与传统的设计方法相比,该技术具有能适应复杂含水地层,满足不同开挖深度的基坑降水设计,具有全面模拟整个降水过程的水量和水位变化的优点。

基于有限元数值模拟的致密储层体积压裂效果影响参数分析

致密储层体积压裂作业后形成复杂裂缝网络,通过改变压裂井附近的渗流场来增大水平井产量,目前针对致密储层体积压裂影响参数对水平井投产后不同时间段压裂效果影响的权重分析较少。通过建立适用于致密油藏油水两相有限元方程对复杂裂缝网络进行有限元数值模拟,采用正交实验与灰色关联法分析水平井体积压裂影响参数在油井投产10年内对产能的影响程度。研究结果表明,裂缝导流能力决定水平井压裂初期产能,裂缝长度决定水平井日产量的递减速度以及后期阶段稳产的产量水平,随着压裂井生产不断进行,裂缝长度取代裂缝导流能力成为影响产量的关键性因素。水平井投产不同时间段内各影响参数对压裂效果的影响程度研究可以为致密储层体积压裂水平井施工设计提供理论支持。

Muglad盆地形成力学机制的有限元数值模拟

以早白垩世Muglad盆地的地质背景为基础,运用弹性力学有限元数值模拟方法,通过对应力场的分析,为盆地形成的力学机制提供依据。模拟结果表明,Muglad盆地是受走滑作用和拉张作用共同影响形成的被动裂谷盆地,指示Muglad盆地是在中生代泛大陆裂解的大背景下,在非洲大陆的走滑作用和拉张作用的共同作用下形成和发育的,而不是简单地由走滑作用诱导形成。

硅钢片轧制过程的有限元数值模拟

利用体积成形软件DEFORM对硅钢片的轧制变形过程进行有限元数值模拟,其中包括对轧件不同厚度层应力及应变分布、整个轧制过程的扭矩分布的数值模拟。结果表明,其模拟结果与传统理论所述结果基本一致。

护舷碰撞过程有限元数值模拟

本文基于ANSYS-LSDYNA,采用有限元方法模拟了钢板撞击橡胶泡沫型护舷的全过程,分析了钢板在撞击过程中速度改变量、不同时刻护舷的应变量以及护舷等效应力,结果表明有限元是一种有效的研究护舷碰撞过程的方法。

基于有限元数值模拟的宽厚板热轧过程宽度控制ANN模型

利用商业有限元分析软件DEFORM-3D,建立宽厚板热轧过程中立辊侧压和随后平轧过程的三维热力耦合刚塑性有限元模型.模拟了不同板坯厚度、板坯宽度、立辊压下量和平辊压下率条件下320种立辊侧压和随后平轧过程,探讨了不同工艺参数对调宽效率的影响规律.模拟结果表明,调宽效率随着板坯宽度和立辊压下量的增加而升高,随着板坯厚度和平辊压下率的增加而降低.基于上述有限元数值模拟结果,借助BP人工神经网络,建立了热轧宽厚板立轧-平轧宽度控制模型,经测试,模型预测的调控效率与有限元数值模拟结果符合很好.将宽度控制模型用于宽厚板实际热轧过程的有限元数值模拟,模拟轧件宽度与实测宽度吻合很好.

剪切带中斑晶旋转的粘弹性有限元数值模拟

采用粘弹性力学分析程序对太行山阜平地区剪切带中斑晶旋转问题进行了有限元数值模拟。在计算模型中假定 :石榴石斑晶为弹性材料 ,在应力作用下只产生很小的弹性变形 ;石英、长石质的基质部分为粘弹性材料 ,它具有随时间变形不断增大的特性。计算结果显示 :沿剪切力方向的最大位移能达到模型宽度的 3～ 4倍 ,相当于简单剪切的角度达到 70°以上 ,斑晶的最大旋转角度能达到 60°左右。这样的变形量级已经十分接近野外实际情况。由于目前地壳物质的粘度值基本上是通过间接推测得到的 ,其可信度较低 ,因此文中还讨论了基质的粘度变化对变形和应力的影响 ,分别取了 7个不同量级的粘度值 (10 17～ 10 2 3Pa·s)进行计算。