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 PDFInfo
- 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
Links
- 239000011435 rock Substances 0.000 title claims abstract description 113
- 238000000034 method Methods 0.000 title claims abstract description 26
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 133
- 239000011707 mineral Substances 0.000 claims abstract description 133
- 239000011159 matrix material Substances 0.000 claims abstract description 60
- 238000004364 calculation method Methods 0.000 claims abstract description 13
- 239000000843 powder Substances 0.000 claims abstract description 12
- 238000000227 grinding Methods 0.000 claims abstract description 10
- 238000001035 drying Methods 0.000 claims abstract description 9
- 238000002441 X-ray diffraction Methods 0.000 claims abstract description 8
- 238000004140 cleaning Methods 0.000 claims abstract description 8
- 238000007639 printing Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims abstract description 7
- 238000012360 testing method Methods 0.000 claims abstract description 7
- 206010017076 Fracture Diseases 0.000 claims description 41
- 208000010392 Bone Fractures Diseases 0.000 claims description 35
- 238000011161 development Methods 0.000 claims description 9
- 238000013507 mapping Methods 0.000 claims description 5
- 239000004744 fabric Substances 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000002050 diffraction method Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims 1
- 238000003756 stirring Methods 0.000 claims 1
- 238000005520 cutting process Methods 0.000 abstract description 46
- 238000005553 drilling Methods 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 2
- 230000015572 biosynthetic process Effects 0.000 description 10
- 238000009533 lab test Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000003086 colorant Substances 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 229910021532 Calcite Inorganic materials 0.000 description 1
- 229910052656 albite Inorganic materials 0.000 description 1
- 229910052612 amphibole Inorganic materials 0.000 description 1
- 229910052626 biotite Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- IUMKBGOLDBCDFK-UHFFFAOYSA-N dialuminum;dicalcium;iron(2+);trisilicate;hydrate Chemical compound O.[Al+3].[Al+3].[Ca+2].[Ca+2].[Fe+2].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-] IUMKBGOLDBCDFK-UHFFFAOYSA-N 0.000 description 1
- YGANSGVIUGARFR-UHFFFAOYSA-N dipotassium dioxosilane oxo(oxoalumanyloxy)alumane oxygen(2-) Chemical compound [O--].[K+].[K+].O=[Si]=O.O=[Al]O[Al]=O YGANSGVIUGARFR-UHFFFAOYSA-N 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 229910052869 epidote Inorganic materials 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- SZVJSHCCFOBDDC-UHFFFAOYSA-N iron(II,III) oxide Inorganic materials O=[Fe]O[Fe]O[Fe]=O SZVJSHCCFOBDDC-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052627 muscovite Inorganic materials 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 229910052655 plagioclase feldspar Inorganic materials 0.000 description 1
- 239000011028 pyrite Substances 0.000 description 1
- 229910052683 pyrite Inorganic materials 0.000 description 1
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 1
- 229910052611 pyroxene Inorganic materials 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 229910052596 spinel Inorganic materials 0.000 description 1
- 239000011029 spinel Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- 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
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification 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射线衍射谱图，确定矿物种类和含量；基质岩石矿物弹性微元模拟器进行矿物弹性微元串联和并联组合，计算弹性参数；测井和实验标定，计算基质岩石弹性参数；计算裂缝开度系数；计算结果随钻实时解释成图并长图打印。解决了潜山储层裂缝预测和随钻实时评价手段缺乏的问题，并有效处理了连续性差和准确性低的问题。
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.
Description
技术领域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:
矿物弹性微元串联模型弹性模量为The elastic modulus of the mineral elastic micro-element series model is
矿物弹性微元串联模型泊松比为The Poisson's ratio of the mineral elastic micro-element series model is
式中：σz为矿物弹性微元串联模型所受总应力；Ez为矿物弹性微元串联模型弹性模量；μz为矿物弹性微元串联模型泊松比；Ei(i＝1,2,…,n)为第i种矿物弹性微元弹性模量；μi(i＝1,2,…,n)为第i种矿物弹性微元泊松比；
进一步的技术方案是步骤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
矿物弹性微元并联模型弹性模量为The elastic modulus of the mineral elastic micro-element parallel model is
矿物弹性微元并联模型泊松比为The Poisson's ratio of the mineral elastic micro-element parallel model is
式中：σ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:
式中：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
式中：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射线衍射谱图，解释不同深度地层岩石基质的矿物种类和含量数据；将岩石基质的矿物种类和含量数据传输给基质岩石矿物弹性微元模拟器，基质岩石矿物弹性微元模拟器依据矿物种类和含量构建矿物弹性微元，并进行矿物弹性微元串联和并联组合，计算串联和并联组合下的弹性参数Ez,μz,Er,μr；将计算好的弹性参数Ez,μz,Er,μr传输给基质岩石弹性参数计算装置，计算基质岩石弹性参数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
步骤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:
矿物弹性微元串联模型弹性模量为The elastic modulus of the mineral elastic micro-element series model is
矿物弹性微元串联模型泊松比为The Poisson's ratio of the mineral elastic micro-element series model is
步骤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
矿物弹性微元并联模型弹性模量为The elastic modulus of the mineral elastic micro-element parallel model is
矿物弹性微元并联模型泊松比为The Poisson's ratio of the mineral elastic micro-element parallel model is
步骤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:
步骤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:
给定驱动应力变化的情况下，裂缝开度系数越高，裂缝越容易张开。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
该装置的工作流程：岩屑采集装置601放在振动筛出砂口下方，岩屑沿筛布斜面落入取样器内；选取新鲜岩屑放入岩屑清洗装置602，在对岩屑充分搅动的情况下对岩屑进行清洗，至岩屑漏出本色；清洗好的岩屑放入岩屑干燥装置603，在温度控制小于85℃的环境下烘干；将干燥的岩屑放入岩屑研磨装置604进行研磨，使岩屑成粉末，粒度小于20um；将研磨好的岩屑粉末放入便携式XRD衍射仪605，测试其X射线衍射谱图，解释不同深度地层岩石基质的矿物种类和含量数据；将岩石基质的矿物种类和含量数据传输给基质岩石矿物弹性微元模拟器606，基质岩石矿物弹性微元模拟器依据矿物种类和含量构建矿物弹性微元，并进行矿物弹性微元串联和并联组合，计算串联和并联组合下的弹性参数Ez,μz,Er,μr；将计算好的弹性参数Ez,μz,Er,μr传输给基质岩石弹性参数计算装置607，计算基质岩石弹性参数EA和μA；将基质岩石弹性参数EA和μA传输给裂缝开度系数计算装置608，计算裂缝开度系数Fw；将基质岩石弹性参数计算装置和裂缝开度系数计算装置计算的基质岩石弹性参数EA和μA及裂缝开度系数Fw传输和存储于解释成图装置609，进行成图，实时预报储层裂缝发育情况；最后通过长图打印装置610进行打印和存档The working process of the device: the
以上所述仅为本发明示意性的具体实施方式，并非用以限定本发明的范围。任何本领域的技术人员，在不脱离本发明的构思和原则的前提下所作出的等同变化与修改，均应属于本发明保护的范围。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)
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)
| 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)
| 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 |
-
2022
- 2022-03-01 CN CN202210199743.0A patent/CN114720494B/en active Active
Patent Citations (3)
| 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)
| Title |
|---|
| 赵丹云;刘修刚;秦可;杜建锋;张宏钧;申瑞臣;: "基于微观弹性模量与矿物组分页岩脆性评价方法研究", 西部探矿工程, no. 05, 15 May 2019 (2019-05-15) * |
Cited By (1)
| 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 |