CN105005079B - 一种井曲线反演方法 - Google Patents

一种井曲线反演方法 Download PDF

Info

Publication number
CN105005079B
CN105005079B CN201510413557.2A CN201510413557A CN105005079B CN 105005079 B CN105005079 B CN 105005079B CN 201510413557 A CN201510413557 A CN 201510413557A CN 105005079 B CN105005079 B CN 105005079B
Authority
CN
China
Prior art keywords
well
curve
seismic properties
inversion
seismic
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
Application number
CN201510413557.2A
Other languages
English (en)
Other versions
CN105005079A (zh
Inventor
姜玉新
雷克辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Pst Service Corp
Original Assignee
Pst Service Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Pst Service Corp filed Critical Pst Service Corp
Priority to CN201510413557.2A priority Critical patent/CN105005079B/zh
Publication of CN105005079A publication Critical patent/CN105005079A/zh
Application granted granted Critical
Publication of CN105005079B publication Critical patent/CN105005079B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明涉及一种井曲线反演方法，所述方法包括(1)选择工区的井和井曲线；(2)井曲线时深转换；(3)沿井轨迹提取地震属性；(4)井曲线地震属性交互关系分析；(5)选择层位；(6)建立地层格架；(7)建立初始模型；(8)神经网络训练；(9)输出训练结果。本发明的有益效果为：本发明由RBF径向基函数神经网络算法进行井曲线反演算法把地质统计和神经网络技术很好地融合到一起，实现了多属性驱动，地质统计为基础的裂缝密度反演。该方法是对现有随机模拟和多属性井曲线反演的一次超越。它即能得到稳定的反演结果，又能体现地震多属性对反演结果的非线性约束，预测精度高于其他单纯依赖地质统计的软件。

Description

一种井曲线反演方法
技术领域
本发明属于石油勘探开发技术，具体涉及一种井曲线反演方法。
背景技术
目前井曲线反演方法有两种：井曲线的随机反演和神经网络井曲线反演。井曲线的随机反演方法采用协模拟方法，用地震数据做“软”约束进行井曲线反演，地震属性在反演中约束作用的大小取决于地震属性和井曲线的相关性，相关系数越大，地震起的约束作用越大。神经网络井曲线反演方法利用神经网络建立井曲线和地震属性的非相关关系，并利用该相关关系进行井曲线反演。上述的两种反演方法都有一定的弊端，井曲线的随机反演方法只能引入一个地震属性做约束，同时约束过程是简单的线性映射，如果地震信息和井曲线的相关性较复杂，地震信息对反演结果的约束作用就大大降低，反演结果基本就是插值结果。神经网络算法不够稳定，如果地震属性和井曲线的相关性不是非常明显，神经网络的井曲线反演结果往往不好，而且多解性很强。
发明内容
为了解决现有技术存在的上述问题，本发明提供了一种井曲线反演方法。采用径向基函数神经网络(RBF)把地质统计和神经网络技术很好地融合到一起，实现了多属性驱动，地质统计为基础的井曲线反演。该方法是对现有随机模拟和多属性井曲线反演的一次超越。它即能得到稳定的反演结果，又能体现地震多属性对反演结果的非线性约束，预测精度比其他单纯依赖地质统计的软件
本发明所采用的技术方案为：
一种井曲线反演方法，其改进之处在于：所述方法包括
(1)选择工区的井和井曲线；
(2)井曲线时深转换；
(3)沿井轨迹提取地震属性；
(4)井曲线地震属性交互关系分析；
(5)选择层位；
(6)建立地层格架；
(7)建立初始模型；
(8)神经网络训练；
(9)输出训练结果。
优选的，所述步骤(1)包括选择曲线目录树中目标井曲线和该曲线的井或在编辑框中输入曲线名直接显示选择的对象。
优选的，所述步骤(2)包括采用数据管理单元将时深转换为时深表；时深转换后生成新的时间域的曲线。
优选的，所述步骤(3)包括用户界面中选择地震测网，该测网中的数据列表中，选择一个或多个地震属性，沿井轨迹提取选择的地震属性。
优选的，所述步骤(4)包括将在目录树中显示选择的井曲线，重采样经过时深转换后，得到的时间域曲线以及沿井轨迹提取的地震属性，进入交会图分析界面，选择X轴和Y轴的井曲线，设定每条曲线的左右刻度，绘制出交会图；在交会图界面中可选择不同的回归公式对交会图进行回归分析。
优选的，所述步骤(5)包括采用层位建立地层格架，在地层格架内进行井曲线反演；选择层位用户界面，选择用户建立地层格架的层位。
优选的，所述步骤(6)包括利用选择的层位建立地层格架，选择地层格架的顶底层位，对于上超选择平行于顶，底部削截；顶超选择平行于底，顶部削截。
优选的，所述步骤(7)包括建立初始模型，初始模型为利用层位做约束，采用插值算法计算；建立初始模型插值方法可采用距离反加权插值方法或克里金插值方法。
优选的，所述步骤(8)包括从列表中选择地震属性，并给定地震属性与井曲线的相关系数，相关系数通过井曲线与地震属性交会图获得，选择地震属性后，进行神经网络训练；利用地震属性与井曲线的相关关系得到由该地震属性重构的新曲线，将重构的新曲线值作为高斯分布的均值，重构曲线与样本曲线的方差作为高斯分布的方差，高斯分布作为神经网络的输入层，经过样本点的网络训练确定输入层到隐层，隐层到输入层的权系数。
优选的，所述步骤(9)包括根据步骤(8)的神经网络训练结果计算井曲线反演结果；选择输出的数据格式和线道范围；输出的数据体为利用地震属性做约束，采用RBF神经网络预测得到的井曲线反演数据体，对于井曲线反演，输入的是裂缝密度曲线，地震属性利用地震几何属性做约束，输出的数据体为裂缝密度数据体。
本发明的有益效果为：
本发明由RBF径向基函数神经网络算法进行井曲线反演算法把地质统计和神经网络技术很好地融合到一起，实现了多属性驱动，地质统计为基础的裂缝密度反演。该方法是对现有随机模拟和多属性井曲线反演的一次超越。它即能得到稳定的反演结果，又能体现地震多属性对反演结果的非线性约束，预测精度高于其他单纯依赖地质统计的软件。
附图说明
图1是本发明提供的的一种井曲线反演方法示意图；
图2是本发明提供的选择井和井曲线界面示意图；
图3是本发明提供的井曲线深时转换界面示意图；
图4是本发明提供的沿井轨迹提取地震属性界面示意图；
图5是本发明提供的时间域井曲线和地震属性界面示意图；
图6是本发明提供的地震属性与井曲线交会分析界面示意图；
图7是本发明提供的建立地层格架界面示意图；
图8是本发明提供的选择层位界面示意图；
图9是本发明提供的地质统计学方法建立初始模型界面示意图；
图10是本发明提供的神经网络训练界面示意图；
图11是本发明提供的神经网格结构示意图。
具体实施方式
下面将结合附图，对本发明做进一步说明：
名词解释：
井曲线反演技术是由井出发，利用地震信息做约束，反演井曲线的三维变化规律。当前的石油勘探开发已经从简单的构造油气藏，转向勘探和开发岩性等隐蔽油气藏。
井曲线反演能很好预测储层和油气的空间分布规律，是进行隐蔽油气探勘开发的关键技术。
如图1所示，本发明提供的一种井曲线反演方法；具体采用下述技术方案实现了RBF(径向基函数神经网络)井曲线反演：
第一步：如图2所示，选择工区的井和井曲线目录树。在底部选择框中选择目标曲线，这时目录树中会自动把包含该曲线的井选择下来。选择井曲线的另外一种方式是在下部的编辑框中输入曲线名，如果同类曲线在不同井中名字不用，可输入曲线名和多个别名，用“，”号分开。例如对时差曲线输入AC，DT。
第二步：井曲线时深转换
如图3所示，由于井曲线反演要用地震属性做约束，需要对井曲线进行时深转换。该页面主要用于对井数据进行时深转换。
时深转换所需的时深表可以通过数据管理输入。
Time Space：时间采样间隔，通常情况选择与输入地震数据体一致。
Resample mode:重采样方式，可选择临近点、算术平均和几何平均三种方式。
时深转换后会生成新的时间域的曲线，新曲线名是原曲线名+“_time”。
第三步：沿井轨迹提取地震属性
如图4所示，为用户界面，首先在Survey栏选择地震测网，Seismic Data列表把该测网中的所有数据列出来，选择一个或多个地震属性，流程会沿井轨迹把选择的地震属性提取出来。
第四步:井曲线地震属性交互关系分析
如图5所示，用户界面，在目录树中会显示选择的井曲线，时深转化重采样得到的时间域曲线以及沿井轨迹提取的地震属性。
如图6所示，用户可以按View Cross plot进入交会图分析界面。选择X轴的Y轴的井曲线，定义每条曲线的左右刻度，就可绘制出交会图。在交会图界面中按右键可选择不同的回归公式对交会图进行回归分析。
第五步：选择层位
如图7所示，井曲线反演需用层位建立地层格架，在地层格架内进行井曲线反演。图7位选择层位用户界面。选择用户建立地层格架的层位，选择层位。
第六步：建立地层格架
如图8所示，利用选择的层位建立地层格架，选择地层格架的顶底层位。地层格架需考虑地层的上超、顶超等地质现象，缺省状态是平行顶底；对于上超选择平行于顶，底部削截。顶超选择平行于底，顶部削截。
第七步：建立初始模型
如图9所示，建立初始模型，初始模型是利用层位做约束，采用插值算法计算得到的。图9为建立初始模型界面，建立初始模型插值方法可以采用距离反加权插值方法或克里金插值方法。
对于克里金算法，各参数意义如下：
Auto Calculation Parameter按钮：按该按钮能根据井信息，采用最小二乘算法自动拟合变差函数参数。
Range1：横向变程1，变差函数长轴的长度。
Range2：横向变程2，变差函数横向短轴的长度。
Angle：变差函数长轴与X轴的夹角。
Max Point nums：最大邻近点数量。
Sill：基台值。
Nuget：块金常数。
Vertical Range：垂向变程。
Vertical Sill：垂向基台值。
第八步：神经网络训练
如图10所示，从列表Select Seismic Attributes中选择地震属性，并给定地震属性与井曲线的相关系数，相关系数通过井曲线与地震属性交会图获得。
Maximum training times:神经网络训练次数。
Network Name：神经网络名称。
如图11所示，为本发明神经网络结构图，与其他网络的主要不同之处在网络的输入端。本发明利用地震属性(attr1,attr2….attri....attrn)与井曲线(curve)的关系能得到由该地震属性重构的新曲线，本发明把重构的新曲线值作为高斯分布的均值，重构曲线与样本曲线的方差作为高斯分布的方差，这样由n个地震属性和样本曲线就能得到n个高斯分布(G1,G2…Gn)。这些高斯分布作为神经网络的输入层，从输入层、隐层到输出层的结构与PNN网络类似。经过样本点的网络训练确定输入层到隐层，隐层到输入层的权系数。
与PNN网络不同的是在输入端的高斯函数，不仅只受地震属性的影响，而且受到了井曲线的影响。这样在地震属性与井曲线关系不好的情况下，仍能得到稳定的反演结果。
结构图中：
其中Ci由地震属性attri重构获得，即根据attri与井曲线的相关关系由地震属性得到的新曲线。
网络训练结束后，利用网络模型对井曲线进行反演，网络所需要的地震属性，可以直接由地震属性体获得。而输入端所需要的“curve”信息，可以由井曲线经过Krig插值获得，插值过程中，用地震层位做约束，沿层分析变差函数。
选择地震属性后，进行神经网络训练。
第九步：输出训练结果
根据第八步的神经网络结果计算井曲线反演结果。如果第八步没有完成神经网络训练，这里将不能输出训练结果。
输出数据可以选择输出的数据格式，线道范围等等。
输出的数据体就是本发明利用地震属性做约束，采用RBF神经网络预测得到的井曲线反演数据体，对于井曲线反演，本发明输入的是裂缝密度曲线，地震属性利用地震几何属性做约束，输出的数据体就是裂缝密度数据体。
本发明不局限于上述最佳实施方式，任何人在本发明的启示下都可得出其他各种形式的产品，但不论在其形状或结构上作任何变化，凡是具有与本申请相同或相近似的技术方案，均落在本发明的保护范围之内。

Claims (2)

1.一种井曲线反演方法，其特征在于：所述方法包括
(1)选择工区的井和井曲线；
(2)井曲线时深转换；
(3)沿井轨迹提取地震属性；
(4)井曲线地震属性交互关系分析；
(5)选择层位；
(6)建立地层格架；
(7)建立初始模型；
(8)神经网络训练；
(9)输出训练结果；
所述步骤(1)包括选择曲线目录树中目标井曲线和该曲线的井或在编辑框中输入曲线名直接显示选择的对象；
所述步骤(2)包括采用数据管理单元将时深转换为时深表；时深转换后生成新的时间域的曲线；
所述步骤(3)包括用户界面中选择地震测网，该测网中的数据列表中，选择一个或多个地震属性，沿井轨迹提取选择的地震属性；
所述步骤(4)包括将在目录树中显示选择的井曲线，重采样经过时深转换后，得到的时间域曲线以及沿井轨迹提取的地震属性，进入交会图分析界面，选择X轴和Y轴的井曲线，设定每条曲线的左右刻度，绘制出交会图；在交会图界面中可选择不同的回归公式对交会图进行回归分析；
所述步骤(5)包括采用层位建立地层格架，在地层格架内进行井曲线反演；选择层位用户界面，选择用户建立地层格架的层位；
所述步骤(6)包括利用选择的层位建立地层格架，选择地层格架的顶底层位，对于上超选择平行于顶，底部削截；顶超选择平行于底，顶部削截；
所述步骤(7)包括建立初始模型，初始模型为利用层位做约束，采用插值算法计算；建立初始模型插值方法可采用距离反加权插值方法或克里金插值方法；
所述步骤(8)包括从列表中选择地震属性，并给定地震属性与井曲线的相关系数，相关系数通过井曲线与地震属性交会图获得，选择地震属性后，进行神经网络训练；利用地震属性与井曲线的相关关系得到由该地震属性重构的新曲线，将重构的新曲线值作为高斯分布的均值，重构曲线与样本曲线的方差作为高斯分布的方差，高斯分布作为神经网络的输入层，经过样本点的网络训练确定输入层到隐层，隐层到输入层的权系数。
2.根据权利要求1所述的一种井曲线反演方法，其特征在于：所述步骤(9)包括根据步骤(8)的神经网络训练结果计算井曲线反演结果；选择输出的数据格式和线道范围；输出的数据体为利用地震属性做约束，采用RBF神经网络预测得到的井曲线反演数据体，对于井曲线反演，输入的是裂缝密度曲线，地震属性利用地震几何属性做约束，输出的数据体为裂缝密度数据体。
CN201510413557.2A 2015-07-14 2015-07-14 一种井曲线反演方法 Active CN105005079B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510413557.2A CN105005079B (zh) 2015-07-14 2015-07-14 一种井曲线反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510413557.2A CN105005079B (zh) 2015-07-14 2015-07-14 一种井曲线反演方法

Publications (2)

Publication Number Publication Date
CN105005079A CN105005079A (zh) 2015-10-28
CN105005079B true CN105005079B (zh) 2017-12-01

Family

ID=54377813

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510413557.2A Active CN105005079B (zh) 2015-07-14 2015-07-14 一种井曲线反演方法

Country Status (1)

Country Link
CN (1) CN105005079B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106873033B (zh) * 2017-03-13 2019-03-22 中国海洋石油集团有限公司 一种基于神经网络的多属性地震信息融合裂缝预测的方法
CN108169817B (zh) * 2017-12-15 2019-11-15 中国石油集团川庆钻探工程有限公司 录井地震联合烃指数反演流体识别方法
CA3103234A1 (en) * 2018-06-10 2019-12-19 Schlumberger Canada Limited Seismic data interpretation system
CN109521470B (zh) * 2018-12-29 2020-04-14 西南石油大学 分析地质构造对地震反演裂缝密度影响的方法
CN110488350B (zh) * 2019-09-20 2021-10-29 西南石油大学 基于卷积神经网络的地震反演大数据生成方法
CN111308549B (zh) * 2020-03-09 2022-01-11 中国石油化工股份有限公司 基于模型反演的变速成图方法
CN111458748B (zh) * 2020-03-30 2021-03-05 青岛理工大学 基于三层数据集神经网络的性能地震动危险性分析方法
CN113589363B (zh) * 2020-04-30 2024-03-19 中国石油化工股份有限公司 融合人工神经网络与地质统计学的油气预测新方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1865341A2 (en) * 2006-06-09 2007-12-12 BHP Billiton Innovation Pty Ltd Method for improving prediction of the viability of potential petroleum reservoirs
CN104502997A (zh) * 2015-01-22 2015-04-08 中国石油化工集团公司 一种利用裂缝密度曲线预测裂缝密度体的方法
CN104614764A (zh) * 2015-01-28 2015-05-13 中国石油大学(华东) 一种基于扰动弹性阻抗反演的裂缝储层流体识别方法
CN104635265A (zh) * 2014-12-31 2015-05-20 中国石油天然气集团公司 一种可控神经网络式预测储层渗透率的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9383465B2 (en) * 2012-03-23 2016-07-05 Schlumberger Technology Corporation Quantitative analysis of time-lapse seismic data

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1865341A2 (en) * 2006-06-09 2007-12-12 BHP Billiton Innovation Pty Ltd Method for improving prediction of the viability of potential petroleum reservoirs
CN104635265A (zh) * 2014-12-31 2015-05-20 中国石油天然气集团公司 一种可控神经网络式预测储层渗透率的方法
CN104502997A (zh) * 2015-01-22 2015-04-08 中国石油化工集团公司 一种利用裂缝密度曲线预测裂缝密度体的方法
CN104614764A (zh) * 2015-01-28 2015-05-13 中国石油大学(华东) 一种基于扰动弹性阻抗反演的裂缝储层流体识别方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Emerge多属性分析理论培训;PST石油技术服务公司;《豆丁网》;20111118;第2、3、26-29、32-33页 *
地震多属性融合曲线重构技术在储层预测中的应用;马光克等;《物探与化探》;20131231;第37卷(第6期);第993-997页 *
地震联合反演技术在储层预测及烃类检测中的应用;赵萍萍;《中国石油和化工标准与质量》;20130420;第181页 *
致密碎屑岩储层地震反演技术方案及应用;刘振峰等;《石油地球物理勘探》;20120430;第47卷(第2期);第298-304,、352页 *

Also Published As

Publication number Publication date
CN105005079A (zh) 2015-10-28

Similar Documents

Publication Publication Date Title
CN105005079B (zh) 一种井曲线反演方法
US8838425B2 (en) Generating facies probablity cubes
Sullivan et al. Architectural analysis of deep-water outcrops: Implications for exploration and development of the Diana sub-basin, western Gulf of Mexico
CN105652323B (zh) 一种储层预测方法
CN104516018B (zh) 一种地球物理勘探中岩性约束下的孔隙度反演方法
CN105093313B (zh) 一种岩溶型油藏单井油气产能预测方法及装置
CN103454686A (zh) 基于地层切片的小尺度沉积相进行储层预测的方法及系统
CN106873033A (zh) 一种基于神经网络的多属性地震信息融合裂缝预测的方法
Sullivan et al. An integrated approach to characterization and modeling of deep-water reservoirs, Diana field, western Gulf of Mexico
CN110056346B (zh) 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法
Hosseini et al. Geostatistical modeling and spatial distribution analysis of porosity and permeability in the Shurijeh-B reservoir of Khangiran gas field in Iran
CN108333624A (zh) 一种基于地球物理信息的虚拟井构建方法
CN106970422A (zh) 一种在三类avo亮点特征储层区识别非亮点油藏的方法
Lee et al. 3D models of Quaternary-aged sedimentary successions within the Canterbury, New Zealand region
Suarez et al. Fracturing-to-Production Simulation Approach for Completion Optimization in the Vaca Muerta Shale
Ahmed et al. Cracking the code of unconventional reservoir simulation
Nguyen et al. Geological characteristics and integrated development plan for giant naturally fractured basement reservoirs
Anna et al. Integrated description and evaluation of reservoirs based on seismic, logging, and geological data: Taking Dongying Formation Member 1 oil reservoir of No. 1 structure, Nanpu Sag as an example
CN115880455A (zh) 基于深度学习的三维智能插值方法
CN115993657A (zh) 一种波阻抗反演方法
Abbaszadeh et al. Development of an integrated reservoir model for a naturally fractured volcanic reservoir in China
Nezhad et al. Simulation of petrophysical parameters of Asmari reservoir using SGS method in Mansuri oil field, southwest of Iran
Pandey et al. Geocellular Model for Tertiary Reservoirs in Manhera Tibba Gas Field, Jaisalmer Basin, Rajasthan, India
Yang et al. Comparison between several multi-parameter seismic inversion methods in identifying plutonic igneous rocks
Ehinola et al. Seismic attributes mapping and 3D static modeling of reservoirs within “OYA” field, offshore Depobelt, Niger delta sedimentary basin, Nigeria

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