WO2018113149A1 - 一种获得动静态弹性参数转换关系的方法 - Google Patents

一种获得动静态弹性参数转换关系的方法 Download PDF

Info

Publication number
WO2018113149A1
WO2018113149A1 PCT/CN2017/081030 CN2017081030W WO2018113149A1 WO 2018113149 A1 WO2018113149 A1 WO 2018113149A1 CN 2017081030 W CN2017081030 W CN 2017081030W WO 2018113149 A1 WO2018113149 A1 WO 2018113149A1
Authority
WO
WIPO (PCT)
Prior art keywords
dynamic
ratio
pressure
static
elastic parameters
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.)
Ceased
Application number
PCT/CN2017/081030
Other languages
English (en)
French (fr)
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 AU2017363224A priority Critical patent/AU2017363224B2/en
Priority to US16/062,309 priority patent/US11175207B2/en
Publication of WO2018113149A1 publication Critical patent/WO2018113149A1/zh
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/40—Investigating hardness or rebound hardness
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00—Measuring force or stress, in general
    • G01L1/24—Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
    • G01N23/207—Diffractometry using detectors, e.g. using a probe in a central position and one or more displaceable detectors in circumferential positions
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/043—Analysing solids in the interior, e.g. by shear waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
    • G01N29/4418—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with a model, e.g. best-fit, regression analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24—Earth materials
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/60—Specific applications or type of materials
    • G01N2223/616—Specific applications or type of materials earth materials
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/028—Material parameters
    • G01N2291/02827—Elastic parameters, strength or force

Definitions

  • the invention relates to a method for obtaining a dynamic-static elastic parameter conversion relationship, and belongs to the technical field of petroleum exploration and logging.
  • the elastic parameter is a parameter used to describe the relationship between the stress applied to the rock and the strain generated, such as Young's modulus, Poisson's ratio, and so on.
  • the static elastic parameter is instructive for the fracturing transformation of the oil and gas layer, but since it can only be obtained by laboratory measurement, it cannot be continuously applied in the well. Therefore, firstly, the log data is used to obtain the dynamic elastic parameters of continuous depth.
  • the dynamic elastic parameters are converted into static elastic parameters by using the conversion rules of dynamic elastic parameters and static elastic parameters obtained based on experimental measurements. Finally, the full-section rock mechanics in the underground is used. Applied during the parameter evaluation process.
  • the existing conversion technology of dynamic and static elastic parameters generally involves placing a series of plunger samples obtained from drilling cores in the laboratory, measuring the dynamic and static elastic parameters under approximate formation conditions, and then establishing dynamic and static elasticity.
  • the linear conversion relationship of the parameters The prior art has obvious drawbacks because it does not consider the influence of the continuous variation of formation stress and lithology that may exist in the same geological horizon on the dynamic and static elastic parameter conversion law. Therefore, relying on the prior art is difficult to accurately and continuously evaluate the rock mechanics parameters of the formation with strong heterogeneous changes, restricting the accurate evaluation of the reservoir completion quality, affecting the optimization of the fracture test oil scheme and the realization of the cost reduction and efficiency improvement goals.
  • the present invention provides a method for obtaining a dynamic-static elastic parameter conversion relationship, which fully considers the influence of formation stress and clay content change on dynamic and static elastic parameter conversion laws, and improves rock mechanical parameters measurement.
  • Well evaluation accuracy is of great significance.
  • the present invention provides a method for obtaining a dynamic-static elastic parameter conversion relationship, the method comprising the following steps:
  • Step S1 acquiring horizontal cores at different depths of the destination layer
  • Step S2 measuring dynamic elastic parameters of the horizontal core under different pressures
  • Step S3 measuring static static parameters of the horizontal core under different pressures
  • Step S4 the clay content of the horizontal core
  • step S5 a relationship between the ratio of the dynamic and static elastic parameters and the formation pressure and the clay content is established; the conversion of the dynamic and static elastic parameters is completed.
  • the technical solution provided by the invention measures the dynamic and static elastic parameters of the core at different depths of the reservoir under different formation pressures and the clay content, and based on the measured results, the ratio of the dynamic and static elastic parameters and the formation pressure and The relationship between the clay content and the conversion law of the dynamic and static elastic parameters can be obtained.
  • the technical solution provided by the present invention is more accurate because it fully considers the effects of multiple influencing factors such as formation pressure and lithological changes.
  • step S1 the obtaining the horizontal core at different depths in the destination layer comprises the following steps:
  • Step S101 drilling m cores in the horizontal direction at the same depth D1 of the formation, which are respectively recorded as D 11 , D 12 , D 13 , ..., D 1m ;
  • Step S102 m cores drilled at the depth D1 of the formation, each piece is divided into two small pieces, and respectively recorded as D 11A , D 11B , D 12A , D 12B , D 13A , D 13B , ..., D 1 mA , D 1mB ;
  • Step S103 according to the operations of step S101 and step S102, drilling m cores at the depth D2 of the formation, and dividing each core of the m cores at the D2 into two small pieces, respectively Is D 21A , D 21B , D 22A , D 22B , D 23A , D 23B , ..., D 2mA , D 2mB ;
  • step S3 measuring the static elastic parameters of the horizontal core under different pressures comprises the following steps:
  • Step S501 calculating the ratio of the dynamic and static elastic parameters of the core samples at the depth D1 at the m pressure values in the P1-Pm, respectively, as E d11 /E s11 , E d12 /E s12 , E d13 /E s13 ,...,E d1m /E s1m , the plot of the ratio of the dynamic and static elastic parameters to the pressure when the formation depth is D1, that is, the ratio of the dynamic and static elastic parameters E d11 /E s11 , E d12 /E s12 , E d13 /E s13 ,...,E d1m /E s1m and their corresponding pressure values P1,P2,P3,...,Pm, the relationship between the two can be obtained by the graph;
  • Step S4 measuring the clay content of the horizontal core:

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • Remote Sensing (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Geology (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Acoustics & Sound (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

一种获得动静态弹性参数转换关系的方法。该方法包括步骤S1，获取目的地层中不同深度处的水平岩芯；步骤S2，在不同的压力下，测量所述水平岩芯的动态弹性参数；步骤S3，在不同的压力下，测量所述水平岩芯的静态弹性参数；步骤S4，测量所述水平岩芯的粘土含量；步骤S5，建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式；完成动静态弹性参数的转换。充分考虑了地层应力和粘土含量对动静态弹性参数转换规律的影响，对提高岩石力学参数测井评价精度具有重要意义。

Description

一种获得动静态弹性参数转换关系的方法 技术领域
本发明涉及一种获得动静态弹性参数转换关系的方法，属于石油勘探测井技术领域。
背景技术
低渗透及致密油气藏日益成为油气勘探的重点和难点领域。建立正确的动静态弹性参数转换方法，对于利用测井资料准确评价地层的岩石力学特征，进而优选压裂试油井段，优化试油完井方案具有重要意义。
弹性参数是用来描述岩石所受应力与所产生应变之间大小关系的参数，如杨氏模量、泊松比等。获得弹性参数的方式有两种：一种是通过实验室超声波及密度测量或者通过井下连续深度的声波及密度测井资料获得岩石的动态弹性参数；另一种是在实验室中通过应力应变测量获得岩石的静态弹性参数。对于油气层的压裂改造具有指导意义的是静态弹性参数，但由于其只能通过实验室测量才能得到，在井下无法连续应用。因此通常首先利用测井资料获得连续深度的动态弹性参数，利用基于实验测量获得的动态弹性参数与静态弹性参数的转换规律，将动态弹性参数转换为静态弹性参数，最后在井下全剖面地层岩石力学参数评价过程中应用。
现有动、静态弹性参数的转换技术，一般是将一系列从钻井取芯获得的柱塞样品置于实验室中，在近似地层条件下对其进行动静态弹性参数测量，然后建立动静态弹性参数的线性转换关系。现有技术存在明显的缺陷，因为它未考虑同一地质层位可能存在的地层应力及岩性的连续变化对动静态弹性参数转换规律的影响。因此依靠现有技术难以用来准确连续评价非均质变化强的地层的岩石力学参数，制约储层完井品质的准确评价，影响压裂试油方案的优化及降本增效目标的实现。
发明内容
为解决上述技术问题，本发明提供了一种获得动静态弹性参数转换关系的方法，该方法充分考虑了地层应力和粘土含量的变化对动静态弹性参数转换规律的影响，对提高岩石力学参数测井评价精度具有重要意义。
为达到上述目的，本发明提供了一种获得动静态弹性参数转换关系的方法，该方法包括以下步骤：
步骤S1，获取目的地层不同深度处的水平岩芯；
步骤S2，在不同的压力下，测量所述水平岩芯的动态弹性参数；
步骤S3，在不同的压力下，测量所述水平岩芯的静态弹性参数；
步骤S4，对所述水平岩芯的粘土含量；
步骤S5，建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式；完成动静态弹性参数的转换。
本发明提供的技术方案通过对储层不同深度水平岩心在不同地层压力下的动、静态弹性参数，以及粘土含量进行测量，根据测量得到的结果，建立了动静态弹性参数的比值与地层压力及粘土含量的函数关系式，由此便可得到动静态弹性参数的转换规律。本发明提供的技术方案由于充分考虑了地层压力及岩性变化等多重影响因素的作用，因而结果更加准确。
在上述方法中，优选地，在步骤S1中，所述获取目的地层中不同深度处的水平岩芯包括以下步骤：
步骤S101，在同一地层深度D1处沿水平方向钻取m块岩芯，分别记为D11，D12，D13，…，D1m；
步骤S102，将地层深度D1处钻取的m块岩芯，每一块切分为两小块，并分别记为D11A，D11B，D12A，D12B，D13A，D13B，…，D1mA，D1mB；
步骤S103，按照步骤S101和步骤S102的操作，钻取地层深度D2处的m块岩芯，并将所述D2处的m块岩芯中的每一块岩芯切分为两小块，分别记为D21A，D21B，D22A，D22B，D23A，D23B，…，D2mA，D2mB；
依此类推，钻取地层深度Dn处的m块岩芯，并将所述Dn处的m块岩芯每一块分为两小块，分别记为Dn1A，Dn1B，Dn2A，Dn2B，Dn3A，Dn3B，…，DnmA，DnmB。
在上述方法中，优选地，在步骤S101中，岩芯的长度至少为8cm，优选为8-10cm；直径为2.5-5cm，优选为2.5cm，但不限于此。
在上述方法中，优选地，在步骤S102中，将每一块岩芯切分为两小块时，两小块中的每一块岩芯的长度至少为4cm，优选为4-5cm；直径为2.5-5cm，优选为2.5cm，但不限于此。
在本发明提供的技术方案中，水平岩心的深度间隔、数量及大小可以根据实际勘测需求和实验精度要求进行合理设置。一般，为保障后续X衍射的测量精度和多地层压力条件的测量需求，获取的水平岩芯数量及重量需要达到一定的设计要求。
在本发明提供的技术方案中，步骤S1中，D11表示在地层深度D1处由第一步获得的岩芯，D12表示在地层深度D1处由第二步获得的岩芯，…，D1m表示在地层深度D1处由第m步获得的岩芯；依此类推，Dnm表示在地层深度Dn处由第m步获得的岩芯。
在上述方法中，优选地，在步骤S2中，在不同的压力下，测量所述水平岩芯的动态弹性参数包括以下步骤：
步骤S201，基于目的地层的压力，从小到大依次设置m个地层压力值，分别记为P1，P2，P3，…，Pm；其中，压力区间[P1，Pm]包含了整个目的地层的压力，即P1-Pm覆盖了整个目的地层的压力；
步骤S202，在压力为P1下，分别对不同地层深度处由第一步获得的水平岩芯D11A，D21A，D31A，…，Dn1A进行测量，计算它们的动态弹性参数，分别记为Ed11，Ed21，Ed31，…，Edn1；在压力为P2的条件下，分别对不同地层深度处由第二步获得的水平岩芯D12A，D22A，D32A，…，Dn2A进行测量，计算它们的动态弹性参数，分别记为Ed12，Ed22，Ed32，…，Edn2；依此类推，在压力为Pm下，分别对不同地层深度处由第m步获得的水平岩芯D1mA，D2mA，D3mA，…，DnmA进行测量，并计算它们的动态弹性参数，分别记为Ed1m，Ed2m，Ed3m，…，Ednm。
在上述方法中，优选地，在步骤S202中，所述测量的项目包括密度、纵波速度和横波速度，所述纵波速度和横波速度的测量方法参照SY/T 6351-2012；更优选地，所述测量的条件为饱和盐水条件。
在上述方法中，地层压力值的个数m与水平方向钻取的岩芯的块数m是相对应的；优选地，m、n的取值均≥3。
在上述方法中，优选地，在步骤S3中，在不同的压力下，测量所述水平岩芯的静态弹性参数包括以下步骤：
步骤S301，参照步骤S201的操作，基于目的地层的压力，从小到大依次设置m个地层压力值，分别记为P1，P2，P3，…，Pm；其中，压力区间[P1，Pm]包含了整个目的地层的压力，即P1-Pm覆盖了整个目的地层的压力；
步骤S302，在压力为P1下，分别对不同地层深度由第一步获得的水平岩芯D11B，D21B，D31B，…，Dn1B进行测量，计算它们的静态弹性参数，分别记为Es11，Es21，Es31，…，Esn1；在压力为P2下，分别对不同地层深度由第二步获得的水平岩芯D12B，D22B，D32B，…，Dn2B进行测量，计算它们的静态弹性参数，分别记为Es12，Es22，Es32，…，Esn2；依此类推，在压力为Pm下，分别对不同地层深度由第m步获得的水平岩芯D1mB，D2mB，D3mB，…，DnmB进行测量，并计算它们的静态弹性参数，分别记为Es1m，Es2m，Es3m，…，Esnm。
在上述方法中，优选地，在步骤S302中，所述测量的项目包括应力和应变，所述应力和应变的测量参照GBT23561.9-2009；更优选地，所述测量的条件为饱和盐水条件。
在上述方法中，优选地，在步骤S4中，测量所述水平岩芯的粘土含量包括以下步骤：
将测量过静态弹性参数的水平岩芯D11B，D12B，D13B，…，D1mB的碎样收集起来，进行X衍射测量，获得地层深度D1处的粘土含量Vcl1；将测量过静态弹性参数的水平岩芯D21B，D22B，D23B，…，D2mB的碎样收集起来，进行X衍射测量，获得地层深度D2处的粘土含量Vcl2；依此类推，将测量过静态弹性参数的水平岩芯Dn1B，Dn2B，Dn3B，…，DnmB的碎样收集起来，进行X衍射测量，获得地层深度Dn处的粘土含量Vcln。所述X衍射测量方法参照行业标准SY/T5163-2010。
在上述方法中，优选地，在步骤S5中，所述建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式包括以下步骤：
步骤S501，计算地层深度D1处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed12/Es12，Ed13/Es13，…，Ed1m/Es1m，构建地层深度为D1时的动静态弹性参数的比值与压力的曲线图，即构建动静态弹性参数的比值Ed11/Es11，Ed12/Es12，Ed13/Es13，…，Ed1m/Es1m与它们各自对应的压力值P1，P2，P3，…，Pm的曲线图，通过曲线图可以获得两者之间的关系；
计算地层深度D2处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Ed21/Es21，Ed22/Es22，Ed23/Es23，…，Ed2m/Es2m，构建地层深度为D2时的动静态弹性参数的比值与压力的曲线图，通过曲线图可以获得两者之间的关系；
依此类推，计算地层深度Dn处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Edn1/Esn1，Edn2/Esn2，Edn3/Esn3，…，Ednm/Esnm，构建地层深度为Dn时的动静态弹性参数的比值与压力的关曲线图，通过曲线图可以获得两者之间的关系；
步骤S502，压力值为P1时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed21/Es21，Ed31/Es31，…，Edn1/Esn1；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为P1时的动静态弹性参数的比值与粘土含量的曲线图，即构建动静态弹性参数的比值Ed11/Es11，Ed21/Es21，Ed31/Es31，…，Edn1/Esn1与它们各自对应的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln的曲线图，通过曲线图可以获得两者之间的关系；
压力值为P2时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed12/Es12，Ed22/Es22，Ed32/Es32，…，Edn2/Esn2；根据获得的D1-Dn中n个 地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为P2时的动静态弹性参数的比值与粘土含量的曲线图；
依此类推，压力值为Pm时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed1m/Es1m，Ed2m/Es2m，Ed3m/Es3m，…，Ednm/Esnm；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为Pm时的动静态弹性参数的比值与粘土含量的曲线图，通过曲线图可以获得两者之间的关系；
步骤S503，根据步骤S501和步骤S502获得的关系，建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式Ed/Es＝f(P,Vcl)。
本发明的有益效果：
1)本发明提供的技术方案通过对储层不同深度水平岩心在不同地层压力下的动、静态弹性参数，以及粘土含量进行测量，根据测量得到的结果，建立了动静态弹性参数的比值与地层压力及粘土含量的函数关系式，由此得到动静态弹性参数转换规律由于充分考虑了地层压力及岩性变化等多重影响因素的作用，因而结果更加准确；
2)本发明提供的技术方案充分考虑了地层应力条件和粘土含量对动静态弹性参数转换规律的影响，对提高岩石力学参数测井评价精度具有重要意义。
附图说明
图1为本发明实施例提供的动静态弹性参数转换方法的流程示意图；
图2为本发明实施例中测量岩芯纵、横波速度的装置的结构示意图；
图3为本发明实施例中动静态弹性参数的比值与地层压力的关系图；
图4为本发明实施例中动静态弹性参数的比值与粘土含量在地层压力为20Mpa时的关系图；
图5为本发明实施例中动静态弹性参数的比值与粘土含量在地层压力为40Mpa时的关系图；
图6为本发明实施例的应用效果对比图；
主要附图标号说明：
1：位移计；2：第一探头；3：第二探头；4：量管；5：孔隙流体控制单元；6：温度控制单元；7：围压控制单元；8：数据采集分析单元。
具体实施方式
为了对本发明的技术特征、目的和有益效果有更加清楚的理解，现对本发明的技术 方案进行以下详细说明，但不能理解为对本发明的可实施范围的限定。
实施例
本实施例提供了一种获得动静态弹性参数转换关系的方法(其流程如图1所示)，其包括以下步骤：
步骤S1，获取目的地层中不同深度处的水平岩芯：
在勘探区域不同深度钻取水平柱塞样岩芯，同一地层深度处钻取m＝4块水平小岩芯，构成一组；一共钻取n＝8组岩芯，分别记为D11，D12，D13，D14；D21，D22，D23，D24；…；D81，D82，D83，D84；其中，每一块水平岩芯的直径为2.5cm，长度为8-10cm。
分别对上述岩芯进行切割，将每一块岩芯切割为两小块，切割面与柱塞样岩芯底面平行，切割点的位置应保证切割后两块岩芯中每一块长度为4-5cm；将切割后满足此条件的岩芯分别记为D11A，D11B，D12A，D12B，…，D14A，D14B；D21A，D21B，D22A，D22B，…，D24A，D24B；…，D81A，D81B，D82A，D82B，…，D84A，D84B；其中，D11A，D11B分别表示：将D11切割为两小块之后，其中一块记为D11A，另一块记为D11B；依此类推。
步骤S2，在不同的压力下，测量所述水平岩芯的动态弹性参数：
基于地层压力，从小到大依次设置4个地层压力值，分别记为P1，P2，P3，P4(单位：MPa)，其中，P1＝20MPa，P2＝30MPa，P3＝40MPa，P4＝60MPa；20MPa-60MPa的压力数值范围包含了整个勘探区域的压力。
在压力为P1＝20MPa的条件下，对由步骤S1获得的水平岩芯D11A，D21A，D31A，…，D81A进行饱和盐水条件下的密度及纵、横波速度测量(纵、横波速度的测量方法参照SY/T 6351-2012)。
测量岩芯纵、横波速度的装置没有特别限定，可以是本领域的常规测量装置；例如可以是图2所示的装置(该装置为本领域常规测量装置中的一种)，该装置包括位移计1、第一探头2、第二探头3、量管4、孔隙流体控制单元5、温度控制单元6、围压控制单元7和数据采集分析单元8。测量时，通过孔隙流体控制单元5、温度控制单元6和围压控制单元7对测试条件进行调节；该装置中位移计1可以进行变形测量，样品两端的探头一方面可以允许流体通过；另一方面可以进行能量的发射和接收，将测得的数据传输给数据采集分析单元8。该测量装置只是本发明实施例的一个例举装置，其他具有相同功能或能够实现相同目的的测量装置均可以用来进行岩芯纵、横波速度的测量。
测量结果如表1所示，计算它们的动态杨氏模量Ed11，Ed21，Ed31，…，Ed81；其中，动态杨氏模量的计算公式如式1所示：
Figure PCTCN2017081030-appb-000001
式1中，E为动态杨氏模量，KPa；ρ为密度，g/cm3；Vp为纵波速度，m/s；Vs为横波速度，m/s；
表1
Figure PCTCN2017081030-appb-000002
Figure PCTCN2017081030-appb-000003
然后，在压力为P2＝30MPa的条件下，对步骤S1获得的水平岩芯D12A，D22A，D32A，…，D82A进行饱和盐水条件下的密度及纵横波速度测量，计算它们的动态杨氏模量Ed12，Ed22，Ed32，…，Ed82；
依此类推，在压力为P4＝60MPa的条件下，对由第一步获得的水平岩芯D14A，D24A，D34A，…，D84A进行饱和盐水条件下的密度及纵横波速度测量，计算它们的动态杨氏模量Ed14，Ed24，Ed34，…，Ed84。
步骤S3，在不同的压力下，测量所述水平岩芯的静态弹性参数：
在压力为P1＝20MPa的条件下，对由步骤S1获得的水平岩芯D11B，D21B，D31B，…，D81B进行饱和盐水条件下的应力应变测量(应力和应变的测量参照GBT23561.9-2009)，并根据测量结果，计算它们的静态杨氏模量(如表2所示)；其中，静态杨氏模量的计算公式如式2所示：
Figure PCTCN2017081030-appb-000004
在式2中：
σab—应力与轴向应变关系曲线上直线段终点与起点的应力差，单位为兆帕(MPa)；
εab—应力与轴向应变关系曲线上直线段终点与起点的应变差，百分数；
依此类推，在压力分别为P2＝30、40、60MPa的条件下，对由第一步获得的水平岩芯D12B、D22B、D32B、…、D82B，…，D14B、D24B、D34B、…、D84B进行饱和盐水条件下的应力应变测量，计算它们的静态杨氏模量。
表2
Figure PCTCN2017081030-appb-000005
Figure PCTCN2017081030-appb-000006
步骤S4，测量所述水平岩芯的粘土含量：
将经过应力应变测量后的岩芯D11B，D12B，…，D14B的碎样收集起来，进行X衍射测量，获得的粘土含量值代表地层深度D1的粘土含量；
依此类推，对经过应力应变测量后的其他各地层深度的水平岩芯碎样，分别收集起来进行X衍射测量，获得的粘土含量值代表各自对应的地层深度的粘土含量，如表3 所示。
表3
样号 粘土总量
D1 3.5
D2 2
D3 19.8
D4 20.8
D5 16.7
D6 12.5
D7 11.8
D8 14.2
步骤S5，建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式：
根据各地层深度在不同压力值下的测量计算结果，建立动静态弹性参数的比值(这里动静态弹性参数的比值即为动静态杨氏模量的比值)与地层压力的关系(如图3所示)；同时，结合各地层深度的粘土含量测量结果，建立固定压力值下的动静态弹性参数的比值与粘土含量的关系(如图4及图5所示，它们揭示了在地层压力分别为20MPa及40MPa时，动静弹性参数的比值与粘土含量的关系)；
根据动静态弹性参数的比值与地层压力的关系及其在不同地层压力值下与粘土含量的关系，建立动静态弹性参数的比值与粘土含量及地层压力的函数关系式，如式2所示：
Figure PCTCN2017081030-appb-000007
式2中，Ed为动态弹性参数(又称之为动态杨氏模量)，Es为静态弹性参数(又称之为静态杨氏模量)，Vclay为粘土含量，Peff为地层压力，A、B、C、D为根据实验结果得到的经验系数，在本实施例中，A＝3.5，B＝-8.14，C＝1.35，D＝-0.3。
图6是在勘探区域一口井中的应用效果对比图。图中共包含11道信息，其中第9道中的Esta为利用本发明实施例提供的方法计算得到的静态杨氏模量，Esta2为利用传统的动静态线性回归方法计算得到的静态杨氏模量，Esta_Core为岩芯测量获得的静态杨氏模量。对比结果表明，利用本发明实施例提供的方法计算的结果与岩芯测量结果最 为接近，三个数据点的平均相对误差为18％，而传统方法的平均相对误差达45％。对该区域3口井7个数据点的统计结果表明，与岩芯测量结果相比，本发明提供的方法的相对误差为13％，而传统方法的相对误差达28％。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (20)

  1. 一种获得动静态弹性参数转换关系的方法，该方法包括以下步骤：
    步骤S1，获取目的地层中不同深度的水平岩芯；
    步骤S2，在不同的压力下，测量所述水平岩芯的动态弹性参数；
    步骤S3，在不同的压力下，测量所述水平岩芯的静态弹性参数；
    步骤S4，测量所述水平岩芯的粘土含量；
    步骤S5，建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式；完成动静态弹性参数的转换。
  2. 根据权利要求1所述的方法，其中，在步骤S1中，所述获取目的地层中不同深度的水平岩芯包括以下步骤：
    步骤S101，在同一地层深度D1处沿水平方向钻取m块岩芯，分别记为D11，D12，D13，…，D1m；
    步骤S102，将地层深度D1处钻取的m块岩芯，每一块切分为两小块，并分别记为D11A，D11B，D12A，D12B，D13A，D13B，…，D1mA，D1mB；
    步骤S103，按照步骤S101和步骤S102的操作，钻取地层深度Dn处的m块岩芯，并将每一块岩芯切分为两小块，分别记为Dn1A，Dn1B，Dn2A，Dn2B，Dn3A，Dn3B，…，DnmA，DnmB。
  3. 根据权利要求2所述的方法，其中，在步骤S101中，岩芯的长度至少为8cm，直径为2.5-5cm。
  4. 根据权利要求2所述的方法，其中，在步骤S102中，将每一块岩芯切分为两小块时，两小块中的每一块岩芯的长度至少为4cm，直径为2.5-5cm。
  5. 根据权利要求2-4任一项所述的方法，其中，在步骤S2中，在不同的压力下，测量所述水平岩芯的动态弹性参数包括以下步骤：
    步骤S201，基于目的地层的压力，从小到大依次设置m个压力值，分别记为P1，P2，P3，…，Pm；其中，压力区间[P1，Pm]包含了整个目的地层的压力；
    步骤S202，在压力为P1下，对获得的水平岩芯D11A，D21A，D31A，…，Dn1A进行测量，计算它们的动态弹性参数，分别记为Ed11，Ed21，Ed31，…，Edn1；
    按照上述操作，在压力为Pm下，对获得的水平岩芯D1mA，D2mA，D3mA，…，DnmA进行测量，并计算它们的动态弹性参数，分别记为Ed1m，Ed2m，Ed3m，…，Ednm。
  6. 根据权利要求5所述的方法，其中，在步骤S202中，所述测量的项目包括密度、纵波速度和横波速度。
  7. 根据权利要求2-4和6任一项所述的方法，其中，在步骤S3中，在不同的压力下，测量所述水平岩芯的静态弹性参数包括以下步骤：
    步骤S301，基于目的地层的压力，从小到大依次设置m个压力值，分别记为P1，P2，P3，…，Pm；其中，压力区间[P1，Pm]包含了整个目的地层的压力；
    步骤S302，在压力为P1下，对获得的水平岩芯D11B，D21B，D31B，…，Dn1B进行测量，计算它们的静态弹性参数，分别记为Es11，Es21，Es31，…，Esn1；
    按照上述操作，在压力为Pm下，对获得的水平岩芯D1mB，D2mB，D3mB，…，DnmB进行测量，并计算它们的静态弹性参数，分别记为Es1m，Es2m，Es3m，…，Esnm。
  8. 根据权利要求5所述的方法，其中，在步骤S3中，在不同的压力下，测量所述水平岩芯的静态弹性参数包括以下步骤：
    步骤S301，基于目的地层的压力，从小到大依次设置m个压力值，分别记为P1，P2，P3，…，Pm；其中，压力区间[P1，Pm]包含了整个目的地层的压力；
    步骤S302，在压力为P1下，对获得的水平岩芯D11B，D21B，D31B，…，Dn1B进行测量，计算它们的静态弹性参数，分别记为Es11，Es21，Es31，…，Esn1；
    按照上述操作，在压力为Pm下，对获得的水平岩芯D1mB，D2mB，D3mB，…，DnmB进行测量，并计算它们的静态弹性参数，分别记为Es1m，Es2m，Es3m，…，Esnm。
  9. 根据权利要求2-4、6和8任一项所述的方法，其中，在步骤S302中，所述测量的项目包括应力和应变。
  10. 根据权利要求5所述的方法，其中，在步骤S302中，所述测量的项目包括应力和应变。
  11. 根据权利要求7所述的方法，其中，在步骤S302中，所述测量的项目包括应力和应变。
  12. 根据权利要求2-4、6、8、10和11任一项所述的方法，其中，在步骤S4中，测量所述水平岩芯的粘土含量包括以下步骤：
    步骤S401，对测量过静态弹性参数的水平岩芯D11B，D12B，D13B，…，D1mB进行X衍射测量，获得地层深度D1处的粘土含量Vcl1；
    步骤S402，按照步骤S401的操作，对测量过静态弹性参数的水平岩芯Dn1B，Dn2B，Dn3B，…，DnmB进行X衍射测量，获得地层深度Dn处的粘土含量Vcln。
  13. 根据权利要求5所述的方法，其中，在步骤S4中，测量所述水平岩芯的粘土含量包括以下步骤：
    步骤S401，对测量过静态弹性参数的水平岩芯D11B，D12B，D13B，…，D1mB进行X 衍射测量，获得地层深度D1处的粘土含量Vcl1；
    步骤S402，按照步骤S401的操作，对测量过静态弹性参数的水平岩芯Dn1B，Dn2B，Dn3B，…，DnmB进行X衍射测量，获得地层深度Dn处的粘土含量Vcln。
  14. 根据权利要求7所述的方法，其中，在步骤S4中，测量所述水平岩芯的粘土含量包括以下步骤：
    步骤S401，对测量过静态弹性参数的水平岩芯D11B，D12B，D13B，…，D1mB进行X衍射测量，获得地层深度D1处的粘土含量Vcl1；
    步骤S402，按照步骤S401的操作，对测量过静态弹性参数的水平岩芯Dn1B，Dn2B，Dn3B，…，DnmB进行X衍射测量，获得地层深度Dn处的粘土含量Vcln。
  15. 根据权利要求9所述的方法，其中，在步骤S4中，测量所述水平岩芯的粘土含量包括以下步骤：
    步骤S401，对测量过静态弹性参数的水平岩芯D11B，D12B，D13B，…，D1mB进行X衍射测量，获得地层深度D1处的粘土含量Vcl1；
    步骤S402，按照步骤S401的操作，对测量过静态弹性参数的水平岩芯Dn1B，Dn2B，Dn3B，…，DnmB进行X衍射测量，获得地层深度Dn处的粘土含量Vcln。
  16. 根据权利要求2-4、6、8、10-11和13-15中任一项所述的方法，其中，在步骤S5中，所述建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式包括以下步骤：
    步骤S501，计算地层深度D1处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed12/Es12，Ed13/Es13，…，Ed1m/Es1m；构建地层深度为D1时的动静态弹性参数的比值与压力的关系；
    按照上述操作，计算地层深度Dn处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Edn1/Esn1，Edn2/Esn2，Edn3/Esn3，…，Ednm/Esnm；构建地层深度为Dn时的动静态弹性参数的比值与压力的关系；
    步骤S502，压力值为P1时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed21/Es21，Ed31/Es31，…，Edn1/Esn1；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为P1时的动静态弹性参数的比值与粘土含量的关系；
    按照上述操作，压力值为Pm时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed1m/Es1m，Ed2m/Es2m，Ed3m/Es3m，…，Ednm/Esnm；根据获 得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为Pm时的动静态弹性参数的比值与粘土含量的关系；
    步骤S503，根据步骤S501和步骤S502获得的关系，建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式Ed/Es＝f(P,Vcl)。
  17. 根据权利要求5所述的方法，其中，在步骤S5中，所述建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式包括以下步骤：
    步骤S501，计算地层深度D1处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed12/Es12，Ed13/Es13，…，Ed1m/Es1m；构建地层深度为D1时的动静态弹性参数的比值与压力的关系；
    按照上述操作，计算地层深度Dn处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Edn1/Esn1，Edn2/Esn2，Edn3/Esn3，…，Ednm/Esnm；构建地层深度为Dn时的动静态弹性参数的比值与压力的关系；
    步骤S502，压力值为P1时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed21/Es21，Ed31/Es31，…，Edn1/Esn1；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为P1时的动静态弹性参数的比值与粘土含量的关系；
    按照上述操作，压力值为Pm时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed1m/Es1m，Ed2m/Es2m，Ed3m/Es3m，…，Ednm/Esnm；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为Pm时的动静态弹性参数的比值与粘土含量的关系；
    步骤S503，根据步骤S501和步骤S502获得的关系，建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式Ed/Es＝f(P,Vcl)。
  18. 根据权利要求7所述的方法，其中，在步骤S5中，所述建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式包括以下步骤：
    步骤S501，计算地层深度D1处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed12/Es12，Ed13/Es13，…，Ed1m/Es1m；构建地层深度为D1时的动静态弹性参数的比值与压力的关系；
    按照上述操作，计算地层深度Dn处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Edn1/Esn1，Edn2/Esn2，Edn3/Esn3，…，Ednm/Esnm；构建地层深度 为Dn时的动静态弹性参数的比值与压力的关系；
    步骤S502，压力值为P1时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed21/Es21，Ed31/Es31，…，Edn1/Esn1；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为P1时的动静态弹性参数的比值与粘土含量的关系；
    按照上述操作，压力值为Pm时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed1m/Es1m，Ed2m/Es2m，Ed3m/Es3m，…，Ednm/Esnm；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为Pm时的动静态弹性参数的比值与粘土含量的关系；
    步骤S503，根据步骤S501和步骤S502获得的关系，建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式Ed/Es＝f(P,Vcl)。
  19. 根据权利要求9所述的方法，其中，在步骤S5中，所述建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式包括以下步骤：
    步骤S501，计算地层深度D1处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed12/Es12，Ed13/Es13，…，Ed1m/Es1m；构建地层深度为D1时的动静态弹性参数的比值与压力的关系；
    按照上述操作，计算地层深度Dn处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Edn1/Esn1，Edn2/Esn2，Edn3/Esn3，…，Ednm/Esnm；构建地层深度为Dn时的动静态弹性参数的比值与压力的关系；
    步骤S502，压力值为P1时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed21/Es21，Ed31/Es31，…，Edn1/Esn1；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为P1时的动静态弹性参数的比值与粘土含量的关系；
    按照上述操作，压力值为Pm时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed1m/Es1m，Ed2m/Es2m，Ed3m/Es3m，…，Ednm/Esnm；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为Pm时的动静态弹性参数的比值与粘土含量的关系；
    步骤S503，根据步骤S501和步骤S502获得的关系，建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式Ed/Es＝f(P,Vcl)。
  20. 根据权利要求12所述的方法，其中，在步骤S5中，所述建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式包括以下步骤：
    步骤S501，计算地层深度D1处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed12/Es12，Ed13/Es13，…，Ed1m/Es1m；构建地层深度为D1时的动静态弹性参数的比值与压力的关系；
    按照上述操作，计算地层深度Dn处的岩芯样品在P1-Pm中m个压力值下的动静态弹性参数的比值，分别记为Edn1/Esn1，Edn2/Esn2，Edn3/Esn3，…，Ednm/Esnm；构建地层深度为Dn时的动静态弹性参数的比值与压力的关系；
    步骤S502，压力值为P1时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed11/Es11，Ed21/Es21，Ed31/Es31，…，Edn1/Esn1；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为P1时的动静态弹性参数的比值与粘土含量的关系；
    按照上述操作，压力值为Pm时，计算岩心样品在D1-Dn中n个地层深度下的动静态弹性参数的比值，分别记为Ed1m/Es1m，Ed2m/Es2m，Ed3m/Es3m，…，Ednm/Esnm；根据获得的D1-Dn中n个地层深度的粘土含量Vcl1，Vcl2，Vcl3，……，Vcln，构建压力值为Pm时的动静态弹性参数的比值与粘土含量的关系；
    步骤S503，根据步骤S501和步骤S502获得的关系，建立动静态弹性参数的比值与地层压力及粘土含量的函数关系式Ed/Es＝f(P,Vcl)。
PCT/CN2017/081030 2016-12-20 2017-04-19 一种获得动静态弹性参数转换关系的方法 Ceased WO2018113149A1 (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
AU2017363224A AU2017363224B2 (en) 2016-12-20 2017-04-19 Method for obtaining conversion relationship between dynamic and static elastic parameters
US16/062,309 US11175207B2 (en) 2016-12-20 2017-04-19 Method for obtaining conversion relationship between dynamic and static elastic parameters

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201611183911.8 2016-12-20
CN201611183911.8A CN106855485B (zh) 2016-12-20 2016-12-20 一种动静态弹性参数的转换方法

Publications (1)

Publication Number Publication Date
WO2018113149A1 true WO2018113149A1 (zh) 2018-06-28

Family

ID=59126411

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2017/081030 Ceased WO2018113149A1 (zh) 2016-12-20 2017-04-19 一种获得动静态弹性参数转换关系的方法

Country Status (4)

Country Link
US (1) US11175207B2 (zh)
CN (1) CN106855485B (zh)
AU (1) AU2017363224B2 (zh)
WO (1) WO2018113149A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111666665A (zh) * 2020-05-26 2020-09-15 中国石油大学(华东) 一种提高碳酸岩动静态杨氏模量转换模型精度的方法
CN117890477A (zh) * 2024-03-13 2024-04-16 西南交通大学 一种基于tsp数据对岩石抗压强度的推算方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107829727A (zh) * 2017-09-12 2018-03-23 中国石油天然气股份有限公司 一种确定地应力的方法及装置
CN109630098A (zh) * 2019-02-15 2019-04-16 中国地质科学院地质力学研究所 一种利用岩芯测量差应力的方法
CN111735696B (zh) * 2020-07-02 2022-09-02 科吉思石油技术咨询(北京)有限公司 一种评价油气藏岩心样品地下原位杨氏模量的方法
CN114112651A (zh) * 2020-08-27 2022-03-01 中国石油化工股份有限公司 一种用于人造岩心的岩石动静态力学参数转换方法及系统
CN119958979A (zh) * 2025-01-27 2025-05-09 中国石油大学(北京) 一种目标岩心的动静态力学参数转换关系确定方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103267678A (zh) * 2013-04-28 2013-08-28 北京大学 一种动静态岩石力学参数同步测量方法及装置
WO2014205248A2 (en) * 2013-06-19 2014-12-24 Conocophillips Company Mechanical characterization of core samples
CN104267435A (zh) * 2014-08-29 2015-01-07 中国石油天然气股份有限公司 一种横观各向同性地层弹性系数的测井计算方法及装置
CN104406849A (zh) * 2014-11-21 2015-03-11 中国石油天然气股份有限公司 一种储层岩石脆性的预测方法及装置
CN104853822A (zh) * 2014-09-19 2015-08-19 杨顺伟 一种评价页岩气储层及寻找甜点区的方法
CN105370267A (zh) * 2014-08-29 2016-03-02 中国石油天然气股份有限公司 一种分析致密砂岩弹性系数应力敏感性的方法及装置
CN105865955A (zh) * 2016-03-23 2016-08-17 成都创源油气技术开发有限公司 页岩脆性测井评价方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8498853B2 (en) * 2009-07-20 2013-07-30 Exxonmobil Upstream Research Company Petrophysical method for predicting plastic mechanical properties in rock formations
US20110246159A1 (en) * 2010-04-02 2011-10-06 Herwanger Jorg V Method and Apparatus to Build a Three-Dimensional Mechanical Earth Model
CN103257081B (zh) 2013-04-28 2015-09-02 北京大学 一种油气藏岩体力学地下原位模型恢复的方法及装置
CN104977618B (zh) 2014-04-09 2017-09-22 中国石油集团东方地球物理勘探有限责任公司 一种评价页岩气储层及寻找甜点区的方法
US10132162B2 (en) * 2014-09-10 2018-11-20 Fracture ID, Inc. Apparatus and method using measurements taken while drilling to map mechanical boundaries and mechanical rock properties along a borehole

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103267678A (zh) * 2013-04-28 2013-08-28 北京大学 一种动静态岩石力学参数同步测量方法及装置
WO2014205248A2 (en) * 2013-06-19 2014-12-24 Conocophillips Company Mechanical characterization of core samples
CN104267435A (zh) * 2014-08-29 2015-01-07 中国石油天然气股份有限公司 一种横观各向同性地层弹性系数的测井计算方法及装置
CN105370267A (zh) * 2014-08-29 2016-03-02 中国石油天然气股份有限公司 一种分析致密砂岩弹性系数应力敏感性的方法及装置
CN104853822A (zh) * 2014-09-19 2015-08-19 杨顺伟 一种评价页岩气储层及寻找甜点区的方法
CN104406849A (zh) * 2014-11-21 2015-03-11 中国石油天然气股份有限公司 一种储层岩石脆性的预测方法及装置
CN105865955A (zh) * 2016-03-23 2016-08-17 成都创源油气技术开发有限公司 页岩脆性测井评价方法

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111666665A (zh) * 2020-05-26 2020-09-15 中国石油大学(华东) 一种提高碳酸岩动静态杨氏模量转换模型精度的方法
CN111666665B (zh) * 2020-05-26 2022-05-27 中国石油大学(华东) 一种提高碳酸岩动静态杨氏模量转换模型精度的方法
CN117890477A (zh) * 2024-03-13 2024-04-16 西南交通大学 一种基于tsp数据对岩石抗压强度的推算方法
CN117890477B (zh) * 2024-03-13 2024-05-17 西南交通大学 一种基于tsp数据对岩石抗压强度的推算方法

Also Published As

Publication number Publication date
CN106855485B (zh) 2019-08-06
AU2017363224B2 (en) 2020-05-07
US11175207B2 (en) 2021-11-16
CN106855485A (zh) 2017-06-16
US20190265142A1 (en) 2019-08-29
AU2017363224A1 (en) 2018-07-05

Similar Documents

Publication Publication Date Title
WO2018113149A1 (zh) 一种获得动静态弹性参数转换关系的方法
CN106368691B (zh) 基于岩石物理地震信息三维异常孔隙压力预测方法
CN107622139B (zh) 裂缝渗透率的计算方法
CN110058323A (zh) 一种致密砂岩地层脆性指数计算方法
CN105931125B (zh) 一种致密油分段多簇体积压裂水平井产量预测方法
CN101413388A (zh) 获得油水同层原始含油饱和度的方法及估算未试油油水同层原始含油饱和度的方法
US10697286B2 (en) Method of identifying a material and/or condition of a material in a borehole
CN105952446B (zh) 一种石油天然气气层成分含量测定方法
CN106368687B (zh) 一种页岩储层脆性评价方法
CN104632201B (zh) 一种页岩油气“甜点区”有机碳含量下限值的确定方法
CN108133086A (zh) 一种应力敏感储层中产水压裂气井裂缝半长反演方法
NO20120953A1 (no) En ny algoritme for noytronporositetsmaling ved bruk av nevralt nettverk
CN108133087A (zh) 一种气-水两相渗流应力敏感储层原始渗透率反演方法
CN109458176A (zh) 碳酸盐岩储层压力的预测方法及其应用
CN115822562A (zh) 一种考虑层内窜流的纵向非均质气藏产能评价方法
CN113777668A (zh) 用于砂泥岩互层致密气藏的地应力计算方法及装置
CN105182421B (zh) 一种定量评价地层脆性特征的方法
CN103306672B (zh) 一种预测页岩地层不同钻进方向的研磨性的方法
WO2018231594A1 (en) Estimation of mechanical properties of transversely isotropic media
CN107247860B (zh) 一种求取页岩储层有机孔隙度的测录井方法
CN105929448B (zh) 一种原状泥页岩地层声波时差响应反演方法
CN102383784A (zh) 存储式漏层位置综合测量仪
CN106250664A (zh) 低孔裂缝性砂岩气藏产能的预测方法和装置
Sanchez et al. Fluid analysis and sampling: the next big step for logging while drilling tools
CN103790581A (zh) 一种确定煤岩工业组分物理参数的方法

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2017363224

Country of ref document: AU

Date of ref document: 20170419

Kind code of ref document: A

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17882536

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 17882536

Country of ref document: EP

Kind code of ref document: A1