成藏三要素的耦合对高效气藏形成的控制作用——以四川盆地川东北飞仙关组鲕滩气藏为例

以川东北地区飞仙关组高效气藏为例，通过能量场的数值模拟，定量．半定量地分析油气成藏关键时刻能量场控制下所发生的烃类生成、运移和聚集及储层等要素演化以及其相互作用、相互影响．指出晚三叠世-早侏罗世是二叠系烃源岩大量成油期，断层沟通油源，油气在飞仙关组古构造圈闭内聚集成大型古油藏；晚侏罗世-白垩纪大巴山前陆盆地巨厚堆积，导致古油藏被深埋，古温度达到170～210℃，古油藏发生裂解成气，古油藏成为高效气源灶；与此同时，高温环境下发生的硫酸盐热还原作用，产生大量H2S等酸性气体，使得古油藏范围内的白云岩储集层发生深部溶蚀作用，储层物性得以改善与保持；喜马拉雅运动，在大巴山前缘产生强烈变形，古气藏发生调整、降温与降压，形成现今气藏分布状态．     

多维虚内键模型(VMIB)及其在岩体数值模拟中的应用

多维虚内键模型(VMIB, virtual multi-dimensional internal bonds)是在虚内键(VIB, virtual internal bond)理论基础上发展起来的一个多尺度力学模型. VIB理论认为固体材料在微观尺度上由随机分布的材料微粒(material particle)组成, 微粒之间由单一法向键连结; 而VMIB模型则在原VIB模型单一法向键基础上又引入了切向键用以约束微粒点对之间的相对转动自由度. 材料的宏观本构方程则由两种虚内键的刚度系数导出, 并在理论上证明了微观两种虚内键的刚度系数与宏观材料常数之间存在一一对应关系, 表明了在VMIB模型中增加切向键的充分性和必要性. 由于材料的断裂准则直接嵌入到本构方程中, 这使得VIB和VMIB模型在模拟材料断裂破坏方面有着较大的优越性. 岩体由于受到分布裂纹的切割作用, 宏观上表现为各向异性力学特征. 在岩体损伤模型中, 岩体裂纹的分布特征通常由损伤张量描述. 岩体在某一方向上的相对刚度大小决定于该方向上损伤值的大小. 在VMIB模型中, 微观虚内键的空间分布密度决定着材料宏观各方向上相对刚度的大小, 这与岩体损伤对岩体力学性能的影响是一致的. 文中在损伤张量与虚内键分布密度之间建立了定量关系, 将岩体等效为VMIB固体. 理论与试验结果表明, 这种等效的VMIB固体可以再现分布裂纹对岩体力学性能的影响, 这为进一步应用VMIB模型对岩体破坏行为进行数值模拟奠定了基础, 同时也为其他多裂纹体的数值计算提供一种新思路.     

城市防汛排水泵站侧向进水集水池水力特性研究

泵站在我国城市防汛排水等方面起着十分重要的作用。城市排水泵站常年运行的流量变化一般很大,配泵机组台数较多,为节省占地面积,常采用布置较紧凑的侧向进水方式。对侧向进水集水池,一般会出现同池各水泵配水不均匀,池内流态紊乱、水泵吸水条件差,继而导致水泵效率下降、机组振动等问题,本文提出了改善侧向进水集水池水力特性的两种技术方案,即配水孔方案与流槽式配水方案。     

一类两步方法的延迟依赖稳定性

讨论实系数延迟微分方程线性多步法的延迟依赖稳定性。重点致力于解析稳定区域和数值方法稳定区域的比较。对一类含有一个自由参数的两步方法，获得了数值稳定区域包含解析稳定区域的条件。数值试验证实了所获理论结果。     

中立型Volterra延迟积分微分方程块θ-方法的稳定性

研究了线性中立型V01terra延迟积分微分方程数值方法的稳定性，给出了块隐式θ-方法保持系统解析解不依赖于延迟的稳定性质的一个充分条件。最后，通过一些数值试验说明了这篇文章的主要结论。     

高阶FD-WENO格式用于Richtmyer-Meshkov不稳定性数值模拟

用双曲守恒律组的高阶加权本质上无振荡有限差分格式(FD-WENO)求解在众多领域有着重要应用的高马赫数二维流体界面Richtmyer-Meshkov不稳定性问题及激波与气泡相互作用问题,获得了较为理想的数值结果。表明高阶FD-WENO格式的确特别适合于求解这类既具有激波有具有复杂流动结构的流体动力学问题,它比通常的二阶Godunov格式(如PPM等)更为优越。     

不确定性结构区间分析的改进Monte Carlo方法

将结构系统中的不确定性参数用区间数来表示，建立了系统的区间有限元控制方程。利用Monte Carlo数值仿真方法对区间有限元控制方程进行求解，得到了结构的位移和应力响应区间。并针对—殷Monte Carlo方法在区间仿真中存在的缺陷，对其进行了改进，使其计算效率大为提高．数值仿真算例表明：经改进的Monte Carlo方法仿真求解计算量大为减少，且能给出较高精度的响应区间，并克服了区间运算中的扩张问题。     

煤矸石淋滤液对地下水污染的数值仿真研究

针对新邱露天煤矿矸石山附近于家沟地区地下水污染的实际情况，分析了矸石山淋滤液在地下水中的、迁移规律。在综合考虑对流、弥散、吸附解吸及存在源汇项的条件下建立煤矸石淋滤液在地下水中迁移的动力学数学模型，采用有限差分的方法对模型进行了数值求解，对不同污染强度、不同时间条件下地下水系统中淋滤液动态变化规律进行了数值仿真，预测了污染物浓度的时空分布特征，并结合实际水质监测数据对模拟结果进行了分析。模拟结果表明，污染物浓度变化趋势及模拟值与监测结果基本吻合，从而验证了模型的可靠性和实用性，为地下水污染的治理提供依据。     

气固两相流风洞内颗粒运动特性研究

对气固两相流风洞试验段的两相三维流场进行了数值模拟。对于两相流采用欧拉双流体模型,计算时选取标准k-ε两方程模型。模拟结果表明：试件迎风面和背风面颗粒浓度较大,但背风面粒子分布梯度随流速变化更明显;背风面区域由于重力作用,粒径越大的颗粒沉积越显著;试件迎风面颗粒沉积总量受粒径和流速的影响较大,但背风面沉积量受到的影响不明显。研究结果有助于分析两相流风洞中的流动特征,并能根据颗粒分布情况合理采取回收措施,为进行砂尘环境试验提供参考。     

污染物传输的两区非平衡吸附数值仿真模型

考虑不可动水和污染物的非平衡吸附特性以及微生物降解等因素,建立了土壤-水环境中污染物运移的两区非平衡吸附模型。基于自行编制的有限差分程序给出了二维仿真算例,模拟了农药等污染物在土壤-水环境中的运移过程,探讨了在不可动水和非平衡吸附影响下土壤中各相浓度分配的规律。模拟结果表明:土壤各相的浓度变化并不同步。可动水相浓度变化最快,对浓度分布起主导作用;不可动水次之;可动区吸附相浓度变化存在明显的滞后现象;而不可动区吸附相浓度变化最为缓慢,滞后现象最为明显。