CN108561127B - 一种基于随机模拟的地层压力预测方法 - Google Patents
一种基于随机模拟的地层压力预测方法 Download PDFInfo
- Publication number
- CN108561127B CN108561127B CN201810253863.8A CN201810253863A CN108561127B CN 108561127 B CN108561127 B CN 108561127B CN 201810253863 A CN201810253863 A CN 201810253863A CN 108561127 B CN108561127 B CN 108561127B
- Authority
- CN
- China
- Prior art keywords
- data
- mode
- formation pressure
- template
- prediction
- 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
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- 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
Abstract
本发明涉及一种基于随机模拟的地层压力预测方法,该方法根据信息熵理论,确定数据模板的最优尺寸,利用数据模板扫描地层压力的训练图像并提取训练图像的模式特征;通过主成分分析或多维尺度变换对该模式降维,并利用k‑means算法进行分类,获取分类数据的均值作为该类的代表均值,再比较已知的地层压力数据与上述均值的差异,随机选取差异最小的分类中的地层压力数据作为预测值。与现有技术相比,本发明可有效实现地层压力数据的预测,且能够提高预测速度、提高预测质量。
Description
技术领域
本发明涉及一种地层压力预测方法,尤其是涉及一种基于随机模拟的地层压力预测方法。
背景技术
地层压力是一个油田在勘探开发过程中所要取得的一个基础数据。在油田的地质勘探过程中,一个区域地层压力的准确预测能为了解该地区油气的分布、运移以及聚集等方面的规律提供重要信息,有助于设计和制定合理的开发方案、提高油气产量;在油田的钻井实践中,地层压力预测的准确程度,关系到所用钻井液密度是否准确,设计的井深结构是否合理,这有助于减少钻井过程中井喷井漏等事故的发生,使得钻井速度提高、钻井周期缩短,减少了钻井施工的成本,直接关系到油田的开发速度和经济效果,并对油田工作人员的生命财产安全有着重要影响,具有极大的经济效益和社会效益。另外,准确地预测地层压力还能够减少泥浆对地层的侵害,有利于保护油气层。
常用的地层压力监测方法主要有:dc指数法、标准化钻速法以及页岩密度法等。随着钻井技术的不断进步,近几年又有不少新的随钻地层压力监测方法相继问世,但这类方法实质上是一种地层压力的监测方法,而并非能实现对某口井或某个区域的地层压力进行可靠的预测。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于随机模拟的地层压力预测方法。
本发明的目的可以通过以下技术方案来实现:
一种基于随机模拟的地层压力预测方法,包括以下步骤:
S1:将已知的地层压力数据指定在三维区域,获取三维地质训练图像;
S2:根据随机模拟的信息熵理论获取最优尺寸的数据模板,利用数据模板扫描三维地质训练图像,获取包含训练图像结构特征的模式,并建立模式库;
S3:对模式库中的模式进行降维;
S4:对降维后的模式分类,获取多个子空间,对每个子空间的所有模式求取平均模式;
S5:在地层压力的预测区域中定义一条随机访问路径,对预测区域内的未知节点进行访问,判断预测区域内是否仍有未知节点,若是,执行下一步,否则,结束访问;
S6:判断待预测区域中是否有条件数据,若是,执行下一步,否则,从模式库随机提取一个模式后,执行步骤八;
S7:对条件数据和子空间中的平均模式求差异,对差异值最小的平均模式所对应的子空间中提取任意一个模式作为模拟结果;
S8:将提取的模式拷贝至待预测区域,获取当前预测区域的预测值;
S9:重复步骤S6到S7,继续对其他节点进行模拟,直到访问路径上的所有节点被模拟完毕。
优选地,所述的步骤S2的具体内容为:
绘制信息熵中的熵值随数据模板尺寸的变化曲线,计算每个模板尺寸对应的曲线斜率,选取斜率小于给定阈值所对应的第一个模板尺寸作为最优模板尺寸,进而确定最优尺寸的数据模板。
优选地,采用主成分分析或多维尺度变换对模式库中的模式进行降维。
优选地,采用k-means算法对降维后的模式分类。
优选地,采用欧氏距离对条件数据和子空间中的平均模式求差异。
与现有技术相比,本发明具有以下优点:
(1)本发明采用信息熵理论获取数据模板,信息熵可用来描述变量的不确定度,不确定度越高,所包含的信息就越多,利用信息熵用来捕获并处理含有地层数据特征的模式,并做特征选择,进而获取最优尺寸的数据模板,提高了后续预测过程中数据的准确性;
(2)本发明采用主成分分析或多维尺度变换对模式库中的模式进行降维处理,使训练图像的维数大大降低,从而加快模拟速度,提高整体的预测速度;
(3)本发明采用k-平均算法对降维后的模式进行了分类,k-平均算法相对可扩展,并且具有较高的效率;因划分获取的条件数据越多,预测的效果越好,进一步可提高地层压力预测的预测质量。
附图说明
图1为本发明方法的流程图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例
将地层压力的已知数据作为训练数据,并使地层压力分布在指定三维区域内,形成三维的地层训练图像。利用数据模板扫描训练图像捕获数据的概率信息。一般采用数据模板扫描训练数据(通常表现为训练图像形式)来捕获模式,这些模式可以视为训练数据的结构特征。如果按照每次移动一个节点(即一个体素),从左到右,从上到下进行扫描,那么就可以获得训练数据全部的模式库。
一般数据模板的大小由用户主观选择。一方面模板尺寸需要尽可能小以降低CPU负荷和内存需求,另一方面需要足够大以捕获训练图像中的特征信息。因此,本发明根据“信息熵”理论来选择合适的模板尺寸。根据“信息熵”理论,在信源中,考虑的不是某一单个符号发生的不确定性,而是要考虑这个信源所有可能发生情况的平均不确定性。若信源符号有n种取值:U1…Ui…Un,其对应的概率为:P1…Pi…Pn,且各种符号的出现彼此独立。这时,信源的平均不确定性应当为单个符号不确定性-logPi的统计平均值,可称为信息熵。即信息熵的表达式为:
其中对数log的底一般取2。
对于本专利来说,Pi表示每种模式Ui的对应概率,假设数据模板的形状为立方体,而各种模式是采用不同尺寸的数据模板扫描训练图像获得的节点的组合方式。熵的值在开始阶段将随着模板尺寸的变化而迅速变化,当模板尺寸达到最优时则变化速度趋于平缓。因此可以先绘制熵的值随数据模板的尺寸变化的变化曲线。通过该变化曲线可以计算每个模板尺寸对应的曲线斜率,并认为斜率小于给定阈值所对应的第一个模板尺寸为最优模板尺寸。
利用最优尺寸的数据模板扫描训练图像,可以获得训练图像的模式库,这些模式库往往维度较高。可采用主成分分析(principal component analysis,PCA)或多维尺度变换(multidimensional scaling,MDS)对训练图像中的历史数据进行降维处理,这样做会使得训练图像的维数大大降低,从而加快模拟速度,获得对应上述训练图像中模式的低维数据。
主成分分析作为一种统计方法,同时也是数学上用来降维的一种方法。主成分分析通过正交变换将一组可能存在相关性的变量转换为一组线性不相关的变量,转换后的这组变量叫主成分。主成分分析是对于原先提出的所有变量,将重复的变量(关系紧密的变量)删去多余,建立尽可能少的新变量,使得这些新变量是两两不相关的,而且这些新变量在反映模型的信息方面尽可能保持原有的信息。多维尺度变换是在低维空间去展示高维多元数据的一种可视化方法。多维尺度变换的基本目标是将原始数据“拟合”到一个低维坐标系中,使得由降维所引起的变形最小。多维尺度变换涉及到的问题可以描述为:当n个项目中的项目之间的相似性(或距离)给定时,求这些项目在低维空间的表示,并使得项目间的接近程度与原先的相似性(或距离)大体匹配。
在利用主成分分析或多维尺度变换完成训练图像中模式的降维后,需要对获取的低维模式进行分类。本发明对这些低维模式采用k-means算法(k-平均算法)进行分类。k-means算法是典型的基于距离的聚类算法,主要思想是:采用距离作为相似性的评价指标,即认为两个对象的距离越近,其相似度就越大。该算法认为簇是由距离靠近的对象组成的,因此把得到紧凑且独立的簇作为最终目标,即分类结果。在利用k-means算法划分模式空间后,获得若干“子空间”,每个子空间包含有若干的模式。对各个子空间的模式求取平均,获得对应每个子空间的平均模式。在模拟过程中,先设定一个模拟路径,利用当前区域内的地层压力条件数据与子空间的平均模式进行比较,选取与条件数据最为接近的平均模式作为当前区域的模拟结果。
如图1所示,本发明基于随机模拟的地层压力预测方法的完整流程如下:
步骤一、将已知的地层压力数据指定在三维区域,获取三维地质训练图像;
步骤二、根据随机模拟的信息熵理论获取最优尺寸的数据模板,利用数据模板扫描三维地质训练图像,获取包含训练图像结构特征的模式,并建立模式库;
步骤三、采用主成分分析或多维尺度变换对模式库中的模式进行降维;
步骤四、利用k-means算法对降维后的模式分类,获取多个子空间,对每个子空间的所有模式求取平均模式;
步骤五、在地层压力的预测区域中定义一条随机访问路径,依据随机访问路径判断预测区域内是否仍有未知节点,若是,执行下一步,否则,结束访问;
步骤六、判断待预测区域中是否有条件数据,若是,执行下一步,否则,从模式库随机提取一个模式后,执行步骤八;
步骤七、采用欧式距离函数对条件数据和分类子空间中的平均模式求差异,对差异值最小的平均模式所对应的子空间中提取任意一个模式作为模拟结果;
步骤八:将提取的模式拷贝至待预测区域,获取当前预测区域的预测值;
步骤九、重复步骤六到步骤七,继续对其他节点进行模拟,直到访问路径上的所有节点被模拟完毕。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的工作人员在本发明揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (3)
1.一种基于随机模拟的地层压力预测方法,其特征在于,包括以下步骤:
S1:将已知的地层压力数据指定在三维区域,获取三维地质训练图像;
S2:根据随机模拟的信息熵理论获取最优尺寸的数据模板,利用数据模板扫描三维地质训练图像,获取包含训练图像结构特征的模式,并建立模式库;
S3:对模式库中的模式进行降维;
S4:对降维后的模式分类,获取多个子空间,对每个子空间的所有模式求取平均模式;
S5:在地层压力的预测区域中定义一条随机访问路径,对预测区域内的未知节点进行访问,判断预测区域内是否仍有未知节点,若是,执行下一步,否则,结束访问;
S6:判断待预测区域中是否有条件数据,若是,执行下一步,否则,从模式库随机提取一个模式后,执行步骤八;
S7:采用欧氏距离对条件数据和子空间中的平均模式求差异,对差异值最小的平均模式所对应的子空间中提取任意一个模式作为模拟结果;
S8:将提取的模式拷贝至待预测区域,获取当前预测区域的预测值;
S9:重复步骤S6到S7,继续对其他节点进行模拟,直到访问路径上的所有节点被模拟完毕;
所述的步骤S2的具体内容为:
绘制信息熵中的熵值随数据模板尺寸的变化曲线,计算每个模板尺寸对应的曲线斜率,选取斜率小于给定阈值所对应的第一个模板尺寸作为最优模板尺寸,进而确定最优尺寸的数据模板。
2.根据权利要求1所述的一种基于随机模拟的地层压力预测方法,其特征在于,所述的步骤S3中,采用主成分分析或多维尺度变换对模式库中的模式进行降维。
3.根据权利要求1所述的一种基于随机模拟的地层压力预测方法,其特征在于,所述的步骤S4中,采用k-means算法对降维后的模式分类。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810253863.8A CN108561127B (zh) | 2018-03-26 | 2018-03-26 | 一种基于随机模拟的地层压力预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810253863.8A CN108561127B (zh) | 2018-03-26 | 2018-03-26 | 一种基于随机模拟的地层压力预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108561127A CN108561127A (zh) | 2018-09-21 |
CN108561127B true CN108561127B (zh) | 2022-04-01 |
Family
ID=63533243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810253863.8A Active CN108561127B (zh) | 2018-03-26 | 2018-03-26 | 一种基于随机模拟的地层压力预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108561127B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110310363A (zh) * | 2019-06-25 | 2019-10-08 | 上海电力学院 | 一种三维建模的表征单元体尺寸获取方法 |
CN112443319B (zh) * | 2019-09-05 | 2024-01-30 | 中国石油化工股份有限公司 | 一种井涌监测方法 |
CN111046341B (zh) * | 2019-12-12 | 2023-06-09 | 重庆地质矿产研究院 | 一种基于主成分分析的非常规天然气压裂效果评价及产能预测方法 |
CN113129426A (zh) * | 2019-12-31 | 2021-07-16 | 中石化石油工程技术服务有限公司 | 基于测井资料的含不确定度地层压力区域三维模型构建方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727578A (zh) * | 2009-11-30 | 2010-06-09 | 上海第二工业大学 | 基于过滤器的连续型信息预测方法 |
CN101726754A (zh) * | 2009-11-30 | 2010-06-09 | 上海第二工业大学 | 基于连续型多点地质统计法和软硬数据的图像模拟方法 |
CN102867301A (zh) * | 2012-08-29 | 2013-01-09 | 西北工业大学 | 一种根据信息熵得到图像显著性特征的方法 |
WO2014062947A3 (en) * | 2012-10-19 | 2014-06-19 | Conocophillips Company | Modeling a reservoir using multiple-point simulations |
CN106407363A (zh) * | 2016-09-08 | 2017-02-15 | 电子科技大学 | 一种基于信息熵的超高维数据降维算法 |
CN106780503A (zh) * | 2016-12-30 | 2017-05-31 | 北京师范大学 | 基于后验概率信息熵的遥感图像最优分割尺度确定方法 |
CN106887040A (zh) * | 2015-12-16 | 2017-06-23 | 中国石油大学(北京) | 多点地质统计学建模方法和装置 |
CN107818348A (zh) * | 2017-11-09 | 2018-03-20 | 上海电力学院 | 基于深度学习的空间数据随机模拟方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2066760A1 (en) * | 1989-09-20 | 1991-03-21 | Robert R. Weakley | Pore pressure prediction method |
US8219366B2 (en) * | 2009-08-26 | 2012-07-10 | Roche Molecular Sytems, Inc. | Determination of elbow values for PCR for parabolic shaped curves |
-
2018
- 2018-03-26 CN CN201810253863.8A patent/CN108561127B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101727578A (zh) * | 2009-11-30 | 2010-06-09 | 上海第二工业大学 | 基于过滤器的连续型信息预测方法 |
CN101726754A (zh) * | 2009-11-30 | 2010-06-09 | 上海第二工业大学 | 基于连续型多点地质统计法和软硬数据的图像模拟方法 |
CN102867301A (zh) * | 2012-08-29 | 2013-01-09 | 西北工业大学 | 一种根据信息熵得到图像显著性特征的方法 |
WO2014062947A3 (en) * | 2012-10-19 | 2014-06-19 | Conocophillips Company | Modeling a reservoir using multiple-point simulations |
CN106887040A (zh) * | 2015-12-16 | 2017-06-23 | 中国石油大学(北京) | 多点地质统计学建模方法和装置 |
CN106407363A (zh) * | 2016-09-08 | 2017-02-15 | 电子科技大学 | 一种基于信息熵的超高维数据降维算法 |
CN106780503A (zh) * | 2016-12-30 | 2017-05-31 | 北京师范大学 | 基于后验概率信息熵的遥感图像最优分割尺度确定方法 |
CN107818348A (zh) * | 2017-11-09 | 2018-03-20 | 上海电力学院 | 基于深度学习的空间数据随机模拟方法 |
Non-Patent Citations (2)
Title |
---|
Integration of Uncertain Data in Geostatistical;Amílcar Soares等;《Mathematical Geosciences》;20171231;第49卷;第253-273页 * |
Stochastic Simulation of Patterns Using Distance-Based;Mehrdad Honarkhah等;《Mathematical Geosciences》;20101231;第42卷;第487-517页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108561127A (zh) | 2018-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108561127B (zh) | 一种基于随机模拟的地层压力预测方法 | |
US20210390458A1 (en) | Data analytics methods for spatial data, and related systems and devices | |
JP6041923B2 (ja) | 情報処理方法及び情報処理装置 | |
CN111475596B (zh) | 一种基于多层级轨迹编码树的子段相似性匹配方法 | |
WO2018081929A1 (zh) | 一种高光谱遥感图像特征提取和分类方法及其系统 | |
CN105069447B (zh) | 一种人脸表情的识别方法 | |
CN108428015B (zh) | 一种基于历史气象数据和随机模拟的风电功率预测方法 | |
CN101930537B (zh) | 基于弯曲不变量相关特征的三维人脸识别方法及系统 | |
JP5567448B2 (ja) | 画像領域分割装置、画像領域分割方法および画像領域分割プログラム | |
WO2014109127A1 (ja) | インデックス生成装置及び方法並びに検索装置及び検索方法 | |
CN111754472A (zh) | 一种肺结节检测方法与检测系统 | |
Bhute et al. | Content based image indexing and retrieval | |
CN105809113A (zh) | 三维人脸识别方法及应用其的数据处理装置 | |
Chen et al. | Mesh convolution: a novel feature extraction method for 3d nonrigid object classification | |
CN114386466B (zh) | 一种用于脉冲星搜寻中候选体信号挖掘的并行的混合聚类方法 | |
CN112966629B (zh) | 基于图像变换和BoF模型的遥感图像场景分类方法 | |
Mercioni et al. | A survey of distance metrics in clustering data mining techniques | |
Luqman et al. | Subgraph spotting through explicit graph embedding: An application to content spotting in graphic document images | |
CN113378021A (zh) | 基于半监督的信息熵主成分分析降维方法 | |
Chen et al. | Automatic spectral representation with improved stacked spectral feature space patch (ISSFSP) for CNN-based hyperspectral image classification | |
Chang et al. | Time series representation learning with contrastive triplet selection | |
CN114359632A (zh) | 一种基于改进PointNet++神经网络的点云目标分类方法 | |
Yang et al. | Adaptive density peak clustering for determinging cluster center | |
Aman et al. | Characterizing colonic detections in CT colonography using curvature-based feature descriptor and bag-of-words model | |
Younso | Nonparametric discrimination of areal functional data |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |