CN114720494A - A method and device for predicting fracture opening coefficient of buried hill reservoir based on XRD whole-rock logging - Google Patents

A method and device for predicting fracture opening coefficient of buried hill reservoir based on XRD whole-rock logging Download PDF

Info

Publication number
CN114720494A
CN114720494A CN202210199743.0A CN202210199743A CN114720494A CN 114720494 A CN114720494 A CN 114720494A CN 202210199743 A CN202210199743 A CN 202210199743A CN 114720494 A CN114720494 A CN 114720494A
Authority
CN
China
Prior art keywords
elastic
mineral
rock
infinitesimal
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210199743.0A
Other languages
Chinese (zh)
Other versions
CN114720494B (en
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.)
CNOOC Energy Technology and Services Ltd
Original Assignee
CNOOC Energy Technology and Services 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 CNOOC Energy Technology and Services Ltd filed Critical CNOOC Energy Technology and Services Ltd
Priority to CN202210199743.0A priority Critical patent/CN114720494B/en
Publication of CN114720494A publication Critical patent/CN114720494A/en
Application granted granted Critical
Publication of CN114720494B publication Critical patent/CN114720494B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00Investigating 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/20Investigating 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V9/00Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/20Identification of molecular entities, parts thereof or of chemical compositions

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Chemical & Material Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Computing Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • General Business, Economics & Management (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Game Theory and Decision Science (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Primary Health Care (AREA)
  • Development Economics (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Mining & Mineral Resources (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开了一种基于XRD全岩录井潜山储层裂缝开度系数预测方法和装置,装置包括:岩屑采集、清洗、干燥和研磨装置,便携式XRD衍射仪,基质岩石矿物弹性微元模拟器,基质岩石弹性参数计算器,裂缝开度系数计算器,解释成图装置和长图打印装置。工作流程为,对上返岩屑取样、清洗、烘干和研磨,制备XRD衍射标准岩屑粉末;测试其X射线衍射谱图,确定矿物种类和含量;基质岩石矿物弹性微元模拟器进行矿物弹性微元串联和并联组合,计算弹性参数;测井和实验标定,计算基质岩石弹性参数;计算裂缝开度系数;计算结果随钻实时解释成图并长图打印。解决了潜山储层裂缝预测和随钻实时评价手段缺乏的问题,并有效处理了连续性差和准确性低的问题。

Figure 202210199743

The invention discloses a method and a device for predicting the fracture opening coefficient of a buried hill reservoir based on XRD whole-rock logging. Simulator, matrix rock elastic parameter calculator, fracture opening coefficient calculator, interpretation map device and long map printing device. The workflow is as follows: sampling, cleaning, drying and grinding the upturned cuttings to prepare XRD diffraction standard cuttings powder; testing the X-ray diffraction pattern to determine the type and content of minerals; The elastic micro-elements are combined in series and in parallel to calculate the elastic parameters; logging and experimental calibration are used to calculate the elastic parameters of the matrix rock; the fracture opening coefficient is calculated; the calculation results are interpreted in real time while drilling into a map and printed as a long map. It solves the problem of lack of means of fracture prediction and real-time evaluation while drilling in buried hill reservoirs, and effectively solves the problems of poor continuity and low accuracy.

Figure 202210199743

Description

一种基于XRD全岩录井潜山储层裂缝开度系数预测方法和 装置A method and device for predicting fracture opening coefficient of buried hill reservoir based on XRD whole-rock logging

技术领域technical field

本发明涉及石油天然气录井领域,尤其涉一种基于XRD全岩录井潜山储层裂缝开度系数预测方法和装置。The invention relates to the field of oil and natural gas logging, in particular to a method and a device for predicting the fracture opening coefficient of buried hill reservoirs based on XRD whole-rock logging.

背景技术Background technique

随着油气勘探的深入,潜山油气藏都有一定程度突破,对复杂潜山裂缝性储层的认识越来越迫切,对储层裂缝的预测和评价,是潜山储层勘探作业的一项核心工作。With the deepening of oil and gas exploration, breakthroughs have been made in buried hill oil and gas reservoirs to a certain extent, and the understanding of complex buried hill fractured reservoirs is becoming more and more urgent. core work.

岩屑录井所获取的岩屑实物信息是地层性质最直接的反映,但随着钻井新技术的使用,井底返出的岩屑通常十分细碎,甚至成粉末状,难于直接对储层的性质判断。XRD衍射全岩录井技术是近几年发展起来的一项主流录井技术,通过对岩屑样品进行X衍射全岩分析,获得矿物组分和含量,目前主要应用于地层岩性命名。实际矿物组分蕴含着丰富的地层信息,是地层的力学性质和裂缝发育难易程度的主控因素,目前这些信息并未得到充分挖掘,深化该信息在岩石力学特征和储层裂缝发育特征方面的应用尤为迫切和重要。The physical information of cuttings obtained by cuttings logging is the most direct reflection of formation properties. However, with the use of new drilling technologies, the cuttings returned from the bottom of the well are usually very fine and even powdery, which is difficult to directly affect the reservoir. judgement of nature. XRD diffraction whole-rock logging technology is a mainstream logging technology developed in recent years. The mineral composition and content are obtained by X-ray diffraction whole-rock analysis of cuttings samples. Currently, it is mainly used in formation lithology naming. The actual mineral composition contains rich stratigraphic information, which is the main controlling factor for the mechanical properties of the stratum and the difficulty of fracture development. At present, this information has not been fully excavated. To deepen this information in terms of rock mechanical characteristics and reservoir fracture development characteristics application is particularly urgent and important.

发明内容SUMMARY OF THE INVENTION

本发明的目的在于提供一种基于XRD全岩录井的潜山储层裂缝开度系数预测方法和装置,该方法和装置能够通过上返岩屑信息,进行XRD衍射全岩测试,随钻实时计算地层裂缝开度系数,预报裂缝发育难易程度。The purpose of the present invention is to provide a method and device for predicting the fracture opening coefficient of buried hill reservoirs based on XRD whole-rock logging, the method and device can perform XRD diffraction whole-rock test through the information of upturned cuttings, and real-time while drilling Calculate the formation fracture opening coefficient to predict the difficulty of fracture development.

第一方面,本发明提供了一种基于XRD全岩录井的潜山储层裂缝开度系数预测方法,按照下述步骤进行:In a first aspect, the present invention provides a method for predicting the fracture opening coefficient of buried hill reservoirs based on XRD whole-rock logging, which is carried out according to the following steps:

步骤1:按照非储层段稀疏储层段加密的原则设计取样间距,然后严格按照迟到时间和取样间距连续在振动筛前挑取真正来自井底的新鲜岩屑。Step 1: Design the sampling interval according to the principle of sparse and intensified reservoir sections in non-reservoir sections, and then continuously pick up fresh cuttings from the bottom of the well in strict accordance with the late arrival time and sampling interval in front of the vibrating screen.

步骤2:将岩屑加工成XRD衍射分析标准样品,测试其X射线衍射谱图,确定不同深度地层岩石基质的矿物含量。Step 2: The cuttings are processed into standard samples for XRD diffraction analysis, and their X-ray diffraction patterns are tested to determine the mineral content of the rock matrix at different depths.

步骤3:构建基质岩石矿物弹性微元模型,将各种矿物理想化为一种弹性微元,弹性微元和各矿物自身的力学属性一致,表征和确定各弹性微元的本构方程和弹性参数。Step 3: Build a matrix rock mineral elastic micro-element model, idealize various minerals into an elastic micro-element, the elastic micro-element is consistent with the mechanical properties of each mineral itself, and characterize and determine the constitutive equation and elasticity of each elastic micro-element parameter.

步骤4:构建矿物弹性微元组合物理模型,基于流变模型理论法,将各种矿物弹性微元进行串联和并联,构建矿物弹性微元组合物理模型。Step 4: Construct a physical model of mineral elastic micro-element combination. Based on the rheological model theory, various mineral elastic micro-elements are connected in series and parallel to build a mineral-elastic micro-element combination physical model.

步骤5:基于矿物弹性微元串联模型总应力和各矿物弹性微元应力相等,矿物弹性微元串联模型总应变等于各矿物弹性微元应变之和的原则,建立矿物弹性微元串联模型的本构方程和弹性参数。Step 5: Based on the principle that the total stress of the mineral-elastic micro-element series model is equal to the stress of each mineral elastic micro-element, and the total strain of the mineral-elastic micro-element series model is equal to the sum of the strains of each mineral elastic micro-element, the basic principle of the mineral-elastic micro-element series model is established. Constitutive equations and elastic parameters.

步骤6:基于矿物弹性微元并联模型总应力等于各矿物弹性微元应力之和,矿物弹性微元并联模型总应变和各矿物弹性微元应变相等的原则,建立矿物弹性微元并联模型的本构方程和弹性参数模型。Step 6: Based on the principle that the total stress of the mineral-elastic micro-element parallel model is equal to the sum of the stress of each mineral-elastic micro-element, and the total strain of the mineral-elastic micro-element parallel model is equal to the strain of each mineral elastic micro-element, the basic principle of the mineral-elastic micro-element parallel model is established. Constitutive equations and elastic parametric models.

步骤7:对矿物弹性微元串联模型和矿物弹性微元并联模型所计算的弹性模量和泊松比,进行测井和室内实验标定,建立适合于区域特征的基质岩石弹性模量和泊松比数学物理方程。Step 7: Perform well logging and laboratory experiment calibration on the elastic modulus and Poisson's ratio calculated by the mineral elastic micro-element series model and the mineral elastic micro-element parallel model, and establish the matrix rock elastic modulus and Poisson's ratio mathematics suitable for regional characteristics physical equations.

步骤8:基于二维平面应变下的均布载荷直裂纹开度增量经典方程,定义裂缝开度系数,并建立裂缝开度系数的弹性模量和泊松比表达形式。Step 8: Based on the classical equation of straight crack opening degree increment under uniform load under two-dimensional plane strain, define the crack opening degree coefficient, and establish the elastic modulus and Poisson's ratio expression form of the crack opening degree coefficient.

步骤9:基于步骤1和步骤2测定的不同深度地层岩石基质的矿物含量结果,运用步骤3和步骤4方法构建基质岩石矿物弹性微元模型和矿物弹性微元组合物理模型,结合步骤5和步骤6原理,推演矿物弹性微元串联模型和矿物弹性微元并联模型的本构方程和弹性参数,再通过步骤7的方法进行测井和室内实验标定,建立特征区域的基质岩石的弹性模量和泊松比数学物理方程,代入步骤8求取基质岩石的裂缝开度系数,预报储层的裂缝发育特征。Step 9: Based on the mineral content results of the formation rock matrix at different depths measured in Steps 1 and 2, use the methods of Steps 3 and 4 to construct the matrix rock mineral-elasticity micro-element model and mineral-elasticity micro-element combined physical model, combining steps 5 and 4. 6. Principle, deduce the constitutive equation and elastic parameters of the mineral-elastic micro-element series model and the mineral-elastic micro-element parallel model, and then use the method of step 7 to perform well logging and laboratory experiment calibration to establish the elastic modulus and poise of the matrix rock in the characteristic area. Substitute into step 8 to obtain the fracture opening coefficient of the matrix rock and predict the fracture development characteristics of the reservoir.

进一步的技术方案是步骤4矿物弹性微元组合物理模型包括矿物弹性微元串联模型和矿物弹性微元并联模型,其中矿物弹性微元串联模型载荷和变形的关系如下:A further technical solution is that the physical model of the mineral-elastic micro-element combination in step 4 includes a mineral-elastic micro-element series model and a mineral-elastic micro-element parallel model, wherein the relationship between the load and the deformation of the mineral-elastic micro-element series model is as follows:

Fz=F1z=F2z=…Fnz F z =F 1z =F 2z =...F nz

ΔLz=ΔL1+ΔL2+…ΔLn ΔL z =ΔL 1 +ΔL 2 +…ΔL n

式中:Fz为矿物弹性微元串联模型所受总载荷,KN;ΔLz为矿物弹性微元串联模型所受总的载荷总变形量;Fiz(i=1,2,…,n)为某个矿物弹性微元所受的载荷;ΔLi(i=1,2,…,n)为某个矿物弹性微元变形量。In the formula: F z is the total load of the mineral-elastic micro-element series model, KN; ΔL z is the total load and deformation of the mineral-elastic micro-element series model; F iz (i=1,2,…,n) is the load on a certain mineral elastic element; ΔL i (i=1,2,...,n) is the deformation amount of a certain mineral elastic element.

矿物弹性微元并联模型的载荷和变形关系如下:The load and deformation relationship of the mineral elastic micro-element parallel model is as follows:

Fr=F1r+F2r+…Fnr F r =F 1r +F 2r +...F nr

ΔL=ΔL1=ΔL2=…ΔLn ΔL=ΔL 1 =ΔL 2 =…ΔL n

式中:Fr为矿物弹性微元并联模型所受总的载荷,KN;ΔLr为矿物弹性微元并联模型所受总的载荷总变形量;Fir(i=1,2,…,n)为某个矿物弹性微元所受的载荷。In the formula: F r is the total load of the mineral elastic micro-element parallel model, KN; ΔL r is the total load and deformation of the mineral-elastic micro-element parallel model; F ir (i=1,2,…,n ) is the load on a mineral elastic element.

进一步的技术方案是步骤5建立的矿物弹性微元串联模型的本构方程和弹性参数表征如下:A further technical solution is that the constitutive equation and elastic parameters of the mineral elastic micro-element series model established in step 5 are characterized as follows:

矿物弹性微元串联模型的本构方程为The constitutive equation of the mineral elastic micro-element series model is:

Figure BDA0003527065200000031
Figure BDA0003527065200000031

矿物弹性微元串联模型弹性模量为The elastic modulus of the mineral elastic micro-element series model is

Figure BDA0003527065200000032
Figure BDA0003527065200000032

矿物弹性微元串联模型泊松比为The Poisson's ratio of the mineral elastic micro-element series model is

Figure BDA0003527065200000033
Figure BDA0003527065200000033

式中:σz为矿物弹性微元串联模型所受总应力;Ez为矿物弹性微元串联模型弹性模量;μz为矿物弹性微元串联模型泊松比;Ei(i=1,2,…,n)为第i种矿物弹性微元弹性模量;μi(i=1,2,…,n)为第i种矿物弹性微元泊松比;

Figure BDA0003527065200000034
为第i种矿物含量。In the formula: σ z is the total stress of the mineral elastic micro-element series model; E z is the elastic modulus of the mineral-elastic micro-element series model; μ z is the Poisson’s ratio of the mineral elastic micro-element series model; E i (i=1, 2,…,n) is the elastic modulus of the i-th mineral elastic element; μ i (i=1,2,…,n) is the Poisson’s ratio of the i-th mineral elastic element;
Figure BDA0003527065200000034
is the i-th mineral content.

进一步的技术方案是步骤6建立的矿物弹性微元并联模型的本构方程和弹性参数表征如下:A further technical solution is that the constitutive equation and elastic parameters of the mineral elastic micro-element parallel model established in step 6 are characterized as follows:

矿物弹性微元并联模型的本构方程The constitutive equation of the parallel model of mineral elastic micro-elements

Figure BDA0003527065200000035
Figure BDA0003527065200000035

矿物弹性微元并联模型弹性模量为The elastic modulus of the mineral elastic micro-element parallel model is

Figure BDA0003527065200000036
Figure BDA0003527065200000036

矿物弹性微元并联模型泊松比为The Poisson's ratio of the mineral elastic micro-element parallel model is

Figure BDA0003527065200000037
Figure BDA0003527065200000037

式中:σr为矿物弹性微元并联模型所受应力;Er为矿物弹性微元并联模型弹性模量;μr为矿物弹性微元串联模型泊松比。In the formula: σ r is the stress of the mineral-elastic micro-element parallel model; E r is the elastic modulus of the mineral-elastic micro-element parallel model; μ r is the Poisson’s ratio of the mineral-elastic micro-element series model.

进一步的技术方案是步骤7对串联模型和并联模型所计算的弹性模量和泊松比测井和室内实验标定,构建适合区域特征的基质弹性参数数学物理方程如下A further technical solution is that in step 7, the elastic modulus and Poisson's ratio calculated by the series model and the parallel model are calibrated by logging and laboratory experiments, and the mathematical and physical equations of matrix elastic parameters suitable for regional characteristics are constructed as follows:

Figure BDA0003527065200000041
Figure BDA0003527065200000041

Figure BDA0003527065200000042
Figure BDA0003527065200000042

式中:EA为岩石基质等效弹性模量;α为岩石基质等效弹性模量修正系数;μA为岩石基质等效泊松比,β为岩石基质等效泊松比修正系数。where E A is the equivalent elastic modulus of the rock matrix; α is the correction coefficient of the equivalent elastic modulus of the rock matrix; μ A is the equivalent Poisson’s ratio of the rock matrix, and β is the correction coefficient of the equivalent Poisson’s ratio of the rock matrix.

进一步的技术方案是步骤8裂缝开度系数的弹性模量和泊松比表达形式如下A further technical solution is that the elastic modulus and Poisson's ratio of the crack opening coefficient in step 8 are expressed as follows

Figure BDA0003527065200000043
Figure BDA0003527065200000043

式中:Fw为裂缝开度系数。In the formula: F w is the crack opening coefficient.

第二方面,本发明提供了一种基于XRD全岩录井的潜山储层裂缝开度系数预测装置,包括:岩屑采集装置,岩屑清洗装置,岩屑干燥装置,岩屑研磨装置,便携式XRD衍射仪,基质岩石矿物弹性微元模拟器,基质岩石弹性参数计算装置,裂缝开度系数计算装置,解释成图装置和长图打印装置。In a second aspect, the present invention provides a device for predicting the fracture opening coefficient of buried hill reservoirs based on XRD whole-rock logging, including: a cuttings collection device, a cuttings cleaning device, a cuttings drying device, and a cuttings grinding device, Portable XRD diffractometer, matrix rock mineral elastic micro-element simulator, matrix rock elastic parameter calculation device, fracture opening coefficient calculation device, interpretation map device and long map printing device.

该装置的工作流程:岩屑采集装置放在振动筛出砂口下方,岩屑沿筛布斜面落入取样器内;选取新鲜岩屑放入岩屑清洗装置,在对岩屑充分搅动的情况下对岩屑进行清洗,至岩屑漏出本色;清洗好的岩屑放入岩屑干燥装置,在温度控制小于85℃的环境下烘干;将烘干好的岩屑放入岩屑研磨装置进行研磨,使岩屑成粉末,粒度小于20um;将研磨好的岩屑粉末放入便携式XRD衍射仪,测试其X射线衍射谱图,解释不同深度地层岩石基质的矿物种类和含量数据;将岩石基质的矿物种类和含量数据传输给基质岩石矿物弹性微元模拟器,基质岩石矿物弹性微元模拟器依据矿物种类和含量构建矿物弹性微元,并进行矿物弹性微元串联和并联组合,计算串联和并联组合下的弹性参数Ezz,Err;将计算好的弹性参数Ezz,Err传输给基质岩石弹性参数计算装置,计算基质岩石弹性参数EA和μA;将基质岩石弹性参数EA和μA传输给裂缝开度系数计算装置,计算裂缝开度系数Fw;将基质岩石弹性参数计算装置和裂缝开度系数计算装置计算的基质岩石弹性参数EA和μA及裂缝开度系数Fw传输和存储于解释成图装置,进行成图,实时预报储层裂缝发育情况;最后通过长图打印装置成图进行打印和存档。The working process of the device: the cuttings collection device is placed under the sand outlet of the vibrating screen, and the cuttings fall into the sampler along the inclined surface of the screen cloth; fresh cuttings are selected and put into the cuttings cleaning device, and the cuttings are fully stirred Clean the cuttings until the cuttings leak out; put the cleaned cuttings into the cuttings drying device, and dry them in an environment where the temperature is controlled less than 85 °C; put the dried cuttings into the cuttings grinding device Grind to make the cuttings into powder with a particle size of less than 20um; put the ground cuttings powder into a portable XRD diffractometer to test the X-ray diffraction pattern, and explain the mineral types and content data of the rock matrix in different depths; The mineral type and content data of the matrix are transmitted to the matrix rock mineral elasticity micro-element simulator. The matrix rock mineral elasticity micro-element simulator constructs the mineral elasticity micro-element according to the mineral type and content, and carries out the series and parallel combination of the mineral elasticity micro-element to calculate the series connection. and the elastic parameters E z , μ z , E r , μ r in parallel combination; transmit the calculated elastic parameters E z , μ z , E r , μ r to the matrix rock elastic parameter calculation device to calculate the matrix rock elastic parameters E A and μ A ; transmit the matrix rock elastic parameters E A and μ A to the fracture opening coefficient calculation device, and calculate the fracture opening coefficient F w ; The rock elastic parameters E A and μ A and the fracture opening coefficient F w are transmitted and stored in the interpretation mapping device, and the mapping is performed to predict the development of reservoir fractures in real time; finally, the long map printing device is used to map and print and archive.

与现有技术相比,本发明的优点在于:Compared with the prior art, the advantages of the present invention are:

本发明通过随钻录井采集的岩屑,进行XRD衍射全岩矿物分析,获得不同深度地层岩石基质的矿物含量,再结合矿物组合物理模型,推演其岩石基质的本构方程和力学参数,预报裂缝开度系数,实现了对储层裂缝的预测和评价,具有成本低,数据采集实时性高和剖面预测连续性等优点。In the present invention, XRD diffraction whole-rock mineral analysis is carried out through the cuttings collected while drilling and logging, to obtain the mineral content of the rock matrix of different depths, and then combined with the physical model of mineral combination, the constitutive equation and mechanical parameters of the rock matrix are deduced, and the prediction is made. The fracture opening coefficient realizes the prediction and evaluation of reservoir fractures, and has the advantages of low cost, high real-time data acquisition and continuity of profile prediction.

附图说明Description of drawings

图1:基于XRD全岩录井的地层裂缝开度系数预测方法流程框图。Figure 1: Flow chart of the prediction method of formation fracture opening coefficient based on XRD whole-rock logging.

图2:矿物弹性微元模型。Figure 2: Mineral elastic microelement model.

图3:矿物弹性微元串联模型。Figure 3: Mineralelastic microelement tandem model.

图4:矿物弹性微元并联模型。Figure 4: Mineralelastic microelement parallel model.

图5:裂缝开度系数应用实例。Figure 5: Example of the application of the crack opening factor.

图6:基于XRD全岩录井的地层裂缝开度系数预测装置示意图。Figure 6: Schematic diagram of the prediction device for formation fracture opening coefficient based on XRD whole-rock logging.

对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,可以根据以上附图获得其他的相关附图。For those of ordinary skill in the art, other related drawings can be obtained from the above drawings without any creative effort.

具体实施方式Detailed ways

以下基于实例对本发明进行描述,但是值得说明的是,本发明并不限于这些实施例。在下文对本发明的细节描述中,详尽描述了一些特定的细节部分。然而,对于没有详尽描述的部分,本领域技术人员也完全可以理解本发明。The present invention is described below based on examples, but it should be noted that the present invention is not limited to these examples. In the following detailed description of the invention, some specific details are described in detail. However, for the parts that are not described in detail, those skilled in the art can fully understand the present invention.

此外,本领域普通技术人员应当理解,所提供的附图只是为说明本发明的目的、特征和优点,附图并不是实际按照比例绘制的。In addition, those of ordinary skill in the art will appreciate that the accompanying drawings are provided only to illustrate the objects, features and advantages of the present invention and are not actually drawn to scale.

实施例Example

图1是本发明基于XRD全岩录井的地层裂缝开度系数预测方法流程框图,对照图1,说明本发明的具体实施方式。Fig. 1 is a flow chart of the method for predicting formation fracture opening degree based on XRD whole-rock logging according to the present invention. With reference to Fig. 1, a specific embodiment of the present invention is described.

步骤1(101):按照非储层段稀疏,储层段加密的原则设计取样间距,然后严格按照迟到时间和取样间距连续在振动筛前挑取岩屑,剔除各种假岩屑以及残留岩屑、坍塌物和掉块,挑选出个体碎小、色调新鲜、棱角明显的真正来自井底的新鲜岩屑。Step 1 (101): Design the sampling interval according to the principle of sparse non-reservoir sections and intensified reservoir sections, and then continuously pick cuttings in front of the vibrating screen in strict accordance with the late arrival time and sampling distance, and remove various false cuttings and residual rocks The cuttings, collapses and falling blocks are selected to select the fresh cuttings from the bottom of the well with small individual fragments, fresh colors and obvious edges and corners.

步骤2(102):把新鲜岩屑研磨成粉末,装入射线衍射仪的样品池中,进行XRD衍射实验,测试样品的X射线衍射图谱,根据每种矿物的特征峰值,确定不同深度地层岩石基质的矿物种类和含量。Step 2 (102): Grind the fresh cuttings into powder, put it into the sample cell of the ray diffractometer, carry out the XRD diffraction experiment, test the X-ray diffraction pattern of the sample, and determine the formation rocks of different depths according to the characteristic peaks of each mineral The type and content of minerals in the matrix.

步骤3(103):构建基质岩石矿物弹性微元模型,将构成岩石各种矿物理想化为一种弹性微元,如附图2所示,弹性微元和各矿物自身的力学属性一致,各矿物自身的本构方程和弹性参数亦为弹性微元的弹性模量的本构方程和弹性参数,常见造岩矿物弹性微元的弹性模量和泊松比,如表1所示。Step 3 (103): constructing a matrix rock mineral elastic micro-element model, idealizing various minerals constituting the rock into an elastic micro-element, as shown in Figure 2, the elastic micro-element is consistent with the mechanical properties of each mineral itself, and each The constitutive equation and elastic parameters of the minerals themselves are also the constitutive equations and elastic parameters of the elastic modulus of the elastic elements. The elastic modulus and Poisson's ratio of the elastic elements of common rock-forming minerals are shown in Table 1.

表1常见造岩矿物弹性微元的弹性模量和泊松比Table 1 Elastic modulus and Poisson's ratio of elastic microelements of common rock-forming minerals

矿物mineral 弹性模量(GPa)Elastic Modulus (GPa) 泊松比Poisson's ratio 石英quartz 96.496.4 0.090.09 方解石Calcite 81.081.0 0.280.28 斜长石plagioclase 74.974.9 0.280.28 钠长石albite 78.578.5 0.290.29 白云石dolomite 121121 0.280.28 磁铁矿magnetite 230.8230.8 0.260.26 黄铁矿Pyrite 299.9299.9 0.160.16 绿帘石Epidote 154.2154.2 0.260.26 尖晶石spinel 293.3293.3 0.260.26 白云母muscovite 78.978.9 0.250.25 黑云母Biotite 69.6669.66 0.250.25 角闪石amphibole 128.8128.8 0.280.28 辉石pyroxene 143.7143.7 0.240.24 粘土矿物Clay minerals 14.214.2 0.300.30

步骤4(104):基于流变模型理论法的矿物弹性微元组合物理模型构建Step 4 (104): Construction of a physical model of mineral elastic micro-element combination based on rheological model theory

基于流变模型理论法,将各种矿物弹性微元进行串联和并联,构建矿物弹性微元组合物理模型,包括矿物弹性微元串联模型和矿物弹性微元并联模型。矿物弹性微元串联模型见图3,矿物弹性微元并联模型见图4。Based on the rheological model theory, various mineral elastic micro-elements are connected in series and parallel to construct a physical model of mineral elastic micro-elements, including the mineral-elastic micro-element series model and the mineral-elastic micro-element parallel model. The mineral-elastic micro-element series model is shown in Figure 3, and the mineral-elastic micro-element parallel model is shown in Figure 4.

其中,矿物弹性微元串联模型的载荷和变形关系遵循如下准则:Among them, the load and deformation relationship of the mineral elastic micro-element series model follows the following criteria:

Fz=F1z=F2z=…Fnz F z =F 1z =F 2z =...F nz

ΔLz=ΔL1+ΔL2+…ΔLn ΔL z =ΔL 1 +ΔL 2 +…ΔL n

其中,矿物弹性微元并联模型的载荷和变形关系遵循如下准则:Among them, the load and deformation relationship of the mineral elastic micro-element parallel model follows the following criteria:

Fr=F1r+F2r+…Fnr F r =F 1r +F 2r +...F nr

ΔL=ΔL1=ΔL2=…ΔLn ΔL=ΔL 1 =ΔL 2 =…ΔL n

步骤5(105):基于矿物弹性微元串联模型总应力和各矿物弹性微元应力相等,矿物弹性微元串联模型总应变等于各矿物弹性微元应变之和的原则,建立矿物弹性微元串联模型的本构方程和弹性参数。Step 5 (105): Based on the principle that the total stress of the mineral elastic micro-element series model is equal to the stress of each mineral elastic micro-element, and the total strain of the mineral-elastic micro-element series model is equal to the sum of the strains of each mineral elastic micro-element, establish the mineral elastic micro-element series The constitutive equations and elastic parameters of the model.

矿物弹性微元串联模型的本构方程为The constitutive equation of the mineral elastic micro-element series model is:

Figure BDA0003527065200000071
Figure BDA0003527065200000071

矿物弹性微元串联模型弹性模量为The elastic modulus of the mineral elastic micro-element series model is

Figure BDA0003527065200000072
Figure BDA0003527065200000072

矿物弹性微元串联模型泊松比为The Poisson's ratio of the mineral elastic micro-element series model is

Figure BDA0003527065200000073
Figure BDA0003527065200000073

步骤6(106):基于矿物弹性微元并联模型总应力等于各矿物弹性微元应力之和,矿物弹性微元并联模型总应变和各矿物弹性微元应变相等的原则,建立矿物弹性微元并联模型的本构方程和弹性参数模型。Step 6 (106): Based on the principle that the total stress of the mineral elastic micro-element parallel model is equal to the sum of the stress of each mineral elastic micro-element, and the total strain of the mineral-elastic micro-element parallel model is equal to the strain of each mineral elastic micro-element, establish the mineral elastic micro-element parallel connection Model constitutive equations and elastic parametric models.

矿物弹性微元并联模型的本构方程The constitutive equation of the parallel model of mineral elastic micro-elements

Figure BDA0003527065200000074
Figure BDA0003527065200000074

矿物弹性微元并联模型弹性模量为The elastic modulus of the mineral elastic micro-element parallel model is

Figure BDA0003527065200000075
Figure BDA0003527065200000075

矿物弹性微元并联模型泊松比为The Poisson's ratio of the mineral elastic micro-element parallel model is

Figure BDA0003527065200000076
Figure BDA0003527065200000076

步骤7(107):对矿物弹性微元串联模型和矿物弹性微元并联模型所计算的弹性模量和泊松比,进行测井和室内实验标定,建立适合与对应区域的基质岩石弹性模量和泊松比数学物理方程,方程形式如下:Step 7 (107): The elastic modulus and Poisson's ratio calculated by the mineral elastic micro-element series model and the mineral elastic micro-element parallel model are calibrated by logging and laboratory experiments, and the matrix rock elastic modulus and Poisson's ratio suitable for the corresponding area are established. Songbi mathematical physics equation, the equation form is as follows:

Figure BDA0003527065200000081
Figure BDA0003527065200000081

Figure BDA0003527065200000082
Figure BDA0003527065200000082

步骤8(108):基于XRD全岩录井的地层裂缝开度系数预测Step 8 (108): Prediction of formation fracture opening coefficient based on XRD whole-rock logging

基于线弹性断裂力学二维平面应变下的均布载荷直裂纹开度增量经典模型,将裂缝开度系数定义为:Based on the classical incremental model of straight crack opening under uniformly distributed load under two-dimensional plane strain of linear elastic fracture mechanics, the crack opening coefficient is defined as:

Figure BDA0003527065200000083
Figure BDA0003527065200000083

给定驱动应力变化的情况下,裂缝开度系数越高,裂缝越容易张开。For a given change in driving stress, the higher the crack opening coefficient, the easier it is for the crack to open.

步骤9(109):基于步骤1和步骤2测定的不同深度地层岩石基质的矿物含量结果,运用步骤3和步骤4方法构建基质岩石矿物弹性微元模型和矿物弹性微元组合物理模型,结合步骤5和步骤6原理,推演响应矿物弹性微元组合物理模型的本构方程和弹性参数,在通过步骤7的方法进行测井和室内实验标定和反算,建立适合区域特征的岩石基质弹性模量和泊松比数学物理方程,代入步骤8求取基质岩石的裂缝开度系数。Step 9 (109): Based on the mineral content results of the formation rock matrix at different depths measured in Steps 1 and 2, use the methods of Steps 3 and 4 to build a matrix rock mineral-elasticity micro-element model and a mineral-elasticity micro-element combined physical model, and combine the steps 5 and the principle of step 6, deduce the constitutive equation and elastic parameters of the physical model in response to the mineral elastic micro-element combination, and perform well logging and laboratory experiment calibration and back calculation by the method of step 7, and establish the elastic modulus of rock matrix suitable for regional characteristics. and Poisson's ratio mathematical and physical equation, and substitute into step 8 to obtain the fracture opening coefficient of the matrix rock.

步骤10(110)将基质岩石弹性参数EA和μA及裂缝开度系数Fw进行解释成图并长图打印,实时预报储层裂缝发育情况。In step 10 (110), the elastic parameters E A and μ A of the matrix rock and the fracture opening coefficient F w are interpreted into a map and printed as a long map, so as to predict the development of fractures in the reservoir in real time.

第二方面,本发明提供了一种基于XRD全岩录井的潜山储层裂缝开度系数预测装置,包括:岩屑采集装置601,岩屑清洗装置602,岩屑干燥装置603,岩屑研磨装置604,便携式XRD衍射仪605,基质岩石矿物弹性微元模拟器606,基质岩石弹性参数计算装置607,裂缝开度系数计算装置608,解释成图装置609和长图打印装置610。In the second aspect, the present invention provides a device for predicting the fracture opening coefficient of buried hill reservoirs based on XRD whole-rock logging, including: a cuttings collection device 601, a cuttings cleaning device 602, a cuttings drying device 603, a cuttings Grinding device 604, portable XRD diffractometer 605, matrix rock mineral elastic micro-element simulator 606, matrix rock elastic parameter calculation device 607, fracture opening coefficient calculation device 608, interpretation map device 609 and long map printing device 610.

该装置的工作流程:岩屑采集装置601放在振动筛出砂口下方,岩屑沿筛布斜面落入取样器内;选取新鲜岩屑放入岩屑清洗装置602,在对岩屑充分搅动的情况下对岩屑进行清洗,至岩屑漏出本色;清洗好的岩屑放入岩屑干燥装置603,在温度控制小于85℃的环境下烘干;将干燥的岩屑放入岩屑研磨装置604进行研磨,使岩屑成粉末,粒度小于20um;将研磨好的岩屑粉末放入便携式XRD衍射仪605,测试其X射线衍射谱图,解释不同深度地层岩石基质的矿物种类和含量数据;将岩石基质的矿物种类和含量数据传输给基质岩石矿物弹性微元模拟器606,基质岩石矿物弹性微元模拟器依据矿物种类和含量构建矿物弹性微元,并进行矿物弹性微元串联和并联组合,计算串联和并联组合下的弹性参数Ezz,Err;将计算好的弹性参数Ezz,Err传输给基质岩石弹性参数计算装置607,计算基质岩石弹性参数EA和μA;将基质岩石弹性参数EA和μA传输给裂缝开度系数计算装置608,计算裂缝开度系数Fw;将基质岩石弹性参数计算装置和裂缝开度系数计算装置计算的基质岩石弹性参数EA和μA及裂缝开度系数Fw传输和存储于解释成图装置609,进行成图,实时预报储层裂缝发育情况;最后通过长图打印装置610进行打印和存档The working process of the device: the cuttings collecting device 601 is placed under the sand outlet of the vibrating screen, and the cuttings fall into the sampler along the inclined surface of the screen cloth; Clean the cuttings until the cuttings leak out their true colors; put the cleaned cuttings into the cuttings drying device 603, and dry them in an environment where the temperature is controlled less than 85°C; put the dried cuttings into the cuttings for grinding The device 604 performs grinding to make the cuttings into powder with a particle size of less than 20um; put the ground cuttings powder into the portable XRD diffractometer 605, test its X-ray diffraction spectrum, and explain the mineral species and content data of the rock matrix in different depths ;Transfer the mineral type and content data of the rock matrix to the matrix rock mineral elasticity micro-element simulator 606, and the matrix rock mineral elasticity micro-element simulator builds the mineral elastic micro-element according to the mineral type and content, and performs the series and parallel connection of the mineral-elastic micro-element. Combine, calculate the elastic parameters E z , μ z , E r , μ r under the series and parallel combination; transmit the calculated elastic parameters E z , μ z , E r , μ r to the matrix rock elastic parameter calculation device 607, Calculate the matrix rock elastic parameters EA and μA ; transmit the matrix rock elastic parameters EA and μA to the fracture opening coefficient calculation device 608, and calculate the fracture opening coefficient Fw ; calculate the matrix rock elastic parameter calculation device and the fracture opening degree The matrix rock elastic parameters E A and μ A and the fracture opening coefficient F w calculated by the coefficient calculation device are transmitted and stored in the interpretation and mapping device 609 for mapping, and the development of reservoir fractures is predicted in real time; finally, through the long map printing device 610 print and archive

以上所述仅为本发明示意性的具体实施方式,并非用以限定本发明的范围。任何本领域的技术人员,在不脱离本发明的构思和原则的前提下所作出的等同变化与修改,均应属于本发明保护的范围。The above descriptions are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention. Equivalent changes and modifications made by any person skilled in the art without departing from the concept and principles of the present invention shall fall within the protection scope of the present invention.

Claims (8)

1. A prediction method for opening coefficient of buried hill reservoir fractures based on XRD full-rock logging is characterized by comprising the following steps:
step 1: designing a sampling interval according to a principle of encryption of a non-reservoir section sparse reservoir section, then continuously picking fresh rock debris really coming from a well bottom in front of a vibrating screen strictly according to late time and the sampling interval, cleaning, drying and grinding the rock debris, and processing the rock debris into XRD diffraction analysis standard rock debris powder;
step 2: placing the standard rock debris powder into a sample pool of a portable XRD diffractometer, testing an X-ray diffraction spectrogram of the standard rock debris powder, and determining the types and the content of minerals of rock matrixes of strata at different depths;
and step 3: constructing a matrix rock mineral elastic infinitesimal model, idealized various minerals into an elastic infinitesimal, the elastic infinitesimal is consistent with the self mechanical properties of each mineral, and the constitutive equation and the elastic parameters of each mineral elastic infinitesimal are represented and determined;
and 4, step 4: constructing a mineral combination physical model, and establishing a mineral elastic infinitesimal combination physical model by connecting various mineral elastic infinitesimal in series and in parallel based on a rheological model theory method, wherein the mineral elastic infinitesimal combination physical model comprises a mineral elastic infinitesimal series model and a mineral elastic infinitesimal parallel model;
and 5: establishing a constitutive equation and an elastic parameter model of the mineral elastic infinitesimal series model based on the principle that the total stress of the mineral elastic infinitesimal series model is equal to the stress of each mineral elastic infinitesimal, and the total strain of the mineral elastic infinitesimal series model is equal to the sum of the stresses of each mineral elastic infinitesimal;
step 6: establishing a constitutive equation and an elastic parameter model of the mineral elastic infinitesimal parallel model based on the principle that the total stress of the mineral elastic infinitesimal parallel model is equal to the sum of the stresses of all mineral elastic infinitesimals and the total strain of the mineral elastic infinitesimal parallel model is equal to the strain of all mineral elastic infinitesimals;
and 7: carrying out logging and indoor experimental calibration on the elastic modulus and the Poisson ratio calculated by the mineral elastic infinitesimal series model and the mineral elastic infinitesimal parallel model, and establishing a matrix rock elastic modulus and Poisson ratio mathematical physical equation suitable for a characteristic region;
and 8: defining a crack opening coefficient based on a classical equation of opening increment of the uniformly distributed load straight cracks under plane strain, and establishing an elastic modulus and Poisson ratio expression form of the crack opening coefficient;
and step 9: based on the mineral content results of the rock matrixes of different-depth strata measured in the steps 1 and 2, a mineral elastic infinitesimal model and a mineral elastic infinitesimal combined physical model of the matrix rock are constructed by using the methods in the steps 3 and 4, an constitutive equation and elastic parameters of a mineral elastic infinitesimal series model and a mineral elastic infinitesimal parallel model are deduced by combining the principles in the steps 5 and 6, logging and indoor experimental calibration are carried out by using the method in the step 7, the elastic modulus and Poisson's ratio mathematical physical equations of the matrix rock in a characteristic region are established, and the fracture opening coefficient of the matrix rock and the fracture development characteristic predicted by a reservoir are obtained by substituting the physical equations in the step 8.
2. The XRD-based whole rock logging buried hill reservoir fracture opening coefficient prediction method is characterized in that: the mineral elastic infinitesimal combination physical model comprises a mineral elastic infinitesimal series model and a mineral elastic infinitesimal parallel model;
the load and deformation relation of the mineral elastic infinitesimal series model is as follows:
Fz=F1z=F2z=…Fnz
ΔLz=ΔL1+ΔL2+…ΔLn
the load and deformation relation of the mineral elastic infinitesimal parallel model is as follows:
Fr=F1r+F2r+…Fnr
ΔL=ΔL1=ΔL2=…ΔLn
3. the prediction method of the opening coefficient of the subsurface reservoir fracture based on the XRD whole rock logging is characterized in that: the constitutive equation and the elastic parameters of the mineral elastic infinitesimal series model are characterized as follows:
constitutive equation of mineral elastic infinitesimal series model
Figure FDA0003527065190000021
Mineral elastic infinitesimal series model elastic modulus
Figure FDA0003527065190000022
Poisson's ratio of mineral elastic infinitesimal series model
Figure FDA0003527065190000023
4. The prediction method of the opening coefficient of the subsurface reservoir fracture based on the XRD whole rock logging is characterized in that: the constitutive equation and the elastic parameter of the mineral elastic infinitesimal parallel model are characterized in that:
constitutive equation of mineral elastic infinitesimal parallel model
Figure FDA0003527065190000024
Elastic modulus of mineral elastic infinitesimal parallel model
Figure FDA0003527065190000025
Poisson's ratio of mineral elastic infinitesimal parallel model
Figure FDA0003527065190000031
5. The prediction method of the opening coefficient of the subsurface reservoir fracture based on the XRD whole rock logging is characterized in that: the matrix rock elastic modulus and Poisson ratio mathematical physical equation is as follows:
Figure FDA0003527065190000032
Figure FDA0003527065190000033
6. the prediction method of the opening coefficient of the subsurface reservoir fracture based on the XRD whole rock logging is characterized in that: the expression forms of the elastic modulus and the Poisson ratio of the crack opening coefficient are as follows
Figure FDA0003527065190000034
In the formula: fwIs the crack opening coefficient.
7. The utility model provides a mine reservoir crack aperture coefficient prediction device dives based on XRD whole rock logging which characterized in that: the device comprises a rock debris collecting device, a rock debris cleaning device, a rock debris drying device, a rock debris grinding device, a portable XRD diffractometer, a matrix rock mineral elastic infinitesimal simulator, a matrix rock elastic parameter calculating device, a crack opening coefficient calculating device, an explanation mapping device and a long map printing device.
8. The XRD-based whole rock logging buried hill reservoir fracture opening coefficient prediction device is characterized in that: the working process is as follows:
the rock debris collecting device is arranged below the sand outlet of the vibrating screen, and rock debris falls into the sampler along the inclined surface of the screen cloth; selecting fresh rock debris, putting the fresh rock debris into a rock debris cleaning device, and cleaning the rock debris under the condition of fully stirring the rock debris until the rock debris leaks out of the natural color; putting the cleaned rock debris into a rock debris drying device, drying the rock debris in an environment with the temperature controlled to be less than 85 ℃, and putting the dried rock debris into a rock debris grinding device for grinding to enable the rock debris to be powder, wherein the granularity is less than 20 mu m; placing the ground rock debris powder into a portable XRD diffractometer, testing an X-ray diffraction spectrogram of the powder, and explaining the mineral types and content data of rock matrixes of strata at different depths; transmitting the mineral type and content data of the rock matrix to a matrix rock mineral elastic infinitesimal simulator, constructing mineral elastic infinitesimal by the matrix rock mineral elastic infinitesimal simulator according to the mineral type and content, combining the mineral elastic infinitesimal in series and in parallel, and calculating the elastic parameter E under the combination of series and in parallelzz,Err(ii) a The calculated elastic parameter Ezz,ErrTransmitting to a matrix rock elasticity parameter calculating device for calculating the matrix rock elasticity parameter EAAnd muA(ii) a The elastic parameter E of the matrix rockAAnd muATransmitting the crack opening coefficient to a crack opening coefficient calculation device to calculate a crack opening coefficient Fw(ii) a The elastic parameter E of the matrix rock calculated by the matrix rock elastic parameter calculating device and the crack opening coefficient calculating deviceAAnd muAAnd crack opening coefficient FwTransmitted and stored in the interpreting device to performThe method comprises the following steps of (1) forecasting the development condition of reservoir fractures in real time; and finally, printing and archiving the images by the long-image printing device.
CN202210199743.0A 2022-03-01 2022-03-01 A method and device for predicting fracture aperture coefficient of buried hill reservoir based on XRD whole rock logging Active CN114720494B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210199743.0A CN114720494B (en) 2022-03-01 2022-03-01 A method and device for predicting fracture aperture coefficient of buried hill reservoir based on XRD whole rock logging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210199743.0A CN114720494B (en) 2022-03-01 2022-03-01 A method and device for predicting fracture aperture coefficient of buried hill reservoir based on XRD whole rock logging

Publications (2)

Publication Number Publication Date
CN114720494A true CN114720494A (en) 2022-07-08
CN114720494B CN114720494B (en) 2024-11-08

Family

ID=82236351

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210199743.0A Active CN114720494B (en) 2022-03-01 2022-03-01 A method and device for predicting fracture aperture coefficient of buried hill reservoir based on XRD whole rock logging

Country Status (1)

Country Link
CN (1) CN114720494B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024087867A1 (en) * 2022-10-27 2024-05-02 中国石油天然气股份有限公司 Characterization method for in-situ openability of deep tectonic fissure, and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016041189A1 (en) * 2014-09-19 2016-03-24 杨顺伟 Method for evaluating shale gas reservoir and seeking desert area
CA2936443A1 (en) * 2015-07-20 2017-01-20 Cgg Services Sa Predicting mechanical and elastic rock properties of the subsurface
CN113138107A (en) * 2021-04-15 2021-07-20 东北石油大学 Rock brittleness evaluation method based on while-drilling rock debris logging information

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016041189A1 (en) * 2014-09-19 2016-03-24 杨顺伟 Method for evaluating shale gas reservoir and seeking desert area
CA2936443A1 (en) * 2015-07-20 2017-01-20 Cgg Services Sa Predicting mechanical and elastic rock properties of the subsurface
CN113138107A (en) * 2021-04-15 2021-07-20 东北石油大学 Rock brittleness evaluation method based on while-drilling rock debris logging information

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵丹云;刘修刚;秦可;杜建锋;张宏钧;申瑞臣;: "基于微观弹性模量与矿物组分页岩脆性评价方法研究", 西部探矿工程, no. 05, 15 May 2019 (2019-05-15) *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024087867A1 (en) * 2022-10-27 2024-05-02 中国石油天然气股份有限公司 Characterization method for in-situ openability of deep tectonic fissure, and storage medium

Also Published As

Publication number Publication date
CN114720494B (en) 2024-11-08

Similar Documents

Publication Publication Date Title
Ceryan et al. Prediction of unconfined compressive strength of carbonate rocks using artificial neural networks
Elkatatny et al. Development of a new correlation to determine the static Young’s modulus
CN105445440B (en) A kind of shale compressibility evaluation method based on landwaste microscopic feature
CN110286421A (en) A modeling method for natural fractures in tight sandstone reservoirs
CN110333200B (en) Method for delineating mineralization center based on short-wave infrared spectrum
Nefeslioglu Evaluation of geo-mechanical properties of very weak and weak rock materials by using non-destructive techniques: Ultrasonic pulse velocity measurements and reflectance spectroscopy
Azimian et al. An empirical correlation of uniaxial compressive strength with P-wave velocity and point load strength index on marly rocks using statistical method
Rezaei et al. Studying the correlation of rock properties with P-wave velocity index in dry and saturated conditions
CN104777035A (en) Shale compressibility comprehensive evaluation method based on uniaxial strength experiment
CN104677734A (en) Multi-scale information integration tight formation brittleness index measuring method
CN104020276B (en) The defining method of transverse isotropy shale reservoir rock mechanics parameter
CN105866835A (en) Fault 3D sealing quantitative evaluating method based on geostress distribution
Rastegarnia et al. Estimation of punch strength index and static properties of sedimentary rocks using neural networks in south west of Iran
CN103867198B (en) Carbonate rock natural gas layer stratum density discrimination method
Kahraman et al. The performance prediction of roadheaders from easy testing methods
CN110163931A (en) Large-scale geological mapping method for complex mining area
Hassan et al. The effect of mineralogical composition on strength and drillability of granitic rocks in Hulu Langat, Selangor Malaysia
CN109870720A (en) A kind of shale gas microcrack Logging Identification Method
CN114720494B (en) A method and device for predicting fracture aperture coefficient of buried hill reservoir based on XRD whole rock logging
Luodes et al. Evaluation of dimension stone in gneissic rocks—a case history from southern Finland
CN116559001A (en) A Quantitative Evaluation Method for Microscopic Compressibility of Tight Sandstone Reservoirs
Wang et al. Application of ultrasonic-rebound method in fast prediction of rock strength
CN102982244A (en) Fuzzy comprehensive evaluation method of geological strength index
CN113917554A (en) A method and electronic device for predicting the density of lamination fractures in shale reservoirs
CN113624603B (en) A method and system for obtaining elastic modulus of rock containing sedimentary rhythm layers

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant