CN117934738A - 区块空间应力场重构的方法、装置、处理器和计算机设备 - Google Patents
区块空间应力场重构的方法、装置、处理器和计算机设备 Download PDFInfo
- Publication number
- CN117934738A CN117934738A CN202311713645.5A CN202311713645A CN117934738A CN 117934738 A CN117934738 A CN 117934738A CN 202311713645 A CN202311713645 A CN 202311713645A CN 117934738 A CN117934738 A CN 117934738A
- Authority
- CN
- China
- Prior art keywords
- target block
- model
- fault
- stress field
- target
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000012512 characterization method Methods 0.000 claims abstract description 29
- 238000003892 spreading Methods 0.000 claims abstract description 19
- 230000007480 spreading Effects 0.000 claims abstract description 18
- 238000010276 construction Methods 0.000 claims description 15
- 238000003860 storage Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 12
- 238000012360 testing method Methods 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 8
- 239000011435 rock Substances 0.000 claims description 8
- 238000009826 distribution Methods 0.000 abstract description 21
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 238000011161 development Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 241001415846 Procellariidae Species 0.000 description 7
- 238000004590 computer program Methods 0.000 description 7
- 238000006073 displacement reaction Methods 0.000 description 6
- 239000003921 oil Substances 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000010801 machine learning Methods 0.000 description 3
- 238000005086 pumping Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 239000010779 crude oil Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000007670 refining Methods 0.000 description 2
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Computer Graphics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请公开了一种区块空间应力场重构的方法、装置、处理器和计算机设备,属于油气田开发技术领域,该方法包括:根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型;根据断层展布模型确定目标区块的全局诱导应力场;构建目标区块的三维地质力学模型;根据全局诱导应力场和三维地质力学模型,得到目标区块的目标空间应力场模型。本申请能够基于区块的断层分布情况对区块的空间应力场进行重构,有利于更加精准地还原出真实地层条件下区块的空间应力场。
Description
技术领域
本申请涉及油气田开发技术领域,具体地涉及一种区块空间应力场重构的方法、装置、处理器和计算机设备。
背景技术
随着常规油气资源的大量减少,现有的常规油气资源产量已经无法满足我们的资源消耗,非常规油气资源开发显得愈发重要。但非常规储层具有物性较差,埋藏较深,非均质性较强,天然裂缝以及断层发育等特点。而传统的压裂设计理论在确定待测区块的应力场时,并未考虑储层物性的上述特征,忽略了天然裂缝和断层对水力裂缝扩展的影响,影响最终得到的待测区块的应力场的准确性。因此,现有技术存在准确性较低的问题。
发明内容
本申请实施例的目的是提供一种区块空间应力场重构的方法、装置、处理器、计算机设备及机器可读存储介质,用以解决现有技术中存在的准确性较低的问题。
为了实现上述目的,本申请第一方面提供一种区块空间应力场重构的方法,该方法包括:
根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型;
根据断层展布模型确定目标区块的全局诱导应力场;
构建目标区块的三维地质力学模型;
根据全局诱导应力场和三维地质力学模型,得到目标区块的目标空间应力场模型。
在本申请实施例中,裂缝表征样本数据库的构建过程包括:获取目标区块的测井数据、压裂测试数据和地震数据;根据目标区块的测井数据、压裂测试数据和地震数据构建目标区块的裂缝表征样本数据库。
在本申请实施例中,根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型包括:基于裂缝表征样本数据库,通过数值模拟软件对目标区块进行断层识别,以得到目标区块的断层展布模型;其中,断层展布模型中包含目标区块中所有断层的断层数据。
在本申请实施例中,根据断层展布模型确定目标区块的全局诱导应力场包括:确定断层展布模型中各条断层的空间坐标;根据各条断层的空间坐标将目标区块的各条断层离散成多个单元网格;基于边界元方法计算出各个单元网格的诱导应力,以得到多个诱导应力;根据多个诱导应力确定目标区块的全局诱导应力场。
在本申请实施例中,构建目标区块的三维地质力学模型包括:获取目标区块的多条测井曲线;基于多条测井曲线确定目标区块的初始地质模型;基于初始地质模型,通过差值函数对初始地质模型中的未知空间区域进行补全处理,以得到目标区块的三维地质力学模型。
在本申请实施例中,基于多条测井曲线确定目标区块的初始地质模型包括:对多条测井曲线进行处理,以得到处理后的测井数据;将处理后的测井数据导入数值模拟软件,以得到目标区块的初始地质模型;其中,处理后的测井数据包括储层物性参数、岩石力学参数、地应力剖面数据。
本申请第二方面提供一种处理器,被配置成执行上述的区块空间应力场重构的方法。
本申请第三方面提供一种区块空间应力场重构的装置,装置包括:
断层展布模型确定模块,用于根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型;
全局诱导应力场确定模块,用于根据断层展布模型确定目标区块的全局诱导应力场;
三维地质力学模型构建模块,用于构建目标区块的三维地质力学模型;
目标空间应力场模型构建模块,用于根据全局诱导应力场和三维地质力学模型,得到目标区块的目标空间应力场模型。
本申请第四方面提供一种计算机设备,包括:
上述的处理器或者上述的区块空间应力场重构的装置。
本申请第五方面提供一种机器可读存储介质,机器可读存储介质上存储程序或指令,程序或指令被处理器执行时实现上述的区块空间应力场重构的方法。
上述技术方案,先根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型,然后根据断层展布模型确定目标区块的全局诱导应力场,同时构建目标区块的三维地质力学模型,最后根据全局诱导应力场和三维地质力学模型,得到目标区块的目标空间应力场模型。本申请能够基于区块的断层分布情况对区块的空间应力场进行重构,有利于更加精准地还原出真实地层条件下区块的空间应力场。
本申请实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本申请实施例的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本申请实施例,但并不构成对本申请实施例的限制。在附图中:
图1为本申请实施例提供的一种区块空间应力场重构的方法的流程示意图;
图2为本申请一具体实施例提供的边界元计算诱导应力场的物理模型示意图;
图3为本申请实施例提高的一种区块空间应力场重构的装置的结构框图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,应当理解的是,此处所描述的具体实施方式仅用于说明和解释本申请实施例,并不用于限制本申请实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。
需要说明,若本申请实施例中有涉及方向性指示(诸如上、下、左、右、前、后……),则该方向性指示仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等,如果该特定姿态发生改变时,则该方向性指示也相应地随之改变。
另外,若本申请实施例中有涉及“第一”、“第二”等的描述,则该“第一”、“第二”等的描述仅用于描述目的,而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外,各个实施例之间的技术方案可以相互结合,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在,也不在本申请要求的保护范围之内。
图1为本申请实施例提供的一种区块空间应力场重构的方法的流程示意图。如图1所示,本申请实施例提供一种区块空间应力场重构的方法,以该方法应用于处理器为例进行说明,该方法可以包括下列步骤。
步骤S101:根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型。
步骤S102:根据断层展布模型确定目标区块的全局诱导应力场。
步骤S103:构建目标区块的三维地质力学模型。
步骤S104:根据全局诱导应力场和三维地质力学模型,得到目标区块的目标空间应力场模型。
在本申请实施例中,目标区块是指待进行应力场重构的区块。预构建的目标区块的裂缝表征样本数据库是指基于技术人员采集的目标区块的现场勘测数据构建的数据库。在对目标区块的空间应力场进行重构之前,处理器可以通过获取到技术人员输入的目标区块的现场勘测数据可以构建目标区块的裂缝表征样本数据库。在进行目标区块的应力场重构时,处理器可以根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型。其中,断层展布模型中包含有目标区块中断层的分布情况的信息,以及目标区块中所有断层的断层数据。
可以理解,考虑到实际断层是三维立体结构,在得到目标区块的断层展布模型后,处理器可以根据断层展布模型确定目标区块的全局诱导应力场。在一个示例中,处理器可以通过断层展布模型中的断层数据确定每一条断层的诱导应力场,再通过每条断层的诱导应力场确定全局诱导应力场。再另一个示例中,处理器可以基于断层展布模型中的断层数据,将断层划分为多个单元网格,确定每个单元网格的诱导应力场,进而通过线性累加确定目标区块的全局诱导应力场。
可以理解,目标区块的全局诱导应力场仅考虑了断层的影响,还需结合目标区块的实际情况,因此,可以构建目标区块的三维地质力学模型。为降低工程量,可以先获取目标区块的多条测井曲线,基于多条测井曲线得到目标区块的三维地质力学模型,再通过插值的方式对未知空间区域进行补全,得到目标区块的三维地质力学模型。
最后,考虑实际地层模型的断层展步和尺度所产生的诱导应力场,将全局诱导应力导入到Petrel软件中的三维地质力学模型,实现模型的修正,明确断层周围地应力场分布特征,以得到在考虑断层影响条件下的三维地质力学模型,即目标区块的目标空间应力场模型,以还原真实地层情况下的三维地质力学模型。
在一个示例中,在得到目标区块的目标空间应力场模型后,可以在petrel软件中输入对应的单井压裂数据,实现单井压裂精细模拟,并由压裂后的裂缝长度与断层之间的距离决定是否需要考虑断层影响,进而对压裂施工设计方案进行修改与优化。在另一个示例中,基于目标区块的目标空间应力场模型,可以利用Petrel软件的kintex模块,完成单井压裂精细模拟,并基于Abaqus建立扩展有限元模型,建立工程尺度模型,明确水力裂缝遇到大尺度断层模型后的水力裂缝扩展规律及形态,并采取相对应的解决措施,实现非常规储层裂缝的均衡改造。
通过上述技术方案,对断层影响条件下的空间应力场进行重构,基于现场压裂施工需要,充分考虑了断层对水力裂缝扩展的影响,引用人工智能方法,自动识别出地下的断层展布,并通过边界元方法计算出的全局诱导应力,对petrel软件的三维地质力学模型进行修正,可实现压裂精细模拟,更为精确的还原出真实地层条件下的裂缝扩展,为后续压裂施工方案设计提出指导意见,以便及时采取对应的措施,实现裂缝均衡扩展,形成复杂裂缝网络结构,提高原油单井产量,最大程度提高非常规储层动用程度。
上述技术方案,先根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型,然后根据断层展布模型确定目标区块的全局诱导应力场,同时构建目标区块的三维地质力学模型,最后根据全局诱导应力场和三维地质力学模型,得到目标区块的目标空间应力场模型。本申请能够基于区块的断层分布情况对区块的空间应力场进行重构,有利于更加精准地还原出真实地层条件下区块的空间应力场。
在本申请实施例中,裂缝表征样本数据库的构建过程可以包括:获取目标区块的测井数据、压裂测试数据和地震数据;根据目标区块的测井数据、压裂测试数据和地震数据构建目标区块的裂缝表征样本数据库。
可以理解,目标区块的测井数据、压裂测试数据和地震数据可以通过技术人员的现场勘测采集得到。处理器获取到技术人员输入的目标区块的测井数压裂测试数据和地震数据并存储,以构成目标区块的裂缝表征样本数据库。在一个实例中,裂缝表征样本数据库中还可以包括技术人员输入的目标区块的动态资料。如此,通过构建目标区块的裂缝表征样本数据库可以为后续计算提供数据基础。
在本申请实施例中,根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型可以包括:基于裂缝表征样本数据库,通过数值模拟软件对目标区块进行断层识别,以得到目标区块的断层展布模型;
其中,断层展布模型中包含目标区块中所有断层的断层数据。
可以理解,目标区块即待进行应力场重构的区块。为了确定目标区块中的断层情况,在构建完成目标区块的裂缝表征样本数据库的情况下,可以基于已有的数值模拟软件对目标区块进行断层识别,进而得到目标区块的断层展布模型。通过断层展布模型可以确定目标区块中断层的分布情况,并且断层展布模型中包含目标区块中所有断层的断层数据。其中,断层数据可以包括目标区块的断层条数、各条断层的断层大小、各条断层的断层倾角、各条断层的断层级别、断层断距以及各条断层的走向等数据。
在一个示例中,可以利用目标区块的裂缝表征样本数据库,基于多专业一体化优化不同尺度断缝体技术组合和精度,利用petrel数值模拟软件对高精度似然体进行识别,通过深度机器学习方法细化断裂解释尺度到四级(5~10米断距),以得到段缝体解释结果,即目标区块的断层展布模型。
如此,通过对目标区块的断缝体进行识别与解释,可以得到真实底层下该目标区块的断层展布模型,有利于提高后续计算的精确性。
在本申请实施例中,根据断层展布模型确定目标区块的全局诱导应力场可以包括:确定断层展布模型中各条断层的空间坐标;根据各条断层的空间坐标将目标区块的各条断层离散成多个单元网格;基于边界元方法计算出各个单元网格的诱导应力,以得到多个诱导应力;根据多个诱导应力确定目标区块的全局诱导应力场。
具体地,考虑到实际断层是三维立体结构,在得到目标区块的断层展布模型后,可以对断层展布模型中划分号的断层读取其空间坐标。然后基于各条断层的空间坐标将目标区块中的各条断层离散成多个单元网格。在一个示例中,离散出来主要针对目标区块中的大断层,处理器可以基于Petrel平台构造建模模块将解释出来的断层采用复杂角点网格的方式嵌入地质模型,并将其导入到matlab软件中,利用该软件自动读取断层所属空间坐标,将大断层划分为一系列的小单元网格,以便后续对单个网格赋值。
进一步地,为确定目标区块的全局诱导应力场,在将断层离散呈多个单元网格后,可以通过边界元计算方法对离散后的单个网格条件下的应力场物理模型进行计算,获得不同基质模量、断层模量、断层长度下的断层软化应力值,进而推断出全局断层诱导应力场。
在一个示例中,各个单元网格的诱导应力均包括法向诱导应力和切向诱导应力,在基于边界元开发应力场计算程序后,可以通过该计算程序对离散化的各个单元网格进行矩阵方程组的计算,得到不同断层软化程度下的法向诱导应力和切向诱导应力,进而判断出断层的影响范围。其中,法向诱导应力和以满足公式(1),切向诱导应力可以满足公式(2):
其中,为第i个单元网格的法向诱导应力,/>为第i个单元网格的切向诱导应力,i为第i个单元网格的编号,j为第j个单元网格的编号,n为法向标记,s为切向标记,Gij为断层高度影响系数,/>为第j个单元网格的切向位移,/>为第j个单元网格的法向位移。和/>均为平面应变,该弹性影响系数矩阵,/>表示由单元网格j的剪切位移不连续引起的单元网格i的正应力,/>表示由于张开位移不连续引起的单元网格j在单元网格i处的正应力,/>表示由于单元网格j的剪切位移不连续引起的单元网格i处的切向应力,/>表示由于单元网格j的张开位移不连续引起的单元网格i处的切向应力。
在本申请实施例中,构建目标区块的三维地质力学模型可以包括:获取目标区块的多条测井曲线;基于多条测井曲线确定目标区块的初始地质模型;基于初始地质模型,通过差值函数对初始地质模型中的未知空间区域进行补全处理,以得到目标区块的三维地质力学模型。
具体地,为了构建目标区块的三维地质力学模型,技术人员需要事先在目标区块上通过打多口井的方式获取目标区块的数据,即目标区块的多条测井曲线。基于现场获得目标区块的多条测井曲线数据,基于多条测井曲线所得到的测井数据,将其导入到techlog软件中,对岩石力学参数和单井地应力剖面进行解释,将解释得到的单井地应力剖面导入到petrel软件,以得到目标区块的初始地质模型。可以理解,由于仅通过目标区块的多条测井曲线构建初始地质模型,初始地质模型中还存在该目标区块的未知空间区域。
进一步地,为了完善初始地质模型,得到目标区块的三维地质力学模型。可以利用techlog软件得到单井地应力剖面,将其导入到petrel软件,通过已知的三维地质力学数据对未知空间区域的参数进行预测,利用变差函数进行克里金差值,得到整个区块的三维地质力学属性,实现目标区块的三维地质力学模型的建立。如此,可以完成杨氏模量、泊松比、最小水平主应力等地质力学参数的模型建立,得到区块的含油饱和度场分布、杨氏模量场分布、泊松比分布、最小水平主应力场分布。
在本申请实施例中,基于多条测井曲线确定目标区块的初始地质模型可以包括:对多条测井曲线进行处理,以得到处理后的测井数据;将处理后的测井数据导入数值模拟软件,以得到目标区块的初始地质模型;其中,处理后的测井数据包括储层物性参数、岩石力学参数、地应力剖面数据。
具体地,基于现场获得目标区块的多条测井曲线数据,基于多条测井曲线所得到的测井数据,将其导入到techlog软件中,对岩石力学参数和单井地应力剖面进行解释,将解释得到的单井地应力剖面导入到petrel数值模拟软件,以得到目标区块的初始地质模型。其中,处理后的测井数据包括储层物性参数、岩石力学参数、地应力剖面数据。需要说明的是,本申请上述实施方式中所采用的techlog软件和petrel数值模拟软件均为举例说明,实际应用中可以根据需求选择其他软件,只要所达到的目的相同即可。
在本申请一具体实施例中,提供了应用本申请所提出的区块空间应力场重构的方法的具体实施例,该具体实施例包括以下步骤:
步骤一:利用高精度似然体及深度机器学习方法细化断层分布。
利用测录井、动态、压裂测试、地震等资料构建裂缝表征样本库,基于多专业一体化优化不同尺度断缝体技术组合和精度,利用petrel数值模拟软件对高精度似然体进行识别,通过深度机器学习方法细化断裂解释尺度到四级(5~10米断距),段缝体解释结果显示出该区块的断层条数、断层级别、断距、倾向、倾角、走向等性质。
步骤二:读取断层空间坐标并完成断层单元离散化。
通过petrel数值模拟软件解释的断层展步二维模型导入到matlab软件中,采用复杂角点网格能够很好的将断层的构造特征展现在地层层位中,方便读取各条断层的空间坐标,并将整个区块的断层划分为有限的单元网格,实现大层的离散化,满足后续模拟断层影响下的应力特征以及裂缝形态的要求。
步骤三:通过边界元方法计算断层软化应力值。
图2为本申请一具体实施例提供的边界元计算诱导应力场的物理模型示意图。如图2所示,利用自主研发的边界元开发应力场计算程序,研究全局断层诱导应力场分布机理,进而计算离散后的断层诱导应力场大小,机理模型中设置基质模量28GPa,断层长度为500m,硬化断层模量为42GPa,软化断层模量为14GPa。得到断层影响范围为120m,峰值大小为5~6MPa。
步骤四:计算单井地应力剖面并建立三维地质力学模型。
基于techlog解释的测井数据,明确储层物性参数、岩石力学参数、地应力剖面,并将处理后的测井数据导入到petrel数值模拟软件,利用插值函数方法完成三维地质力学建模,得到区块的含油饱和度场分布、杨氏模量场分布、泊松比分布、最小水平主应力场分布。
步骤五:将断层诱导应力场导入地质模型实现三维应力场重构。
开展断层影响下的地应力反演,计算真实断层影响下应力场分布,结合Petrel软件地质力学模型,更新地应力,断层造成地应力降低(垂直断层5-6MPa,沿断层4-5MPa)左右,影响距离约100-150m,超过150m后几乎不受断层影响。
步骤六:开展单井压裂精细压裂模拟,明确断层影响下裂缝扩展形态。
基于Petrel软件的kintex模块,输入射孔位置,段簇间距、泵注参数、液体参数等压裂施工参数,实现单井压裂精细化模拟。在此基础上,基于Abaqus软件,建立尺寸100×100m,断层60m的工程尺度大小的断层模型,研究大断层对水力裂缝扩展影响,硬化断层排斥水力裂缝扩展,导致水力裂缝转向与泵注压力抬升,需考虑高泵压风险。断层软化吸引截断水力裂缝扩展,产生迂曲缝,易出现砂堵和漏失风险。
这样,通过上述技术方案可以实现目标区块的压裂精细模拟,更为精确的还原出真实地层条件下的裂缝扩展。为后续压裂施工方案设计提出指导意见,以便及时采取对应的措施。实现裂缝均衡扩展,形成复杂裂缝网络结构,提高原油单井产量,最大程度提高非常规储层动用程度。
本申请实施例还提供一种处理器,被配置成执行上述实施方式中的区块空间应力场重构的方法。
图3为本申请实施例提高的一种区块空间应力场重构的装置的结构框图。如图3所示,本申请实施例还提供一种区块空间应力场重构的装置300,装置300包括:
断层展布模型确定模块310,用于根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型。
全局诱导应力场确定模块320,用于根据断层展布模型确定目标区块的全局诱导应力场。
三维地质力学模型构建模块330,用于构建目标区块的三维地质力学模型。
目标空间应力场模型构建模块340,用于根据全局诱导应力场和三维地质力学模型,得到目标区块的目标空间应力场模型。
上述区块空间应力场重构的装置300,先根据预构建的目标区块的裂缝表征样本数据库,确定目标区块的断层展布模型,然后根据断层展布模型确定目标区块的全局诱导应力场,同时构建目标区块的三维地质力学模型,最后根据全局诱导应力场和三维地质力学模型,得到目标区块的目标空间应力场模型。本申请能够基于区块的断层分布情况对区块的空间应力场进行重构,有利于更加精准地还原出真实地层条件下区块的空间应力场。
在一个实例中,裂缝表征样本数据库的构建过程包括:获取目标区块的测井数据、压裂测试数据和地震数据;根据目标区块的测井数据、压裂测试数据和地震数据构建目标区块的裂缝表征样本数据库。
在一个实例中,断层展布模型确定模块310还用于:基于裂缝表征样本数据库,通过数值模拟软件对目标区块进行断层识别,以得到目标区块的断层展布模型;其中,断层展布模型中包含目标区块中所有断层的断层数据。
在一个实例中,全局诱导应力场确定模块320还用于:确定断层展布模型中各条断层的空间坐标;根据各条断层的空间坐标将目标区块的各条断层离散成多个单元网格;基于边界元方法计算出各个单元网格的诱导应力,以得到多个诱导应力;根据多个诱导应力确定目标区块的全局诱导应力场。
在一个实例中,三维地质力学模型构建模块330还用于:获取目标区块的多条测井曲线;基于多条测井曲线确定目标区块的初始地质模型;基于初始地质模型,通过差值函数对初始地质模型中的未知空间区域进行补全处理,以得到目标区块的三维地质力学模型。
在一个实例中,三维地质力学模型构建模块330还用于:对多条测井曲线进行处理,以得到处理后的测井数据;将处理后的测井数据导入数值模拟软件,以得到目标区块的初始地质模型;其中,处理后的测井数据包括储层物性参数、岩石力学参数、地应力剖面数据。
本申请实施例还提供一种计算机设备,包括:上述实施方式中的处理器或者上述实施方式中的区块空间应力场重构的装置。
本申请实施例还提供一种机器可读存储介质,机器可读存储介质上存储程序或指令,程序或指令被处理器执行时实现上述实施方式中的区块空间应力场重构的方法。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (10)
1.一种区块空间应力场重构的方法,其特征在于,所述方法包括:
根据预构建的目标区块的裂缝表征样本数据库,确定所述目标区块的断层展布模型;
根据所述断层展布模型确定所述目标区块的全局诱导应力场;
构建所述目标区块的三维地质力学模型;
根据所述全局诱导应力场和所述三维地质力学模型,得到所述目标区块的目标空间应力场模型。
2.根据权利要求1所述的方法,其特征在于,所述裂缝表征样本数据库的构建过程包括:
获取所述目标区块的测井数据、压裂测试数据和地震数据;
根据所述目标区块的测井数据、所述压裂测试数据和所述地震数据构建所述目标区块的裂缝表征样本数据库。
3.根据权利要求1所述的方法,其特征在于,所述根据预构建的目标区块的裂缝表征样本数据库,确定所述目标区块的断层展布模型包括:
基于所述裂缝表征样本数据库,通过数值模拟软件对所述目标区块进行断层识别,以得到所述目标区块的断层展布模型;
其中,所述断层展布模型中包含所述目标区块中所有断层的断层数据。
4.根据权利要求1所述的方法,其特征在于,所述根据所述断层展布模型确定所述目标区块的全局诱导应力场包括:
确定所述断层展布模型中各条断层的空间坐标;
根据所述各条断层的空间坐标将所述目标区块的各条断层离散成多个单元网格;
基于边界元方法计算出各个单元网格的诱导应力,以得到多个诱导应力;
根据所述多个诱导应力确定所述目标区块的全局诱导应力场。
5.根据权利要求1所述的方法,其特征在于,所述构建所述目标区块的三维地质力学模型包括:
获取所述目标区块的多条测井曲线;
基于所述多条测井曲线确定所述目标区块的初始地质模型;
基于所述初始地质模型,通过差值函数对所述初始地质模型中的未知空间区域进行补全处理,以得到所述目标区块的三维地质力学模型。
6.根据权利要求5所述的方法,其特征在于,所述基于所述多条测井曲线确定所述目标区块的初始地质模型包括:
对所述多条测井曲线进行处理,以得到处理后的测井数据;
将所述处理后的测井数据导入数值模拟软件,以得到所述目标区块的初始地质模型;
其中,所述处理后的测井数据包括储层物性参数、岩石力学参数、地应力剖面数据。
7.一种处理器,其特征在于,被配置成执行根据权利要求1至6中任意一项所述的区块空间应力场重构的方法。
8.一种区块空间应力场重构的装置,其特征在于,所述装置包括:
断层展布模型确定模块,用于根据预构建的目标区块的裂缝表征样本数据库,确定所述目标区块的断层展布模型;
全局诱导应力场确定模块,用于根据所述断层展布模型确定所述目标区块的全局诱导应力场;
三维地质力学模型构建模块,用于构建所述目标区块的三维地质力学模型;
目标空间应力场模型构建模块,用于根据所述全局诱导应力场和所述三维地质力学模型,得到所述目标区块的目标空间应力场模型。
9.一种计算机设备,其特征在于,包括:
根据权利要求7所述的处理器或者根据权利要求8所述的区块空间应力场重构的装置。
10.一种机器可读存储介质,所述机器可读存储介质上存储程序或指令,其特征在于,所述程序或所述指令被处理器执行时实现根据权利要求1至6中任意一项所述的区块空间应力场重构的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311713645.5A CN117934738A (zh) | 2023-12-13 | 2023-12-13 | 区块空间应力场重构的方法、装置、处理器和计算机设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311713645.5A CN117934738A (zh) | 2023-12-13 | 2023-12-13 | 区块空间应力场重构的方法、装置、处理器和计算机设备 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117934738A true CN117934738A (zh) | 2024-04-26 |
Family
ID=90769198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311713645.5A Pending CN117934738A (zh) | 2023-12-13 | 2023-12-13 | 区块空间应力场重构的方法、装置、处理器和计算机设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117934738A (zh) |
-
2023
- 2023-12-13 CN CN202311713645.5A patent/CN117934738A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112381937B (zh) | 一种基于钻孔和复杂地质剖面的多源地质数据耦合建模方法 | |
AU2011283192B2 (en) | Methods and systems for machine-learning based simulation of flow | |
US10087721B2 (en) | Methods and systems for machine—learning based simulation of flow | |
AU2014404414B2 (en) | Optimizing multistage hydraulic fracturing design based on three-dimensional (3D) continuum damage mechanics | |
AU2011283193B2 (en) | Methods and systems for machine-learning based simulation of flow | |
EP3018502A2 (en) | Modeling fluid-conducting fractures in reservoir simulation grids | |
US20130096899A1 (en) | Methods And Systems For Machine - Learning Based Simulation of Flow | |
CN106014399B (zh) | 一种非均质地层高精度三维地应力模型建立方法 | |
EP3102968B1 (en) | Method to generate the in-situ state of stress in a domain omega in a geological structure | |
BRPI1101623A2 (pt) | método para modelagem de fraturas de um campo, mìdia legìvel por computador armazenando instruções para modelagem de fraturas de um campo, e sistema para determinar fraturas em um campo | |
US11353622B2 (en) | Systems and methods for hydrocarbon reservoir three dimensional unstructured grid generation and development | |
Zuo et al. | Investigation on failure behavior of collapse column in China’s coal mine based on discontinuous deformation numerical method | |
CN117934738A (zh) | 区块空间应力场重构的方法、装置、处理器和计算机设备 | |
Omagbon et al. | Experiences in developing a dual porosity model of the Leyte geothermal production field | |
AU2020250261A1 (en) | Automatic cartesian gridding with logarithmic refinement at arbitrary locations | |
CN112419493A (zh) | 页岩储层三维属性模型建立方法及装置 | |
CN113568047B (zh) | 一种叠前成像道集生成方法和装置 | |
US10417354B2 (en) | Model order reduction technique for discrete fractured network simulation | |
Lin et al. | The development of GIS-PMWIN and its application for mine-water modelling in the Far West Rand, South Africa | |
Székely | Three-dimensional mesh resolution control in finite difference groundwater flow models through boxed spatial zooming | |
CN118364703A (zh) | 用于确定油藏中co2有效埋存量的方法及装置 | |
CN118442044A (zh) | 用于确定多井平台的压裂方案的方法、装置及处理器 | |
CN113449407A (zh) | 油藏地质模型的优化方法、装置、电子设备及存储介质 | |
CN109190082A (zh) | 航空重力数据向下延拓时抑制边界效应的余弦函数扩边法 | |
CN115482344A (zh) | 基于钻孔数据的地层力学参数有限元建模方法及终端设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |