[发明专利]一种大尺度裂缝性油藏数值模拟方法

摘要

本发明提供了一种大尺度裂缝性油藏数值模拟方法，属于油气田开发领域。所述方法将油藏中的大尺度裂缝作为油藏的内边界进行降维处理，建立油藏几何模型，然后采用三角形网格对油藏几何模型进行几何剖分形成离散单元，同时形成油藏几何信息描述文件、几何离散信息文件及井信息描述文件，在离散单元基础上构建数值计算单元，建立表征此类油藏流体流动的数学模型，采用积分法在每个数值计算单元上建立其数值计算格式，形成大型线性代数方程组，并进行求解，从而实现大尺度裂缝性油藏数值模拟。本发明能够精确表征大尺度裂缝性油藏，物理意义明确，对大尺度裂缝油藏数值模拟具有很好的适应性。

主权项

一种大尺度裂缝性油藏数值模拟方法，其特征在于：所述方法将油藏中的大尺度裂缝作为油藏的内边界进行降维处理，建立油藏几何模型，然后采用三角形网格对油藏几何模型进行几何剖分形成离散单元，同时形成油藏几何信息描述文件、几何离散信息文件及井信息描述文件，在离散单元基础上构建数值计算单元，建立表征此类油藏流体流动的数学模型，采用积分法在每个数值计算单元上建立其数值计算格式，形成线性代数方程组，并进行求解，从而实现大尺度裂缝性油藏数值模拟，所述方法具体包括以下步骤：(1)根据实际油藏特征及岩心、测井及流体的测试资料确定油藏储层物性参数，建立表征大尺度裂缝性油藏流体流动的数学模型；(2)依据油藏实际情况及油藏描述成果，确定油藏几何边界范围，建立油藏几何模型，形成油藏几何信息描述文件；(3)对步骤(2)建立的油藏模型进行几何离散，生成油藏模型的几何离散信息文件：平面上基岩部分采用具有厚度的Delaunay三角形单元进行划分，纵向上所有三角形厚度相加即为储层厚度，大尺度裂缝平面上用线单元代表，所述线单元具有代表裂缝直实高度的厚度，生成油藏模型的几何离散信息文件；(4)基于步骤(3)的几何离散信息文件，连接所有三角形的重心与三边的中点形成联络线，以三角形的顶点为中心，其周围的所有联络线组成的多边形即为数值计算单元；(5)计算数值计算单元节点参数值：设定饱和度变量(Sw，So)在每个数值计算单元内恒定不变，其数值通过测井数据取得，而数值计算单元的流动势变量(Φw，Φo，Φc)由以该数值计算单元的节点为顶点的相邻三角网格的值通过线性逼近计算得到：Φ(x)=Σi=1nNi(x)Φi---(7)]]>其中，w,o分别表示水相和油相；Sw为油藏含水饱和度；So为油藏含油饱和度；Φw为油藏水流动势；Φo为油藏原油流动势；Φc为毛管压力势；n为三角形单元的顶点数；Φi代表节点i的任意流动势变量；x为节点i的坐标；Ni(x)为形函数，为坐标的独立函数；由式(7)可得在一个三角形内任意变量的梯度恒为：▿Φ=Σi=1nΦi▿Ni(x)---(8)]]>(6)对步骤(1)建立的数学模型在步骤(4)建立的数值计算单元上进行积分，建立流动控制方程的数值计算方程，数值计算单元节点参数值是由步骤(5)计算得到的；(7)确定油藏的边界条件与初始条件，然后根据步骤(6)建立的每一个数值计算单元上的数值计算方程，构建线性代数方程组，并进行求解，完成大尺度裂缝性油藏的数值模拟，所述步骤(2)包括：(21)建立外边界轨迹文件，该文件的内容包括：描述油藏边界轨迹的点的X、Y坐标值及储层深度值，其中，所述外边界轨迹文件为油藏边界范围文件；(22)将大尺度裂缝作为油藏内边界进行简化降维处理，根据油藏描述得到的大尺度裂缝具体位置，建立大尺度裂缝的具体信息文件，具体信息文件的内容包括；裂缝轨迹点对应的裂缝编号、X、Y坐标值、裂缝高度、裂缝开度；(23)根据油藏中井的情况，建立井信息文件，该井信息文件的内容包括井号、X、Y坐标值、井所在的组别及井类别。