[go: up one dir, main page]

CN116136606A - Fracture logging identification method and system under oil-based mud conditions - Google Patents

Fracture logging identification method and system under oil-based mud conditions Download PDF

Info

Publication number
CN116136606A
CN116136606A CN202111357624.5A CN202111357624A CN116136606A CN 116136606 A CN116136606 A CN 116136606A CN 202111357624 A CN202111357624 A CN 202111357624A CN 116136606 A CN116136606 A CN 116136606A
Authority
CN
China
Prior art keywords
wave impedance
impedance value
crack
rock
corrected
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
CN202111357624.5A
Other languages
Chinese (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 CN202111357624.5A priority Critical patent/CN116136606A/en
Publication of CN116136606A publication Critical patent/CN116136606A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/53Statics correction, e.g. weathering layer or transformation to a datum
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6226Impedance
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/64Geostructures, e.g. in 3D data cubes
    • G01V2210/642Faults
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

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)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a crack logging identification method and a crack logging identification system under an oil-based mud condition, which can improve the precision of crack logging identification under the oil-based mud condition, effectively evaluate the development degree of stratum cracks and judge the filling property of the cracks. S1, calculating a rock wave impedance value according to a density curve and a longitudinal wave acoustic wave time difference curve based on conventional logging data of a stratum to be analyzed; s2, correcting the calculated rock wave impedance value according to a clay content curve in conventional logging data, calculating according to a preset correction formula to obtain a corrected wave impedance value, and drawing a corrected wave impedance curve graph; and S3, judging whether the crack of the stratum to be analyzed is developed or not based on the calculated corrected wave impedance value, if the corrected wave impedance value is not greater than a preset wave impedance threshold value, indicating that the crack of the corresponding stratum is a developed crack, otherwise, indicating that the crack is an undeveloped crack.

Description

一种油基泥浆条件下的裂缝测井识别方法及系统Fracture logging identification method and system under oil-based mud conditions

技术领域technical field

本发明涉及致密砂岩裂缝测井评价方法技术领域,具体为一种油基泥浆条件下的裂缝测井识别方法及系统。The invention relates to the technical field of tight sandstone fracture logging evaluation methods, in particular to a fracture logging identification method and system under oil-based mud conditions.

背景技术Background technique

目前,针对致密油气储集层,裂缝的发育与否决定了致密油气的富集与高产。但由于取心和薄片资料较少,利用测井资料尤其是成像测井资料实现对裂缝的测井判别显得尤为重要。At present, for tight oil and gas reservoirs, the development of fractures determines the enrichment and high production of tight oil and gas. However, due to the lack of coring and thin section data, it is particularly important to use logging data, especially imaging logging data, to realize logging discrimination of fractures.

然而,油基泥浆条件下,微电阻率成像测井由于地层受油基钻井液影响,导电性失灵,很难对裂缝进行精细识别和评价,目前针对油基泥浆条件下的裂缝的测井识别主要存在以下问题:However, under the condition of oil-based mud, micro-resistivity imaging logging is affected by oil-based drilling fluid and the electrical conductivity of the formation fails, so it is difficult to finely identify and evaluate fractures. Currently, the logging identification of fractures under oil-based mud conditions There are mainly the following problems:

1、精度限制,油基泥浆条件下微电阻率成像测井难以准确识别出裂缝,与水基泥浆条件下相比,裂缝形迹识别较困难;1. Due to the limitation of accuracy, it is difficult to accurately identify fractures in micro-resistivity imaging logging under oil-based mud conditions. Compared with water-based mud conditions, it is more difficult to identify fracture traces;

2、由于受油基泥浆影响,微电阻率成像测井上难以识别裂缝的充填性。2. Due to the influence of oil-based mud, it is difficult to identify the filling properties of fractures on micro-resistivity imaging logging.

发明内容Contents of the invention

为了解决现有技术中存在的问题,本发明提供一种油基泥浆条件下的裂缝测井识别方法及系统,可以提高油基泥浆条件下裂缝测井识别的精度,对地层裂缝发育程度进行有效评价,同时对于裂缝的充填性进行判断。In order to solve the problems existing in the prior art, the present invention provides a fracture logging identification method and system under oil-based mud conditions, which can improve the accuracy of fracture logging identification under oil-based mud conditions, and effectively monitor the development degree of formation fractures. Evaluation, and at the same time judge the filling performance of the fracture.

为实现上述目的,本发明提供如下技术方案:To achieve the above object, the present invention provides the following technical solutions:

一种油基泥浆条件下的裂缝测井识别方法,包括如下步骤,A fracture logging identification method under oil-based mud conditions, comprising the following steps,

S1、基于待分析地层的常规测井资料,根据密度曲线及纵波声波时差曲线计算出岩石波阻抗值;S1. Based on the conventional logging data of the formation to be analyzed, the rock wave impedance value is calculated according to the density curve and the compressional wave acoustic wave time difference curve;

S2、根据常规测井资料中的泥质含量曲线对计算得到的岩石波阻抗值进行校正,根据预设的校正公式计算获得校正波阻抗值,并绘制校正波阻抗曲线图;S2. Correct the calculated rock wave impedance value according to the shale content curve in the conventional logging data, calculate and obtain the corrected wave impedance value according to the preset correction formula, and draw the corrected wave impedance curve;

S3、基于计算得到的校正波阻抗值判断待分析地层的裂缝是否发育,若校正波阻抗值不大于预设的波阻抗阈值,则表示对应地层的裂缝为发育裂缝,否则为未发育裂缝。S3. Based on the calculated corrected wave impedance value, it is judged whether the fractures in the formation to be analyzed are developed. If the corrected wave impedance value is not greater than the preset wave impedance threshold, it means that the fractures in the corresponding formation are developed fractures, otherwise they are not developed fractures.

优选地,在S1中还包括对于计算得到的异常的岩石波阻抗值进行剥离;Preferably, S1 also includes stripping the calculated abnormal rock wave impedance value;

其中,异常的判断标准为计算的岩石波阻抗值偏离正常岩石物理性质对应岩石波阻抗值的预设范围。Among them, the abnormal judgment standard is that the calculated rock wave impedance value deviates from the preset range of the rock wave impedance value corresponding to the normal rock physical properties.

优选地,在S2中,所述校正包括,Preferably, in S2, the correction includes,

选取岩心中不同泥质含量位置处的多个未发育裂缝,计算泥质含量变化时对应的岩石波阻抗值的变化值;Select multiple undeveloped fractures at different shale content positions in the core, and calculate the change value of rock wave impedance corresponding to the change of shale content;

对比泥质含量变化与岩石波阻抗值变化的关系,得到校正公式为,Comparing the relationship between the change of shale content and the change of rock wave impedance value, the correction formula is obtained as,

当埋深≤5000m时,校正波阻抗值Z’=Z+1.3VSH,其中,VSH为泥质含量,Z为岩石波阻抗值;When the buried depth is ≤5000m, the corrected wave impedance value Z’=Z+1.3VSH, where VSH is the mud content and Z is the rock wave impedance value;

当埋深>5000m时,校正波阻抗值Z’=Z-1.3VSH。When the buried depth is >5000m, the corrected wave impedance value Z'=Z-1.3VSH.

优选地,基于多井资料,设置预设的波阻抗阈值为1300Kg/cm3·s。Preferably, based on multi-well data, the preset acoustic impedance threshold is set to 1300Kg/cm 3 ·s.

优选地,还包括根据校正波阻抗值定性判别待分析井段的裂缝发育程度;Preferably, it also includes qualitatively judging the fracture development degree of the well section to be analyzed according to the corrected wave impedance value;

根据校正波阻抗曲线图划分井段的裂缝发育层段;According to the corrected wave impedance curve, the fracture-developed layers of the well section are divided;

裂缝发育程度为裂缝发育层段的厚度与整个井段的深度的比,比值越大,则该井段的裂缝发育程度越高。The degree of fracture development is the ratio of the thickness of the fracture-developed interval to the depth of the entire well section. The larger the ratio, the higher the fracture development degree of the well section.

优选地,还包括获取油基泥浆成像测井图像,并与校正波阻抗曲线图进行对比建立裂缝识别图版;Preferably, it also includes acquiring the oil-based mud imaging logging image, and comparing it with the corrected wave impedance curve to establish a fracture identification plate;

其中,对于在油基泥浆成像测井图像上能够直接识别的裂缝,基于计算得到的校正波阻抗值与预设的波阻抗阈值对比,判断待分析地层的裂缝是否被充填,若校正波阻抗值大于波阻抗阈值,则表示对应地层的裂缝被充填;Among them, for the fractures that can be directly identified on the oil-based mud imaging logging image, based on the comparison between the calculated corrected wave impedance value and the preset wave impedance threshold, it is judged whether the fractures in the formation to be analyzed are filled. If the corrected wave impedance value If it is greater than the wave impedance threshold, it means that the fractures in the corresponding formation are filled;

对于在油基泥浆成像测井图像上无法直接识别出的裂缝,继续进行S3步骤判断裂缝是否发育。For the fractures that cannot be directly identified on the oil-based mud imaging logging image, proceed to step S3 to determine whether the fractures are developed.

优选地,在S1中岩石波阻抗值的计算公式为,Preferably, the calculation formula of the rock wave impedance value in S1 is,

岩石波阻抗值Z=105×DEN÷3.28DTC,Rock wave impedance value Z=10 5 ×DEN÷3.28DTC,

其中,DEN为密度值,DTC为纵波声波时差。Among them, DEN is the density value, and DTC is the longitudinal wave sound wave time difference.

一种油基泥浆条件下的裂缝测井识别系统,包括,A fracture logging identification system under oil-based mud conditions, comprising:

计算模块,用于基于待分析地层的常规测井资料,根据密度曲线及纵波声波时差曲线计算出岩石波阻抗值;The calculation module is used to calculate the rock wave impedance value based on the conventional logging data of the stratum to be analyzed according to the density curve and the compressional wave acoustic time difference curve;

校正模块,用于根据常规测井资料中的泥质含量曲线和预设的校正公式计算获得校正波阻抗值,并绘制校正波阻抗曲线图;The correction module is used to calculate and obtain the corrected wave impedance value according to the shale content curve in the conventional logging data and the preset correction formula, and draw the corrected wave impedance curve;

比较模块,用于基于计算得到的校正波阻抗值判断待分析地层的裂缝是否发育,若校正波阻抗值小于预设的波阻抗阈值,则表示对应地层的裂缝为发育裂缝。The comparison module is used to determine whether the fractures in the formation to be analyzed are developed based on the calculated corrected wave impedance value. If the corrected wave impedance value is less than the preset wave impedance threshold, it indicates that the fractures in the corresponding formation are developed fractures.

优选地,还包括成像模块,用于获取油基泥浆成像测井图像,并与校正波阻抗曲线图进行对比建立裂缝识别图版。Preferably, an imaging module is also included, which is used to obtain the oil-based mud imaging logging image, and compare it with the corrected wave impedance curve to establish a fracture identification plate.

与现有技术相比,本发明具有以下有益效果:Compared with the prior art, the present invention has the following beneficial effects:

本发明提供一种油基泥浆条件下的裂缝测井识别方法,包括基于常规测井资料中的密度曲线及纵波声波时差曲线计算出岩石波阻抗值;根据地层埋深不同,利用测井计算得到的地层泥质含量来校正所计算出的波阻抗值,从而能够排除砂泥岩岩性对波阻抗值的影响,并通过校正波阻抗值与波阻抗阈值进行对比判断待分析地层的裂缝发育状况,可以提高油基泥浆条件下裂缝测井识别的精度,能够在缺乏成像资料的情况下,对地层裂缝发育程度进行有效评价。The invention provides a fracture logging identification method under the condition of oil-based mud, including calculating the rock wave impedance value based on the density curve and the longitudinal wave acoustic wave time difference curve in the conventional logging data; The shale content of the formation is used to correct the calculated wave impedance value, so that the influence of the sandstone lithology on the wave impedance value can be eliminated, and the fracture development status of the formation to be analyzed can be judged by comparing the corrected wave impedance value with the wave impedance threshold value. It can improve the accuracy of fracture logging identification under oil-based mud conditions, and can effectively evaluate the development degree of formation fractures in the absence of imaging data.

本发明能够对油基泥浆背景下的裂缝进行更精确的测井识别。一是本发明利用声电结合的方法结合成像测井图像以及岩心资料能实现油基泥浆背景下的高精度的裂缝测井识别;二是本发明也能适用于水基泥浆条件下的裂缝识别,能够在缺乏成像资料的情况下,对地层裂缝发育程度进行有效评价;三是本发明可以在一定程度上识别出油基泥浆背景下裂缝充填性。The invention can perform more accurate logging identification on fractures under the background of oil-based mud. One is that the present invention can realize high-precision fracture logging identification under the background of oil-based mud by combining acoustic and electric methods with imaging logging images and core data; the other is that the present invention can also be applied to fracture identification under water-based mud conditions , in the absence of imaging data, the development degree of formation fractures can be effectively evaluated; third, the present invention can identify the fracture filling property in the background of oil-based mud to a certain extent.

附图说明Description of drawings

图1是本发明识别方法流程图;Fig. 1 is the flow chart of identification method of the present invention;

图2是本发明实施例中致密砂岩裂缝岩心识别图版;Fig. 2 is a tight sandstone fracture core identification plate in the embodiment of the present invention;

图3是本发明实施例中单井上岩石波阻抗值曲线图;Fig. 3 is a curve diagram of rock wave impedance value on a single well in an embodiment of the present invention;

图4是本发明实施例中对应地层发育裂缝的校正波阻抗值曲线图;Fig. 4 is a curve diagram of the corrected wave impedance value corresponding to formation fractures in the embodiment of the present invention;

图5是本发明实施例中油基泥浆背景下裂缝识别图版;Fig. 5 is a chart for identifying cracks under the background of oil-based mud in an embodiment of the present invention;

图6是本发明实施例中成像图像上识别出的裂缝形迹图;Fig. 6 is a diagram of crack traces identified on the imaging image in the embodiment of the present invention;

图7是本发明实施例中油基泥浆单井裂缝测井识别结果图。Fig. 7 is a diagram showing the results of oil-based mud single well fracture logging identification in the embodiment of the present invention.

图中,a、水平张裂缝;b、低角度裂缝;c、高角度裂缝;d、多组高角度裂缝;e、网状裂缝;f、石膏充填低角度裂缝;g、方解石充填高角度裂缝;h、裂缝溶蚀;i、方解石充填与裂缝溶蚀。In the figure, a, horizontal tensile fractures; b, low-angle fractures; c, high-angle fractures; d, multiple groups of high-angle fractures; e, network fractures; f, low-angle fractures filled with gypsum; g, high-angle fractures filled with calcite ; h, fracture dissolution; i, calcite filling and fracture dissolution.

具体实施方式Detailed ways

下面结合具体的实施例对本发明做进一步的详细说明,所述是对本发明的解释而不是限定。The present invention will be further described in detail below in conjunction with specific embodiments, which are explanations of the present invention rather than limitations.

本发明一种油基泥浆条件下的裂缝测井识别方法,如图1所示,包括如下步骤,A fracture logging identification method under the oil-based mud condition of the present invention, as shown in Figure 1, comprises the following steps,

S1、基于待分析地层的常规测井资料,根据密度曲线及纵波声波时差曲线计算出岩石波阻抗值;S1. Based on the conventional logging data of the formation to be analyzed, the rock wave impedance value is calculated according to the density curve and the compressional wave acoustic wave time difference curve;

S2、根据常规测井资料中的泥质含量曲线对计算得到的岩石波阻抗值进行校正,根据预设的校正公式计算获得校正波阻抗值,并绘制校正波阻抗曲线图;S2. Correct the calculated rock wave impedance value according to the shale content curve in the conventional logging data, calculate and obtain the corrected wave impedance value according to the preset correction formula, and draw the corrected wave impedance curve;

S3、基于计算得到的校正波阻抗值判断待分析地层的裂缝是否发育,若校正波阻抗值小于预设的波阻抗阈值,则表示对应地层的裂缝为发育裂缝,否则为未发育裂缝。S3. Based on the calculated corrected wave impedance value, it is judged whether the fractures in the formation to be analyzed are developed. If the corrected wave impedance value is less than the preset wave impedance threshold, it means that the fractures in the corresponding formation are developed fractures, otherwise they are not developed fractures.

本发明提供一种油基泥浆条件下的裂缝测井识别方法,包括基于常规测井资料中的密度曲线及纵波声波时差曲线计算出岩石波阻抗值;根据地层埋深不同,建立以5000m为界限的利用测井计算得到的地层泥质含量来校正所计算出的波阻抗值,从而能够排除砂泥岩岩性对波阻抗值的影响,并通过校正波阻抗值与波阻抗阈值进行对比判断待分析地层的裂缝发育状况,可以提高油基泥浆条件下裂缝测井识别的精度,能够在缺乏成像资料的情况下,对地层裂缝发育程度进行有效评价。The invention provides a fracture logging identification method under the condition of oil-based mud, which includes calculating the rock wave impedance value based on the density curve and the longitudinal wave acoustic wave time difference curve in the conventional logging data; according to the difference in the buried depth of the formation, a limit of 5000m is established Correct the calculated wave impedance value by using the formation shale content calculated by well logging, so as to eliminate the influence of sandstone and mudstone lithology on the wave impedance value, and compare the corrected wave impedance value with the wave impedance threshold value to be analyzed The development status of formation fractures can improve the accuracy of fracture logging identification under oil-based mud conditions, and can effectively evaluate the development degree of formation fractures in the absence of imaging data.

进一步地,在S1中还包括对于计算得到的异常的岩石波阻抗值进行剥离,通过剥离出测井值异常段,可以使得判断结果更加准确,减小误差;Further, S1 also includes stripping the calculated abnormal rock wave impedance value, and by stripping out the abnormal section of the logging value, the judgment result can be made more accurate and the error can be reduced;

其中,异常的判断标准为计算的岩石波阻抗值偏离正常岩石物理性质对应岩石波阻抗值的预设范围,一般为与正常范围相差一两个数量级。本实施例中,在S2中,所述校正包括,Among them, the criterion for judging abnormality is that the calculated rock wave impedance value deviates from the preset range of rock wave impedance value corresponding to the normal rock physical properties, which is generally one or two orders of magnitude different from the normal range. In this embodiment, in S2, the correction includes,

选取岩心中不同泥质含量位置处的多个未发育裂缝,计算泥质含量变化时对应的岩石波阻抗值的变化值;Select multiple undeveloped fractures at different shale content positions in the core, and calculate the change value of rock wave impedance corresponding to the change of shale content;

对比泥质含量变化与岩石波阻抗值变化的关系,得到校正公式为,Comparing the relationship between the change of shale content and the change of rock wave impedance value, the correction formula is obtained as,

当埋深≤5000m时,校正波阻抗值Z’=Z+1.3VSH,其中,VSH为泥质含量,Z为岩石波阻抗值;When the buried depth is ≤5000m, the corrected wave impedance value Z’=Z+1.3VSH, where VSH is the mud content and Z is the rock wave impedance value;

当埋深>5000m时,校正波阻抗值Z’=Z-1.3VSH。When the buried depth is >5000m, the corrected wave impedance value Z'=Z-1.3VSH.

本发明通过泥质含量获得校正波阻抗值,排除砂泥岩岩性对波阻抗值的影响。The invention obtains the corrected wave impedance value through the shale content, and eliminates the influence of sand and mudstone lithology on the wave impedance value.

本实施例中,基于多井资料,设置预设的波阻抗阈值为1300Kg/cm3·s。In this embodiment, based on multi-well data, the preset acoustic impedance threshold is set to 1300Kg/cm 3 ·s.

进一步地,还包括根据校正波阻抗值定性判别待分析井段的裂缝发育程度;Further, it also includes qualitatively judging the fracture development degree of the well section to be analyzed according to the corrected wave impedance value;

根据校正波阻抗曲线图划分井段的裂缝发育层段;According to the corrected wave impedance curve, the fracture-developed layers of the well section are divided;

裂缝发育程度为裂缝发育层段的厚度与整个井段的深度的比,比值越大,则该井段的裂缝发育程度越高。The degree of fracture development is the ratio of the thickness of the fracture-developed interval to the depth of the entire well section. The larger the ratio, the higher the fracture development degree of the well section.

进一步地,还包括获取油基泥浆成像测井图像,并与校正波阻抗曲线图进行对比建立裂缝识别图版;Further, it also includes obtaining the oil-based mud imaging logging image, and comparing it with the corrected wave impedance curve to establish a fracture identification plate;

其中,对于在油基泥浆成像测井图像上能够直接识别的裂缝,基于计算得到的校正波阻抗值与预设的波阻抗阈值对比,判断待分析地层的裂缝是否被充填,若校正波阻抗值大于波阻抗阈值,则表示对应地层的裂缝被充填;Among them, for the fractures that can be directly identified on the oil-based mud imaging logging image, based on the comparison between the calculated corrected wave impedance value and the preset wave impedance threshold, it is judged whether the fractures in the formation to be analyzed are filled. If the corrected wave impedance value If it is greater than the wave impedance threshold, it means that the fractures in the corresponding formation are filled;

对于在油基泥浆成像测井图像上无法直接识别出的裂缝,继续进行S3步骤判断裂缝是否发育。For the fractures that cannot be directly identified on the oil-based mud imaging logging image, proceed to step S3 to determine whether the fractures are developed.

本发明还提供一种综合利用岩石波阻抗计算与油基泥浆成像测井图像对比建立油基泥浆背景下裂缝识别图版的方法,通过声电结合的方法能够更精确识别油基泥浆背景下的裂缝等步骤。其中,基于常规测井获取岩石密度值及纵波声波时差值数据能计算出岩石波阻抗值(岩石波阻抗=岩石密度x纵波速度);由于波阻抗值会受到岩性影响,因此针对致密砂岩地层,需要根据泥质含量值校正泥质对波阻抗的影响,由此得出校正后的波阻抗值,并结合多井资料利用校正波阻抗曲线图划分裂缝发育层带;最后通过常规测井声电结合—波阻抗法与成像测井图像以及岩心对比即可实现油基泥浆背景下高精度裂缝测井识别。通过上述步骤,本发明能够对油基泥浆背景下的裂缝进行更精确的测井识别。一是本发明利用声电结合的方法结合成像测井图像以及岩心资料能实现油基泥浆背景下的高精度的裂缝测井识别;二是本发明也能适用于水基泥浆条件下的裂缝识别,能够在缺乏成像资料的情况下,对地层裂缝发育程度进行有效评价;三是本发明可以在一定程度上识别出油基泥浆背景下裂缝充填性。The present invention also provides a method for comprehensively utilizing the calculation of rock wave impedance and the comparison of oil-based mud imaging and logging images to establish a fracture identification plate under the background of oil-based mud, and the method of combining sound and electricity can more accurately identify cracks under the background of oil-based mud and other steps. Among them, the rock impedance value can be calculated based on the rock density value and P-wave acoustic time difference data obtained by conventional logging (rock wave impedance = rock density x P-wave velocity); since the wave impedance value will be affected by lithology, the tight sandstone In the formation, it is necessary to correct the influence of shale on wave impedance according to the shale content value, so as to obtain the corrected wave impedance value, and combine multi-well data with the corrected wave impedance curve to divide the fracture-developed zone; finally, through conventional logging Acoustic-electric combination—wave impedance method, imaging logging image and core comparison can realize high-precision fracture logging identification in the background of oil-based mud. Through the above steps, the present invention can perform more accurate logging identification on fractures in the background of oil-based mud. One is that the present invention can realize high-precision fracture logging identification under the background of oil-based mud by combining acoustic and electric methods with imaging logging images and core data; the other is that the present invention can also be applied to fracture identification under water-based mud conditions , in the absence of imaging data, the development degree of formation fractures can be effectively evaluated; third, the present invention can identify the fracture filling property in the background of oil-based mud to a certain extent.

本实施例中,在S1中岩石波阻抗值的计算公式为,In this embodiment, the calculation formula of the rock wave impedance value in S1 is:

岩石波阻抗值Z=105×DEN÷3.28DTC,Rock wave impedance value Z=10 5 ×DEN÷3.28DTC,

其中,DEN为密度值,DTC为纵波声波时差。Among them, DEN is the density value, and DTC is the longitudinal wave sound wave time difference.

一种油基泥浆条件下的裂缝测井识别系统,包括,A fracture logging identification system under oil-based mud conditions, comprising:

计算模块,用于基于待分析地层的常规测井资料,根据密度曲线及纵波声波时差曲线计算出岩石波阻抗值;The calculation module is used to calculate the rock wave impedance value based on the conventional logging data of the stratum to be analyzed according to the density curve and the compressional wave acoustic time difference curve;

校正模块,用于根据常规测井资料中的泥质含量曲线和预设的校正公式计算获得校正波阻抗值,并绘制校正波阻抗曲线图;The correction module is used to calculate and obtain the corrected wave impedance value according to the shale content curve in the conventional logging data and the preset correction formula, and draw the corrected wave impedance curve;

比较模块,用于基于计算得到的校正波阻抗值判断待分析地层的裂缝是否发育,若校正波阻抗值小于预设的波阻抗阈值,则表示对应地层的裂缝为发育裂缝。The comparison module is used to determine whether the fractures in the formation to be analyzed are developed based on the calculated corrected wave impedance value. If the corrected wave impedance value is less than the preset wave impedance threshold, it indicates that the fractures in the corresponding formation are developed fractures.

其中,还包括成像模块,用于获取油基泥浆成像测井图像,并与校正波阻抗曲线图进行对比建立裂缝识别图版。Among them, an imaging module is also included, which is used to obtain the oil-based mud imaging logging image, and compare it with the corrected wave impedance curve to establish a fracture identification plate.

进一步地,本发明所提供的识别方法可以在Techlog或者Resform软件平台中实现,但不限于此。Further, the identification method provided by the present invention can be implemented in Techlog or Reform software platform, but not limited thereto.

实施例Example

本实施例中,通过在Resform中实现,其具体实施方式如下:In this embodiment, by implementing in Reform, its specific implementation is as follows:

(1)通过岩心观察获取致密砂岩裂缝发育特征,包括裂缝产状、裂缝形迹、裂缝充填性等。(1) Obtain the fracture development characteristics of tight sandstone through core observation, including fracture occurrence, fracture trace, fracture filling, etc.

①如图2所示,裂缝从产状上可以分为高角度裂缝、斜交缝、低角度裂缝及近水平裂缝;① As shown in Fig. 2, fractures can be divided into high-angle fractures, oblique fractures, low-angle fractures and near-horizontal fractures in terms of occurrence;

②由于受不同充填物充填,裂缝呈现出不同的特征,按照充填性裂缝可以分为未充填裂缝和充填裂缝。②Because of being filled with different fillers, the fractures present different characteristics. According to the filling properties, fractures can be divided into unfilled fractures and filled fractures.

(2)利用常规测井资料中的密度曲线及纵波声波时差曲线计算出岩石波阻抗值,岩石波阻抗Z=岩石密度×纵波速度。(2) Calculate the rock wave impedance value by using the density curve and P-wave acoustic time-difference curve in the conventional logging data, and the rock wave impedance Z=rock density×P-wave velocity.

其中,利用常规测井资料的实际计算公式为Z=105×DEN÷3.28DTC;Among them, the actual calculation formula using conventional logging data is Z=10 5 ×DEN÷3.28DTC;

Z为波阻抗,单位是Kg/cm3·s,DEN为密度值,单位是g/cm3,DTC为纵波声波时差,单位是μs/ft。Z is the wave impedance, the unit is Kg/cm 3 ·s, DEN is the density value, the unit is g/cm 3 , DTC is the longitudinal sound wave time difference, the unit is μs/ft.

其中,对于测井值异常段,即测井曲线值偏离正常岩石物理性质对应测井值,通常与正常范围相差一两个数量级,应进行剥离,如图3所示。Among them, for the section with abnormal logging value, that is, the logging curve value deviates from the corresponding logging value of normal petrophysical properties, and usually differs from the normal range by one or two orders of magnitude, it should be stripped, as shown in Figure 3.

(3)基于常规测井资料得到的泥质含量曲线VSH对计算出的波阻抗进行校正,获取校正波阻抗值。(3) Correct the calculated wave impedance based on the shale content curve VSH obtained from conventional logging data to obtain the corrected wave impedance value.

①选取成像测井及岩心中多个未发育裂缝且泥质含量不同处,得到泥质含量变化下,波阻抗值变化值。① Select the imaging logging and several places in the core without fractures and with different shale content to obtain the change value of wave impedance under the change of shale content.

②通过对比多井未发育裂缝地层泥质含量与波阻抗关系发现,泥质含量每增加1%,波阻抗值变化约为1.3Kg/cm3·s;②By comparing the relationship between shale content and wave impedance in unfractured formations in multiple wells, it is found that for every 1% increase in shale content, the change of wave impedance is about 1.3Kg/cm 3 ·s;

受埋深影响,当埋深<5000m时,砂岩波阻抗>泥岩波阻抗,其校正公式:校正波阻抗Z’=Z+1.3VSH;Affected by the burial depth, when the burial depth is <5000m, the wave impedance of sandstone > the wave impedance of mudstone, and the correction formula is: corrected wave impedance Z’=Z+1.3VSH;

当埋深>5000m时,砂岩波阻抗<泥岩波阻抗,其校正公式:校正波阻抗Z’=Z-1.3VSH。When the buried depth is >5000m, sandstone wave impedance < mudstone wave impedance, the correction formula is: corrected wave impedance Z’=Z-1.3VSH.

③结合多井资料计算波阻抗对比,认为未发育裂缝时,地层波阻抗值>1300Kg/cm3·s,因此设置以1300Kg/cm3·s为界,当波阻抗值<1300Kg/cm3·s时,判定为地层发育裂缝,如图4所示。③Combined with the multi-well data to calculate the wave impedance comparison, it is considered that when no fractures are developed, the formation wave impedance value is >1300Kg/cm 3 ·s, so it is set to 1300Kg/cm 3 ·s as the boundary, when the wave impedance value is <1300Kg/cm 3 · s, it is determined that fractures are developed in the formation, as shown in Fig. 4.

(4)利用校正波阻抗曲线图与获取的油基泥浆成像测井图像以及岩心对比建立油基泥浆背景下裂缝识别图版。(4) Using the corrected wave impedance curve and the obtained oil-based mud imaging logging image and core comparison to establish a fracture identification chart under the oil-based mud background.

①根据校正波阻抗曲线图所划分出的裂缝发育层带与油基泥浆成像测井图像进行对比建立图版,如图5所示;① Comparing the fracture-developed zone zone divided by the corrected wave impedance curve with the oil-based mud imaging logging image to establish a chart, as shown in Fig. 5;

②油基泥浆成像测井图像上发现明显的裂缝形迹,可以直接在成像上拾取,并判断该裂缝的校正波阻抗值是否大于1300Kg/cm3·s,若大于则说明其为充填裂缝,由此可在一定程度上识别裂缝充填性,如图6所示。②Obvious fracture traces are found on the oil-based mud imaging logging image, which can be picked up directly on the image, and judge whether the corrected wave impedance value of the fracture is greater than 1300Kg/cm 3 ·s. If it is greater, it indicates that it is a filled fracture. This can identify the fracture filling property to a certain extent, as shown in Fig. 6.

③油基泥浆成像测井图像上的裂缝形迹模糊,可能出现漏识的现象,,或是对于缺少成像资料的井通过对比校正波阻抗曲线来判定是否能拾取该裂缝,若其校正波阻抗值小于1300Kg/cm3·s,则表明其为发育裂缝,并通过校正波阻抗曲线图划分井段的裂缝发育层段定性判别出该井段的裂缝发育程度,裂缝发育层段占整个井段的比值越大,裂缝越发育,如图7所示。③The traces of fractures on the oil-based mud imaging logging image are blurred, and there may be a phenomenon of missing recognition, or for wells lacking imaging data, it is judged whether the fracture can be picked by comparing the corrected wave impedance curve, if the corrected wave impedance value If it is less than 1300Kg/cm 3 ·s, it indicates that it is developed fractures, and the fracture development degree of the well interval can be qualitatively judged by dividing the fracture development interval of the well interval through the corrected wave impedance curve. The larger the ratio, the more developed the cracks, as shown in Figure 7.

本发明提供的识别方法,通过波阻抗法与成像测井图像对比,反映其较好的识别效果,同时也可以看出波阻抗法对于成像测井图中漏识的裂缝有很好的校正判别作用。The identification method provided by the present invention, through the comparison of the wave impedance method with the imaging logging image, reflects its better identification effect, and it can also be seen that the wave impedance method has a good correction and discrimination for the missed fractures in the imaging logging image effect.

Claims (9)

1. A crack logging identification method under the condition of oil-based mud is characterized by comprising the following steps,
s1, calculating a rock wave impedance value according to a density curve and a longitudinal wave acoustic wave time difference curve based on conventional logging data of a stratum to be analyzed;
s2, correcting the calculated rock wave impedance value according to a clay content curve in conventional logging data, calculating according to a preset correction formula to obtain a corrected wave impedance value, and drawing a corrected wave impedance curve graph;
and S3, judging whether the crack of the stratum to be analyzed is developed or not based on the calculated corrected wave impedance value, if the corrected wave impedance value is not greater than a preset wave impedance threshold value, indicating that the crack of the corresponding stratum is a developed crack, otherwise, indicating that the crack is an undeveloped crack.
2. The method for identifying a fracture log under oil-based mud conditions according to claim 1, wherein in S1, the method further comprises stripping the calculated abnormal rock wave impedance value;
the abnormal judgment standard is that the calculated rock wave impedance value deviates from a preset range of the rock wave impedance value corresponding to the normal rock physical property.
3. A method of fracture log identification in oil based mud conditions as set forth in claim 1 wherein, in S2, said correcting comprises,
selecting a plurality of undeveloped cracks at different argillaceous content positions in the rock core, and calculating a change value of a rock wave impedance value corresponding to the change of the argillaceous content;
comparing the relation between the change of the clay content and the change of the rock wave impedance value to obtain a correction formula,
when the burial depth is less than or equal to 5000m, correcting the wave impedance value Z' =Z+1.3 VSH, wherein VSH is the clay content, and Z is the rock wave impedance value;
when the burial depth is >5000m, the corrected wave impedance value Z' =z-1.3 VSH.
4. The method for identifying a fracture log under oil-based mud conditions according to claim 1, wherein the preset wave impedance threshold is set to 1300Kg/cm based on multi-well data 3 ·s。
5. The method for identifying a fracture log under oil-based mud conditions according to claim 1, further comprising qualitatively determining a fracture development level of the well section to be analyzed based on the corrected wave impedance value;
dividing a crack development interval of the well section according to the correction wave impedance curve chart;
the crack development degree is the ratio of the thickness of the crack development interval to the depth of the whole well section, and the greater the ratio is, the higher the crack development degree of the well section is.
6. The method for identifying a fracture log under oil-based mud conditions of claim 1, further comprising acquiring an oil-based mud imaging log image and comparing with a corrected wave impedance profile to create a fracture identification template;
for the cracks which can be directly identified on the oil-based mud imaging logging image, judging whether the cracks of the stratum to be analyzed are filled or not based on the comparison between the calculated corrected wave impedance value and a preset wave impedance threshold value, and if the corrected wave impedance value is larger than the wave impedance threshold value, indicating that the cracks of the corresponding stratum are filled;
and (3) continuing to carry out the step (S3) to judge whether the crack is developed or not for the crack which cannot be directly identified on the oil-based mud imaging logging image.
7. The method for identifying a fracture log under oil-based mud conditions according to claim 1, wherein the calculation formula of the rock wave impedance value in S1 is,
rock wave impedance z=10 5 ×DEN÷3.28DTC,
Wherein DEN is a density value, and DTC is a longitudinal wave acoustic wave time difference.
8. A fracture log identification system under oil-based mud conditions, characterized in that it comprises, based on the fracture log identification method of any one of claims 1-7,
the calculation module is used for calculating a rock wave impedance value according to a density curve and a longitudinal wave sound wave time difference curve based on conventional logging data of the stratum to be analyzed;
the correction module is used for calculating to obtain a correction wave impedance value according to a clay content curve in conventional logging data and a preset correction formula, and drawing a correction wave impedance curve graph;
and the comparison module is used for judging whether the crack of the stratum to be analyzed is developed or not based on the calculated corrected wave impedance value, and if the corrected wave impedance value is smaller than a preset wave impedance threshold value, the crack of the corresponding stratum is a development crack.
9. The system of claim 8, further comprising an imaging module for acquiring an oil-based mud imaging log image and comparing with the corrected wave impedance profile to create a fracture identification map.
CN202111357624.5A 2021-11-16 2021-11-16 Fracture logging identification method and system under oil-based mud conditions Pending CN116136606A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111357624.5A CN116136606A (en) 2021-11-16 2021-11-16 Fracture logging identification method and system under oil-based mud conditions

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111357624.5A CN116136606A (en) 2021-11-16 2021-11-16 Fracture logging identification method and system under oil-based mud conditions

Publications (1)

Publication Number Publication Date
CN116136606A true CN116136606A (en) 2023-05-19

Family

ID=86332942

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111357624.5A Pending CN116136606A (en) 2021-11-16 2021-11-16 Fracture logging identification method and system under oil-based mud conditions

Country Status (1)

Country Link
CN (1) CN116136606A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014169499A1 (en) * 2013-04-19 2014-10-23 中国石油大学(华东) Method for identifying and interpreting three-dimensional structure of ancient karst reservoir layer of carbonate rock
CN108897066A (en) * 2018-08-23 2018-11-27 中国石油天然气股份有限公司 Quantitative prediction method and device for fracture density of carbonate rock
US20190235106A1 (en) * 2016-10-04 2019-08-01 Landmark Graphics Corporation Multivariate Analysis Of Seismic Data, Microseismic Data, And Petrophysical Properties In Fracture Modeling
CN112731556A (en) * 2019-10-28 2021-04-30 中国石油化工股份有限公司 Crack development region prediction method and computer storage medium for predicting crack development region

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014169499A1 (en) * 2013-04-19 2014-10-23 中国石油大学(华东) Method for identifying and interpreting three-dimensional structure of ancient karst reservoir layer of carbonate rock
US20190235106A1 (en) * 2016-10-04 2019-08-01 Landmark Graphics Corporation Multivariate Analysis Of Seismic Data, Microseismic Data, And Petrophysical Properties In Fracture Modeling
CN108897066A (en) * 2018-08-23 2018-11-27 中国石油天然气股份有限公司 Quantitative prediction method and device for fracture density of carbonate rock
CN112731556A (en) * 2019-10-28 2021-04-30 中国石油化工股份有限公司 Crack development region prediction method and computer storage medium for predicting crack development region

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
林会喜: "河流相地层波阻抗特征及其控制因素", 石油与天然气地质, 31 October 2004 (2004-10-31), pages 524 - 527 *
程四洪 等: "变质岩潜山内幕裂缝表征及储层预测研究", 地球物理学进展, 31 December 2017 (2017-12-31), pages 598 - 600 *

Similar Documents

Publication Publication Date Title
CN105986815B (en) Method for identifying geological dessert in shale stratum
CN110320562B (en) Method for correcting acoustic time difference in horizontal well of shale gas reservoir
CN101950034B (en) A Baseline Drift Correction Method for Spontaneous Potential Curves Used in Seismic Inversion
WO2017127058A1 (en) Fracture interpretation with resistivity and sonic logs in biaxial anisotropic formations
CN104747163A (en) A method and device for identifying reservoir fractures in tight sandstone
WO2025092278A1 (en) Model well-based electric imaging logging fracture width calculation method
CN115327630B (en) Acoustic logging correction method, apparatus, device, storage medium and program product
CN110939428A (en) Identification method for reservoir fractures of tight sandstone oil and gas reservoir
CN116136606A (en) Fracture logging identification method and system under oil-based mud conditions
CN112698424B (en) Method for obtaining fault rock displacement pressure based on rock debris
CN114580233A (en) Method and system for identifying natural fractures of tight oil reservoir based on R/S analysis and finite difference method
CN116010789B (en) Carbonate reservoir type identification method, device, equipment and application
CN116359925B (en) Shallow profile position coordinate correction method based on water depth matching
CN109254330B (en) Method and device for identifying fracture segments of igneous rock reservoirs
CN115012903B (en) Logging evaluation method for judging shale bedding structure development
CN118346267A (en) Formation pressure prediction method, device, equipment and medium based on resistivity while drilling
CN105178950B (en) A method for determining total organic carbon in fine-grained sedimentary rocks
CN114065097B (en) A reservoir permeability determination method, system, device and readable storage medium
CN116931112A (en) Method for evaluating logging while drilling of tidal plateau phase dolomite reservoir
CN113971351A (en) Method and device for determining porosity of crack
CN118091786B (en) Method and device for correcting linear density of cracks based on underground logging data
CN113379209A (en) Shale oil saturation evaluation method based on geochemical pyrolysis correction
CN116263852A (en) Method and device for predicting total organic carbon content of shale and storage medium
CN112100254A (en) Method and system for docking logging data with logging data
CN114086939B (en) Oil and gas quantitative identification method and system for complex fault-block oil and gas reservoirs

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