CN113919119A - 一种方位距离加权测井建模方法 - Google Patents

一种方位距离加权测井建模方法 Download PDF

Info

Publication number
CN113919119A
CN113919119A CN202010666198.2A CN202010666198A CN113919119A CN 113919119 A CN113919119 A CN 113919119A CN 202010666198 A CN202010666198 A CN 202010666198A CN 113919119 A CN113919119 A CN 113919119A
Authority
CN
China
Prior art keywords
logging
azimuth
distance
point
data
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
Application number
CN202010666198.2A
Other languages
English (en)
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.)
China Petroleum and Chemical Corp
China Petrochemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Geophysical Research Institute of Sinopec Shengli Oilfield Co
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 China Petroleum and Chemical Corp, Geophysical Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN202010666198.2A priority Critical patent/CN113919119A/zh
Publication of CN113919119A publication Critical patent/CN113919119A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • 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
    • G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48—Processing data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Geophysics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及数据建模技术领域，具体涉及一种方位距离加权测井建模方法所述方法包括：步骤1，测井坐标和数据加载；步骤2，测井方位距离加权系数求取：求取插值点到已知测井点的距离，计算反距离加权系数，计算综合非方位遮挡系数，进而可得测井方位距离加权系数；步骤3，利用方位距离加权数据插值公式进行模型构建。本发明方法针对测井数据空间相关性和方向性的特点，采用数据插值建模不仅考虑了距离因素，而且考虑了测井间的相互方位遮挡关系，能够提高测井数据建模精度，对提高地震反演约束能力和油藏描述能力具有重要作用。

Description

一种方位距离加权测井建模方法
技术领域
本发明涉及数据建模技术领域，具体涉及一种方位距离加权测井建模方法。
背景技术
通过数据插值建模能够利用采样点数据预测未知点数据，构建连续数据面，从而实现整个空间区域的数据分析应用。数据插值建模方法众多，像反距离加权法、克里金插值法、多元回归法、三角网法等，这些方法需要根据数据类型、分布特征及数据相关性等影响因素进行方法的选择。
其中反距离加权插值方法是一种常用的数据插值建模方法，适用于空间相关性强的数据建模。测井数据反映了地质沉积岩性空间变化特征，具有强空间相关性，因此反距离加权法在测井数值建模中得到广泛应用。但实际应用中发现，测井数据分布受地质构造和沉积物源方向影响，测井数据分布不仅具有强空间相关性，而且具有强方向性，距离预测点最近测井会对较远测井产生方位遮挡，由于反距离加权插值建模方法只考虑距离而没考虑这种方位遮挡关系，影响数据插值建模精度。
发明内容
本发明主要目的是提供一种方位距离加权测井建模方法，使得具有空间关性和方向性的数据插值建模精度更高。
为实现上述目的，本发明采用以下技术方案：
提供一种方位距离加权测井建模方法，所述方法包括：
步骤1，测井坐标和数据加载；
步骤2，测井方位距离加权系数求取：求取插值点到已知测井点的距离，计算反距离加权系数，计算综合非方位遮挡系数，进而可得测井方位距离加权系数；
步骤3，利用方位距离加权数据插值公式进行模型构建。
优选地，步骤2中方位距离加权系数为：
Figure BDA0002578620840000021
其中，βi为每个测井的综合非方位遮挡系数；di为插值点到已知测井点的距离，p和q称为加权幂指数。
优选地，在步骤2中，利用欧式距离公式，求取插值点到已知测井点的距离di：
Figure BDA0002578620840000022
反距离加权系数为
Figure BDA0002578620840000023
利用得到距离计算每口测井的反距离加权系数
Figure BDA0002578620840000024
该距离大小反映已知点对插值点影响的强弱，距离插值点越远的已知点影响越小。
优选地，综合非方位遮挡系数计算方法为：
求取已知测井点zk(xi,yi)与其它已知测井点zk(xj,yj)遮挡角度
Figure BDA0002578620840000025
比较遮挡角度与最大遮挡角度ω，同时比较两点与插值点距离di和dj的大小，已知点zk(xi,yi)被已知点zk(xj,yj)的两两非方位遮挡权系数计算公式为：
Figure BDA0002578620840000026
其中，
Figure BDA0002578620840000027
为遮挡角度，
Figure BDA0002578620840000028
ω为根据实际数据分析确定的最大遮挡角度，0°≤ω≤90°；i和j分别为测井编号。
通过对所有的已知测井点两两组合，分别计算两两测点间的非方位遮挡权系数，将所有的非方位遮挡权系数相乘，得每个测井的综合非方位遮挡权系数。
所述遮挡角度就是以插值点为顶点，以任意两个已知测井点zk(xi,yi)和zk(xj,yj)到插值点z(x,y)连线为边组成的夹角。当一个已知点比另一个已知点更靠近插值点，并且遮挡角度小时，就会产生遮挡。
进一步优选地，测井综合非方位遮挡系数计算公式为：
Figure BDA0002578620840000031
优选地，方位距离数据插值公式为：
Figure BDA0002578620840000032
其中，βi为每个测井的综合非方位遮挡系数；di为插值点到已知测井点的距离，p和q称为加权幂指数，zi为已知测井数据。
进一步优选地，当插值点为已知点时用已知点数据代替插值点，当插值点不是任何已知点时，利用公式(4)对未知点进行插值。利用该方位数据插值公式，完成所有未知点的数据插值。
优选地，所述数据为测井数据。
与现有技术相比，本发明具有以下优势：
本发明方法针对测井数据空间相关性和方向性的特点，采用数据插值建模不仅考虑了距离因素，而且考虑了测井间的相互方位遮挡关系，能够提高测井数据建模精度，对提高地震反演约束能力和油藏描述能力具有重要作用。
附图说明
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
图1为本发明一种方位距离加权测井建模方法的流程图；
图2为本发明一个实施例中的方位距离加权系数计算示意图；
图3为采用常规反距离加权测井建模方法得到测井属性平面图；
图4为采用本发明方位距离加权测井建模方法得到测井属性平面图。
具体实施方式
应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作和/或它们的组合。
为了使得本领域技术人员能够更加清楚地了解本发明的技术方案，以下将结合具体的实施例详细说明本发明的技术方案。
实施例1
如图1所示，图1为本发明所述一种方位距离加权测井建模方法流程图。所述方法包括：
步骤1，测井坐标和数据加载需要记录测井的个数，用n表示。测井坐标为(xi,yi)，其中i＝1,2…,n。测井数据表示为zi；
步骤2，测井方位距离加权系数求取：求取插值点到已知测井点的距离，计算反距离加权系数，计算综合非方位遮挡系数，进而可得测井方位距离加权系数；
利用欧式距离定义公式(1)，求取插值点到已知测井点的距离di：
Figure BDA0002578620840000051
利用得到距离计算每口测井的反距离加权系数
Figure BDA0002578620840000052
求取已知测井点zk(xi,yi)与其它已知测井点zk(xj,yj)遮挡角度，也就是以插值点为顶点，以任意两个已知测井点zk(xi,yi)和zk(xj,yj)到插值点z(x,y)连线为边组成的夹角，表示为
Figure BDA0002578620840000053
比较遮挡角度与最大遮挡角度ω，同时比较两者与插值点距离di和dj的大小，当一个已知点比另一个已知点更靠近插值点，并且遮挡角度小时，就会产生遮挡，可用非方位遮挡权系数来描述这种遮挡关系，已知点zk(xi,yi)被已知点zk(xj,yj)的两两非方位遮挡权系数计算公式为(2)：
Figure BDA0002578620840000054
Figure BDA0002578620840000055
为遮挡角度，
Figure BDA0002578620840000056
ω为根据实际数据分析确定的最大遮挡角度，0°≤ω≤90°。
通过对所有的已知测井点两两组合，就可以得到相互的间的非方位遮挡权系数，然后对所有的非方位遮挡权系数相乘，就可以得到每个测井已知点对插值点的综合非方位遮挡权系数，每个测井的综合非方位遮挡系数计算公式表示为：
Figure BDA0002578620840000057
通过上步骤得到每个已知测井点对插值点的反距离加权系数
Figure BDA0002578620840000058
和综合非方位遮挡系数βi，进而可以得到测井方位距离加权系数
Figure BDA0002578620840000059
其中，p和q称为加权幂指数。
步骤3，利用方位距离加权数据插值公式进行模型构建。在得到测井点方位距离加权系数的基础上，构建方位距离加权数据插值公式(4)：
Figure BDA0002578620840000061
当插值点为已知点时用已知点数据代替插值点，当插值点不是任何已知点时，利用公式(4)对未知点进行插值。利用该方位数据插值公式，完成所有未知点的数据插值。
图2为采用本发明实施例中方位距离加权系数计算示意图，具体示意了测井1被测井2方位遮挡下的距离方位加权系数求取。示意图中显示了已经加载的60口测井的坐标分布。图2中插值点到其中两个测井点1和测井点2的距离d1和d2，以及夹角
Figure BDA0002578620840000062
然后利用公式(2)求取测井1被测井2遮挡而具有的非方位遮挡系数β12，将测井1与其他测井两两组合，就可以得到测井1被其他测井非方位遮挡系数列β12,β13…β1i，然后公式(3)求取测井1的综合非方位遮挡系数β1。在得到测井1到插值点距离d1和综合非方位遮挡系数β1后，就可以得到测井1的加权系数
Figure BDA0002578620840000063
这里q和p分别取1，测井1的距离方位加权系数为
Figure BDA0002578620840000064
利用图2中所示方位距离加权系数计算方法，依次求取所有测井的距离方位加权系数
Figure BDA0002578620840000065
同样q和p分别取1。
分别利用原有的反距离加权方法和本发明方法提供的数据插值公式进行了数据插值建模。其中反距离加权公式为：
Figure BDA0002578620840000071
图3为利用反距离权系数插值建模公式(5)，其中p取1，进行了数据插值建模后提取的属性平面图。图4为利用方位距离加权数据插值公式(4)，进行数据插值建模后提取的属性平面。可以看出本发明能够较好的考虑测井数据间的的相关性和方向性，消除已知点附近局部异常效应，测井插值模型更能体现地质构造沉积自然变化。
上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (8)

1.一种方位距离加权测井建模方法，其特征在于，所述方法包括：
步骤1，测井坐标和数据加载；
步骤2，测井方位距离加权系数求取：求取插值点到已知测井点的距离，计算反距离加权系数，计算综合非方位遮挡系数，进而可得测井方位距离加权系数；
步骤3，利用方位距离加权数据插值公式进行模型构建。
2.根据权利要求1所述方位距离加权测井建模方法，其特征在于，步骤2中方位距离加权系数为：
Figure FDA0002578620830000011
其中，βi为每个测井的综合非方位遮挡系数；di为插值点到已知测井点的距离，p和q称为加权幂指数。
3.根据权利要求1或2所述方位距离加权测井建模方法，其特征在于，在步骤2中，利用欧式距离公式，求取插值点到已知测井点的距离di：
Figure FDA0002578620830000012
反距离加权系数为
Figure FDA0002578620830000013
4.根据权利要求1或2所述方位距离加权测井建模方法，其特征在于，综合非方位遮挡系数计算方法为：
求取已知测井点zk(xi,yi)与其它已知测井点zk(xj,yj)遮挡角度
Figure FDA0002578620830000015
比较遮挡角度与最大遮挡角度ω，同时比较两点与插值点距离di和dj的大小，已知点zk(xi,yi)被已知点zk(xj,yj)的两两非方位遮挡权系数计算公式为：
Figure FDA0002578620830000014
其中，
Figure FDA0002578620830000023
为遮挡角度，
Figure FDA0002578620830000024
ω为根据实际数据分析确定的最大遮挡角度，0°≤ω≤90°；
通过对所有的已知测井点两两组合，分别计算两两测点间的非方位遮挡权系数，将所有的非方位遮挡权系数相乘，得每个测井的综合非方位遮挡权系数。
5.根据权利要求4所述方位距离加权测井建模方法，其特征在于，测井综合非方位遮挡系数计算公式为：
Figure FDA0002578620830000021
6.根据权利要求1所述方位距离加权测井建模方法，其特征在于，方位距离数据插值公式为：
Figure FDA0002578620830000022
其中，βi为每个测井的综合非方位遮挡系数；di为插值点到已知测井点的距离，p和q称为加权幂指数，zi为已知测井数据。
7.根据权利要求6所述方位距离加权测井建模方法，其特征在于，当插值点为已知点时用已知点数据代替插值点，当插值点不是任何已知点时，利用公式(4)对未知点进行插值。
8.根据权利要求1所述方位距离加权测井建模方法，其特征在于，所述数据为测井数据。
CN202010666198.2A 2020-07-10 2020-07-10 一种方位距离加权测井建模方法 Pending CN113919119A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010666198.2A CN113919119A (zh) 2020-07-10 2020-07-10 一种方位距离加权测井建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010666198.2A CN113919119A (zh) 2020-07-10 2020-07-10 一种方位距离加权测井建模方法

Publications (1)

Publication Number Publication Date
CN113919119A true CN113919119A (zh) 2022-01-11

Family

ID=79231408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010666198.2A Pending CN113919119A (zh) 2020-07-10 2020-07-10 一种方位距离加权测井建模方法

Country Status (1)

Country Link
CN (1) CN113919119A (zh)

Similar Documents

Publication Publication Date Title
DK2356493T3 (en) PROCEDURE FOR GEOLOGICAL MODELING OF SEISMIC DATA BY TRACK CORRELATION
CN108318921B (zh) 一种基于横向约束的快速地震随机反演方法
US7117093B2 (en) Method, media, and signals for processing seismic data to obtain a velocity field
US10795039B2 (en) Generating pseudo pressure wavefields utilizing a warping attribute
CN110031898B (zh) 数据优化方法及积分法叠前深度偏移方法
CN113093272A (zh) 基于卷积编码的时间域全波形反演方法
WO2017136133A1 (en) Efficient seismic attribute gather generation with data synthesis and expectation method
WO2021127382A1 (en) Full waveform inversion in the midpoint-offset domain
CN114609669B (zh) 基于方位弹性阻抗的hti型裂缝储层参数预测方法及系统
MXPA04009215A (es) Metodo para migracion sismica que utiliza operadores de extrapolacion de profundidad explicita con longitud de operador dinamicamente variable.
CN112099082B (zh) 一种共面元共方位角道集的地震回折波走时反演方法
CN113919119A (zh) 一种方位距离加权测井建模方法
CN109188522B (zh) 速度场构建方法及装置
CN109387872B (zh) 表面多次波预测方法
Xie et al. 3D wavefront attribute determination and conflicting dip processing
US12013508B2 (en) Method and system for determining seismic processing parameters using machine learning
US12000971B2 (en) Method and system for seismic processing using virtual trace bins based on offset attributes and azimuthal attributes
CN111538080B (zh) 地震成像的方法
US10082588B2 (en) Adaptive structure-oriented operator
Waldeland et al. Fast estimation of prestack common reflection surface parameters
US10067254B2 (en) Removal of an estimated acquisition effect from a marine survey measurement
CN112083490A (zh) 地震数据噪声衰减方法及装置
US20230266492A1 (en) Enhancement of single sensor seismic data
CN114002741B (zh) 叠前深度偏移方法及装置、计算机可读存储介质
CN113126163B (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