CN104895550B - 一种致密气压裂水平井数值试井模型建立求解方法 - Google Patents
一种致密气压裂水平井数值试井模型建立求解方法 Download PDFInfo
- Publication number
- CN104895550B CN104895550B CN201510308810.8A CN201510308810A CN104895550B CN 104895550 B CN104895550 B CN 104895550B CN 201510308810 A CN201510308810 A CN 201510308810A CN 104895550 B CN104895550 B CN 104895550B
- Authority
- CN
- China
- Prior art keywords
- pressure
- dimensionless
- fracture
- well
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012360 testing method Methods 0.000 title claims abstract description 57
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000008878 coupling Effects 0.000 claims abstract description 20
- 238000010168 coupling process Methods 0.000 claims abstract description 20
- 238000005859 coupling reaction Methods 0.000 claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims description 28
- 230000035699 permeability Effects 0.000 claims description 21
- 230000015572 biosynthetic process Effects 0.000 claims description 16
- 239000007788 liquid Substances 0.000 claims description 14
- 230000009466 transformation Effects 0.000 claims description 13
- 239000006185 dispersion Substances 0.000 claims description 9
- 239000000203 mixture Substances 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 239000002131 composite material Substances 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000005251 gamma ray Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000007789 sealing Methods 0.000 claims description 3
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 230000035945 sensitivity Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 3
- 210000002706 plastid Anatomy 0.000 abstract 1
- 239000007789 gas Substances 0.000 description 41
- 238000005755 formation reaction Methods 0.000 description 11
- 230000007246 mechanism Effects 0.000 description 6
- 230000005514 two-phase flow Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明提供了一种致密气压裂水平井数值试井模型建立求解方法，包括以下步骤：步骤一：致密气藏压裂水平井的二维地质体和三维地质体的生成；步骤二：对生成的致密气藏压裂水平井二维地质体和三维地质体进行网格离散；步骤三：水平井筒无压差的渗流模型计算；步骤四：建立耦合模型，并对建立的耦合模型进行求解，并把得到的解生成试井理论曲线；步骤五：将步骤四中得到的理论曲线与实测曲线进行拟合，得到试井解释的参数；具有计算速度快、曲线拟合好、解释结果准确。
Description
技术领域
本发明涉及一种致密气压裂水平井数值试井模型建立求解的方法，属于石油工业油气井试井领域。
背景技术
致密气作为三大非常规天然气之一，资源量丰富，开发潜力大。致密气藏具有低渗、低压、低丰度等特点，气井自然生产能力低，需要经过储层改造措施后才具有工业开采价值。水力压裂技术和水平井技术是提高致密气藏产能的有效方法。目前致密气藏的开采普遍采用多段压裂水平井技术。
致密气藏压裂水平井试井解释是获取压后裂缝参数和储层参数的重要手段，也是对渗流机理进行直接验证的有效方法。由于致密气藏多段压裂水平井渗流机理和井型的复杂性，目前国内还没有专门针对致密气藏压裂水平井的试井解释模型，主要采用常规商业软件Saphir和EPS软件所提供的多段压裂水平井解析模型来进行解释分析，这严重影响到致密气压裂水平井试井资料的正确解释。
目前致密气压裂水平井试井解释主要存在的问题：
1)储层渗流机理方面，目前采用的试井解释模型没有考虑致密气藏的应力敏感、启动压力梯度等非线性渗流机理。基于达西线性渗流机理的常规试井模型不适用于致密气藏。倘若直接采用常规压裂水平井试井模型对致密气藏压裂水平井试井资料进行解释分析，将会给试井解释的拟合带来困难，拟合得到的结果也可能会存在很大误差。
2)储层非均质方面，常规试井解析模型假设储层为均匀介质，无法考虑储层的非均质性，然而实际储层具有明显的非均质性，储层非均质势必会对致密气渗流及井底压力响应产生重要影响，因此目前没考虑储层非均质性也会对试井曲线拟合及解释结果带来一定影响。
3)裂缝分布方面，常规压裂水平井试井解析模型目前无法处理裂缝非等间距的情况，而实际情况的压裂的位置几乎都是不等间距的，这也给实际应用带了不小的误差。
4)井筒多相流及井眼轨迹方面，目前在试井解释模型中几乎都没有考虑到井筒多相流和井眼轨迹的影响。致密气的开采通常伴随着水的产出，造成井筒内为气液两相流，气液两相流会加大井筒内流体的流动阻力，引起井筒内较大的压差，从而影响到试井解释结果。此外，目前水平井试井测试中压力计通常下入在造斜点以上10～20m的位置，距离水平井段有500m以上的距离。而在试井解释中，通常将压力计所测得的压力当作水平井筒段的压力，这势必会引起试井解释的误差。解决这一问题的最好方法就是建立考虑井筒多相流以及真实的井眼轨迹的压裂水平井试井模型。
发明内容
为了克服上述现有技术的缺点，本发明的目的在于提供一种计算速度快、曲线拟合好、解释结果准确的致密气压裂水平井数值试井模型建立求解的方法。
为了达到上述目的，本发明采取的技术方案为：一种致密气压裂水平井数值试井模型建立求解方法，包括以下步骤：
步骤一：致密气藏压裂水平井的二维地质体和三维地质体的生成；
步骤二：对生成的致密气藏压裂水平井二维地质体和三维地质体进行网格离散；
步骤三：井筒无压差的渗流模型计算；
步骤四：建立耦合模型，并对建立的耦合模型进行求解，并把得到的解生成试井理论曲线；
步骤五：讲步骤四中得到的理论曲线与实测曲线进行拟合，得到试井解释的参数。
所述的步骤一中生成致密气藏压裂水平井的二维地质体和三维地质体，具体步骤如下：
1)根据压裂水平井所处地质体的外边界，井筒内边界、裂缝以及复合区，再通过设置内外边界以及裂缝属性确定地质体的具体大小及形状，绘制建立二维地质体；
2)根据建立的二维体质和井眼轨迹以及储层上下边界的位置，利用几何体布尔运算生成三维地质体。
所述的步骤二中的二维地质体和三维地质体进行网格离散的具体步骤如下：
首先，将Netgen开源软件包成功进行编译，搭建运行环境；
然后，将二维地质体和三维地质体中的内外边界属分别按照Netgen网格文件格式的要求形成二维网格文件和三维网格文件，再按照Netgen设置的网格离散步骤进行网格离散。
所述的水平井筒无压差的渗流模型计算具体方法为：
1)考虑应力敏感的地层和裂缝渗流方程
地层渗流方程：
裂缝渗流方程：
初始条件：
pD(x,y,z,0)＝0 (3)
内边界条件：
外边界条件：
封闭：
定压：
式中符号含义：
pDR为地层区域的无量纲压力；pDf为裂缝区域的无量纲压力；tD为无量纲时间；CDL为无量纲井筒储存系数；KxD为x方向无量纲渗透率；KyD为y方向无量纲渗透率；KzD为z方向无量纲渗透率；KfD为无量纲裂缝渗透率；γD为无量纲渗透率模量；pwD为第一条裂缝与井筒交点处的无量纲压力；MPjD为j点无量纲压力与pwD之间的差值；Aj为内边界三角形无量纲面积；hD为无量纲储层厚度；St为井筒表皮系数；
2)方程求解
首先引入变换，将非线性渗流方程线性化，再采用有限元方法进行求解，变换公式为：
采用混合有限元法将变换后的地层和裂缝渗流方程联立进行求解,地层和裂缝***的有限元方程分解为地层区域的有限元方程(式9右边第一项)和代表裂缝***的有限元方程(式9右边第二项)。
A.地层区域三维有限元方程为：
B.裂缝面二维有限元方程为：
将有限元方程(10)～(15)联立组成***刚度矩阵，利用并行化的SuperLU数值求解器对大型线性方程组进行求解，可以得到整个储层的压力场分布以及内边界法向压力梯度，再由此计算各裂缝生产流量：
式中符号含义：
η为线性化变换参数；ηw为无量纲井底压力值所对应的变换参数；wf为裂缝宽度，m；wfD为无量纲裂缝宽度；LjD为裂缝内边界单元线长度；V为四面体体积；b，c，d为有限元系数；i，j，k，m为有限元四面体四个顶点序号；Qfi为第i条裂缝的流量，m3/d；Qsc为标况下气井的流量，m3/d。
所述的步骤四中的建立耦合模型，并对建立的耦合模型进行求解按照以下步骤进行的：
A，根据步骤四计算出的井筒压力pWD和裂缝流量Qfi，采用井筒多相流计算公式进行计算，得到井筒各点的无量纲压差MPiD，具体的公式如下：
井筒多相流计算的基础方程为：
式中符号含义：
ρL为液体密度，kg/m3；ρg为气体密度，kg/m3；G为气液混合物质量流量，kg/s；vm为混合物流动速度，m/s；vsg为气体表观流速，m/s；A为井筒油管截面积，m2；D为油管内径，m。
其中持液率HL和摩擦阻力系数λ采用Beggs-Brill方法计算，根据公式(17)得到井筒压力梯度，再根据井筒内边界相对距离累加得到各点与井底标准点间的压差值：
将得到的各点压差MPi进行无量纲化，即可获得各点无量纲压差MPiD。
B，把步骤A中求出的内边界各点压差MPiD，带入到渗流模型中计算出井筒压力p’WD和裂缝流量Q’fi，井筒多相流模型与渗流模型的耦合条件如下：
地层与井筒交接面上：
裂缝与井筒交接面上：
C，把前后迭代步计算得到的p’WD和pWD进行相减，当两者之间差的绝对值小于ε，则继续下一时间步的计算，当两者之差的绝对值大于或等于ε，则把当前得到的p’WD和Q’fi与前一步的值取平均后得到新的井筒压力pWD和裂缝流量Qfi，并带入到步骤A中进行迭代，直到得到的p’WD与pWD进行相减之差的绝对值小于ε，ε一般取值为10-4；
D，把步骤C确定的p’WD和对应的Q’fi进行记录，若此时时间步k<总时间步n，则继续下一时间步的计算，根据当前的井筒压力pWD，裂缝流量Qfi，从步骤A开始计算得到新一轮的p’WD和Q’fi；
E，把步骤D中n次计算中所得到p’WD和Q’fi的值生成试井理论曲线。
所述步骤五主要是理论曲线与实测曲线在双对数曲线图、半对数曲线图和全历史压力曲线图中对比拟合，根据曲线的拟合程度，可进行拟合参数的调整，最终使理论曲线与实测曲线在双对数曲线图、半对数曲线图和全历史压力曲线图中都能够得到较好的吻合，曲线拟合才完成；曲线拟合完成后，即可以得到试井解释的参数，包括裂缝参数和储层参数。
本发明采用以上技术方案，具有以下优点，具有计算速度快、曲线拟合好、解释结果准确。
附图说明
图1考虑井眼轨迹及储层非均质性的压裂水平井三维地质体及网格离散；
图2耦合模型建立及求解的流程图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
实施例1
如图2所示一种致密气压裂水平井数值试井模型建立求解方法，包括以下步骤：
步骤一：致密气藏压裂水平井的二维地质体和三维地质体的生成；主要
1)根据压裂水平井所处地质体的外边界，井筒内边界、裂缝以及复合区，再通过设置内外边界以及裂缝属性确定地质体的具体大小及形状，绘制建立二维地质体；
2)根据建立的二维体质和井眼轨迹以及储层上下边界的位置，利用几何体布尔运算生成三维地质体。
步骤二：如图1所示对生成的致密气藏压裂水平井二维地质体和三维地质体进行网格离散；首先，将Netgen开源软件包成功进行编译，搭建运行环境；
然后，将二维地质体和三维地质体中的内外边界属分别按照Netgen网格文件格式的要求形成二维网格文件和三维网格文件，再按照Netgen设置的网格离散步骤进行网格离散。
步骤三：水平井筒无压差的渗流模型计算；
1)考虑应力敏感的地层和裂缝渗流方程
地层渗流方程：
裂缝渗流方程：
初始条件：
pD(x,y,z,0)＝0 (3)
内边界条件：
外边界条件：
封闭：
定压：
式中符号含义：
pDR为地层区域的无量纲压力；pDf为裂缝区域的无量纲压力；tD为无量纲时间；CDL为无量纲井筒储存系数；KxD为x方向无量纲渗透率；KyD为y方向无量纲渗透率；KzD为z方向无量纲渗透率；KfD为无量纲裂缝渗透率；γD为无量纲渗透率模量；pwD为第一条裂缝与井筒交点处的无量纲压力；MPjD为j点无量纲压力与pwD之间的差值；Aj为内边界三角形无量纲面积；hD为无量纲储层厚度；St为井筒表皮系数；
2)方程求解
首先引入变换，将非线性渗流方程线性化，再采用有限元方法进行求解，变换公式为：
采用混合有限元法将变换后的地层和裂缝渗流方程联立进行求解,地层和裂缝***的有限元方程分解为地层区域的有限元方程(式9右边第一项)和代表裂缝***的有限元方程(式9右边第二项)。
A.地层区域三维有限元方程为：
B.裂缝面二维有限元方程为：
将有限元方程(10)～(15)联立组成***刚度矩阵，利用并行化的SuperLU数值求解器对大型线性方程组进行求解，可以得到整个储层的压力场分布以及内边界法向压力梯度，再由此计算各裂缝生产流量：
式中符号含义：
η为线性化变换参数；ηw为无量纲井底压力值所对应的变换参数；wf为裂缝宽度，m；wfD为无量纲裂缝宽度；LjD为裂缝内边界单元线长度；V为四面体体积；b，c，d为有限元系数；i，j，k，m为有限元四面体四个顶点序号；Qfi为第i条裂缝的流量，m3/d；Qsc为标况下气井的流量，m3/d。
步骤四：建立耦合模型，并对建立的耦合模型进行求解，并把得到的解生成试井理论曲线；并对建立的耦合模型进行求解按照以下步骤进行的：
A，根据步骤四计算出的井筒压力pWD和裂缝流量Qfi，采用井筒多相流计算公式进行计算，得到井筒各点的无量纲压差MPiD，具体的公式如下：
井筒多相流计算的基础方程为：
式中符号含义：
ρL为液体密度，kg/m3；ρg为气体密度，kg/m3；G为气液混合物质量流量，kg/s；vm为混合物流动速度，m/s；vsg为气体表观流速，m/s；A为井筒油管截面积，m2；D为油管内径，m。
其中持液率HL和摩擦阻力系数λ采用Beggs-Brill方法计算，根据公式(17)得到井筒压力梯度，再根据井筒内边界相对距离累加得到各点与井底标准点间的压差值：
将得到的各点压差MPi进行无量纲化，即可获得各点无量纲压差MPiD。
B，把步骤A中求出的内边界各点压差MPiD，带入到渗流模型中计算出井筒压力p’WD和裂缝流量Q’fi，井筒多相流模型与渗流模型的耦合条件如下：
地层与井筒交接面上：
裂缝与井筒交接面上：
C，把前后迭代步计算得到的p’WD和pWD进行相减，当两者之间差的绝对值小于ε，则继续下一时间步的计算，当两者之差的绝对值大于或等于ε，则把当前得到的p’WD和Q’fi与前一步的值取平均后得到新的井筒压力pWD和裂缝流量Qfi，并带入到步骤A中进行迭代，直到得到的p’WD与pWD进行相减之差的绝对值小于ε，ε一般取值为10-4；
D，把步骤C确定的p’WD和对应的Q’fi进行记录，若此时时间步k<总时间步n，则继续下一时间步的计算，根据当前的井筒压力pWD，裂缝流量Qfi，从步骤A开始计算得到新一轮的p’WD和Q’fi；
E，把步骤D中n次计算中所得到p’WD和Q’fi的值生成试井理论曲线。步骤五：讲步骤四中得到的理论曲线与实测曲线进行拟合，得到试井解释的参数；
所述步骤五主要是理论曲线与实测曲线在双对数曲线图、半对数曲线图和全历史压力曲线图中对比拟合，根据曲线的拟合程度，可进行拟合参数的调整，最终使理论曲线与实测曲线在双对数曲线图、半对数曲线图和全历史压力曲线图中都能够得到较好的吻合，曲线拟合才完成；曲线拟合完成后，即可以得到试井解释的参数，包括裂缝参数和储层参数。
实施例2
1.实现致密气藏压裂水平井地质体快速生成
首先建立二维地质体，再由二维地质体转换为三维地质体。通过编写软件几何图形绘制的代码，实现绘制压裂水平井所处地质体的外边界，井筒内边界、裂缝以及复合区等功能，再通过设置内外边界以及裂缝属性确定地质体的具体大小及形状，通过这种方式可快速建立二维地质体。由二维地质体，结合井眼轨迹以及储层上下边界的位置，利用几何体布尔运算可直接生成三维地质体，采用OpenCasCade工具进行三维视图显示，如图1所示的三维地质体。
2.实现二维地质体和三维地质体网格自动离散功能
首先将Netgen开源软件包成功进行编译，搭建运行环境，然后将二维地质体中的内外边界属性按照Netgen网格文件格式的要求形成二维网格文件，再按照Netgen设置的网格离散步骤进行网格离散。三维地质体的离散过程与二维地质体的离散过程类似，不同的是三维网格离散首先需将三维地质体输出为常用格式(如STEP格式)，再采用Netgen对该输出格式的文件进行网格离散。二维地质体的离散时间比较短，所需时间一般在10s～30s之间，三维地质体的离散时间较长，需根据网格离散的密度而定，通常离散时间不超过5分钟。网格离散后自动显示离散的网格结点及网格线，如图1所示。
3.建立试井解释模型，划分模型参数的类型
根据压裂水平井所涉及到的流动方式，建立考虑多因素耦合的试井解释模型。将模型中所涉及到的参数分为已知参数和未知参数(待解释参数)，已知参数与未知参数应该设置不同的输入接口，以免混淆。已知参数根据实际情况输入，而未知参数可在曲线拟合完成后确定。
4.编制模型数值求解的代码，实现模型快速求解
建立完试井解释模型后，根据模型特点，设计模型的数值求解算法，编制模型数值求解的计算机代码，如图2所示的模型求解流程图。首先输入基础数据，建立地质体模型，将地质体模型进行网格离散。初始步按照井筒无压差的方式进行渗流模型计算，计算得到初始的井底流压和各裂缝流量分布，根据此结果利用井筒多相流计算模型计算井筒的压差分布，再由此压差进行渗流模型计算，获取新的井底压力和裂缝流量分布，与之前的计算结果对比判断该时间步的迭代是否结束。若两者之差小于小量ε则进行下个时间步的计算，反之则继续迭代，直到稳定为止后再进行下一时间步的计算。利用此方法最终完成所有设定时间步的计算，即完成对耦合模型的求解。
5.进行实测曲线与理论曲线拟合
模型求解完成后，自动生成试井理论曲线，将理论曲线与实测曲线在双对数曲线图、半对数曲线图和全历史压力曲线图中对比拟合。双对数曲线图在主视图中，而半对数曲线图和全历史压力曲线图在辅助视图区中，根据曲线的拟合程度，可进行拟合参数的调整。最终使理论曲线与实测曲线在双对数曲线图、半对数曲线图和全历史压力曲线图中都能够得到较好的吻合，曲线拟合才完成，曲线拟合完成后，即可以得到试井解释的参数，包括裂缝参数和储层参数。
实施例3
1)致密气藏压裂水平井二维和三维地质体快速生成及显示方法；
首先建立二维地质体，再由二维地质体转换为三维地质体。根据水平井测井解释结果设置储层的非均质特征，包括孔隙度、含气饱和度以及初始渗透率分布，在地质体中划分不同的区域，每个区域代表一个均质体，每个均质体具有不同的孔隙度、含水饱和度和初始渗透率，由许多个不同均质体组成整个地质***的非均质性。均质体数目越多，参数差异越大，则地质体的非均质性越强；由二维地质体，再结合井眼轨迹以及储层上下边界的位置，可直接生成三维地质体；为了更好地进行三维地质体的显示，基于OpenCasCade工具进行三维视图的全方位显示，显示内容包括三维地质体、网格离散图、计算云图等。
2)致密气藏压裂水平井网格自动离散的实现；
将生成的地质体转换为STEP格式文件、IGES格式文件或者BREP格式文件，然后采用Netgen开源工具实现二维地质体和三维地质体的网格自动离散功能，二维网格为三角形网格，三维网格为四面体网格。针对压裂水平井试井问题的特点，在网格离散过程中通过设置裂缝及井筒离散结点的数量来自动加密井筒及裂缝附近的网格，通过设置网格离散参数来控制网格离散的质量。网格越密集，网格量越多，网格离散的时间越长，计算所需时间也越长，而网格量过小会影响到计算精度，因此实际情况应该选择适中的网格密度。通常情况下二维网格量为1万～3万之间，三维网格量为10万～50万之间。
3)致密气藏压裂水平井地层渗流及井筒多相流耦合模型的建立；
将致密气藏压裂水平井划分为三个流动区域：井筒、裂缝和地层。不同区域的流动机理不同，井筒为气液两相流，裂缝为高速非达西渗流，地层为具有应力敏感的非线性渗流。首先基于致密气压裂水平井渗流特征，将裂缝区看作为高渗透高速非达西渗流区域，由于裂缝宽度很小，裂缝区流体的流动视为二维流动，建立一种压裂水平井二维数值试井模型和三维数值试井模型；二维模型主要用于水平井筒以套管形式完井的情况，而三维模型用于水平井筒以裸眼方式完井的情况；在该渗流模型的基础上再考虑井筒气液两相流，采用Beggs-Brill方法建立井筒气液两相流的流动阻力计算方法；渗流模型与井筒多相流模型根据地层与井筒接触面上压力连续的方式耦合。
4)致密气藏压裂水平井耦合模型的快速求解方法。
采用混合有限元方法对模型进行数值求解，将渗流模型与井筒气液两相流计算模型进行耦合迭代计算。根据试井渗流模型的特点，将计算时间节点设置为对数分布，每个对数周期内设置10～20个计算点，计算的起始时间和终止时间可以根据实际情况调整。每个时间步迭代结束的条件是计算结果既满足渗流模型又满足井筒多相流计算模型。利用SuperLU求解器对二维模型的线性方程组进行求解，利用并行化的SuperLU求解器对三维模型的线性方程组进行求解，并行计算采用的线程数根据计算的网格量来确定，通常情况下并行线程数设置为6～10个即可满足计算需求。
5.进行实测曲线与理论曲线拟合
模型求解完成后，自动生成试井理论曲线，将理论曲线与实测曲线在双对数曲线图、半对数曲线图和全历史压力曲线图中对比拟合。双对数曲线图在主视图中，而半对数曲线图和全历史压力曲线图在辅助视图区中，根据曲线的拟合程度，可进行拟合参数的调整。最终使理论曲线与实测曲线在双对数曲线图、半对数曲线图和全历史压力曲线图中都能够得到较好的吻合，曲线拟合才完成，曲线拟合完成后，即可以得到试井解释的参数，包括裂缝参数和储层参数。
以上例举仅仅是对本发明的举例说明，并不构成对本发明的保护范围的限制，凡是与本发明相同或相似的设计均属于本发明的保护范围之内。
Claims (5)
1.一种致密气压裂水平井数值试井模型建立求解方法，其特征在于，包括以下步骤：
步骤一：致密气藏压裂水平井的二维地质体和三维地质体的生成；
步骤二：对生成的致密气藏压裂水平井二维地质体和三维地质体进行网格离散；
步骤三：井筒无压差的渗流模型计算；
步骤四：建立耦合模型，并对建立的耦合模型进行求解，并把得到的解生成试井理论曲线；
步骤五：将步骤四中得到的试井理论曲线与实测曲线进行拟合，得到试井解释的参数；
所述的步骤三中井筒无压差的渗流模型计算具体方法为：
1)考虑应力敏感的地层和裂缝渗流方程
地层渗流方程：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>x</mi>
<mi>D</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>x</mi>
<mi>D</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>y</mi>
<mi>D</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>y</mi>
<mi>D</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>z</mi>
<mi>D</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>z</mi>
<mi>D</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>x</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>&gamma;</mi>
<mi>D</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>y</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>&gamma;</mi>
<mi>D</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>z</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>&gamma;</mi>
<mi>D</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>z</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>=</mo>
<mfrac>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>&gamma;</mi>
<mi>D</mi>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
</mrow>
</msub>
</mrow>
</msup>
<msub>
<mi>K</mi>
<mrow>
<mi>l</mi>
<mi>D</mi>
</mrow>
</msub>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
裂缝渗流方程：
<mrow>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>y</mi>
<mi>D</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>z</mi>
<mi>D</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>K</mi>
<mrow>
<mi>f</mi>
<mi>D</mi>
</mrow>
</msub>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>f</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
初始条件：
pD(x,y,z,0)＝0 (3)
内边界条件：
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>j</mi>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&gamma;</mi>
<mi>D</mi>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>&pi;h</mi>
<mi>D</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>L</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msub>
<mi>dp</mi>
<mrow>
<mi>w</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>dt</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>w</mi>
<mi>D</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>D</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>A</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&gamma;</mi>
<mi>D</mi>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;h</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>MP</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
外边界条件：
封闭：
定压：
式中符号含义：
pDR为地层区域的无量纲压力；pDf为裂缝区域的无量纲压力；tD为无量纲时间；CDL为无量纲井筒储存系数；KxD为x方向无量纲渗透率；KyD为y方向无量纲渗透率；KzD为z方向无量纲渗透率；KfD为无量纲裂缝渗透率；γD为无量纲渗透率模量；pwD为第一条裂缝与井筒交点处的无量纲压力；MPjD为j点无量纲压力与pwD之间的差值；Aj为内边界三角形无量纲面积；hD为无量纲储层厚度；St为井筒表皮系数；
2)方程求解
首先引入变换，将非线性渗流方程线性化，再采用有限元方法进行求解，变换公式为：
<mrow>
<msub>
<mi>p</mi>
<mi>D</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&gamma;</mi>
<mi>D</mi>
</msub>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&gamma;</mi>
<mi>D</mi>
</msub>
<mi>&eta;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
采用混合有限元法将变换后的地层和裂缝渗流方程联立进行求解,地层和裂缝***的有限元方程分解为地层区域的有限元方程式9右边第一项和代表裂缝***的有限元方程式9右边第二项；
<mrow>
<mo>&Integral;</mo>
<mo>&Integral;</mo>
<msub>
<mo>&Integral;</mo>
<mi>&Omega;</mi>
</msub>
<mi>F</mi>
<mi>E</mi>
<mi>Q</mi>
<mi>d</mi>
<mi>&Omega;</mi>
<mo>=</mo>
<mo>&Integral;</mo>
<mo>&Integral;</mo>
<msub>
<mo>&Integral;</mo>
<msub>
<mi>&Omega;</mi>
<mi>m</mi>
</msub>
</msub>
<msub>
<mi>FEQd&Omega;</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>w</mi>
<mi>f</mi>
</msub>
<mo>&CenterDot;</mo>
<mo>&Integral;</mo>
<msub>
<mo>&Integral;</mo>
<mover>
<msub>
<mi>&Omega;</mi>
<mi>f</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
</msub>
<mi>F</mi>
<mi>E</mi>
<mi>Q</mi>
<mi>d</mi>
<mover>
<msub>
<mi>&Omega;</mi>
<mi>f</mi>
</msub>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
A.地层区域三维有限元方程为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>VK</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>x</mi>
<mi>D</mi>
</mrow>
</msub>
<msubsup>
<mi>b</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>y</mi>
<mi>D</mi>
</mrow>
</msub>
<msubsup>
<mi>c</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>z</mi>
<mi>D</mi>
</mrow>
</msub>
<msubsup>
<mi>d</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>10</mn>
<msub>
<mi>K</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msub>
<mi>VK</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>x</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>y</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<msub>
<mi>c</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>z</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<msub>
<mi>d</mi>
<mi>j</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>20</mn>
<msub>
<mi>K</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<msubsup>
<mi>&eta;</mi>
<mi>j</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>VK</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>x</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>y</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>z</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<msub>
<mi>d</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>20</mn>
<msub>
<mi>K</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&eta;</mi>
<mi>k</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msub>
<mi>VK</mi>
<mi>L</mi>
</msub>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>x</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>b</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>y</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<msub>
<mi>c</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>K</mi>
<mrow>
<mi>z</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<msub>
<mi>d</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>20</mn>
<msub>
<mi>K</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<msubsup>
<mi>&eta;</mi>
<mi>m</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>AK</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mn>6</mn>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>AK</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mn>12</mn>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>&eta;</mi>
<mrow>
<mi>j</mi>
<mo>/</mo>
<mi>k</mi>
<mo>/</mo>
<mi>m</mi>
</mrow>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>AK</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mn>12</mn>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>&eta;</mi>
<mrow>
<mi>k</mi>
<mo>/</mo>
<mi>m</mi>
<mo>/</mo>
<mi>j</mi>
</mrow>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mi>V</mi>
<mrow>
<mn>10</mn>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>+</mo>
<mfrac>
<mi>V</mi>
<mrow>
<mn>20</mn>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>j</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>+</mo>
<mfrac>
<mi>V</mi>
<mrow>
<mn>20</mn>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>k</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>+</mo>
<mfrac>
<mi>V</mi>
<mrow>
<mn>20</mn>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>m</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>A</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;h</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>w</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>w</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&eta;</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>D</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>&eta;</mi>
<msub>
<mo>|</mo>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>A</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;h</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&eta;</mi>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
B.裂缝面二维有限元方程为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>b</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>c</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>6</mn>
<msub>
<mi>K</mi>
<mrow>
<mi>f</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>b</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<msub>
<mi>c</mi>
<mi>j</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>12</mn>
<msub>
<mi>K</mi>
<mrow>
<mi>f</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&eta;</mi>
<mi>j</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<msub>
<mi>c</mi>
<mi>i</mi>
</msub>
<msub>
<mi>c</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>12</mn>
<msub>
<mi>K</mi>
<mrow>
<mi>f</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&eta;</mi>
<mi>k</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<mfrac>
<mi>l</mi>
<mn>3</mn>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mi>l</mi>
<mn>6</mn>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>&eta;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfrac>
<mi>A</mi>
<mrow>
<mn>6</mn>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>i</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>+</mo>
<mfrac>
<mi>A</mi>
<mrow>
<mn>12</mn>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>j</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>+</mo>
<mfrac>
<mi>A</mi>
<mrow>
<mn>12</mn>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>k</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>L</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&eta;</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<mi>&pi;</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>w</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mn>1</mn>
<mo>+</mo>
<mfrac>
<msub>
<mi>C</mi>
<mrow>
<mi>D</mi>
<mi>L</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>&Delta;t</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<msubsup>
<mi>&eta;</mi>
<mi>w</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>n</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&eta;</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>D</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&eta;</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>t</mi>
</msub>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>L</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>L</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;h</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>&eta;</mi>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
将有限元方程(10)～(15)联立组成***刚度矩阵，利用并行化的SuperLU数值求解器对大型线性方程组进行求解，得到整个储层的压力场分布以及内边界法向压力梯度，再由此计算各裂缝生产流量：
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>f</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msub>
<mi>w</mi>
<mrow>
<mi>f</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>L</mi>
<mrow>
<mi>i</mi>
<mi>D</mi>
</mrow>
</msub>
<msub>
<mi>K</mi>
<mrow>
<mi>f</mi>
<mi>D</mi>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&gamma;</mi>
<mi>D</mi>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;h</mi>
<mi>D</mi>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>i</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>n</mi>
<mo>&prime;</mo>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mo>|</mo>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
式中符号含义：
η为线性化变换参数；ηw为无量纲井底压力值所对应的变换参数；wf为裂缝宽度，m；wfD为无量纲裂缝宽度；LjD为裂缝内边界单元线长度；V为四面体体积；b，c，d为有限元系数；i，j，k，m为有限元四面体四个顶点序号；Qfi为第i条裂缝的流量，m3/d；Qsc为标况下气井的流量，m3/d。
2.根据权利要求1所述的一种致密气压裂水平井数值试井模型建立求解方法，其特征在于，所述的步骤一中生成致密气藏压裂水平井的二维地质体和三维地质体，具体步骤如下：
1)根据压裂水平井所处地质体的外边界，井筒内边界、裂缝以及复合区，再通过设置内外边界以及裂缝属性确定地质体的具体大小及形状，绘制建立二维地质体；
2)根据建立的二维体质和井眼轨迹以及储层上下边界的位置，利用几何体布尔运算生成三维地质体。
3.根据权利要求1所述的一种致密气压裂水平井数值试井模型建立求解方法，其特征在于，所述的步骤二中的二维地质体和三维地质体进行网格离散的具体步骤如下：
首先，将Netgen开源软件包成功进行编译，搭建运行环境；
然后，将二维地质体和三维地质体中的内外边界属性分别按照Netgen网格文件格式的要求形成二维网格文件和三维网格文件，再按照Netgen设置的网格离散步骤进行网格离散。
4.根据权利要求1所述的一种致密气压裂水平井数值试井模型建立求解方法，其特征在于，所述的步骤四中的建立耦合模型，并对建立的耦合模型进行求解按照以下步骤进行的：
A，根据步骤四计算出的井筒压力pWD和裂缝流量Qfi，采用井筒多相流计算公式进行计算，得到井筒各点的无量纲压差MPiD，具体的公式如下：
井筒多相流计算的基础方程为：
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>dp</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>Z</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>&rho;</mi>
<mi>l</mi>
</msub>
<msub>
<mi>H</mi>
<mi>L</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&rho;</mi>
<mi>g</mi>
</msub>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>H</mi>
<mi>L</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>&lambda;Gv</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<mn>2</mn>
<mi>A</mi>
<mi>D</mi>
</mrow>
</mfrac>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mo>{</mo>
<mo>&lsqb;</mo>
<msub>
<mi>&rho;</mi>
<mi>l</mi>
</msub>
<msub>
<mi>H</mi>
<mi>L</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&rho;</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>H</mi>
<mi>L</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msub>
<mi>v</mi>
<mi>m</mi>
</msub>
<msub>
<mi>v</mi>
<mrow>
<mi>s</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>}</mo>
<mo>/</mo>
<mi>p</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
式中符号含义：
ρL为液体密度，kg/m3；ρg为气体密度，kg/m3；G为气液混合物质量流量，kg/s；vm为混合物流动速度，m/s；vsg为气体表观流速，m/s；A为井筒油管截面积，m2；D为油管内径，m；
持液率HL和摩擦阻力系数λ采用Beggs-Brill方法计算，根据公式(17)得到井筒压力梯度，再根据井筒内边界相对距离累加得到各点与井底标准点间的压差值：
<mrow>
<msub>
<mi>MP</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>i</mi>
</munderover>
<msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>dp</mi>
<mi>w</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>Z</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>i</mi>
</msub>
<msub>
<mi>&Delta;Z</mi>
<mi>j</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
将得到的各点压差MPi进行无量纲化，即获得各点无量纲压差MPiD；
B，把步骤A中求出的内边界各点压差MPiD，带入到渗流模型中计算出井筒压力p’WD和裂缝流量Q’fi，井筒多相流模型与渗流模型的耦合条件如下：
地层与井筒交接面上：
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>R</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mo>|</mo>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>w</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>MP</mi>
<mrow>
<mi>i</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
裂缝与井筒交接面上：
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>D</mi>
<mi>f</mi>
<mi>j</mi>
</mrow>
</msub>
<msub>
<mo>|</mo>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>w</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>MP</mi>
<mrow>
<mi>j</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>20</mn>
<mo>)</mo>
</mrow>
</mrow>
C，把前后迭代步计算得到的p’WD和pWD进行相减，当两者之间差的绝对值小于ε，则继续下一时间步的计算，当两者之差的绝对值大于或等于ε，则把当前得到的p’WD和Q’fi与前一步的值取平均后得到新的井筒压力pWD和裂缝流量Qfi，并带入到步骤A中进行迭代，直到得到的p’WD与pWD进行相减之差的绝对值小于ε，ε一般取值为10-4；
D，把步骤C确定的p’WD和对应的Q’fi进行记录，若此时时间步k<总时间步n，则继续下一时间步的计算，根据当前的井筒压力pWD，裂缝流量Qfi，从步骤A开始计算得到新一轮的p’WD和Q’fi；
E，把步骤D中n次计算中所得到p’WD和Q’ fi的值生成试井理论曲线。
5.根据权利要求1所述的一种致密气压裂水平井数值试井模型建立求解方法，其特征在于，所述步骤五主要是理论曲线与实测曲线在双对数曲线图、半对数曲线图和全历史压力曲线图中对比拟合，根据曲线的拟合程度，进行拟合参数的调整，最终使理论曲线与实测曲线在双对数曲线图、半对数曲线图和全历史压力曲线图中都能够得到较好的吻合，曲线拟合才完成；曲线拟合完成后，即得到试井解释的参数，包括裂缝参数和储层参数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510308810.8A CN104895550B (zh) | 2015-06-04 | 2015-06-04 | 一种致密气压裂水平井数值试井模型建立求解方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510308810.8A CN104895550B (zh) | 2015-06-04 | 2015-06-04 | 一种致密气压裂水平井数值试井模型建立求解方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104895550A CN104895550A (zh) | 2015-09-09 |
| CN104895550B true CN104895550B (zh) | 2018-03-13 |
Family
ID=54028465
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510308810.8A Active CN104895550B (zh) | 2015-06-04 | 2015-06-04 | 一种致密气压裂水平井数值试井模型建立求解方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104895550B (zh) |
Families Citing this family (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106547930A (zh) * | 2015-09-16 | 2017-03-29 | 中国石油化工股份有限公司 | 考虑致密气藏渗流机理的泄气半径计算方法 |
| CN105334293B (zh) * | 2015-10-19 | 2017-03-29 | 中国石油天然气股份有限公司 | 多重介质中致密气的流态识别方法及其装置 |
| CN107237626B (zh) * | 2016-03-29 | 2020-11-13 | 中国石油化工股份有限公司 | 一种井钻遇溶洞体积的确定方法 |
| CN107313759B (zh) * | 2016-04-26 | 2019-08-30 | 中国石油化工股份有限公司 | 低渗稠油油藏直井热采压力分布预测方法及*** |
| CN106326549B (zh) * | 2016-08-23 | 2019-11-08 | 中国石油天然气股份有限公司 | 确定直井页岩气藏地层参数的方法和*** |
| CN106408208B (zh) * | 2016-10-14 | 2019-05-28 | 东北石油大学 | 体积压裂改造效果评价方法 |
| CN107143317B (zh) * | 2017-05-31 | 2020-09-22 | 西南石油大学 | 低渗透油藏空气驱数值模拟方法及装置 |
| CN107480383A (zh) * | 2017-08-21 | 2017-12-15 | 中国石油大学(北京) | 一种通过测压数据监测注水动态裂缝的方法 |
| CN107462936B (zh) * | 2017-08-28 | 2019-06-18 | 中国石油大学(北京) | 利用压力监测资料反演低渗透储层非达西渗流规律的方法 |
| EP3679221A1 (en) * | 2017-09-08 | 2020-07-15 | Roxar Software Solutions AS | Well fracture modelling |
| CN108019196B (zh) * | 2017-11-22 | 2019-11-05 | 中国石油大学(北京) | 一种体积压裂水平井试井解释处理方法及装置 |
| CN108756834B (zh) * | 2018-05-28 | 2020-03-31 | 中国科学技术大学 | 页岩气生产数据分析方法 |
| CN109446649B (zh) * | 2018-10-29 | 2019-07-26 | 西安石油大学 | 致密油藏体积压裂水平井三维渗流模型的建立方法 |
| CN109522634B (zh) * | 2018-11-09 | 2022-08-19 | 中国石油天然气集团有限公司 | 一种致密气多段体积压裂水平井数值分析方法 |
| CN109670220B (zh) * | 2018-12-05 | 2019-08-27 | 西南石油大学 | 一种基于非结构网格的水平井气水两相数值模拟方法 |
| CN110486008B (zh) * | 2019-09-02 | 2023-11-03 | 合肥鼎盛锦业科技有限公司 | 一种径向复合油藏的参数解释方法及*** |
| GB2595062B (en) * | 2020-04-24 | 2022-05-04 | Univ China Mining | Experimental method and system for simulating evolution of reservoir fracture stress field |
| CN111764890B (zh) * | 2020-06-08 | 2021-12-28 | 中海石油(中国)有限公司深圳分公司 | 一种特低渗储层的试井分析方法 |
| CN113944458B (zh) * | 2020-06-30 | 2023-06-30 | 中国石油天然气股份有限公司 | 单井试井分析方法 |
| CN114459911B (zh) * | 2022-01-24 | 2023-06-20 | 湖南继善高科技有限公司 | 一种油气压裂裂缝的方位解释方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6549879B1 (en) * | 1999-09-21 | 2003-04-15 | Mobil Oil Corporation | Determining optimal well locations from a 3D reservoir model |
| CN101446196A (zh) * | 2008-04-14 | 2009-06-03 | 中国石油大学(北京) | 三重介质油藏分支水平井的试井分析方法及装置 |
| CN102243680A (zh) * | 2011-07-21 | 2011-11-16 | 中国科学技术大学 | 一种网格划分方法及*** |
| CN103266881A (zh) * | 2013-05-22 | 2013-08-28 | 中国石化集团华北石油局 | 一种致密低渗气田多级压裂水平井产量预测方法 |
| CN104594872A (zh) * | 2015-01-04 | 2015-05-06 | 西南石油大学 | 一种优化致密气藏压裂水平井裂缝导流能力的方法 |
-
2015
- 2015-06-04 CN CN201510308810.8A patent/CN104895550B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6549879B1 (en) * | 1999-09-21 | 2003-04-15 | Mobil Oil Corporation | Determining optimal well locations from a 3D reservoir model |
| CN101446196A (zh) * | 2008-04-14 | 2009-06-03 | 中国石油大学(北京) | 三重介质油藏分支水平井的试井分析方法及装置 |
| CN102243680A (zh) * | 2011-07-21 | 2011-11-16 | 中国科学技术大学 | 一种网格划分方法及*** |
| CN103266881A (zh) * | 2013-05-22 | 2013-08-28 | 中国石化集团华北石油局 | 一种致密低渗气田多级压裂水平井产量预测方法 |
| CN104594872A (zh) * | 2015-01-04 | 2015-05-06 | 西南石油大学 | 一种优化致密气藏压裂水平井裂缝导流能力的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 水平井分段压裂在特低渗透油藏开发中的应用;李春芹;《西南石油大学学报（自然科学版）》;20111231;85-86 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104895550A (zh) | 2015-09-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104895550B (zh) | 一种致密气压裂水平井数值试井模型建立求解方法 | |
| WO2016192077A1 (zh) | 一种致密气压裂水平井数值试井模型建立求解方法 | |
| CN109441422B (zh) | 一种页岩气井间距优化开采方法 | |
| CN109241588B (zh) | 一种基于拟连续地质力学模型的单裂缝扩展的模拟方法 | |
| KR101671983B1 (ko) | 기하등가투과도를 이용한 다공성 매질의 투과도 산출방법 | |
| CN104533370B (zh) | 压裂水平井油藏、裂缝、井筒全耦合模拟方法 | |
| Xu et al. | Radial permeability of fractured porous media by Monte Carlo simulations | |
| CN104535472B (zh) | 一种煤岩动态渗透率检测方法及装置 | |
| CN104750896A (zh) | 一种缝洞型碳酸盐岩油藏数值模拟方法 | |
| Liu et al. | Investigation of the influence of natural cavities on hydraulic fracturing using phase field method | |
| CN114201932B (zh) | 一种复杂情况下的致密油藏压裂井试井模拟方法 | |
| CN102339325A (zh) | 一种分析离散裂缝性油藏流体流动的方法 | |
| CN108661631A (zh) | 一种页岩气井产量预测方法 | |
| NO345808B1 (no) | Flersegmentsprekker | |
| CN104463345B (zh) | 一种预测地层油水界面位置的方法及装置 | |
| Zhang et al. | High-order streamline simulation and macro-scale visualization experimental studies on waterflooding under given pressure boundaries | |
| Wang et al. | Investigation of imbibition areas during well shut-in based on mercury injection experiment and BP neural network | |
| Liu et al. | Rock crack propagation mechanism of oriented perforation hydraulic fracture under different perforation parameters | |
| Mei et al. | Fractal analysis of shape factor for matrix-fracture transfer function in fractured reservoirs | |
| Lei et al. | Characterization of multiphase flow in shale oil reservoirs considering multiscale porous media by high-resolution numerical simulation | |
| Wang et al. | Comprehensive evaluation of waterflooding front in low‐permeability reservoir | |
| Emanuel et al. | Application of streamtube techniques to full-field waterflood simulation | |
| CN116384267A (zh) | 致密储层压裂水平井的最终可采储量的确定方法及设备 | |
| CN108664678B (zh) | 一种产量预测方法 | |
| CN115828785A (zh) | 一种深部煤层气藏嵌套介质气水两相生产动态预测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20201027 Address after: 100007 Dongcheng District, Dongzhimen, China, North Street, No. 9 Oil Mansion, No. Patentee after: CHINA NATIONAL PETROLEUM Corp. Patentee after: CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. Address before: 710018 Shaanxi province Xi'an Weiyang Changqing Changqing Park District building 1207. Patentee before: CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |