CN108051860A - 一种利用趋势面形态生成误差校正面的方法 - Google Patents

一种利用趋势面形态生成误差校正面的方法 Download PDF

Info

Publication number
CN108051860A
CN108051860A CN201711211609.3A CN201711211609A CN108051860A CN 108051860 A CN108051860 A CN 108051860A CN 201711211609 A CN201711211609 A CN 201711211609A CN 108051860 A CN108051860 A CN 108051860A
Authority
CN
China
Prior art keywords
trend surface
well
trend
face
value
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.)
Granted
Application number
CN201711211609.3A
Other languages
English (en)
Other versions
CN108051860B (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.)
Beijing Gdf Oil And Gas Tech Inc
Original Assignee
Beijing Gdf Oil And Gas Tech Inc
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 Beijing Gdf Oil And Gas Tech Inc filed Critical Beijing Gdf Oil And Gas Tech Inc
Priority to CN201711211609.3A priority Critical patent/CN108051860B/zh
Publication of CN108051860A publication Critical patent/CN108051860A/zh
Application granted granted Critical
Publication of CN108051860B publication Critical patent/CN108051860B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Orthopedics, Nursing, And Contraception (AREA)
  • Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)

Abstract

本发明提出了一种利用趋势面形态生成误差校正面的方法，包括：获取每口井的井点校正值和地质趋势面；计算地质趋势面的趋势面数值，并建立所述井点校正值和趋势面数值之间对应关系；将每口井的井点校正值和趋势面数值做交会图，进行曲线拟合，拟合多项式趋势线以得到拟合多项式；根据拟合得到的多项式生成初始校正面；将所述初始校正面中的误差点进行常规网格化生成作为剩余误差，实现对剩余误差校正；将剩余误差校正结果加到初始校正面，生成最终校正面。本发明加入了地质趋势面做约束，没井控制的地方也可以有新的地质依据，因此生成的结果更加接近地质预期，符合生产者的实际需求。

Description

一种利用趋势面形态生成误差校正面的方法
技术领域
本发明涉及石油地质绘图技术领域，特别涉及一种利用趋势面形态生成误差校正面的方法。
背景技术
石油地质生产中，构造等值线图需要钻井数据校正。方法是用井点处的误差生成一个误差校正曲面，再用这个曲面来校正原构造图，得到和钻井数据相符的构造图。
目前主要的做法就是将井点直接网格化生成曲面，然后用于校正。使用该方法存在一些不合理的地方。如下图1所示，图左为误差点分布位置图，图右是主流方法得到的误差校正量曲面。在没有井数据控制的地方，由于存在多解性，无任何地质意义，光靠数学算法生成的校正面，达不到地质工作者的需求。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此，本发明的目的在于提出一种利用趋势面形态生成误差校正面的方法。
为了实现上述目的，本发明的实施例提供一种利用趋势面形态生成误差校正面的方法，包括如下步骤：
步骤S1，获取每口井的井点校正值Wi，得校正值集合{W1，W2，…，Wn}，其中，n表示有n口井和地质趋势面；
步骤S2，计算每口井在地质趋势面的趋势面数值Zi，得趋势值集合{Z1，Z2，…，Zn}，并建立所述井点校正值W和趋势面数值Z之间对应关系；
步骤S3，将每口井井点校正值W和趋势面数值Z做交会图，进行曲线拟合，拟合多项式趋势线以得到拟合多项式；
步骤S4，根据拟合得到的多项式生成初始校正面；
步骤S5，将所述初始校正面中的误差点进行常规网格化生成作为剩余误差，实现对剩余误差校正；
步骤S6，将所述步骤S5中的剩余误差校正结果加到步骤S4中的初始校正面，生成最终校正面。
进一步，在所述步骤S1中，根据井点构造图对应时间等值图数据，获取地质趋势面。
进一步，在所述步骤S2中，所述计算所述井点校正值和地质趋势面之间对应关系，包括：将每口井的井点坐标在所述地址趋势面相应的位置进行插值，得到趋势面数值Zn，将所述每口井的井点校正值Wn和趋势面数值Zn建立对应关系，形成数据对。
进一步，在所述步骤S3中，每口井点拟合得到的多项式为：
W＝Cn*Zm+Cn-1*Zm-1+…+C1*Z1+C0，
其中，W为校正量，Cn到C0为多项式常量系数，Z为趋势面的值，m表示是m次项。
进一步，每口井点拟合得到的多项式为分段函数，随Z的值域变化的W＝f(Z)曲线。
进一步，在所述步骤S4中，将地质趋势控制曲面，趋势值Z代入到对应的拟合函数W＝f(Z)，运算后生成对应的校正值W，并写入新的曲面，得到初始校正面。
根据本发明实施例的利用趋势面形态生成误差校正面的方法，引入地质趋势面，通过井的校正量值与井在趋势面上的取值和做数值拟合，得到拟合函数，再通过拟合函数系将趋势面换算成校正曲面。本发明通过趋势面约束将校正量数据生成有地质意义的平面网格的过程，并采用校正量与趋势面之间的数值拟合思路，由拟合关系将趋势面计算得到校正面的过程。本发明由于加入了地质趋势面做约束，没井控制的地方也可以有新的地质依据，因此生成的结果更加接近地质预期，符合生产者的实际需求。
本发明实施例的利用趋势面形态生成误差校正面的方法，具有以下有益效果：
1)传统方法生成的校正面，在无井约束的地方变形失真情况出现机率较大，而本方案不会造成这一现象；
2)基础数据获取容易，建立成果耗时少，效率高；
3)生成的校正面，有井控制的位置很准确，没有井控制的位置取值合理，效果好。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1为现有的原始误差点的示意图；
图2为常规网格化得到的校正曲图；
图3为根据本发明实施例的利用趋势面形态生成误差校正面的方法的流程图；
图4为根据本发明实施例的利用趋势面形态生成误差校正面的方法的示意图；
图5为根据本发明实施例的用于约束校正面的地址趋势面的示意图；
图6为根据本发明实施例的拟合曲线的示意图；
图7为根据本发明实施例的计算层速度界面图；
图8为根据本发明实施例的生成最终校正面的示意图；
图9a和图9b分别为常规方法和本发明方法生成的校正曲面的示意图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
如图3和图4所示，本发明实施例的利用趋势面形态生成误差校正面的方法，包括如下步骤：
步骤S1，获取每口井的井点校正值Wi，得校正值集合{W1，W2，…，Wn}和地质趋势面，其中，n表示有n口井，下同。
在本发明的一个实施例中，每口井的井点校正值Wi，如图1左图所示，井点校正值做的平面分布图。
并且，根据井点构造图对应时间等值图数据，获取地质趋势面，如图5所示。
步骤S2，计算每口井在地质趋势面的趋势面数值Zi，得趋势值集合{Z1，Z2，…，Zn}，并建立井点校正值W和趋势面数值Z之间对应关系。
在本发明的一个实施例中，计算井点校正值和地质趋势面之间对应关系，包括：将每口井的井点坐标在地址趋势面相应的位置进行插值，得到趋势面数值Zi，将每口井的井点校正值Wi和趋势面数值Zi建立对应关系，形成数据对，全部记录下来。
步骤S3，将每口井井点校正值W和趋势面数值Z做交会图，进行曲线拟合，拟合多项式趋势线以得到拟合多项式。
具体的，将所有井的校正量{W1，W2，…，Wn},趋势面数据{Z1，Z2，…，Zn}，做交会图，拟合多项式趋势线。即，用一条尽可以通过多数有效数据点的曲线，把这些点连接起来，得到误差W与数值Z的函数关系W＝f(Z)。
该函数关系一般为多项式W＝Cn*Zm+Cn-1*Zm-1+…+C1*Z1+C0(式中W为校正量，Cn到C0为多项式常量系数，Z为趋势面的值，m表示是m次项)。
W＝f(Z)也可以分段函数，随Z的值域变化的W＝f(Z)曲线。
图6为根据本发明实施例的拟合曲线的示意图。
步骤S4，根据拟合得到的多项式生成初始校正面。
在本步骤中，将地质趋势控制曲面，每个网格节点处的趋势值Z代入到对应的拟合函数W＝f(Z)，运算后生成对应的校正值W，并将所有计算好的值W写入新的曲面，得到初始校正面，如图7所示。
在本发明的一个实施例中，如果对生成的误差面效果不满意，可以调整步骤S3中的拟合曲线，重新生成。
步骤S5，将初始校正面中的误差点进行常规网格化生成作为剩余误差，实现对剩余误差校正。
具体的，因为拟合曲线是尽可能过所有的有效点得来的，所以生成的初始曲面在井点处还会有少量的剩余误差。由于这个误差已经很小，只需要将这些误差点做常规网格化生成即可作为剩余误差。
步骤S6，将步骤S5中的剩余误差校正结果加到步骤S4中的初始校正面，生成最终校正面。
在本步骤中，将上一步生成的剩余误差加到初始校正面上，就能得到最终校正曲面，此校正量曲面，满足在有井的位置，和校正量值相同，没有井的地方，和趋势背景相同。图8为根据本发明实施例的生成最终校正面的示意图。
如果对得到的结果效果不满意，可以返回步骤S3重新曲线拟合，直至得到满意的结果，流程到此结束。
参考图9a和图9b进行效果对比：图9a为常规方法网格化生成的校正曲面；图9b为本发明在地质趋势面约束后生成的校正曲面。通过对比可知，采用本发明提供的方法生成的校正曲面效果更好，结合了地质认识，更符合实际生产需要。
根据本发明实施例的利用趋势面形态生成误差校正面的方法，引入地质趋势面，通过井的校正量值与井在趋势面上的取值和做数值拟合，得到拟合函数，再通过拟合函数系将趋势面换算成校正曲面。本发明通过趋势面约束将校正量数据生成有地质意义的平面网格的过程，并采用校正量与趋势面之间的数值拟合思路，由拟合关系将趋势面计算得到校正面的过程。本发明由于加入了地质趋势面做约束，没井控制的地方也可以有新的地质依据，因此生成的结果更加接近地质预期，符合生产者的实际需求。
本发明实施例的利用趋势面形态生成误差校正面的方法，具有以下有益效果：
1)传统方法生成的校正面，在无井约束的地方变形失真情况出现机率较大，而本方案不会造成这一现象；
2)基础数据获取容易，建立成果耗时少，效率高；
3)生成的校正面，有井控制的位置很准确，没有井控制的位置取值合理，效果好。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。

Claims (6)

1.一种利用趋势面形态生成误差校正面的方法，其特征在于，包括如下步骤：
步骤S1，获取每口井的井点校正值Wi，得校正值集合{W1，W2，…，Wn}和地质趋势面，其中，n表示有n口井；
步骤S2，计算每口井在地质趋势面的趋势面数值Zi，得趋势值集合{Z1，Z2，…，Zn}，并建立所述井点校正值W和趋势面数值Z之间对应关系；
步骤S3，将每口井的井点校正值W和趋势面数值Z做交会图，进行曲线拟合，拟合多项式趋势线以得到拟合多项式；
步骤S4，根据拟合得到的多项式生成初始校正面；
步骤S5，将所述初始校正面中的误差点进行常规网格化生成作为剩余误差，实现对剩余误差校正；
步骤S6，将所述步骤S5中的剩余误差校正结果加到步骤S4中的初始校正面，生成最终校正面。
2.如权利要求1所述的利用趋势面形态生成误差校正面的方法，其特征在于，在所述步骤S1中，根据井点构造图对应时间等值图数据，获取地质趋势面。
3.如权利要求1所述的利用趋势面形态生成误差校正面的方法，其特征在于，在所述步骤S2中，所述计算所述井点校正值和地质趋势面之间对应关系，包括：将每口井的井点坐标在所述地址趋势面相应的位置进行插值，得到趋势面数值Zn，将所述每口井的井点校正值Wn和趋势面数值Zn建立对应关系，形成数据对。
4.如权利要求1所述的利用趋势面形态生成误差校正面的方法，其特征在于，在所述步骤S3中，每口井点拟合得到的多项式为：
W＝Cn*Zm+Cn-1*Zm-1+…+C1*Z1+C0，
其中，W为校正量，Cn到C0为多项式常量系数，Z为趋势面的值，m表示是m次项。
5.如权利要求1所述的利用趋势面形态生成误差校正面的方法，其特征在于，每口井点拟合得到的多项式为分段函数，随Z的值域变化的W＝f(Z)曲线。
6.如权利要求1所述的利用趋势面形态生成误差校正面的方法，其特征在于，在所述步骤S4中，将地质趋势控制曲面，趋势值Z代入到对应的拟合函数W＝f(Z)，运算后生成对应的校正值W，并写入新的曲面，得到初始校正面。
CN201711211609.3A 2017-11-28 2017-11-28 一种利用趋势面形态生成误差校正面的方法 Active CN108051860B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711211609.3A CN108051860B (zh) 2017-11-28 2017-11-28 一种利用趋势面形态生成误差校正面的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711211609.3A CN108051860B (zh) 2017-11-28 2017-11-28 一种利用趋势面形态生成误差校正面的方法

Publications (2)

Publication Number Publication Date
CN108051860A true CN108051860A (zh) 2018-05-18
CN108051860B CN108051860B (zh) 2021-08-03

Family

ID=62120727

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711211609.3A Active CN108051860B (zh) 2017-11-28 2017-11-28 一种利用趋势面形态生成误差校正面的方法

Country Status (1)

Country Link
CN (1) CN108051860B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110703347A (zh) * 2019-10-24 2020-01-17 中国石油化工股份有限公司 基于构造背景的重力断裂影像识别方法
CN112037306A (zh) * 2020-07-14 2020-12-04 中国石油天然气集团有限公司 砂体厚度图的绘制方法、装置、计算机设备及存储介质

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070021951A1 (en) * 2003-07-23 2007-01-25 Lee Seislink Wook B 3D velocity modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved for curved-ray prestack time migration and for such migration followed by prestack depth migration
CN102089677A (zh) * 2008-07-09 2011-06-08 雪佛龙美国公司 测井记录归一化
CN103226210A (zh) * 2013-04-01 2013-07-31 中国石油天然气股份有限公司 一种圈闭落实方法
CN103901478A (zh) * 2012-12-28 2014-07-02 中国石油天然气集团公司 一种井震信息联合确定储层沉积特征和分布的方法
CN104678440A (zh) * 2015-02-15 2015-06-03 山东科技大学 井约束二维地震变速度场非线性误差校正方法
CN105005077A (zh) * 2015-07-06 2015-10-28 成都理工大学 稀井条件下实钻井与虚拟井联合的薄层厚度预测方法
CN105093311A (zh) * 2015-06-29 2015-11-25 成都理工大学 叠合盆地多期差异隆升区地层剥蚀厚度测量方法
CN105604483A (zh) * 2015-12-29 2016-05-25 中国石油天然气股份有限公司 一种基于等时界面推演的水平井靶点设计方法
CN105717540A (zh) * 2016-03-14 2016-06-29 中国海洋石油总公司 一种微幅度构造精细预测方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070021951A1 (en) * 2003-07-23 2007-01-25 Lee Seislink Wook B 3D velocity modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved for curved-ray prestack time migration and for such migration followed by prestack depth migration
CN102089677A (zh) * 2008-07-09 2011-06-08 雪佛龙美国公司 测井记录归一化
CN103901478A (zh) * 2012-12-28 2014-07-02 中国石油天然气集团公司 一种井震信息联合确定储层沉积特征和分布的方法
CN103226210A (zh) * 2013-04-01 2013-07-31 中国石油天然气股份有限公司 一种圈闭落实方法
CN104678440A (zh) * 2015-02-15 2015-06-03 山东科技大学 井约束二维地震变速度场非线性误差校正方法
CN105093311A (zh) * 2015-06-29 2015-11-25 成都理工大学 叠合盆地多期差异隆升区地层剥蚀厚度测量方法
CN105005077A (zh) * 2015-07-06 2015-10-28 成都理工大学 稀井条件下实钻井与虚拟井联合的薄层厚度预测方法
CN105604483A (zh) * 2015-12-29 2016-05-25 中国石油天然气股份有限公司 一种基于等时界面推演的水平井靶点设计方法
CN105717540A (zh) * 2016-03-14 2016-06-29 中国海洋石油总公司 一种微幅度构造精细预测方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
T ARMSTRONG 等: "Removal of Overburden Velocity Anomaly Effects for Depth Conversion", 《GEOPHYSICAL PROSPECTING》 *
刘建芳 等: "相关分析构造图误差校正方法应用及效果", 《勘探地球物理进展》 *
满益志 等: "复杂逆掩构造及低幅度构造变速成图技术与应用", 《天然气工业》 *
王兴军 等: "低幅度构造变速成图技术", 《石油地球物理勘探》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110703347A (zh) * 2019-10-24 2020-01-17 中国石油化工股份有限公司 基于构造背景的重力断裂影像识别方法
CN110703347B (zh) * 2019-10-24 2021-04-16 中国石油化工股份有限公司 基于构造背景的重力断裂影像识别方法
CN112037306A (zh) * 2020-07-14 2020-12-04 中国石油天然气集团有限公司 砂体厚度图的绘制方法、装置、计算机设备及存储介质

Also Published As

Publication number Publication date
CN108051860B (zh) 2021-08-03

Similar Documents

Publication Publication Date Title
CN103939066B (zh) 一种一注多采井组定注水量确定油井产液量的方法
CN108051860A (zh) 一种利用趋势面形态生成误差校正面的方法
CN104265279B (zh) 断层条件下随钻测井曲线预测方法
CN104678440A (zh) 井约束二维地震变速度场非线性误差校正方法
CN109101745A (zh) 一种斜拉桥跨中挠度随温度变化的估算方法
CN107194104A (zh) 一种采用测井曲线计算页岩储层中无机质孔隙度、总孔隙度以及组分的方法
CN111367307A (zh) 一种用校正网络代替角加速度计的飞行器侧向过载跟踪方法
CN106802164A (zh) 一种适用于粮情监测的绝对水势计算方法及云图生成方法
CN103323031B (zh) 一种基于星敏的地平仪系统误差在线补偿的方法
CN105093280A (zh) 表层模型对地震数据影响的低频与高频成分的分解方法
CN113935253A (zh) 一种基于数据赋权的页岩气井经验产量递减模型拟合方法
CN112096374B (zh) 一种分流法测量持油率的动态测量误差补偿方法
CN109299402A (zh) 基于要素预分段的地址匹配方法
US20130013277A1 (en) Ghost Region Approaches for Solving Fluid Property Re-Distribution
CN105301657A (zh) 一种基于岩石物理意义的曲线校正方法
CN111209646B (zh) 一种致密低渗砂岩气藏多层合采井产气量劈分方法和装置
CN108074282B (zh) 基于三维地质建模软件的物性参数的厚度校正方法
CN116384267A (zh) 致密储层压裂水平井的最终可采储量的确定方法及设备
CN113219542B (zh) 一种基于改进的阻尼最小二乘法的频率域电磁反演方法
CN110244385B (zh) 一种大斜度井测井曲线的校正方法
CN110761780A (zh) 一种基于井震结合的三维地质导向方法
CN112922590B (zh) 基于异常产能试井资料进行指数校正评价气井产能的方法
CN111550239B (zh) 一种分段变参数异常井-模型数据耦合校正方法
CN113482596A (zh) 一种基于lwd随钻仪器实时三维伽马成像数据处理方法
CN108168508B (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
GR01 Patent grant
GR01 Patent grant