CN107807386A - 一种预测点坝砂体位置的方法及装置 - Google Patents

一种预测点坝砂体位置的方法及装置 Download PDF

Info

Publication number
CN107807386A
CN107807386A CN201710889005.8A CN201710889005A CN107807386A CN 107807386 A CN107807386 A CN 107807386A CN 201710889005 A CN201710889005 A CN 201710889005A CN 107807386 A CN107807386 A CN 107807386A
Authority
CN
China
Prior art keywords
mrow
sampled point
variation
full
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
CN201710889005.8A
Other languages
English (en)
Other versions
CN107807386B (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201710889005.8A priority Critical patent/CN107807386B/zh
Publication of CN107807386A publication Critical patent/CN107807386A/zh
Application granted granted Critical
Publication of CN107807386B publication Critical patent/CN107807386B/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/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616—Data from specific type of measurement
    • G01V2210/6161—Seismic or acoustic, e.g. land or sea measurements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
    • G01V2210/632—Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种预测点坝砂体位置的方法及装置，该方法包括：获取曲流河的振幅属性平面图，在振幅属性平面图上确定振幅属性值大于预设阈值的第一区域；根据第一区域内每一个采样点的振幅属性值，计算每一个采样点的变异系数值，并根据所有采样点的变异系数值，得到振幅属性变异系数平面图；在振幅属性变异系数平面图上确定曲流河的河道砂体的边界；根据第一区域内每一个采样点的倾角大小和方位角大小，计算每一个采样点的全角变化梯度值，并得到全角变化梯度平面图；在全角变化梯度平面图上确定曲流河的迁移方向；根据曲流河的河道砂体的边界及曲流河的迁移方向，预测点坝砂体的位置。本发明可比较精确的预测点坝砂体的位置。

Description

一种预测点坝砂体位置的方法及装置
技术领域
本发明涉及油气开采技术领域，特别涉及一种预测点坝砂体位置的方法及装置。
背景技术
点坝砂体区域是曲流河主要的储油区，确定曲流河的点坝砂体区域对于后续开采具有非常重要的指导意义。曲流河河道迁移摆动频繁，点坝砂体复杂展布多变，因此，提供一种方法预测曲流河的点坝砂体的位置十分必要。
现有技术在通常采用地震属性分析法来确定点坝砂体的位置，将振幅平面图中地震波的振幅值较高的区域对应的曲流河的位置确定为点坝砂体的位置。
在实现本发明的过程中，本发明人发现现有技术中至少存在以下问题：
现有技术的方法确定的是曲流河当前的点坝砂体的位置，随着曲流河河道迁移摆动，点坝砂体的位置发生变化，现有技术的方法并不能精确预测点坝砂体的位置。
发明内容
有鉴于此，本发明提供一种预测点坝砂体位置的方法及装置，以比较精确的预测点坝砂体的位置。
具体而言，包括以下的技术方案：
第一方面，本发明提供了一种预测点坝砂体位置的方法，包括：
获取曲流河的振幅属性平面图，在所述振幅属性平面图上确定振幅属性值大于预设阈值的第一区域；
根据所述第一区域内多个采样点中每一个采样点的振幅属性值，计算每一个所述采样点的变异系数值，并根据所述多个采样点的变异系数值，得到振幅属性变异系数平面图；
在所述振幅属性变异系数平面图上，根据所述多个采样点的位置及变异系数值，确定所述曲流河的河道砂体的边界；
根据所述第一区域内每一个所述采样点的倾角大小和方位角大小，计算每一个所述采样点的全角变化梯度值，并根据所述多个采样点的全角变化梯度值，得到全角变化梯度平面图；
在所述全角变化梯度平面图上，根据所述多个采样点的位置及全角变化梯度值，确定所述曲流河的迁移方向；
根据所述曲流河的河道砂体的边界及所述曲流河的迁移方向，预测所述点坝砂体的位置。
可选择地，所述在所述振幅属性变异系数平面图上，根据所述多个采样点的位置及变异系数值，确定所述曲流河的河道砂体的边界，包括：
在所述振幅属性变异系数平面图上，将变异系数大于第一预设阈值的所述采样点所在的第二区域确定为所述曲流河的河道砂体的边界。
可选择地，所述在所述全角变化梯度图上，根据所述多个采样点的位置及全角变化梯度值，确定所述曲流河的迁移方向，包括：
在所述全角变化梯度平面图上，全角变化梯度大于第二预设阈值的所述采样点所在的区域为所述曲流河的凹岸，全角变化梯度小于第三预设阈值的所述采样点所在的区域为所述曲流河的凸岸；
所述曲流河的迁移方向由所述凸岸指向所述凹岸。
可选择地，所述采样点的变异系数的计算公式为：
其中，
Vi——采样点i的变异系数；
Ai——第i个样品的振幅属性值；
——第一区域内所有采样点的振幅属性平均值；
n——采样范围内采样点的个数。
可选择地，所述采样点的全角变化梯度的计算公式为：
其中，
Kab——A、B两点的全角变化梯度；
φa、φb——A、B两点的倾角；
θa、θb——A、B两点的方位角。
第二方面，本发明还提供了一种预测点坝砂体位置的装置，包括：
第一确定模块，用于获取曲流河的振幅属性平面图，在所述振幅属性平面图上确定振幅属性值大于预设阈值的第一区域；
第一获取模块，用于根据所述第一区域内多个采样点中每一个采样点的振幅属性值，计算每一个所述采样点的变异系数值，并根据所述多个采样点的变异系数值，得到振幅属性变异系数平面图；
第二确定模块，用于在所述振幅属性变异系数平面图上，根据所述多个采样点的位置及变异系数值，确定所述曲流河的河道砂体的边界；
第二获取模块，根据所述第一区域内每一个所述采样点的倾角大小和方位角大小，计算每一个所述采样点的全角变化梯度值，并根据所述多个采样点的全角变化梯度值，得到全角变化梯度平面图；
第三确定模块，用于在所述全角变化梯度平面图上，根据所述多个采样点的位置及全角变化梯度值，确定所述曲流河的迁移方向；
预测模块，用于根据所述曲流河的河道砂体的边界及所述曲流河的迁移方向，预测所述点坝砂体的位置。
可选择地，所述第二确定模块具体用于：
在所述振幅属性变异系数平面图上，将变异系数大于第一预设阈值的所述采样点所在的第二区域确定为所述曲流河的河道砂体的边界。
可选择地，所述第三确定模块具体用于：
在所述全角变化梯度平面图上，全角变化梯度大于第二预设阈值的所述采样点所在的区域为所述曲流河的凹岸，全角变化梯度小于第三预设阈值的所述采样点所在的区域为所述曲流河的凸岸；
所述曲流河的迁移方向由所述凸岸指向所述凹岸。
可选择地，所述采样点的变异系数的计算公式为：
其中，
Vi——采样点i的变异系数；
Ai——第i个样品的振幅属性值；
——第一区域内所有采样点的振幅属性平均值；
n——采样范围内采样点的个数。
可选择地，所述采样点的全角变化梯度的计算公式为：
其中，
Kab——A、B两点的全角变化梯度；
φa、φb——A、B两点的倾角；
θa、θb——A、B两点的方位角。
本发明实施例提供的技术方案的有益效果：
本发明实施例提供了一种预测点坝砂体位置的方法及装置，首先在振幅属性平面图上确定振幅属性值大于预设阈值的第一区域，并计算第一区域内所有采样点的变异系数值，得到振幅属性变异系数平面图，进而在振幅属性变异系数平面图上确定曲流河的河道砂体的边界；计算每一个采样点的全角变化梯度值，得到全角变化梯度平面图，确定曲流河的河道的迁移方向；最后根据曲流河河道砂体的边界及曲流河的河道的迁移方向，确定点坝砂体的位置。可知，本发明实施例通过确定曲流河的河道砂体的边界及曲流河的河道侧向摆动，从而预测点坝砂体的位置，预测结果比较精确。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例中一种预测点坝砂体位置的方法的流程图；
图2是本发明实施例中点坝砂体位置的示意图；
图3是本发明实施例中一种预测点坝砂体位置的装置的框图。
具体实施方式
为使本发明的技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
本发明实施例提供了一种预测点坝砂体位置的方法，如图1所示，包括步骤101、102、103、104、105和106。下面将进行具体说明。
101、获取曲流河的振幅属性平面图，在振幅属性平面图上确定振幅属性值大于预设阈值的第一区域。
在振幅属性平面图上，点坝砂体所在区域的振幅属性值较高。当曲流河附近的点坝砂体分布发生变化后，曲流河所在区域的振幅属性平面图就会发生变化。因此可在振幅属性平面图上确定振幅属性值大于预设阈值的第一区域，初步确定点坝砂体的分布位置。
其中，预设阈值根据曲流河的实际情况确定，一般为6000～8000。
102、根据第一区域内多个采样点中每一个采样点的振幅属性值，计算每一个采样点的变异系数值，并根据多个采样点的变异系数值，得到振幅属性变异系数平面图。
在计算采样点的振幅变异系数时，首先以该采样点为中心确定该采样点的采样范围，根据该采样点及该采样点的采样范围内的所有采样点的振幅属性值，求取该采样点的变异系数。采样范围的选取是计算变异系数的关键，一般采样范围的半径不大于河道宽度的一半，不小于75m，并且采样范围内的采样点的数量不少于9个。
采样点的变异系数的计算公式为：
其中，
Vi——采样点i的变异系数；
Ai——第i个样品的振幅属性值；
——第一区域内所有采样点的振幅属性平均值；
n——采样范围内采样点的个数。
103、在振幅属性变异系数平面图上，根据所述多个采样点的位置及变异系数值，确定曲流河的河道砂体的边界。
具体地，在振幅属性变异系数平面图上，将变异系数大于第一预设阈值的采样点所在的第二区域确定为曲流河的河道砂体的边界。其中，第一预设阈值可为0.6。在振幅属性变异系数平面图上，采样点的变异系数小于0.4时，反映采样点的振幅属性值变化程度较小，变异系数小于0.4的采样点所在的区域为砂泥岩内部区域；采样点的变异系数在0.4～0.6之间时，反映采样点的振幅属性值变化程度中等，变异系数在0.4～0.6之间的采样点所在的区域为砂泥岩变化区域；采样点的变异系数大于0.6时，反映该采样点的振幅属性值变化程度较大，变异系数大于0.6的采样点所在的区域为砂泥岩区域的边界，也即为曲流河的河道砂体的边界。
104、根据第一区域内每一个采样点的倾角大小和方位角大小，计算每一个采样点的全角变化梯度值，并根据多个采样点的全角变化梯度值，得到全角变化梯度平面图。
地震全角变化梯度主要是在三维空间内衡量地震波形变化情况，指示地震波的不连续性。可根据地震全角变化梯度预测河道迁移方向。
具体地，采样点的全角变化梯度的计算公式为：
其中，
Kab——A、B两点的全角变化梯度；
φa、φb——A、B两点的倾角；
θa、θb——A、B两点的方位角。
A、B两点为相邻最近的两个采样点。
其中，倾角为河道表面与水平面所成的角，方位角为目标地震道在水平面的投影线与正北方向所成的角。根据河道底界面的倾角属性平面图和方位角属性平面图，可得到所有采样点的倾角大小和方位角大小，进而计算出每一个采样点的全角变化梯度。当曲流河的河道摆动时，倾角和方位角也会发生变化。
105、在全角变化梯度平面图上，根据所述多个采样点的位置及全角变化梯度值，确定曲流河的迁移方向。
在全角变化梯度平面图上，全角变化梯度大于第二预设阈值的区域被认为是曲流河的凹岸，全角变化梯度小于第三预设阈值的区域被认为是曲流河的凸岸，全角变化梯度在第二预设阈值与第三预设阈值之间的区域被认为是曲流河的河道。曲流河的河道的迁移方向由凸岸指向凹岸。其中，第二预设阈值可为40，第三预设阈值可为25。
106、根据曲流河的河道砂体的边界以及曲流河的迁移方向，预测点坝砂体的位置。
具体地，可对全角变化梯度平面图进行处理，将全角变化梯度小于40的区域去掉，则处理后的全角变化梯度平面图很好的反映了曲流河的河道的迁移方向，也即河道砂体的侧向迁移方向。将处理后的全角变化梯度平面图与曲流河河道砂体的边界相结合，即可确定点坝砂体的具体位置，从而精确地预测曲流河的储层。如图2所示，可确定点坝砂体发育于曲流河的凸岸，向曲流河的凹岸迁移。
本发明实施例提供的预测点坝砂体位置的方法，首先在振幅属性平面图上确定振幅属性值大于预设阈值的第一区域，并计算第一区域内所有采样点的变异系数值，得到振幅属性变异系数平面图，进而在振幅属性变异系数平面图上确定曲流河河道砂体的边界；计算每一个采样点的全角变化梯度值，得到全角变化梯度平面图，确定曲流河的河道的迁移方向；最后根据曲流河的河道砂体的边界及曲流河的河道的迁移方向，确定点坝砂体的位置。可知，本发明实施例通过确定曲流河河道砂体的边界及曲流河的河道侧向摆动，从而预测点坝砂体的位置，预测结果比较精确。
本发明实施例还提供了一种预测点坝砂体位置的装置，如图3所示，包括第一确定模块301、第一获取模块302、第二确定模块303、第二获取模块304、第三确定模块305和预测模块306。下面将对各模块进行具体介绍。
第一确定模块301，用于获取曲流河的振幅属性平面图，在振幅属性平面图上确定振幅属性值大于预设阈值的第一区域；
第一获取模块302，用于根据第一区域内每一个采样点的振幅属性值，计算每一个采样点的变异系数值，并根据多个采样点的变异系数值，得到振幅属性变异系数平面图；
第二确定模块303，用于在振幅属性变异系数平面图上，根据所述多个采样点的位置及变异系数值，确定曲流河的河道砂体的边界；
第二获取模块304，根据第一区域内多个采样点中的每一个采样点的倾角大小和方位角大小，计算每一个采样点的全角变化梯度值，并根据多个采样点的全角变化梯度值，得到全角变化梯度平面图；
第三确定模块305，用于在全角变化梯度平面图上，根据所述多个采样点的位置及全角变化梯度值，确定曲流河的迁移方向；
预测模块306，用于根据曲流河河道砂体的边界及曲流河的迁移方向，预测点坝砂体的位置。
作为一种可能的实现方式，第二确定模块303具体用于：
在振幅属性变异系数平面图上，将变异系数大于第一预设阈值的采样点所在的第二区域确定为曲流河的河道砂体的边界。
作为一种可能的实现方式，第三确定模块305具体用于：
在全角变化梯度平面图上，全角变化梯度大于第二预设阈值的采样点所在的区域为曲流河的凹岸，全角变化梯度小于第三预设阈值的采样点所在的区域为曲流河的凸岸；
曲流河的迁移方向由凸岸指向凹岸。
作为一种可能的实现方式，采样点的变异系数的计算公式为：
其中，
Vi——采样点i的变异系数；
Ai——第i个样品的振幅属性值；
——第一区域内所有采样点的振幅属性平均值；
n——采样范围内采样点的个数。
作为一种可能的实现方式，采样点的全角变化梯度的计算公式为：
其中，
Kab——A、B两点的全角变化梯度；
φa、φb——A、B两点的倾角；
θa、θb——A、B两点的方位角。
本发明实施例提供的预测点坝砂体位置的装置，首先在振幅属性平面图上确定振幅属性值大于预设阈值的第一区域，并计算第一区域内所有采样点的变异系数值，得到振幅属性变异系数平面图，进而在振幅属性变异系数平面图上确定曲流河的河道砂体的边界；计算每一个采样点的全角变化梯度值，得到全角变化梯度平面图，确定曲流河的河道的迁移方向；最后根据曲流河的河道砂体的边界及曲流河的河道的迁移方向，确定点坝砂体的位置。可知，本发明实施例通过确定曲流河河道砂体的边界及曲流河的河道侧向摆动，从而预测点坝砂体的位置，预测结果比较精确。
以上所述仅是为了便于本领域的技术人员理解本发明的技术方案，并不用以限制本发明。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种预测点坝砂体位置的方法，其特征在于，包括：
获取曲流河的振幅属性平面图，在所述振幅属性平面图上确定振幅属性值大于预设阈值的第一区域；
根据所述第一区域内多个采样点中每一个采样点的振幅属性值，计算每一个所述采样点的变异系数值，并根据所述多个采样点的变异系数值，得到振幅属性变异系数平面图；
在所述振幅属性变异系数平面图上，根据所述多个采样点的位置及变异系数值，确定所述曲流河的河道砂体的边界；
根据所述第一区域内每一个所述采样点的倾角大小和方位角大小，计算每一个所述采样点的全角变化梯度值，并根据所述多个采样点的全角变化梯度值，得到全角变化梯度平面图；
在所述全角变化梯度平面图上，根据所述多个采样点的位置及全角变化梯度值，确定所述曲流河的迁移方向；
根据所述曲流河的河道砂体的边界及所述曲流河的迁移方向，预测所述点坝砂体的位置。
2.根据权利要求1所述的预测点坝砂体位置的方法，其特征在于，所述在所述振幅属性变异系数平面图上，根据所述多个采样点的位置及变异系数值，确定所述曲流河的河道砂体的边界，包括：
在所述振幅属性变异系数平面图上，将变异系数大于第一预设阈值的所述采样点所在的第二区域确定为所述曲流河的河道砂体的边界。
3.根据权利要求1所述的预测点坝砂体位置的方法，其特征在于，所述在所述全角变化梯度图上，根据所述多个采样点的位置及全角变化梯度值，确定所述曲流河的迁移方向，包括：
在所述全角变化梯度平面图上，全角变化梯度大于第二预设阈值的所述采样点所在的区域为所述曲流河的凹岸，全角变化梯度小于第三预设阈值的所述采样点所在的区域为所述曲流河的凸岸；
所述曲流河的迁移方向由所述凸岸指向所述凹岸。
4.根据权利要求1所述的预测点坝砂体位置的方法，其特征在于，所述采样点的变异系数的计算公式为：
<mrow> <msub> <mi>V</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <msqrt> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>/</mo> <mi>n</mi> </mrow> </msqrt> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> </mfrac> </mrow>
其中，
Vi——采样点i的变异系数；
Ai——第i个样品的振幅属性值；
——第一区域内所有采样点的振幅属性平均值；
n——采样范围内采样点的个数。
5.根据权利要求1所述的预测点坝砂体位置的方法，其特征在于，所述采样点的全角变化梯度的计算公式为：
<mrow> <msub> <mi>K</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mi>a</mi> </msub> <mo>-</mo> <msub> <mi>&amp;phi;</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>a</mi> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow>
其中，
Kab——A、B两点的全角变化梯度；
φa、φb——A、B两点的倾角；
θa、θb——A、B两点的方位角。
6.一种预测点坝砂体位置的装置，其特征在于，包括：
第一确定模块，用于获取曲流河的振幅属性平面图，在所述振幅属性平面图上确定振幅属性值大于预设阈值的第一区域；
第一获取模块，用于根据所述第一区域内多个采样点中每一个采样点的振幅属性值，计算每一个所述采样点的变异系数值，并根据所述多个采样点的变异系数值，得到振幅属性变异系数平面图；
第二确定模块，用于在所述振幅属性变异系数平面图上，根据所述多个采样点的位置及变异系数值，确定所述曲流河的河道砂体的边界；
第二获取模块，根据所述第一区域内每一个所述采样点的倾角大小和方位角大小，计算每一个所述采样点的全角变化梯度值，并根据所述多个采样点的全角变化梯度值，得到全角变化梯度平面图；
第三确定模块，用于在所述全角变化梯度平面图上，根据所述多个采样点的位置及全角变化梯度值，确定所述曲流河的迁移方向；
预测模块，用于根据所述曲流河的河道砂体的边界及所述曲流河的迁移方向，预测所述点坝砂体的位置。
7.根据权利要求1所述的预测点坝砂体位置的装置，其特征在于，所述第二确定模块具体用于：
在所述振幅属性变异系数平面图上，将变异系数大于第一预设阈值的所述采样点所在的第二区域确定为所述曲流河的河道砂体的边界。
8.根据权利要求1所述的预测点坝砂体位置的装置，其特征在于，所述第三确定模块具体用于：
在所述全角变化梯度平面图上，全角变化梯度大于第二预设阈值的所述采样点所在的区域为所述曲流河的凹岸，全角变化梯度小于第三预设阈值的所述采样点所在的区域为所述曲流河的凸岸；
所述曲流河的迁移方向由所述凸岸指向所述凹岸。
9.根据权利要求1所述的预测点坝砂体位置的装置，其特征在于，所述采样点的变异系数的计算公式为：
<mrow> <msub> <mi>V</mi> <mi>i</mi> </msub> <mo>=</mo> <mfrac> <msqrt> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msup> <mrow> <mo>(</mo> <msub> <mi>A</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>/</mo> <mi>n</mi> </mrow> </msqrt> <mover> <mi>A</mi> <mo>&amp;OverBar;</mo> </mover> </mfrac> </mrow>
其中，
Vi——采样点i的变异系数；
Ai——第i个样品的振幅属性值；
——第一区域内所有采样点的振幅属性平均值；
n——采样范围内采样点的个数。
10.根据权利要求1所述的预测点坝砂体位置的装置，其特征在于，所述采样点的全角变化梯度的计算公式为：
<mrow> <msub> <mi>K</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>&amp;phi;</mi> <mi>a</mi> </msub> <mo>-</mo> <msub> <mi>&amp;phi;</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>&amp;theta;</mi> <mi>a</mi> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow>
其中，
Kab——A、B两点的全角变化梯度；
φa、φb——A、B两点的倾角；
θa、θb——A、B两点的方位角。
CN201710889005.8A 2017-09-27 2017-09-27 一种预测点坝砂体位置的方法及装置 Active CN107807386B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710889005.8A CN107807386B (zh) 2017-09-27 2017-09-27 一种预测点坝砂体位置的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710889005.8A CN107807386B (zh) 2017-09-27 2017-09-27 一种预测点坝砂体位置的方法及装置

Publications (2)

Publication Number Publication Date
CN107807386A true CN107807386A (zh) 2018-03-16
CN107807386B CN107807386B (zh) 2019-04-09

Family

ID=61591780

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710889005.8A Active CN107807386B (zh) 2017-09-27 2017-09-27 一种预测点坝砂体位置的方法及装置

Country Status (1)

Country Link
CN (1) CN107807386B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109541683A (zh) * 2018-12-26 2019-03-29 中国海洋石油集团有限公司 一种地下点坝倾斜地层模型网格的构建方法及系统
CN110319919A (zh) * 2018-03-29 2019-10-11 上海梅山钢铁股份有限公司 一种应用于风机混泥土基础的松动判别方法
CN111460680A (zh) * 2020-04-17 2020-07-28 中海石油(中国)有限公司 一种曲流河相复合砂体构型单元几何特征的描述方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1664616A (zh) * 2005-04-07 2005-09-07 西北大学 等速砂泥岩地层的储层预测油气的方法
US20140334262A1 (en) * 2013-04-22 2014-11-13 Robert H. Brune Method and Apparatus for Active Seismic Shear Wave Monitoring of Hydro-Fracturing of Oil and Gas Reservoirs Using Arrays of Multi-Component Sensors and Controlled Seismic Sources
CN105607146A (zh) * 2015-09-10 2016-05-25 中国海洋石油总公司 一种曲流河砂体规模的定量表征方法
CN105717552A (zh) * 2016-02-04 2016-06-29 中国石油大学(华东) 基于测井曲线高频旋回预测湖相滩坝砂体空间展布形态的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1664616A (zh) * 2005-04-07 2005-09-07 西北大学 等速砂泥岩地层的储层预测油气的方法
US20140334262A1 (en) * 2013-04-22 2014-11-13 Robert H. Brune Method and Apparatus for Active Seismic Shear Wave Monitoring of Hydro-Fracturing of Oil and Gas Reservoirs Using Arrays of Multi-Component Sensors and Controlled Seismic Sources
CN105607146A (zh) * 2015-09-10 2016-05-25 中国海洋石油总公司 一种曲流河砂体规模的定量表征方法
CN105717552A (zh) * 2016-02-04 2016-06-29 中国石油大学(华东) 基于测井曲线高频旋回预测湖相滩坝砂体空间展布形态的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王谦: "复合曲流砂体中单一点坝精细描述方法研究", 《内蒙古石油化工》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110319919A (zh) * 2018-03-29 2019-10-11 上海梅山钢铁股份有限公司 一种应用于风机混泥土基础的松动判别方法
CN109541683A (zh) * 2018-12-26 2019-03-29 中国海洋石油集团有限公司 一种地下点坝倾斜地层模型网格的构建方法及系统
CN109541683B (zh) * 2018-12-26 2020-04-07 中国海洋石油集团有限公司 一种地下点坝倾斜地层模型网格的构建方法及系统
CN111460680A (zh) * 2020-04-17 2020-07-28 中海石油(中国)有限公司 一种曲流河相复合砂体构型单元几何特征的描述方法
CN111460680B (zh) * 2020-04-17 2024-03-15 中海石油(中国)有限公司 一种曲流河相复合砂体构型单元几何特征的描述方法

Also Published As

Publication number Publication date
CN107807386B (zh) 2019-04-09

Similar Documents

Publication Publication Date Title
CN103645516B (zh) 基于岩石物理相控油气作用确定油气产能的方法
Zhao et al. An approach to predict the height of fractured water-conducting zone of coal roof strata using random forest regression
CN107807386A (zh) 一种预测点坝砂体位置的方法及装置
Hewett et al. Theory for the semi-analytical calculation of oil recovery and effective relative permeabilities using streamtubes
CN105425292A (zh) 一种油气预测方法及装置
CN103150421A (zh) 利用位移监测同时确定抗滑桩桩位和临界深度的方法
CN106202813B (zh) 一种计算地层原油粘度的方法及装置
CN107808032A (zh) 页岩气储层含水饱和度计算方法及系统
Berlin et al. Steepened channels upstream of knickpoints: Controls on relict landscape response
CN104564069B (zh) 一种基于方格网法的地面动态沉陷预测与复垦方法
Siddiqui et al. DEM-based drainage network analysis using steepness and Hack SL indices to identify areas of differential uplift in Emilia–Romagna Apennines, northern Italy
CN104459800A (zh) 一种预测砂体尖灭的方法和装置
Zakirov et al. Optimal control of field development in a closed loop
CN106483562A (zh) 一种基于地震振幅预测地层厚度变化的方法
Jayne et al. Permeability correlation structure of the Columbia River Plateau and implications for fluid system architecture in continental large igneous provinces
CN104060624A (zh) 复杂软土地基中桩基的设置方法
CN110807573B (zh) 一种海相碳酸盐岩油气有利勘探区的确定方法及系统
Pei et al. Determination of infill drilling time window based on depletion-induced stress evolution of shale reservoirs with complex natural fractures
Sedghi et al. Groundwater flow to a general well configuration in an unconfined aquifer overlying a fractured bedrock
Tumlinson et al. Numerical evaluation of pumping well transmissivity estimates in laterally heterogeneous formations
Guo et al. A fractal simulation method for simulating the resource abundance of oil and gas and its application
Vargas-Silva et al. Integration of Improved Asymmetric Frac Design Using Strain Derived From Geomechanical Modeling in Reservoir Simulation
CN107329169A (zh) 一种角度道集的提取方法及装置
Klemetsdal et al. Implicit high-resolution compositional simulation with optimal ordering of unknowns and adaptive spatial refinement
Movahedian et al. Management of an aquifer with an emphasis on the interaction of the aquifer and river using MODFLOW model in Gotvand–Aghili plain, Khuzestan, Iran

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