CN110749928B - 钙质砂岩识别方法和装置 - Google Patents

钙质砂岩识别方法和装置 Download PDF

Info

Publication number
CN110749928B
CN110749928B CN201810821414.9A CN201810821414A CN110749928B CN 110749928 B CN110749928 B CN 110749928B CN 201810821414 A CN201810821414 A CN 201810821414A CN 110749928 B CN110749928 B CN 110749928B
Authority
CN
China
Prior art keywords
wave impedance
longitudinal wave
seismic data
data
curve
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
CN201810821414.9A
Other languages
English (en)
Other versions
CN110749928A (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
Daqing Oilfield Co Ltd
Original Assignee
Petrochina Co Ltd
Daqing Oilfield 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, Daqing Oilfield Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201810821414.9A priority Critical patent/CN110749928B/zh
Publication of CN110749928A publication Critical patent/CN110749928A/zh
Application granted granted Critical
Publication of CN110749928B publication Critical patent/CN110749928B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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
    • G01V1/50—Analysing data
    • 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/6169—Data from specific type of measurement using 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)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种钙质砂岩识别方法和装置，属于地质勘探技术领域。该方法包括对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据；将纵波阻抗数据与地震子波进行褶积，以得到合成地震数据；计算合成地震数据与实际地震数据的相关系数；当相关系数小于第一设定值时，再次依次获取纵波阻抗数据和合成地震数据，直至相关系数不小于第一设定值；将纵波阻抗数据中纵波阻抗不低于第二设定值的区域解释为钙质砂岩区域，相关系数可以反映出合成地震数据与实际地震数据的相关程度，相关程度越大，则表明合成地震数据越接近实际地震数据，在纵波阻抗数据中纵波阻抗不低于第二设定值的区域则可以被解释为钙质砂岩区域。

Description

钙质砂岩识别方法和装置
技术领域
本发明涉及地震勘探技术领域，特别涉及一种钙质砂岩识别方法和装置。
背景技术
钙质砂岩指的是含有钙质胶结物的致密砂岩。钙质砂岩是碎屑岩储层中广泛存在的岩石类型。
钙质砂岩的分布对剩余油分布有着较大的影响，在油田开发中，识别出钙质砂岩，弄清钙质砂岩的分布有利于提高开发效率，提高产量。
目前常用的识别钙质砂岩的方法为测井识别和井间建模识别，这种方法在距离井较近的区域可以较好的识别钙质砂岩，但是在距离井较远的地方则识别的准确度较低，且随着距离的增加，准确度会越来越低。
发明内容
本发明实施例提供了一种钙质砂岩识别方法和装置，以解决现有方法难以准确识别钙质砂岩的问题。所述技术方案如下：
一方面，本发明实施例提供了一种钙质砂岩识别方法，所述方法包括：
对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据；
将所述纵波阻抗数据与地震子波进行褶积，以得到合成地震数据；
计算所述合成地震数据与实际地震数据的相关系数；
当所述相关系数小于第一设定值时，再次依次获取所述纵波阻抗数据和所述合成地震数据，直至所述相关系数不小于第一设定值；
当所述相关系数不小于第一设定值时，将所述纵波阻抗数据中纵波阻抗不低于第二设定值的区域解释为钙质砂岩区域。
可选地，在所述对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟之前，所述方法还包括：
获取研究区目标层的井点深度域纵波阻抗曲线；
将所述井点深度域纵波阻抗曲线转换成所述井点时间域纵波阻抗曲线。
可选地，所述获取研究区目标层的井点深度域纵波阻抗曲线，包括：
获取研究区目标层的声波速度曲线和密度曲线；
根据所述声波速度曲线和所述密度曲线得到井点深度域纵波阻抗曲线。
可选地，所述计算所述合成地震数据与实际地震数据的相关系数，包括：
根据以下公式计算所述相关系数：
Figure BDA0001741479530000021
其中，ρxy为相关系数，x为所述合成地震数据，y为所述实际地震数据，cov(x，y)为x和y的协方差，D(x)为x的方差，D(y)为y的方差。
可选地，在所述将所述纵波阻抗数据中不低于第二设定值的区域解释为钙质砂岩区域之前，所述方法还包括：
根据取芯井处的所述纵波阻抗曲线与所述取芯井中获取的岩心样本确定所述第二设定值。
另一方面，本发明实施例提供了一种钙质砂岩识别装置，所述钙质砂岩识别装置包括：
获取模块，用于对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据，并将所述纵波阻抗数据与地震子波进行褶积，以得到合成地震数据；
处理模块，用于计算所述合成地震数据与实际地震数据的相关系数，并且当所述相关系数小于第一设定值时，控制所述获取模块再次依次获取所述纵波阻抗数据和所述合成地震数据，直至所述相关系数不小于第一设定值；当所述相关系数不小于所述第一设定值时，将所述纵波阻抗数据中纵波阻抗不低于第二设定值的区域解释为钙质砂岩区域。
可选地，所述获取模块还用于获取研究区目标层的井点深度域纵波阻抗曲线，并将所述井点深度域纵波阻抗曲线转换成所述井点时间域纵波阻抗曲线。
可选地，所述获取模块还用于获取研究区目标层的声波速度曲线和密度曲线，并根据所述声波速度曲线和所述密度曲线得到井点深度域纵波阻抗曲线。
可选地，所述处理模块被配置为根据以下公式计算所述相关系数：
Figure BDA0001741479530000031
其中，ρxy为相关系数，x为所述合成地震数据，y为所述实际地震数据，cov(x，y)为x和y的协方差，D(x)为x的方差，D(y)为y的方差。
可选地，所述处理模块还用于根据取芯井处的所述纵波阻抗曲线与所述取芯井中获取的岩心样本确定所述第二设定值。
本发明实施例提供的技术方案带来的有益效果至少包括：通过对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据，虽然在井点与井点之间并没有直接获取纵波阻抗曲线，但是通过序贯高斯模拟可以模拟出纵波阻抗数据，纵波阻抗数据中包括了研究区目标层的各个位置的纵波阻抗曲线，将纵波阻抗数据与地震子波进行褶积，得到合成地震数据，通过计算合成地震数据与实际地震数据的相关系数，相关系数可以反映出合成地震数据与实际地震数据的相关程度，相关程度越大，则表明合成地震数据越接近实际地震数据，模拟出的纵波阻抗数据越准确，当相关系数不小于第一设定值时认为模拟出的纵波阻抗数据可以准确反映出不同位置真实的纵波阻抗，由于钙质砂岩与泥岩、砂岩等岩石的纵波阻抗有着较大的区别，因此在纵波阻抗数据中纵波阻抗不低于第二设定值的区域则可以被解释为钙质砂岩区域。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种钙质砂岩识别方法流程图；
图2是本发明实施例提供的另一种钙质砂岩识别方法流程图；
图3是本发明实施例提供的一个井点深度域纵波阻抗曲线；
图4是本发明实施例提供的一条井点时间域纵波阻抗曲线；
图5是本发明实施例提供的一种纵波阻抗与地震子波的褶积示意图；
图6是本发明实施例提供的一种纵波阻抗与岩心样品的统计示意图；
图7是本发明实施例提供的一种钙质砂岩区域分布示意图；
图8是本发明实施例提供的一种钙质砂岩识别装置的结构框图；
图9示出了本发明一个示例性实施例提供的计算机设备的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的一种钙质砂岩识别方法流程图。如图1所示，该方法包括：
S11：对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据。
S12：将纵波阻抗数据与地震子波进行褶积，以得到合成地震数据。
S13：计算合成地震数据与实际地震数据的相关系数。
S14：比较相关系数与第一设定值。
当相关系数小于第一设定值时，返回步骤S11。
当相关系数不小于第一设定值时，执行步骤S15。
S15：将纵波阻抗数据中纵波阻抗不低于第二设定值的区域解释为钙质砂岩区域。
通过对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据，虽然在井点与井点之间并没有直接获取纵波阻抗曲线，但是通过序贯高斯模拟可以模拟出纵波阻抗数据，纵波阻抗数据中包括了各个位置的纵波阻抗曲线，将纵波阻抗数据与地震子波进行褶积，得到合成地震数据，通过计算合成地震数据与实际地震数据的相关系数，相关系数可以反映出合成地震数据与实际地震数据的相关程度，相关程度越大，则表面合成地震数据越接近实际地震数据，模拟出的纵波阻抗数据越准确，当相关系数不小于第一设定值时认为模拟出的纵波阻抗数据可以准确反映出不同位置真实的纵波阻抗，由于钙质砂岩与泥岩、砂岩等岩石的纵波阻抗有着较大的区别，因此在纵波阻抗数据中纵波阻抗不低于第二设定值的区域则可以被解释为钙质砂岩区域。
图2是本发明实施例提供的另一种钙质砂岩识别方法流程图。如图2所示，该方法包括：
S21：获取研究区目标层的井点深度域纵波阻抗曲线。
研究区目标层为需要确定钙质砂岩分布的油藏区域，井点深度域纵波阻抗曲线可以从已有的测井资料中获取。
此外，井点深度域纵波阻抗曲线也可以采用以下方式获取：
步骤一：获取研究区目标层的声波速度曲线和密度曲线。
声波速度曲线和密度曲线可以通过测井获取。
步骤二：根据声波速度曲线和密度曲线得到井点深度域纵波阻抗曲线。
声波速度曲线可以反映出在不同深度的岩层中纵波的传播速度，密度曲线可以反映出不同深度的岩层的密度，由于声波速度与岩层的密度相乘可以得到岩层的纵波阻抗，因此根据声波速度曲线和密度曲线可以获取到井点深度域纵波阻抗曲线。
图3是本发明实施例提供的一个井点深度域纵波阻抗曲线，在图3中，纵坐标表示深度，横坐标表示纵波阻抗，由该井点深度域纵波阻抗曲线可以知道，对于地层中不同的岩石结构，纵波阻抗也不同，因此井点深度域纵波阻抗曲线可以反映出井点处的地层结构。同时由图3可以知道钙质砂岩具有较高的纵波阻抗。
S22：将井点深度域纵波阻抗曲线转换成井点时间域纵波阻抗曲线。
图4是本发明实施例提供的一条井点时间域纵波阻抗曲线，在图4中，纵坐标表示时间，横坐标表示纵波阻抗，该井点时间域纵波阻抗曲线可以通过井震标定的方法将井点深度域纵波阻抗曲线转换到时间域得到。
此外，若存在记载有井点时间域纵波阻抗曲线的地质勘探资料，也可以从地质勘探资料中直接获取。
S23：对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据。
序贯高斯模拟得到的纵波阻抗数据为三维数据体，可以显示出研究区内各位置的纵波阻抗。
对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟后，获取的纵波阻抗数据是研究区目标层的所有位置的纵波阻抗曲线的集合，即包括了井所在位置的纵波阻抗曲线，也包括了井间位置的纵波阻抗曲线。由于是根据井点时间域纵波阻抗曲线进行的序贯高斯模拟，因此在井所在位置的纵波阻抗曲线最接近实际的地质结构，而井间位置的纵波阻抗曲线则与实际的地质结构存在一定差别。
S24：将纵波阻抗数据与地震子波进行褶积，以得到合成地震数据。
图5是本发明实施例提供的一种纵波阻抗与地震子波的褶积示意图，图5示出了图4所示的纵波阻抗曲线51与地震子波52褶积后得到的合成记录53，图5中还示出了部分实际地震数据54。
将纵波阻抗数据与地震子波进行褶积时，可以选取多个不同位置，以多个不同位置的纵波阻抗曲线与地震子波进行褶积，这样就可以得到合成地震数据。合成地震数据为三维数据体。
在选取多个不同位置时，可以将研究区的地表划分成网格状，选取网格点处的纵波阻抗曲线与地震子波进行褶积，在划分研究区时，网格的大小可以与研究区的实际地震数据的网格大小相同，这样使得合成地震数据与实际地震数据具有相同的分辨率。例如实际地震数据的网格为10m×10m，则可以将研究区划分成10m×10m的网格状。
S25：计算合成地震数据与实际地震数据的相关系数。
相关系数可以体现出合成地震数据与实际地震数据的相关程度，相关系数越高，则相关程度越高，表明合成地震数据与实际地震数据越接近，这就说明在S23中通过序贯高斯模拟获取的纵波阻抗数据越准确，越能够反映出实际的地层结构。
实际地震数据可以通过地震勘探技术获取，也可以是已有资料中所记载的地震数据。
可以根据以下公式计算相关系数ρxy：
Figure BDA0001741479530000061
其中，x为合成地震数据，y为实际地震数据，cov(x，y)为x和y的协方差，D(x)为x的方差，D(y)为y的方差。
S26：比较相关系数与第一设定值。
当相关系数小于第一设定值时，返回步骤S23。
当相关系数不小于第一设定值时，执行步骤S27。
示例性地，第一设定值可以设置为80％。当相关系数大于或等于第一设定值时可以认为合成地震数据与实际地震数据足够接近，通过序贯高斯模拟获取的纵波阻抗数据已经可以较好的反应实际的地层结构，当相关系数小于第一设定值时，则认为通过序贯高斯模拟获取的纵波阻抗数据与实际的地层结构有较大的偏差，因此需要重新进行序贯高斯模拟，以再次获取纵波阻抗数据。
S27：将纵波阻抗数据中纵波阻抗不低于第二设定值的区域解释为钙质砂岩区域。
在执行S27之前，可以根据取芯井处的纵波阻抗曲线与取芯井中获取的岩心样本确定第二设定值。
图6是本发明实施例提供的一种纵波阻抗与岩心样品的统计示意图，如图6所示，地层中不同的岩石的纵波阻抗不同，通常钙质砂岩的纵波阻抗比砂岩和泥岩要高，通过将取芯井处的纵波阻抗曲线与取芯井中获取的岩心样本对比，可以得到钙质砂岩的纵波阻抗的下限值，该下限值可以作为第二设定值。纵波阻抗高于第二设定值的区域则可以解释为钙质砂岩区域。由图6可以得到第二设定值为8.2×106kg/m3*m/s。
需要说明的是，第二设定值的实际取值对于不同的研究区可能会不同，由图6所确定的第二设定值的数值仅为举例。
此外，若存在记载有第二设定值的地质勘探资料，也可以从地质勘探资料中直接获取。
图7是本发明实施例提供的一种钙质砂岩区域分布示意图，图7是在序贯高斯模拟得到的纵波阻抗数据的基础上提取纵波阻抗不低于8.2×106kg/m3*m/s的区域得到，根据图7可以准确获知钙质砂岩的分布。
图8是本发明实施例提供的一种钙质砂岩识别装置的结构框图，如图8所示，该钙质砂岩识别装置包括获取模块10和处理模块20。获取模块10和处理模块20可以根据前述的步骤S11～S15或步骤S21～S27识别钙质砂岩。
其中，获取模块10用于对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据，并将纵波阻抗数据与地震子波进行褶积，以得到合成地震数据。
处理模块20用于计算合成地震数据与实际地震数据的相关系数，并且当相关系数小于第一设定值时，控制获取模块再次依次获取纵波阻抗数据和合成地震数据，直至相关系数不小于第一设定值；当相关系数不小于第一设定值时，将纵波阻抗数据中纵波阻抗不低于第二设定值的区域解释为钙质砂岩区域。
通过对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据，虽然在井点与井点之间并没有直接获取纵波阻抗曲线，但是通过序贯高斯模拟可以模拟出纵波阻抗数据，纵波阻抗数据中包括了各个位置的纵波阻抗曲线，将纵波阻抗数据与地震子波进行褶积，得到合成地震数据，通过计算合成地震数据与实际地震数据的相关系数，相关系数可以反映出合成地震数据与实际地震数据的相关程度，相关程度越大，则表面合成地震数据越接近实际地震数据，模拟出的纵波阻抗数据越准确，当相关系数不小于第一设定值时认为模拟出的纵波阻抗数据可以准确反映出不同位置真实的纵波阻抗，由于钙质砂岩与泥岩、砂岩等岩石的纵波阻抗有着较大的区别，因此在纵波阻抗数据中纵波阻抗不低于第二设定值的区域则可以被解释为钙质砂岩区域。
示例性地，第一设定值可以设置为80％。当相关系数大于或等于第一设定值时可以认为合成地震数据与实际地震数据足够接近，通过序贯高斯模拟获取的纵波阻抗数据已经可以较好的反应实际的地层结构，当相关系数小于第一设定值时，则认为通过序贯高斯模拟获取的纵波阻抗数据与实际的地层结构有较大的偏差，因此需要重新进行序贯高斯模拟，以再次获取纵波阻抗数据。
获取模块10还可以用于获取研究区目标层的井点深度域纵波阻抗曲线，并将井点深度域纵波阻抗曲线转换成井点时间域纵波阻抗曲线。研究区目标层为需要确定钙质砂岩分布的油藏区域，井点深度域纵波阻抗曲线可以从已有的测井资料中获取。获取模块具体可以根据前述的步骤S21获取研究区目标层的井点深度域纵波阻抗曲线，并将井点深度域纵波阻抗曲线转换成井点时间域纵波阻抗曲线。
此外，获取模块10还可以用于获取研究区目标层的声波速度曲线和密度曲线，并根据声波速度曲线和密度曲线得到井点深度域纵波阻抗曲线。声波速度曲线和密度曲线可以通过测井获取。声波速度曲线可以反映出在不同深度的岩层中纵波的传播速度，密度曲线可以反映出不同深度的岩层的密度，由于声波速度与岩层的密度相乘可以得到岩层的纵波阻抗，因此根据声波速度曲线和密度曲线可以得到到井点深度域纵波阻抗曲线。
在计算合成地震数据与实际地震数据的相关系数时，处理模块20被配置为根据以下公式计算相关系数ρxy：
Figure BDA0001741479530000091
其中，x为合成地震数据，y为实际地震数据，cov(x，y)为x和y的协方差，D(x)为x的方差，D(y)为y的方差。相关系数可以体现出合成地震数据与实际地震数据的相关程度，相关系数越高，则相关程度越高，表明合成地震数据与实际地震数据越接近，这就说明通过序贯高斯模拟得到的纵波阻抗数据越准确，越能够反映出实际的地层结构。
处理模块20还可以用于根据取芯井处的纵波阻抗曲线与取芯井中获取的岩心样本确定第二设定值。地层中不同的岩石的纵波阻抗不同，通常钙质砂岩的纵波阻抗比砂岩和泥岩要高，通过将取芯井处的纵波阻抗曲线与取芯井中获取的岩心样本对比，可以得到钙质砂岩的纵波阻抗的下限值，该下限值可以作为第二设定值。纵波阻抗高于第二设定值的区域则可以解释为钙质砂岩区域。本实施例中第二设定值为8.2×106kg/m3*m/s。第二设定值的实际取值对于不同的研究区可能会不同，本实施例中的第二设定值的具体数值仅为举例。
图9示出了本发明一个示例性实施例提供的计算机设备的结构框图。如图9所示，计算机设备300包括中央处理单元(CPU)301、包括随机存取存储器(RAM)302和只读存储器(ROM)303的系统存储器304，以及连接系统存储器304和中央处理单元301的系统总线305。计算机设备300还包括帮助计算机内的各个器件之间传输信息的基本输入/输出系统(I/O系统)306，和用于存储操作系统313、应用程序314和其他程序模块315的大容量存储设备307。
基本输入/输出系统306包括有用于显示信息的显示器308和用于用户输入信息的诸如鼠标、键盘之类的输入设备309。其中显示器308和输入设备309都通过连接到系统总线305的输入输出控制器310连接到中央处理单元301。基本输入/输出系统306还可以包括输入输出控制器310以用于接收和处理来自键盘、鼠标、或电子触控笔等多个其他设备的输入。类似地，输入输出控制器310还提供输出到显示屏、打印机或其他类型的输出设备。
大容量存储设备307通过连接到系统总线305的大容量存储控制器(未示出)连接到中央处理单元301。大容量存储设备307及其相关联的计算机可读介质为计算机设备300提供非易失性存储。也就是说，大容量存储设备307可以包括诸如硬盘或者CD-ROM驱动器之类的计算机可读介质(未示出)。
不失一般性，计算机可读介质可以包括计算机存储介质和通信介质。计算机存储13介质包括以用于存储诸如计算机可读指令、数据结构、程序模块或其他数据等信息的任何方法或技术实现的易失性和非易失性、可移动和不可移动介质。计算机存储介质包括RAM、ROM、EPROM、EEPROM、闪存或其他固态存储其技术，CD-ROM、DVD或其他光学存储、磁带盒、磁带、磁盘存储或其他磁性存储设备。当然，本领域技术人员可知计算机存储介质不局限于上述几种。上述的系统存储器304和大容量存储设备307可以统称为存储器。
根据本发明的各种实施例，计算机设备300还可以通过诸如因特网等网络连接到网络上的远程计算机运行。也即计算机设备300可以通过连接在系统总线305上的网络接口单元311连接到网络312，或者说，也可以使用网络接口单元311来连接到其他类型的网络或远程计算机系统(未示出)。
上述存储器还包括一个或者一个以上的程序，一个或者一个以上程序存储于存储器中，被配置由CPU执行。所述一个或者一个以上程序包含用于进行本发明实施例提供的如前所述的地震能量均衡方法的指令。
本发明实施例还提供了一种非临时性计算机可读存储介质，当该存储介质中的指令由钙质砂岩识别装置的处理器执行时，使得钙质砂岩识别装置能够执行图1～2中任一项提供的钙质砂岩识别方法。
一种包含指令的计算机程序产品，当其在计算机上运行时，使得计算机执行图1～2其中任一项提供的钙质砂岩识别方法的指令。
需要说明的是：上述实施例提供的钙质砂岩识别装置在进行钙质砂岩识别时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的钙质砂岩识别装置与钙质砂岩识别方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种钙质砂岩识别方法，其特征在于，所述方法包括：
对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据；
将所述纵波阻抗数据与地震子波进行褶积，以得到合成地震数据；
计算所述合成地震数据与实际地震数据的相关系数；
当所述相关系数小于第一设定值时，对所述研究区目标层的所述井点时间域纵波阻抗曲线再次进行序贯高斯模拟，以得到新的纵波阻抗数据和新的合成地震数据，计算所述新的合成地震数据与实际地震数据的相关系数，直至所述相关系数不小于第一设定值；
当所述相关系数不小于第一设定值时，将所述纵波阻抗数据中纵波阻抗不低于第二设定值的区域解释为钙质砂岩区域。
2.根据权利要求1所述的钙质砂岩识别方法，其特征在于，所述计算所述合成地震数据与实际地震数据的相关系数，包括：
根据以下公式计算所述相关系数：
Figure FDA0003042398850000011
其中，ρxy为相关系数，x为所述合成地震数据，y为所述实际地震数据，cov(x，y)为x和y的协方差，D(x)为x的方差，D(y)为y的方差。
3.根据权利要求1或2所述的钙质砂岩识别方法，其特征在于，在所述对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟之前，所述方法还包括：
获取研究区目标层的井点深度域纵波阻抗曲线；
将所述井点深度域纵波阻抗曲线转换成所述井点时间域纵波阻抗曲线。
4.根据权利要求3所述的钙质砂岩识别方法，其特征在于，所述获取研究区目标层的井点深度域纵波阻抗曲线，包括：
获取研究区目标层的声波速度曲线和密度曲线；
根据所述声波速度曲线和所述密度曲线得到井点深度域纵波阻抗曲线。
5.根据权利要求1或2所述的钙质砂岩识别方法，其特征在于，在所述将所述纵波阻抗数据中不低于第二设定值的区域解释为钙质砂岩区域之前，所述方法还包括：
根据取芯井处的所述纵波阻抗曲线与所述取芯井中获取的岩心样本确定所述第二设定值。
6.一种钙质砂岩识别装置，其特征在于，所述钙质砂岩识别装置包括：
获取模块，用于对研究区目标层的井点时间域纵波阻抗曲线进行序贯高斯模拟，以得到纵波阻抗数据，并将所述纵波阻抗数据与地震子波进行褶积，以得到合成地震数据；
处理模块，用于计算所述合成地震数据与实际地震数据的相关系数，并且当所述相关系数小于第一设定值时，控制所述获取模块对所述研究区目标层的所述井点时间域纵波阻抗曲线再次进行序贯高斯模拟，以得到新的纵波阻抗数据和新的合成地震数据，计算所述新的合成地震数据与实际地震数据的相关系数，直至所述相关系数不小于第一设定值；当所述相关系数不小于所述第一设定值时，将所述纵波阻抗数据中纵波阻抗不低于第二设定值的区域解释为钙质砂岩区域。
7.根据权利要求6所述的钙质砂岩识别装置，其特征在于，所述获取模块还用于获取研究区目标层的井点深度域纵波阻抗曲线，并将所述井点深度域纵波阻抗曲线转换成所述井点时间域纵波阻抗曲线。
8.根据权利要求7所述的钙质砂岩识别装置，其特征在于，所述获取模块还用于获取研究区目标层的声波速度曲线和密度曲线，并根据所述声波速度曲线和所述密度曲线得到井点深度域纵波阻抗曲线。
9.根据权利要求6～8任一项所述的钙质砂岩识别装置，其特征在于，所述处理模块被配置为根据以下公式计算所述相关系数：
Figure FDA0003042398850000031
其中，ρxy为相关系数，x为所述合成地震数据，y为所述实际地震数据，cov(x，y)为x和y的协方差，D(x)为x的方差，D(y)为y的方差。
10.一种计算机可读存储介质，其特征在于，所述存储介质中存储有至少一条指令，所述指令由处理器加载并执行以实现如权利要求1至5任一项所述的钙质砂岩识别方法。
CN201810821414.9A 2018-07-24 2018-07-24 钙质砂岩识别方法和装置 Active CN110749928B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810821414.9A CN110749928B (zh) 2018-07-24 2018-07-24 钙质砂岩识别方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810821414.9A CN110749928B (zh) 2018-07-24 2018-07-24 钙质砂岩识别方法和装置

Publications (2)

Publication Number Publication Date
CN110749928A CN110749928A (zh) 2020-02-04
CN110749928B true CN110749928B (zh) 2021-09-28

Family

ID=69275497

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810821414.9A Active CN110749928B (zh) 2018-07-24 2018-07-24 钙质砂岩识别方法和装置

Country Status (1)

Country Link
CN (1) CN110749928B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101126815A (zh) * 2006-08-17 2008-02-20 中国石油天然气股份有限公司 一种利用地震岩性因子和岩性阻抗进行油气检测的方法
WO2014042866A1 (en) * 2012-09-13 2014-03-20 Chevron U.S.A. Inc. System and method for performing simultaneous petrophysical analysis of composition and texture of rock formations
CN103954995A (zh) * 2014-04-22 2014-07-30 核工业北京地质研究院 一种砂岩型铀矿勘查中的砂体识别方法
CN105046097A (zh) * 2015-09-09 2015-11-11 电子科技大学 基于非参数各向异性变差函数的随机建模方法
CN106842359A (zh) * 2015-12-07 2017-06-13 中国石油化工股份有限公司 利用波阻抗定量识别复杂砂砾岩体岩性的方法
CN107728202A (zh) * 2016-08-10 2018-02-23 中国石油化工股份有限公司 致密砂岩储层精细描述方法及系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101126815A (zh) * 2006-08-17 2008-02-20 中国石油天然气股份有限公司 一种利用地震岩性因子和岩性阻抗进行油气检测的方法
WO2014042866A1 (en) * 2012-09-13 2014-03-20 Chevron U.S.A. Inc. System and method for performing simultaneous petrophysical analysis of composition and texture of rock formations
CN103954995A (zh) * 2014-04-22 2014-07-30 核工业北京地质研究院 一种砂岩型铀矿勘查中的砂体识别方法
CN105046097A (zh) * 2015-09-09 2015-11-11 电子科技大学 基于非参数各向异性变差函数的随机建模方法
CN106842359A (zh) * 2015-12-07 2017-06-13 中国石油化工股份有限公司 利用波阻抗定量识别复杂砂砾岩体岩性的方法
CN107728202A (zh) * 2016-08-10 2018-02-23 中国石油化工股份有限公司 致密砂岩储层精细描述方法及系统

Also Published As

Publication number Publication date
CN110749928A (zh) 2020-02-04

Similar Documents

Publication Publication Date Title
RU2669948C2 (ru) Оптимизация многоступенчатого проекта нефтяного месторождения в условиях неопределенности
RU2602409C2 (ru) Инверсия анизотропии многоскважинной системы
CA2823710C (en) Methods and systems regarding models of underground formations
RU2601232C2 (ru) Способы и системы для включения мест выделенных точек псевдоповерхности в модели скорости сейсмических волн
RU2631407C1 (ru) Способ и устройство для обработки сейсмических сигналов
MX2007006657A (es) Estimado de modelo de baja frecuencia.
US10087733B2 (en) Fracture mapping using vertical seismic profiling wave data
US20170175492A1 (en) Methodology for building realistic numerical forward stratigraphic models in data sparse environment
US20220299476A1 (en) Systems and Methods for Ultrasonic Characterization of Permafrost, Frozen Soil and Saturated Soil Samples
CN112051609B (zh) 地震波成像方法、系统、存储介质、计算机程序、终端
EA028337B1 (ru) Моделирование методом конечных элементов скважинных сейсмических сигналов из слоистых анизотропных формаций и использование его при упругой инверсии
CN109709602B (zh) 一种远探测声波偏移成像方法、装置及系统
WO2016115437A1 (en) Approaches to tying well logs to seismic data
CN110749928B (zh) 钙质砂岩识别方法和装置
CN115877464A (zh) 一种岩性识别方法、装置、计算机设备及存储介质
CN113534263B (zh) 一种不依赖测井资料的含油气饱和度预测方法
US20170336523A1 (en) Seismic signal processing method, apparatus and system
AU2014394076B2 (en) Methods and systems for identifying and plugging subterranean conduits
CN111624662B (zh) 速度场校正方法、装置和设备
CN108983288B (zh) 基于时频谱图像特征分析的油水识别方法
CN114594513A (zh) 岩石物理模板建立方法、含气饱和度定量预测方法及装置
CN114063169B (zh) 一种波阻抗反演方法、系统、设备和存储介质
CN112888970A (zh) 确定声波慢度的方法和装置
WO2020236153A1 (en) Enhanced-resolution sonic data processing for formation body wave slowness with full offset waveform data
CN116774285B (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