In the SiChuan gas fields of China, because of the complex geologic conditions, almost all the active water invasions appear to be water breakthrough along fractures.These seal the path of gas flowing in fractured gas reservoirs, as a result largely reducing the efficiency ofgas production-Currently, it is very difficult to simulate such fracture water breakthrough using known fracture models (e.g. dual porosity, single permeability and dual porosity, dual permeability). In this paper, we established a fracture water breakthrough numerical simulation model of water drive gas reservoir after we extensively investigated typical water breakthrough in SiChuan gas fields. Next, we applied this model to simulate the bottom water gas field of WeiYuan of SiChuan province and suggested effective water preventing measures. Then we gave four kinds of water invasion models: slow water coning, vertical water breakthrough, horizontal water breakthrough, vertical water breakthrough and horizontal breakthrough. Our simulation increased our knowledge of the mechanism of water invasion. Our models have been widely used in gas production in SiChuan fractured gas reservoirs.

液体在煤层中的传输机理包括：气体在煤内表面解吸，并通过基岩和微孔隙扩散进入裂缝网络中。若岩块表面烷气体的解入速度比气、水粗在煤层割理中的流动速度快得多，那么在模拟煤层开果过程时，解吸动能是可以不考虑的。这个假设允许吸附在煤层表面上的甲烷气可以作为溶解在非流动油中的气体来模拟；煤层中的朗格缪尔等温曲线可视为常规油藏中的溶解气油比曲线。可用常规油藏模描述煤层气，而不需要对模型源码做任凭修改，基于上述思路，用热果模型模拟煤层气开采过程，并与用煤层气模拟软件（COMETPC）的计算结果进行了比较，趋势非常接近。

基于L.K.THOMAS等人的双也单渗模型的基本思想和对理论渗流问题的研究，本文提出了一个新型的全面的全隐式三维三相裂缝性油藏水平井模型，模型中对正确考虑油水、油气情况、毛管吸收、重力掺泄作用作了合理的处理。

In China, reservoirs with low permeability and heterogeneity play an important role in recoverable reserves and production rates. These types of reservoirs have complex geological features, poor reserve employment, production rate decline, and unbalanced water-intake capacity between layers, low water sweep efficiency and low recovery. These phenomena are more pronounced especially when the reservoir reaches high water cut. Different types of well stimulation such as fracturing, acidizing and recompletion etc. are used in these reservoirs recently, but these remedies just temporarily slow down production decline, while they cannot change the fact of low recovery in essence. Therefore, how to improve EOR has a remarkably practical meaning for these types of reservoirs. Artificial seismic recovery process is a very effective physical process that can improve oil recovery. Artificial seismic recovery method is a low frequency vibration recovery enhancement process. It is to apply low frequent wave produced by the surface high power extra-low frequency controllable seismic source on the formation, so as to improve the physical characteristics of groundfluid and the percolation condition of the formation, and achieve the purpose of enhancing the ultimate oil recovery. Field tests such as Xinjiang, Jiling oilfield in China are very successful and encouraging. In order to get a further understanding of vibration recovery, this paper uses a coupled geomechanic black oil simulator to simulate flow behavior of a reservoir with vibration recovery −Block 1764 of Xinjiang under surface vibration and have successfully matched the production history of this technology of the reservoir. This paper provides field test, simulation results and some in-sights to the mechanisms of the process.

气体以泡沫的形式注入油藏是一种极有效的提高采收率方法。泡沫可作为流度控制剂，也可作为封堵剂或被用来控制生产气油化。到目前为止，已有许多泡沫驱模型用来模拟泡沫驱的动态，如机理的、经验的和半经验泡沫模型。介绍了两种最常用的泡沫驱模型；PBM模型。其中，经验模型由于需要的实验数据量少，只不过需增加一个额外的表面性剂流动方程等， 所以被广泛采用。为此，开发了一个简易的泡沫马铃薯经验模型。

矩形局部网格加密方法，由于在模拟区域的任何地方均使用相同的坐标体系，因而反映井眼周围的径向流动特征不精细，同时细分网格后造成节点数增大，也限制了模拟的规模和精度。为此，本文采用了一种环绕井眼的混合网格加密方法来模拟水平井井区内的流体流动，即在水平井眼附近采用柱坐标网格，而在远离井眼的模拟区域采用直角坐标网格，从而在模拟上实现了多种坐标体系的结合，既能够快为准确地反映井眼周围液体流动特征，又能大大地减小网格数目，克服常规矩形网格在模拟过程中的不足，文中推导了混合网格的差分，介绍了网格划分及网格的推序，并将这种技术应用于三级三相水平井的模拟程序中，取得了良好的效果。

通常模拟混相驱动态采用组分模型，由于分组分模型求解工作量大，计算速度慢，要求输入的参数较多，而这些参数很难取全取准，同时组分模型不易被油藏工程师所掌握，因此，介绍了一种模拟混相驱动态的模拟模型，该模型是在黑油模型基础上发展起来的拟四级分模型，它兼有黑油模型和多组分模型的某些优点，如果油模型的稳写性好和计算速度快、组分模型的可模拟凝析气藏和混相驱等；给出了数学模型及各种物理参数以及混相流体质的处理方法。该模型还以SPE为例进行计算，结果与SPE的结果一致，从而验证了本模型的可靠性。

对达西行为描述出发点是基于1901年Forehheimer提出一个经典方程，即Forehheimer方程，利用该方程来弥补流体高速流动时达西定律的不足。在该方程中增加了一个非达西项到达西方程中以考虑达方程的缺陷。非达西项是非达西系数、流体密度和流速平方的乘积。在确定非达西效应时一个最重要方面自是心可以能准确地估计非达西系数。为此，在大量文献研究的基础上，在对高速气体渗流时非遮西行为的产生原因和目前各种非达西效应的数学插述进行分析的基础上，建立了非达西系数的求取通式，并且该通式适用储集层类型、渗透率范围较大。进一步建立了一个新型的全隐式三维气藏渗流数学模型，该模型能有效地考虑高速气体渗流对非达西现象。模型采用极小化方法求解代方程组，速度快，稳定性好。

本文介绍了4种块不完全分解ORTHOMIN加速方法其中自然同格推序上有三种，即一、二、三级块不完全分解和D4同格排序上的一级块不完全分解，并将此方法用于水平井的油藏模拟，取得了满意的效果。

对于低渗透气藏，压力变化引起的介质变形导致储层孔隙度、渗透率变化明显，对气藏开发动态有较大的影响，并且低渗透气藏的流-固耦合效应较中高渗透率气藏强得多。因此，在低渗透气藏模拟中考虑流-固耦合效应可以合理地预测气藏开发动态。此外，由于低渗透气渗汉机制的特殊性，其渗流曲线偏离达西定律，且存在启动压力梯度，造成低渗气藏的耦合模拟更加复杂。文章基于对低透气藏储层特征及渗流机理的认识，将渗流力学和弹塑性力学相结合，建立了低渗气藏流-固耦合渗流数学模型、数值模型和模拟模型，并利用该模型对四川某低渗气藏开发过程中储层层参数、开发指标的变化进行了模拟计算，还比较了考虑介质变形的耦合型和不考虑介质高精尖形的刚性模型的计算结果。