CN106291714A - 获得地层吸收衰减品质因子的方法及其地层划分方法 - Google Patents
获得地层吸收衰减品质因子的方法及其地层划分方法 Download PDFInfo
- Publication number
- CN106291714A CN106291714A CN201610785617.8A CN201610785617A CN106291714A CN 106291714 A CN106291714 A CN 106291714A CN 201610785617 A CN201610785617 A CN 201610785617A CN 106291714 A CN106291714 A CN 106291714A
- Authority
- CN
- China
- Prior art keywords
- stratum
- micro
- formation
- logging
- thickness
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 35
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 54
- 238000004445 quantitative analysis Methods 0.000 claims abstract description 9
- 235000019738 Limestone Nutrition 0.000 claims description 9
- 239000006028 limestone Substances 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 6
- 239000011435 rock Substances 0.000 claims description 2
- 239000004575 stone Substances 0.000 description 31
- 238000004458 analytical method Methods 0.000 description 10
- 238000013517 stratification Methods 0.000 description 8
- 239000002245 particle Substances 0.000 description 5
- 239000004576 sand Substances 0.000 description 5
- 238000005267 amalgamation Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000003475 lamination Methods 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
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
-
- 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
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)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提供了一种获得地层吸收衰减品质因子的方法及其地层划分方法。所述地层划分方法包括：对目标地层进行岩性录井分层和微测井分层；综合分析岩性录井分层结果和微测井分层结果，得综合分层结果。其中，微测井分层包括：定量分析微测井资料，得微测井资料中每个深度点所对应的属性值；找出所述属性值内的拐点位置及其对应的地层深度值，并据此对目标地层进行分层；岩性录井分层包括：分析岩性录井资料，统计出各岩性的岩性特征及其对应的深度值，并据此对目标地层进行分层。本发明能够获得复杂地层(尤其适合于近地表地层、或者具有岩性复杂变化特征的地层等)的吸收衰减品质因子，且所获得的吸收衰减品质因子的准确度得以提高。
Description
技术领域
本发明涉及工程地球物理勘探技术领域，具体来讲，涉及一种用于获得地层吸收衰减品质因子的地层划分方法、以及获得地层吸收衰减品质因子的方法，本发明的方法尤其适合于近地表地层、或者具有岩性复杂变化特征的地层等，并且能够提高计算其吸收衰减品质因子的精确度。
背景技术
一般来讲，地震波在地层中传播时，因能量被介质吸收，从而产生了波的衰减和频散。在地球物理中，一般用吸收衰减品质因子(亦可简称为Q或Q值)来描述地层对地震波的衰减程度。吸收衰减品质因子Q越小的地层，对地震波的衰减越快。
不同的岩性地层的吸收衰减品质因子Q不同，具体来讲，一般是砂岩地层的吸收衰减品质因子Q大于泥岩地层；泥岩地层的吸收衰减品质因子Q大于灰岩地层；疏松地层的吸收衰减品质因子Q最小，对地震波的衰减最强。通常，在地震勘探中，需要尽可能减少地层吸收衰减对资料的影响，故而需要获得较精确的吸收衰减品质因子Q值，特别是获得高精度近地表松散地层衰减品质因子Q值。
现有技术在计算地层衰减品质因子Q值时，通常根据地层的岩性划分地层，这种方法在地层较厚(例如，地层在几十米或更大厚度内岩性变化不大的情况)时，其所划分地层计算得到的Q值精度较高。
然而，在地球近地表区域，地层岩性变化很快(例如，通常在几～30米范围内岩性多次变化)，Q值变化也较大，如果仍然根据岩性分层计算Q值，则由于岩层太薄，会导致同一地层内采集样点不足，从而使计算过程中统计效应变差，最终导致所计算的Q值不准。
发明内容
本发明的目的在于解决现有技术存在的上述不足中的至少一项。
例如，本发明的目的之一在于提供一种适合于复杂地层，尤其适合于近地表地层、或者具有岩性复杂变化特征(例如，岩性在小于40米内(或称30几米内变化2次或3次以上))的地层等的获得地层吸收衰减品质因子的方法及其地层划分方法。
为了实现上述目的，本发明的一方面提供了一种用于获得地层吸收衰减品质因子的地层划分方法。所述方法包括以下步骤：利用岩性录井资料对目标地层进行分层，以得到岩性录井分层结果，利用微测井资料对所述目标地层进行分层，以得到微测井分层结果；综合分析岩性录井分层结果和微测井分层结果，以得到综合分层结果。其中，所述利用微测井资料对目标地层进行分层的步骤包括：定量分析微测井资料，以得到所述微测井资料中每个深度点所对应的属性值；找出所述属性值内的拐点位置以及该拐点位置对应的地层深度值，并依据该地层深度值对目标地层进行分层；所述利用岩性录井资料对目标地层进行分层的步骤包括：分析目标地层的岩性录井资料，统计出各岩性的岩性特征及其对应的深度值，并依据所述对应的深度值对目标地层进行分层。
在本发明的一个示例性实施例中，所述综合分析岩性录井分层结果和微测井分层结果的步骤可包括：根据微测井分层结果对岩性录井分层结果中的各子地层进行合并，并且岩性录井分层结果中的厚度大于5米且为风化层、泥岩、砂岩或灰岩的子地层不应被合并，而且合并所获得的地层的厚度不应大于10米。
在本发明的一个示例性实施例中，所述属性值可包括频率、波形特征、能量和速度中的一种或两种以上。
在本发明的一个示例性实施例中，所述目标地层可以为近地表地层，或者可以为具有岩性复杂变化特征的地层，所述岩性复杂变化特征可以指岩性在小于40米内变化2次或3次以上。
在本发明的一个示例性实施例中，所述微测井资料可以由单井微测井方法或双井微测井方法采集得到。
在本发明的一个示例性实施例中，所述微测井资料可以通过井中激发地表接收或井中激发井中接收的微测井方法得到。
本发明的另一方面提供了一种获得地层吸收衰减品质因子的方法。所述方法通过如上所述的地层划分方法的综合分层结果来计算目标地层的吸收衰减品质因子。
与现有技术相比，本发明的有益效果包括：能够获得复杂地层(尤其适合于近地表地层、或者具有岩性复杂变化特征的地层等)的吸收衰减品质因子，且所获得的吸收衰减品质因子的准确度得以提高。
附图说明
图1示出了根据本发明获得地层吸收衰减品质因子的方法的一个示例性实施例的过程示意图。
具体实施方式
在下文中，将结合示例性实施例来详细说明本发明的用于获得地层吸收衰减品质因子的地层划分方法以及获得地层吸收衰减品质因子的方法。
在本发明的一个示例性实施例中，用于获得地层吸收衰减品质因子的地层划分方法包括以下步骤：
(1)岩性录井分层和微测井分层
岩性录井分层步骤可通过利用岩性录井资料对目标地层进行分层，从而得到岩性录井分层结果。具体来讲，利用岩性录井资料对目标地层进行分层的步骤可包括：分析目标地层的岩性录井资料，统计出各岩性的岩性特征及其(或这些岩性特征)对应的深度值，并依据该(或这些)对应的深度值对目标地层进行分层，从而得到岩性录井分层结果。
微测井分层步骤可通过利用微测井资料对所述目标地层进行分层，从而得到微测井分层结果。具体来讲，利用微测井资料对目标地层进行分层的步骤可包括：定量分析微测井资料，以得到微测井资料中每个深度点所对应的属性值；找出这些属性值内的拐点位置以及该(或这些)拐点位置对应的地层深度值，并依据该(或这些)地层深度值对目标地层进行分层，从而得到微测井分层结果。
在本示例性实施例中，微测井资料可由单井微测井方法或双井微测井方法采集得到。此外，微测井资料也可通过井中激发地表接收或井中激发井中接收等方式得到。
在本示例性实施例中，所述属性值包括频率、波形特征、能量和速度值中的一种或两种以上。
需要说明的是，在本发明的方法中，岩性录井分层步骤和微测井分层步骤无先后顺序的要求。也就是说，岩性录井分层步骤和微测井分层步骤可先后进行，也可同时进行。
(2)综合分层步骤
综合分层步骤通过对岩性录井分层结果和微测井分层结果进行综合分析，从而得到综合分层结果。
具体来讲，综合分析岩性录井分层结果和微测井分层结果的步骤可以为：根据微测井分层结果对岩性录井分层结果中的各子地层进行合并，并且岩性录井分层结果中的厚度大于5米且为风化层、泥岩、砂岩或灰岩的子地层不应被合并，而且合并所获得的地层的厚度不应大于10米。
一般情况下，微测井属性分层(例如，频率属性分层、波形特征属性分层、能量属性分层和速度值属性分层)的分层结果中的子层数少于岩性录井分层的子层数。
本发明的方法尤其适合以近地表地层作为目标地层，或者以具有岩性复杂变化特征(岩性在小于40米内变化2次或3次以上)的地层作为目标地层。然而，本发明的目标地层不限于此。
在本发明的另一示例性实施例中，获得地层吸收衰减品质因子的方法，可根据前述示例性实施例中的地层划分方法所得到的综合分层结果，来计算目标地层的吸收衰减品质因子。
下面将结合具体示例来进一步详细描述本发明的示例性实施例。
示例1
目标地层为某一近地表地层，其地层总厚度为30m(即，由地表面的0米至地下30米)。
采集并获取该近地表地层的岩性录井资料，经过分析和统计，获得按地层由浅到深方向计风化层(A1子地层)、泥岩(A2子地层)、泥砂岩性(A3子地层)、泥岩(A4子地层)、砂泥岩(A5子地层)、砂岩(A6子地层)、灰岩(A7子地层)等共7层岩性，所对应底界深度值分别为0～5m、5～8m、8～12m、12～14m、14～18m、18～24m、24～30m。也就是说，A1子地层的厚度为5m，A2子地层的厚度为3m，A3子地层的厚度为4m，A4子地层的厚度为2m，A5子地层的厚度为4m，A6子地层的厚度为6m，A7子地层的厚度为6m。
采集并获取该近地表地层的微测井资料，定量分析微测井资料，以得到微测井资料中每个深度点所对应的频率。找出这些频率的拐点位置以及这些拐点位置对应的地层深度值，获得频率属性的分层结果。频率属性的各子地层及其地层深度值分别为：B1子地层0～4m、B2子地层4～8m、B3子地层8～17m、B4子地层17～25m、B5子地层25～30m。也就是说，B1子地层的厚度为4m，B2子地层的厚度为4m，B3子地层的厚度为9m、B4子地层的厚度为8m、B5子地层的厚度为5m。
对上述岩性录井分层结果所获得的A1子地层至A7子地层、以及上述微测井频率属性分层结果所获得的B1子地层至B5子地层进行综合分析，可根据微测井频率属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据B3子地层的地层深度值8～17m，考虑将A3子地层至A5子地层进行合并。A3子地层、A4子地层、A5子地层中只有A4子地层为泥岩，且A4子地层的地层厚度为2m，小于5m，故而A4子地层可以被合并而不用单独分层。另外，A3子地层至A5子地层合并后地层总厚度为10m，不大于10m。因此，可将A3子地层至A5子地层合并。从而，得到合并后的分层结果为：0～5m、5～8m、8～17m、17～25m、25～30m。
示例2
目标地层为某一近地表地层，其地层总厚度为30m。
采集并获取该近地表地层的岩性录井资料，经过分析和统计，获得按地层由浅到深方向计风化层(A1子地层)、泥岩(A2子地层)、泥砂岩性(A3子地层)、泥岩(A4子地层)、砂泥岩(A5子地层)、砂岩(A6子地层)、灰岩(A7子地层)等共7层岩性，所对应底界深度值分别为0～5m、5～8m、8～12m、12～14m、14～18m、18～24m、24～30m。也就是说，A1子地层的厚度为5m，A2子地层的厚度为3m，A3子地层的厚度为4m，A4子地层的厚度为2m，A5子地层的厚度为4m，A6子地层的厚度为6m，A7子地层的厚度为6m。
采集并获取该近地表地层的微测井资料，定量分析微测井资料，以得到微测井资料中每个深度点所对应的能量。找出这些能量的拐点位置以及这些拐点位置对应的地层深度值，获得能量属性的分层结果。微测井能量属性的各子地层及其地层深度值分别为：B1子地层0～5m、B2子地层5～14m、B3子地层14～17m、B4子地层17～26m、B5子地层26～30m。也就是说，B1子地层的厚度为5m，B2子地层的厚度为9m，B3子地层的厚度为3m、B4子地层的厚度为9m、B5子地层的厚度为4m。
对上述岩性录井分层结果所获得的A1子地层至A7子地层、以及上述微测井能量属性分层结果所获得的B1子地层至B5子地层进行综合分析，可根据微测井能量属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据B2子地层的地层深度值5～14m，考虑将A2子地层至A4子地层进行合并。A2子地层、A3子地层、A4子地层中A2子地层为泥岩、A4子地层为泥岩，且A2子地层、A4子地层的地层厚度分别为3m、2m，都小于5m，故而A2子地层、A4子地层不用单独分层。另外，考虑到A2子地层至A4子地层合并后地层总厚度为9m，小于10m。因此，可将A2子地层至A4子地层合并。从而，得到合并后的分层结果为：0～5m、5～14m、14～17m、17～26m、26～30m。
示例3
目标地层为某一近地表地层，其地层总厚度为30m。
采集并获取该近地表地层的岩性录井资料，经过分析和统计，获得按地层由浅到深方向计风化层(A1子地层)、泥岩(A2子地层)、泥砂岩性(A3子地层)、泥岩(A4子地层)、砂泥岩(A5子地层)、砂岩(A6子地层)、灰岩(A7子地层)等共7层岩性，所对应底界深度值分别为0～5m、5～8m、8～12m、12～14m、14～18m、18～24m、24～30m。也就是说，A1子地层的厚度为5m，A2子地层的厚度为3m，A3子地层的厚度为4m，A4子地层的厚度为2m，A5子地层的厚度为4m，A6子地层的厚度为6m，A7子地层的厚度为6m。
采集并获取该近地表地层的微测井资料，定量分析微测井资料，以得到微测井资料中每个深度点所对应的频率和速度值。找出这些频率的拐点位置以及这些拐点位置对应的地层深度值，同时，找出这些速度值的拐点位置以及这些拐点位置对应的地层深度值，获得频率属性和速度属性的分层结果。
其中，频率属性的各子地层及其地层深度值分别为：B1子地层0～4m、B2子地层4～8m、B3子地层8～17m、B4子地层17～25m、B5子地层25～30m。也就是说，B1子地层的厚度为4m，B2子地层的厚度为4m，B3子地层的厚度为9m、B4子地层的厚度为8m、B5子地层的厚度为5m。速度属性的各子地层及其地层深度值分别为：C1子地层0～5m、C2子地层5～17m、C3子地层17～24m、C4子地层24～30m。也就是说，C1子地层的厚度为5m，C2子地层的厚度为12m，C3子地层的厚度为7m、C4子地层的厚度为6m。
对上述岩性录井分层结果所获得的A1子地层至A7子地层、上述微测井频率属性分层结果所获得的B1子地层至B5子地层、以及上述微测井速度属性分层结果所获得的C1子地层至C4子地层进行综合分析。
根据微测井频率属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据B3子地层的地层深度值8～17m，考虑将A3子地层至A5子地层进行合并。A3子地层、A4子地层、A5子地层中只有A4子地层为泥岩，且A4子地层的地层厚度为2m，小于5m，故而A4子地层可以被合并而不用单独分层。另外，A3子地层至A5子地层合并后地层总厚度为10m，不大于10m。因此，可将A3子地层至A5子地层合并。
根据微测井速度属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据C2子地层的地层深度值5～17m，考虑将A2子地层至A5子地层进行合并。A2子地层、A3子地层、A4子地层、A5子地层中A2子地层、A4子地层为泥岩，且A2子地层、A4子地层的地层厚度都小于5m，故而A2子地层、A4子地层不用单独分层。然而，考虑到A2子地层至A5子地层合并后地层总厚度为12m，大于10m。因此，不建议将A2子地层至A5子地层进行合并。然而，可将A2子地层至A4子地层合并，或者将A3子地层至A5子地层合并。
综合考虑上述合并结果，可将A3子地层至A5子地层合并，从而得到最终的分层结果为：0～5m、5～8m、8～17m、17～25m、25～30m。
示例4
目标地层为某一近地表地层，其地层总厚度为30m。
采集并获取该近地表地层的岩性录井资料，经过分析和统计，获得按地层由浅到深方向计风化层(A1子地层)、泥岩(A2子地层)、泥砂岩性(A3子地层)、泥岩(A4子地层)、砂泥岩(A5子地层)、砂岩(A6子地层)、灰岩(A7子地层)等共7层岩性，所对应底界深度值分别为0～5m、5～8m、8～12m、12～14m、14～18m、18～24m、24～30m。也就是说，A1子地层的厚度为5m，A2子地层的厚度为3m，A3子地层的厚度为4m，A4子地层的厚度为2m，A5子地层的厚度为4m，A6子地层的厚度为6m，A7子地层的厚度为6m。
采集并获取该近地表地层的微测井资料，定量分析微测井资料，以得到微测井资料中每个深度点所对应的频率、速度值和能量。找出这些频率的拐点位置以及这些拐点位置对应的地层深度值，找出这些速度值的拐点位置以及这些拐点位置对应的地层深度值，同时，找出这些能量的拐点位置以及这些拐点位置对应的地层深度值，获得频率属性、速度属性和能量属性的分层结果。
其中，频率属性的各子地层及其地层深度值分别为：B1子地层0～4m、B2子地层4～8m、B3子地层8～17m、B4子地层17～25m、B5子地层25～30m。也就是说，B1子地层的厚度为4m，B2子地层的厚度为4m，B3子地层的厚度为9m、B4子地层的厚度为8m、B5子地层的厚度为5m。速度属性的各子地层及其地层深度值分别为：C1子地层0～5m、C2子地层5～17m、C3子地层17～24m、C4子地层24～30m。也就是说，C1子地层的厚度为5m，C2子地层的厚度为12m，C3子地层的厚度为7m、C4子地层的厚度为6m。能量属性的各子地层及其地层深度值分别为：D1子地层0～5m、D2子地层5～14m、D3子地层14～17m、D4子地层17～26m、D5子地层26～30m。也就是说，D1子地层的厚度为5m，D2子地层的厚度为9m，D3子地层的厚度为3m、D4子地层的厚度为9m、D5子地层的厚度为4m。
对上述岩性录井分层结果所获得的A1子地层至A7子地层、上述微测井频率属性分层结果所获得的B1子地层至B5子地层、上述微测井速度属性分层结果所获得的C1子地层至C4子地层、以及上述微测井能量属性分层结果所获得的D1子地层至D5子地层进行综合分析。
根据微测井频率属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据B3子地层的地层深度值8～17m，考虑将A3子地层至A5子地层进行合并。A3子地层、A4子地层、A5子地层中只有A4子地层为泥岩，且A4子地层的地层厚度为2m，小于5m，故而A4子地层可以被合并而不用单独分层。另外，A3子地层至A5子地层合并后地层总厚度为10m，不大于10m。因此，可将A3子地层至A5子地层合并。
根据微测井速度属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据C2子地层的地层深度值5～17m，考虑将A2子地层至A5子地层进行合并。A2子地层、A3子地层、A4子地层、A5子地层中A2子地层、A4子地层为泥岩，且A2子地层、A4子地层的地层厚度都小于5m，故而A2子地层、A4子地层不用单独分层。然而，考虑到A2子地层至A5子地层合并后地层总厚度为12m，大于10m。因此，不建议将A2子地层至A5子地层进行合并。然而，可将A2子地层至A4子地层合并，或者将A3子地层至A5子地层合并。
根据微测井能量属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据D2子地层的地层深度值5～14m，考虑将A2子地层至A4子地层进行合并。A2子地层、A3子地层、A4子地层中A2子地层为泥岩、A4子地层为泥岩，且A2子地层、A4子地层的地层厚度分别为3m、2m，都小于5m，故而A2子地层、A4子地层不用单独分层。另外，考虑到A2子地层至A4子地层合并后地层总厚度为9m，小于10m。因此，可将A2子地层至A4子地层合并。
综合考虑上述合并结果，可将A3子地层至A5子地层合并，从而得到最终的分层结果为：0～5m、5～8m、8～17m、17～25m、25～30m。或者，也可将A2子地层至A4子地层合并，从而得到最终的分层结果为：0～5m、5～14m、14～17m、17～26m、26～30m。
示例5
目标地层为某一近地表地层，其地层总厚度为30m。
采集并获取该近地表地层的岩性录井资料，经过分析和统计，获得按地层由浅到深方向计风化层(A1子地层)、泥岩(A2子地层)、泥砂岩性(A3子地层)、泥岩(A4子地层)、砂泥岩(A5子地层)、砂岩(A6子地层)、灰岩(A7子地层)等共7层岩性，所对应底界深度值分别为0～5m、5～8m、8～12m、12～14m、14～18m、18～24m、24～30m。也就是说，A1子地层的厚度为5m，A2子地层的厚度为3m，A3子地层的厚度为4m，A4子地层的厚度为2m，A5子地层的厚度为4m，A6子地层的厚度为6m，A7子地层的厚度为6m。
采集并获取该近地表地层的微测井资料，定量分析微测井资料，以得到微测井资料中每个深度点所对应的频率、速度值、能量和波形特征。找出这些频率的拐点位置以及这些拐点位置对应的地层深度值；找出这些速度值的拐点位置以及这些拐点位置对应的地层深度值；找出这些能量的拐点位置以及这些拐点位置对应的地层深度值；同时，找出这些波形特征的拐点位置以及这些拐点位置对应的地层深度值，获得频率属性、速度属性、能量属性和波形特征属性的分层结果。
其中，频率属性的各子地层及其地层深度值分别为：B1子地层0～4m、B2子地层4～8m、B3子地层8～17m、B4子地层17～25m、B5子地层25～30m。也就是说，B1子地层的厚度为4m，B2子地层的厚度为4m，B3子地层的厚度为9m、B4子地层的厚度为8m、B5子地层的厚度为5m。速度属性的各子地层及其地层深度值分别为：C1子地层0～5m、C2子地层5～17m、C3子地层17～24m、C4子地层24～30m。也就是说，C1子地层的厚度为5m，C2子地层的厚度为12m，C3子地层的厚度为7m、C4子地层的厚度为6m。能量属性的各子地层及其地层深度值分别为：D1子地层0～5m、D2子地层5～14m、D3子地层14～17m、D4子地层17～26m、D5子地层26～30m。也就是说，D1子地层的厚度为5m，D2子地层的厚度为9m，D3子地层的厚度为3m、D4子地层的厚度为9m、D5子地层的厚度为4m。波形特征属性的各子地层及其地层深度值分别为：E1子地层0～4m、E2子地层4～16m、E3子地层16～24m、E4子地层24～30m。也就是说，E1子地层的厚度为4m，E2子地层的厚度为12m，E3子地层的厚度为8m、E4子地层的厚度为6m。
对上述岩性录井分层结果所获得的A1子地层至A7子地层、上述微测井频率属性分层结果所获得的B1子地层至B5子地层、上述微测井速度属性分层结果所获得的C1子地层至C4子地层、以及上述微测井能量属性分层结果所获得的D1子地层至D5子地层进行综合分析。
根据微测井频率属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据B3子地层的地层深度值8～17m，考虑将A3子地层至A5子地层进行合并。A3子地层、A4子地层、A5子地层中只有A4子地层为泥岩，且A4子地层的地层厚度为2m，小于5m，故而A4子地层可以被合并而不用单独分层。另外，A3子地层至A5子地层合并后地层总厚度为10m，不大于10m。因此，可将A3子地层至A5子地层合并。
根据微测井速度属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据C2子地层的地层深度值5～17m，考虑将A2子地层至A5子地层进行合并。A2子地层、A3子地层、A4子地层、A5子地层中A2子地层、A4子地层为泥岩，且A2子地层、A4子地层的地层厚度都小于5m，故而A2子地层、A4子地层不用单独分层。然而，考虑到A2子地层至A5子地层合并后地层总厚度为12m，大于10m。因此，不建议将A2子地层至A5子地层进行合并。然而，可将A2子地层至A4子地层合并，或者将A3子地层至A5子地层合并。
根据微测井能量属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据D2子地层的地层深度值5～14m，考虑将A2子地层至A4子地层进行合并。A2子地层、A3子地层、A4子地层中A2子地层为泥岩、A4子地层为泥岩，且A2子地层、A4子地层的地层厚度分别为3m、2m，都小于5m，故而A2子地层、A4子地层不用单独分层。另外，考虑到A2子地层至A4子地层合并后地层总厚度为9m，小于10m。因此，可将A2子地层至A4子地层合并。
根据微测井波形特征属性分层结果对岩性录井分层结果中的子层进行合并。例如，根据E2子地层的地层深度值4～16m，考虑将A2子地层至A4子地层进行合并。A2子地层、A3子地层、A4子地层中A2子地层、A4子地层为泥岩，且A2子地层、A4子地层的地层厚度都小于5m，故而A2子地层、A4子地层不用单独分层。而且，A2子地层至A4子地层合并后地层总厚度为9m，不大于10m。因此，可将A2子地层至A4子地层进行合并。根据E3子地层的地层深度值16～24m，考虑将A5子地层和A6子地层进行合并。A5子地层、A6子地层中A6子地层为砂岩，然而，A6子地层的地层厚度并不小于5m，故而A6子地层需单独分层。
综合考虑上述合并结果，优选，将A2子地层至A4子地层合并，从而得到最终的分层结果为：0～5m、5～14m、14～17m、17～26m、26～30m。另外，也可将A3子地层至A5子地层合并，从而得到最终的分层结果为：0～5m、5～8m、8～17m、17～25m、25～30m。
综上所述，本发明能够增加计算Q值地层的厚度，避免薄地层计算Q值难的问题；同时，由于是将地震波频率、波形特征、能量、速度和岩性性质差别不大的多个薄地层划分为一层，各薄层的Q值差异也较小，有利于保证综合分层后计算Q值的精度。
尽管上面已经结合附图和示例性实施例描述了本发明，但是本领域普通技术人员应该清楚，在不脱离权利要求的精神和范围的情况下，可以对上述实施例进行各种修改。
Claims (7)
1.一种用于获得地层吸收衰减品质因子的地层划分方法，其特征在于，所述方法包括以下步骤：
利用岩性录井资料对目标地层进行分层，以得到岩性录井分层结果，利用微测井资料对所述目标地层进行分层，以得到微测井分层结果；
综合分析岩性录井分层结果和微测井分层结果，以得到综合分层结果，
其中，所述利用微测井资料对目标地层进行分层的步骤包括：定量分析微测井资料，以得到所述微测井资料中每个深度点所对应的属性值；找出所述属性值内的拐点位置以及该拐点位置对应的地层深度值，并依据该地层深度值对目标地层进行分层；
所述利用岩性录井资料对目标地层进行分层的步骤包括：分析目标地层的岩性录井资料，统计出各岩性的岩性特征及其对应的深度值，并依据所述对应的深度值对目标地层进行分层。
2.根据权利要求1所述的用于获得地层吸收衰减品质因子的地层划分方法，所述综合分析岩性录井分层结果和微测井分层结果的步骤包括：根据微测井分层结果对岩性录井分层结果中的各子地层进行合并，并且岩性录井分层结果中的厚度大于5米且为风化层、泥岩、砂岩或灰岩的子地层不应被合并，而且合并所获得的地层的厚度不应大于10米。
3.根据权利要求1所述的用于获得地层吸收衰减品质因子的地层划分方法，其中，所述属性值包括频率、波形特征、能量和速度中的一种或两种以上。
4.根据权利要求1所述的用于获得地层吸收衰减品质因子的地层划分方法，其中，所述目标地层为近地表地层，或者为具有岩性复杂变化特征的地层，所述岩性复杂变化特征是指岩性在小于40米内变化2次或3次以上。
5.根据权利要求1所述的用于获得地层吸收衰减品质因子的地层划分方法，其中，所述微测井资料由单井微测井方法或双井微测井方法采集得到。
6.根据权利要求1所述的用于获得地层吸收衰减品质因子的地层划分方法，其中，所述微测井资料通过井中激发地表接收或井中激发井中接收的微测井方法得到。
7.一种获得地层吸收衰减品质因子的方法，其特征在于，所述方法根据权利要求1至6中任意一项所述的地层划分方法的综合分层结果来计算目标地层的吸收衰减品质因子。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610785617.8A CN106291714A (zh) | 2016-08-31 | 2016-08-31 | 获得地层吸收衰减品质因子的方法及其地层划分方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610785617.8A CN106291714A (zh) | 2016-08-31 | 2016-08-31 | 获得地层吸收衰减品质因子的方法及其地层划分方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106291714A true CN106291714A (zh) | 2017-01-04 |
Family
ID=57673061
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610785617.8A Pending CN106291714A (zh) | 2016-08-31 | 2016-08-31 | 获得地层吸收衰减品质因子的方法及其地层划分方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106291714A (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106646616A (zh) * | 2017-01-19 | 2017-05-10 | 中国石油大学(北京) | 一种补偿地层吸收衰减的方法及装置 |
| CN106908838A (zh) * | 2017-03-15 | 2017-06-30 | 徐诗薇 | 构建目标区域地层非弹性衰减品质因子立体模型的方法 |
| CN107389794A (zh) * | 2017-08-14 | 2017-11-24 | 中国石油大学(北京) | 测量岩石衰减系数的方法与系统 |
| CN110673211A (zh) * | 2019-10-13 | 2020-01-10 | 东北石油大学 | 一种基于测井与地震数据的品质因子建模方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2196620A1 (en) * | 2008-12-15 | 2010-06-16 | Services Pétroliers Schlumberger | A micro-logging system and method |
| CN102109617A (zh) * | 2010-12-15 | 2011-06-29 | 大庆油田有限责任公司 | 用双井微测井技术测量近地表地层q值的方法 |
| CN104793247A (zh) * | 2015-03-12 | 2015-07-22 | 中国石油天然气集团公司 | 近地表速度建模方法及建模装置 |
-
2016
- 2016-08-31 CN CN201610785617.8A patent/CN106291714A/zh active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2196620A1 (en) * | 2008-12-15 | 2010-06-16 | Services Pétroliers Schlumberger | A micro-logging system and method |
| CN102109617A (zh) * | 2010-12-15 | 2011-06-29 | 大庆油田有限责任公司 | 用双井微测井技术测量近地表地层q值的方法 |
| CN104793247A (zh) * | 2015-03-12 | 2015-07-22 | 中国石油天然气集团公司 | 近地表速度建模方法及建模装置 |
Non-Patent Citations (2)
| Title |
|---|
| 尚新民: "准中沙漠区近地表特征及针对性勘探技术", 《西南石油大学学报（自然科学版）》 * |
| 李敏等: "近地表Q值模型建立实例", 《SPG/SEG北京2016国际地球物理会议》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106646616A (zh) * | 2017-01-19 | 2017-05-10 | 中国石油大学(北京) | 一种补偿地层吸收衰减的方法及装置 |
| CN106646616B (zh) * | 2017-01-19 | 2018-10-30 | 中国石油大学(北京) | 一种补偿地层吸收衰减的方法及装置 |
| CN106908838A (zh) * | 2017-03-15 | 2017-06-30 | 徐诗薇 | 构建目标区域地层非弹性衰减品质因子立体模型的方法 |
| CN107389794A (zh) * | 2017-08-14 | 2017-11-24 | 中国石油大学(北京) | 测量岩石衰减系数的方法与系统 |
| CN107389794B (zh) * | 2017-08-14 | 2020-01-10 | 中国石油大学(北京) | 测量岩石衰减系数的方法与系统 |
| CN110673211A (zh) * | 2019-10-13 | 2020-01-10 | 东北石油大学 | 一种基于测井与地震数据的品质因子建模方法 |
| CN110673211B (zh) * | 2019-10-13 | 2021-06-04 | 东北石油大学 | 一种基于测井与地震数据的品质因子建模方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101329405B (zh) | 一种简单的多参数地震反演方法 | |
| CN101551466B (zh) | 一种利用与偏移距有关的地震属性提高油气储层预测精度的方法 | |
| Yasin et al. | Estimation of petrophysical parameters from seismic inversion by combining particle swarm optimization and multilayer linear calculator | |
| CN101446645B (zh) | 一种利用地震流体阻抗进行流体确定的方法 | |
| US7069150B2 (en) | Method for optimizing migration fields using time slice analysis | |
| CN104237945B (zh) | 一种地震资料自适应高分辨处理方法 | |
| CN105938203B (zh) | 一种储层特性的检测方法及装置 | |
| CN104166161A (zh) | 一种基于各向异性的椭圆速度反演的裂缝预测方法及装置 | |
| CN106291714A (zh) | 获得地层吸收衰减品质因子的方法及其地层划分方法 | |
| US7668042B2 (en) | Method for obtaining a seismic wavelet using seismic data and well logs in a deviated well | |
| CN114609675A (zh) | 基于高频旋回对碳酸盐岩地层沉积微地貌的定量恢复方法 | |
| CN105510965B (zh) | 流体识别方法和装置 | |
| CN111983683B (zh) | 一种少井条件下盐下湖相灰岩储层预测方法和系统 | |
| CN104422955B (zh) | 一种利用旅行时变化量进行各向异性参数提取的方法 | |
| CN106443791B (zh) | 求取倾斜地层或各向异性地层横波剩余静校正量的方法 | |
| CN103278852A (zh) | 利用地震数据体波形结构特征模型预测油气的方法 | |
| CN111142159B (zh) | 中深层异常体地震速度建模方法及系统 | |
| CN111077577B (zh) | 井地联合储层描述方法及装置 | |
| CN109143368B (zh) | 山前带地震去噪方法 | |
| Cappetti et al. | A new deep exploration program and preliminary results of a 3D seismic survey in the Larderello-Travale geothermal field (Italy) | |
| CN102478665B (zh) | 一种确定地震波入射角和振幅的方法 | |
| Paris et al. | Predicting reservoir quality in the Bakken Formation, North Dakota, using petrophysics and 3C seismic data | |
| CN112305594B (zh) | 非均质储层的油气分布确定方法及系统 | |
| CN114462703A (zh) | 声波参数曲线预测、测井曲线预测方法及电子设备 | |
| CN106353796A (zh) | 一种提高地面地震资料分辨率方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan Applicant after: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180402 Address after: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei Applicant after: BGP INC., CHINA NATIONAL PETROLEUM Corp. Address before: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd. |
|
| TA01 | Transfer of patent application right | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170104 |
|
| WD01 | Invention patent application deemed withdrawn after publication |