CN105675635B - Tight rock component relative amount and brittleness index determine method and apparatus - Google Patents

Tight rock component relative amount and brittleness index determine method and apparatus Download PDF

Info

Publication number
CN105675635B
CN105675635B CN201511031926.8A CN201511031926A CN105675635B CN 105675635 B CN105675635 B CN 105675635B CN 201511031926 A CN201511031926 A CN 201511031926A CN 105675635 B CN105675635 B CN 105675635B
Authority
CN
China
Prior art keywords
relative amount
indicatrix
brittleness index
rock
mineral constituent
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.)
Active
Application number
CN201511031926.8A
Other languages
Chinese (zh)
Other versions
CN105675635A (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.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201511031926.8A priority Critical patent/CN105675635B/en
Publication of CN105675635A publication Critical patent/CN105675635A/en
Application granted granted Critical
Publication of CN105675635B publication Critical patent/CN105675635B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00Investigating materials by wave or particle radiation
    • G01N2223/05Investigating materials by wave or particle radiation by diffraction, scatter or reflection
    • G01N2223/056Investigating materials by wave or particle radiation by diffraction, scatter or reflection diffraction

Landscapes

  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The present invention provides a kind of tight rock component relative amounts and brittleness index to determine method and apparatus, wherein, this method includes:It is analyzed based on X diffraction total rock, establishes each mineral constituent relative amount indicatrix in tight rock;Reconstruct each mineral constituent relative amount indicatrix and brittleness index indicatrix in tight rock;Each mineral constituent relative amount indicatrix and brittleness index indicatrix carry out seismic waveform instruction inverting, obtain each mineral constituent relative amount and the inverting data volume of brittleness index in tight rock as constraints using in the tight rock of reconstruct.The present invention solves each mineral constituent relative amount and brittleness index relation unobvious in compact reservoir in the prior art and causes be difficult on seismic profile the technical issues of, has reached the simple technique effect for accurately determining relation between each mineral constituent relative amount and brittleness index.

Description

Tight rock component relative amount and brittleness index determine method and apparatus
Technical field
The present invention relates to technical field of geological exploration, more particularly to a kind of tight rock component relative amount and brittleness index Determine method and apparatus.
Background technology
With the increasing of conventional oil and gas resource exploration development difficulty, unconventional petroleum resources will play important role. In recent years, not only fine and close petroleum resources is found that in large-scale terrestrial lake basin, also have breakthrough in succession in small-sized fault depressed lacustrine basin.In view of This, under the background of compact reservoir, the problem of prediction of " dessert " reservoir becomes extensive concern.Especially compact reservoir rock ore deposit Object constituent content and brittleness index are both one of object that unconventionaloil pool exploration and development is paid close attention to and unconventionaloil pool Explore the important parameter of " dessert " area selection.
The content of brittle mineral determines the effect of later stage fracturing reform in compact reservoir rock forming mineral component, and directly Affect oil and gas production.Rock brittleness with high quartz content or high-carbon Carbonate Rocks content is relatively high, in fracturing process Failure by shear can be generated and advantageously form complicated netted seam, by volume fracturing technology, well yield can be improved.In depth Hollow band or sloped region have the rock of high potassium feldspar content, secondary pore growing, rich accumulation of oil and gas.
Therefore, how effectively to predict that brittle mineral constituent content and brittleness index spatial distribution are directly closed in compact reservoir It is the exploration deployment work to next step.
At present, less to the research of rock forming mineral component and brittleness index, rock forming mineral component is mainly captured using element Well logging (ECS) calculates mineral content, expensive, and the content of each mineral constituent can only can be just calculated in specific software.It is crisp Sex index can indirectly be predicted using 2 elastic parameters of Young's modulus and Poisson's ratio, but cannot be quantified and be asked in rock respectively Mineral constituent content can not disclose rock brittleness Producing reason.
How rock forming mineral component and brittleness index to be studied as a whole, not yet propose effective solution party at present Case.
The content of the invention
An embodiment of the present invention provides a kind of tight rock component relative amounts and brittleness index to determine method, this method bag It includes:
It is analyzed based on X diffraction total rock, establishes each mineral constituent relative amount indicatrix in tight rock;The fine and close rock of reconstruct Each mineral constituent relative amount indicatrix and brittleness index indicatrix in stone;With each mineral constituent in the tight rock of reconstruct Relative amount indicatrix and brittleness index indicatrix carry out seismic waveform instruction inverting, obtain densification as constraints Each mineral constituent relative amount and the inverting data volume of brittleness index in rock.
In one embodiment, each mineral constituent relative amount indicatrix in tight rock is reconstructed, including:Pass through survey Using multiple linear regression analysis method, it is bent to reconstruct each mineral constituent relative amount feature in tight rock for well curve sensitivity analysis Line.
In one embodiment, each mineral constituent relative amount indicatrix is expressed as in the tight rock of reconstruct:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep side To RS represents shallow lateral, and a, b, c, d, e, f, g are the parameter after reconstruct, is constant.
In one embodiment, the brittleness index indicatrix in tight rock is reconstructed, including:According to the densification of reconstruct Each mineral constituent relative amount indicatrix in rock reconstructs the brittleness index indicatrix in tight rock, wherein, the cause The brittleness index indicatrix of close rock is expressed as:
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, CARBONATE Represent the relative amount of carbonate rock, CLAY represents the relative amount of clay.
In one embodiment, referred to each mineral constituent relative amount indicatrix and brittleness in the tight rock of reconstruct Number indicatrix carries out seismic waveform instruction inverting, obtains each mineral constituent relative amount in tight rock as constraints And the inverting data volume of brittleness index, including:With each mineral constituent relative amount indicatrix in the tight rock of reconstruct and crisp Sex index indicatrix is as constraints, by analyzing prospect pit, determine the sample number of model foundation and high frequency into Point;According to the sample number and radio-frequency component determined, initial model is established;Based on the initial model, pass through seismic wave Shape indicates inverting, is finally inversed by the spatial distribution and brittleness index of the relative amount of each mineral constituent on well profile and even well profile Spatial distribution;The spatial distribution for crossing the relative amount of each mineral constituent on well profile and even well profile that inverting is obtained It is bent with each mineral constituent relative amount indicatrix and brittleness index feature of the reconstruct with the spatial distribution of brittleness index Line is compared;If difference is less than predetermined error threshold, it is anti-that the initial model is applied to entire seismic data cube Drill to obtain each mineral constituent relative amount and the inverting data volume of brittleness index in rock.
The embodiment of the present invention additionally provides a kind of tight rock component relative amount and brittleness index determining device, the device Including:
Module is established, for being based on the analysis of X diffraction total rock, it is bent to establish each mineral constituent relative amount feature in tight rock Line;Reconstructed module, for reconstructing each mineral constituent relative amount indicatrix and brittleness index indicatrix in tight rock;Instead Drill module, for each mineral constituent relative amount indicatrix in the tight rock using reconstruct and brittleness index indicatrix as Constraints carries out seismic waveform instruction inverting, obtains the anti-of each mineral constituent relative amount and brittleness index in tight rock Drill data volume.
In one embodiment, the reconstructed module is specifically used for passing through logging sensitivity analy-sis, and utilization is polynary Linear regression method reconstructs each mineral constituent relative amount indicatrix in tight rock.
In one embodiment, each mineral constituent relative amount feature is bent in the tight rock of the reconstructed module reconstruct Line is expressed as:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep side To RS represents shallow lateral, and a, b, c, d, e, f, g are the parameter after reconstruct, is constant.
In one embodiment, the reconstructed module is specifically used for each mineral constituent phase in the tight rock according to reconstruct To content characteristics curve, the brittleness index indicatrix in tight rock is reconstructed, wherein, the brittleness index of the tight rock is special Sign curve is expressed as:
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, CARBONATE Represent the relative amount of carbonate rock, CLAY represents the relative amount of clay.
In one embodiment, the inverting module includes:Determination unit, for each ore deposit in the tight rock with reconstruct Object component relative amount indicatrix and brittleness index indicatrix, by analyzing prospect pit, are determined as constraints Go out the sample number and radio-frequency component of model foundation;Unit is established, for according to the sample number and radio-frequency component determined, establishing just Beginning model;First inverting unit, for based on the initial model, indicating inverting by seismic waveform, being finally inversed by well The spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index on section and even well profile;Comparison unit, For the spatial distribution and brittleness for crossing the relative amount of each mineral constituent on well profile and even well profile for obtaining inverting The spatial distribution of index, each mineral constituent relative amount indicatrix and brittleness index indicatrix with the reconstruct carry out Comparison;Second inverting unit, in the case of being less than predetermined error threshold in definite difference, by the initial model application Each mineral constituent relative amount and the inverting data volume of brittleness index in rock are obtained in entire seismic data cube inverting.
In embodiments of the present invention, it is special that each mineral constituent relative amount indicatrix and brittleness index in tight rock are reconstructed Levy curve, and using in the tight rock of reconstruct each mineral constituent relative amount indicatrix and brittleness index indicatrix as about Beam condition carries out seismic waveform instruction inverting, obtains the inverting of each mineral constituent relative amount and brittleness index in tight rock Data volume solves each mineral constituent relative amount and brittleness index relation unobvious in compact reservoir in the prior art and causes Be difficult on seismic profile the technical issues of, reach simple and accurately determined each mineral constituent relative amount and brittleness index Between relation technique effect.
Description of the drawings
Attached drawing described herein is used for providing a further understanding of the present invention, forms the part of the application, not Form limitation of the invention.In the accompanying drawings:
Fig. 1 is the method flow that tight rock constituent content according to embodiments of the present invention and brittleness index determine method Figure;
Fig. 2 is each mineral group in utilization multiple linear regression analysis method reconstruct L1 well tight rocks according to embodiments of the present invention Divide relative amount and brittleness index indicatrix schematic diagram;
Fig. 3 is that multiple linear regression analysis method according to embodiments of the present invention predicts mineral constituent content and measured value correlation Analyze schematic diagram;
Fig. 4 is that quartz content (Q) inverting according to embodiments of the present invention connects well profile schematic diagram;
Fig. 5 is that feldspar content (F) inverting according to embodiments of the present invention connects well profile schematic diagram;
Fig. 6 is that brittleness index (BI) inverting according to embodiments of the present invention connects well profile schematic diagram;
Fig. 7 is the plane isopleth distribution schematic diagram of quartz content (Q) according to embodiments of the present invention;
Fig. 8 is the plane isopleth distribution schematic diagram of feldspar content (F) according to embodiments of the present invention;
Fig. 9 is the plane isopleth distribution schematic diagram of brittleness index (BI) according to embodiments of the present invention;
Figure 10 is tight rock constituent content according to embodiments of the present invention and the structure diagram of brittleness index determining device.
Specific embodiment
Understand to make the object, technical solutions and advantages of the present invention clearer, it is right with reference to embodiment and attached drawing The present invention is described in further details.Here, the exemplary embodiment and its explanation of the present invention be for explaining the present invention, but simultaneously It is not as a limitation of the invention.
In embodiments of the present invention, provide a kind of tight rock constituent content and brittleness index determines method, such as Fig. 1 institutes Show, comprise the following steps:
Step 101:It is analyzed based on X diffraction total rock, establishes each mineral constituent relative amount indicatrix in tight rock;
Step 102:Reconstruct each mineral constituent relative amount indicatrix and brittleness index indicatrix in tight rock;
Specifically, using multiple linear regression analysis method, can be reconstructed by logging sensitivity analy-sis in tight rock Each mineral constituent relative amount indicatrix, each mineral constituent relative amount indicatrix can represent in the tight rock of reconstruct For:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep side To RS represents shallow lateral, and a, b, c, d, e, f, g are the parameter after reconstruct, is constant.
Correspondingly, can each mineral constituent relative amount indicatrix in the tight rock according to reconstruct, reconstruct fine and close rock Brittleness index indicatrix in stone, wherein, the brittleness index indicatrix of tight rock can be expressed as:
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, CARBONATE Represent the relative amount of carbonate rock, CLAY represents the relative amount of clay.
When specific implementation, it can be analyzed by logging sensitivity analy-sis with reference to different methods, with Reconstruct each mineral constituent relative amount indicatrix in tight rock.
Step 103:It is bent with each mineral constituent relative amount indicatrix in the tight rock of reconstruct and brittleness index feature Line carries out seismic waveform instruction inverting, obtains each mineral constituent relative amount and brittleness in tight rock and refer to as constraints Several inverting data volumes.
I.e., it is possible to each mineral constituent relative amount indicatrix and brittleness index indicatrix in the tight rock of reconstruct As constraints, first to crossing well profile and even the instruction inverting of well profile progress seismic wave, each mineral group in tight rock is obtained Divide relative amount and the inverting data volume of brittleness index.
Specifically, seismic waveform can be carried out by following steps and indicates inverting:
S1:Using in the tight rock of reconstruct each mineral constituent relative amount indicatrix and brittleness index indicatrix as Constraints by analyzing prospect pit, determines the sample number and radio-frequency component of model foundation;
S2:According to the sample number and radio-frequency component determined, initial model is established;
S3:Based on the initial model, inverting is indicated by seismic waveform, was finally inversed by well profile and even well profile The upper spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index;
S4:The spatial distribution for crossing the relative amount of each mineral constituent on well profile and even well profile that inverting is obtained It is bent with each mineral constituent relative amount indicatrix and brittleness index feature of the reconstruct with the spatial distribution of brittleness index Line is compared;
S5:If difference is less than predetermined error threshold, it is anti-that the initial model is applied to entire seismic data cube Drill to obtain each mineral constituent relative amount and the inverting data volume of brittleness index in rock.
Each component content in above-mentioned tight rock and brittleness index are determined with reference to a specific embodiment method into Row explanation, however, it is noted that, which is not formed to the present invention not merely to the present invention is better described Work as restriction.
In this example primarily to solving how " dessert " storage is effectively predicted under compact reservoir background in the prior art The problem of layer rock forming mineral component and brittleness index spatial distribution, is proposed one kind and is combined based on experiment, well logging, seismic facies therefore Compact reservoir in rock forming mineral component and brittleness index spatial prediction techniques.With experimental analysis (analysis of X diffraction total rock) for base Plinth, it is effective using seismic waveform instruction inversion technique different from conventional wave impedance inversion on Logging Curves analysis foundation Rock forming mineral constituent content and brittleness index distribution and development degree in compact reservoir are predicted, has promoted unconventional reservoir The development of middle rock forming mineral constituent content and brittleness index Predicting Technique.
Specifically, each component content and brittleness index determine that method mainly includes the following aspects in the tight rock:
1) brittleness index is calculated using each mineral constituent relative amount of experimental analysis in rock, can included:
S1:Rock core to typical prospect pit or typical outcrop samples systematic sampling in sampling process, can be set between sampling 0.5-1m is divided into, preferably has sample spot distribution in the different depth section of interval of interest, it so could utility well-log information The indicatrix precision of forward modeling is higher;
S2:X diffraction total rock experimental analyses are carried out to the sample of acquisition;
S3:X diffraction total rock analysis results are classified, you can be divided into quartzy (Q), feldspar (F), carbonate rock (CARBONATE), clay (CLAY) and other (OTHERS) five major class mineral, if conditions permit, can also segment feldspar For potassium feldspar (KF) and plagioclase (XF) etc., the measured value for dividing these mineral and brittleness index (BI) is finely playbacked in lithology On the premise of, it demarcates on well, for carrying out logging sensitivity analy-sis.
Specifically, the relative amount of each mineral constituent in tight rock can be utilized to obtain the expression of brittleness index (BI) Formula:
Wherein, BI represents rock brittleness index, and V represents mineral volume, if V is well logging calculated value, then V is represented just It is a volumn concentration, that is, represents the relative amount of each ore deposit component.
2) each mineral constituent relative amount and brittleness index in well logging quantitative judge tight rock, can include:
S1:Each mineral constituent relative amount and brittleness index log parameter sensitivity analysis in tight rock:
Each mineral constituent of rock core or outcrop samples in tight rock at different depth can be determined by above-mentioned analysis Then relative amount and brittleness index size, can sum up different minerals component log response according to various conventional log Feature.For example, it can sum up:Quartzy (Q) relative amount electrical property feature is more apparent, and resistivity value is higher, natural gamma value, Middle subvalue and density value are relatively low.More than one high three low logging response character be quartzy relative amount main log response it is special Sign;Clay mineral (CLAY) relative amount is more sensitive to natural gamma, density, interval transit time and resistivity, shows as nature Gamma, density and interval transit time are higher, and resistivity is relatively low.The above are quartz and clay mineral relative amount logging response character, Other mineral constituent relative amounts and brittleness index logging response character are different, the log type sensitive to its Difference, so-called sensitive reference is exactly height of the log to the influence degree of each mineral constituent relative amount.
S2:Using multiple linear regression analysis method, each mineral constituent relative amount and brittleness index in quantitative forecast rock.
Since single factor test logging method multi-solution is larger, log parameter is multifactor concentrated expression.Therefore, can combine Log parameter sensitivity analysis by analyzing distinct methods, preferably goes out each in multiple linear regression analysis method reconstruct tight rock Mineral constituent relative amount and brittleness index indicatrix.For example, reconstruct quartzy (Q), feldspar (F), carbonate rock (CARBONATE), clay (CLAY) and other (OTHERS) mineralogical composition relative amounts and the indicatrix of brittleness index can be with It is expressed as:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep side To RS represents shallow lateral, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, and CARBONATE represents carbonate rock The relative amount of (main including mineral such as calcite and dolomites), CLAY represent the relative amount of clay, a, b, c, d, e, f, g It is constant for the parameter that multiple linear regression analysis method is fitted.
It will be compared by the result of each mineral constituent relative amount and experiment test in the tight rock for calculating of logging well, And ensureing its correlation more than 0.7, correlation is higher, illustrates that fitting effect is preferable.
3) each mineral constituent relative amount and brittleness index in seismic inversion identification tight rock.
Using each mineral constituent relative amount in above-mentioned reconstruct tight rock and the indicatrix of brittleness index as constraints, By analyzing prospect pit, preferably go out rational sample number and suitable radio-frequency component, rational initial model is established, with introductory die Based on type, inverting is indicated with seismic waveform, first the company's of being finally inversed by well profile different minerals component relative amount and brittleness index Spatial distribution, by inversion result and reconstruct indicatrix compare, if inversion result and match result on well are preferable, Indicate that inverting finally obtains each mineral constituent relative amount and the inverting data volume of brittleness index in rock by seismic waveform.
Above-mentioned seismic waveform instruction inverting takes full advantage of the intensive seismic waveform information of spatial distribution in the horizontal, vertical High-resolution well-log information is taken full advantage of upwards.Seismic waveform is densely distributed space structure data, and it is heavy to reflect The spatial variations of product environment and lithology combination, seismic waveform instruction inverting utilize the preferably related well sample sheet of seismic waveform similitude, Sample for reference spatial distribution distance and curve distribution feature establish initial model, instead of variogram analysis space variant structure, Unbiased optimal estimation is carried out to radio-frequency component.
Specifically, which can include following three steps:
S1:Loading experiment analysis of data, well-log information and seismic data, and by the experimental analysis data in seismic data It is demarcated on well-log information;
S2:It is preferred that related well sample sheet, sample for reference spatial distribution distance and curve distribution feature establish initial model;
S3:Preferably go out near well point with the most similar seismic waveform of indicatrix, using global optimization approach, promote To no wellblock, inverting certainty greatly enhances, from completely random to progressively determining.By company's well profile of inverting and indicatrix value Comparison is repeated, until establishing suitable geological model, geological model is applied to entire seismic data cube, it is anti-so as to reach Drill the purpose of each mineral constituent relative amount and brittleness index data volume in tight rock.
Method is determined to each component content in above-mentioned tight rock and brittleness index with a specific experimental implementation process It illustrates, comprises the following steps:
Step 1:Each mineral constituent relative amount and brittleness index in well logging quantitative judge tight rock.
First, log parameter sensitivity analysis is done to each mineral constituent relative amount and brittleness index in tight rock, It is secondary by analyzing distinct methods, preferably go out multiple linear regression analysis method reconstruct tight rock in each mineral constituent relative amount and Brittleness index indicatrix.As shown in Fig. 2, for using multiple linear regression analysis method reconstruct L1 well tight rock Mineralss and Brittleness index indicatrix can sum up each mineral constituent relative amount logging response character in tight rock by the figure, That is, quartzy (Q) relative amount electrical property feature is more apparent, and resistivity value is higher, natural gamma value, middle subvalue and density value compared with It is low.More than one high three low logging response character be quartzy relative amount main logging response character;Clay mineral (CLAY) Relative amount is more sensitive to natural gamma value, density value, interval transit time and resistivity, shows as natural gamma value, density value Higher with interval transit time, resistivity value is relatively low.The above are quartz and clay mineral relative amount logging response character, other ore deposits Object component relative amount and brittleness index logging response character are different.
Using the method for multiple linear regression, each mineral constituent relative amount in tight rock can be reconstructed, for example, stone English relative amount (Q), feldspar relative amount (F), clay mineral relative amount (CLAY) and brittleness index (BI) and main well logging The relation of response characteristic curve reconstructs indicatrix, and basis is provided for subsequent inverting work.
Quartz content (Q), feldspar content (F), clay mineral content (CLAY) and brittleness index (BI) feature of reconstruct are bent Line is as follows:
Q=0.31*AC-1.23*CN-66.52*DEN+0.02*GR-0.04*RD+0.14*RS+136. 61;
F=-0.06*AC+1.44*CN-22.28*DEN-0.34*GR+0.002*RD+0.06*RS+97 .55;
CLAY=0.10*AC+3.28*CN+40.33*DEN-0.26*GR-0.08*RD+0.234*RS- 141.63;
It is analyzed by cross plot as shown in Figure 3, it can be with comparing calculation value (the fine and close rock calculated using log Each mineral constituent relative amount in stone) with measured value (passing through X diffraction total rocks analysis result) correlation, correlation is higher, shows Precision of prediction and effect are better.As shown in figure 3, quartz, feldspar, clay mineral and other mineral facies to content correlation all compared with Height, substantially more than 0.75, in this way, reconstruct quartz and feldspar relative amount indicatrix calculated value are related to measured value Property is more preferable, and related coefficient is up to more than 0.9.
Step 2:Each mineral constituent relative amount and brittleness index in seismic inversion identification rock:
Analyzed with reference to the logging character of step 1, with each mineral constituent relative amount in the tight rock that is reconstructed in step 1 and The indicatrixes such as brittleness index are constraints, preferably go out rational sample number and suitable radio-frequency component, are established rational initial Model crosses each mineral constituent relative amount and brittleness index on well profile and even well profile with seismic waveform instruction inversion prediction Spatial distribution, wherein, spatial distribution refers to the geology distribution rule of each mineral constituent relative amount and brittleness index in the plane Rule.If efficiency of inverse process is preferable with coincideing on well, indicate that inverting finally obtains each mineral group in tight rock by seismic waveform Divide relative amount and the inverting data volume of brittleness index, wherein, inverting data volume is to have passed through seismic wave to original earthquake data The data volume that shape inverting obtains, this data volume can have geologic rule that can also not have geologic rule.Divide as shown in Figs. 4-6 Not Wei in the inverting data volume quartzy relative amount (Q), feldspar relative amount (F), brittleness index (BI) it is anti-by seismic waveform The company's well profile drilled, the quartzy relative amount (Q) of the reconstruct such as the log in Fig. 4 to 6 respectively Jing Guo squared processing, The indicatrix of feldspar relative amount (F), brittleness index (BI).It can be seen that by company's well profile in Fig. 4 to 6 and pass through seismic wave Quartzy relative amount (Q), feldspar relative amount (F) and brittleness index (BI) obtained by shape inverting and statistical result (well column on well The signature logging curve on sub- side) it is consistent substantially, this explanation is with each in multiple linear regression analysis method reconstruct tight rock Mineral constituent relative amount and brittleness index indicatrix are indicated inverting by seismic waveform, fine and close rock can be gone out with Accurate Prediction Each mineral constituent relative amount and brittleness index distribution characteristics and development degree in stone.
However, due to each mineral constituent relative amount in rock and brittleness index and wave impedance relation unobvious, in addition It is difficult on seismic profile, it is opposite that each mineral constituent in rock cannot be preferably predicted on conventional wave impedance inversion data volume Content and brittleness index distribution characteristics and development degree, and use each mineral group in multiple linear regression analysis method reconstruct tight rock Divide relative amount and brittleness index indicatrix, inverting is indicated by waveform, each mineral in tight rock can be gone out with Accurate Prediction Component relative amount and brittleness index distribution characteristics and development degree, this is also the advantage place of the application.
According to each mineral constituent relative amount and brittleness index in the rock as shown in Fig. 7 to 9 on inversion result Plane distribution situation and development degree carry out " dessert " reservoir prediction under compact reservoir background system evaluation, and brittleness are produced Raw reason has carried out rational geologic interpretation, and the exploration portion of fine and close oil gas and shale oil gas is instructed by these achievements of acquirement Administration, wherein, Fig. 7 show the plane isopleth distribution schematic diagram of quartz content (Q), and Fig. 8 show the plane of feldspar content (F) Distribution of contours schematic diagram, Fig. 9 show the plane isopleth distribution schematic diagram of brittleness index (BI).Because by quartz content or Based on carbonate rock content, the Brittleness of generation can generate failure by shear in fracturing process and advantageously form complexity Netted seam by volume fracturing technology, can improve well yield.It is special by feldspar, particularly brittleness caused by potassium feldspar content Sign, in deep hollow band or sloped region, secondary pore growing, secondary pore growing zone is lower step exploration targets.
In upper example, experimental analysis technology (analysis of X diffraction total rock), characteristic curve construction technology and earthquake mainly make use of Waveform inversion technology, based on X diffraction total rock analysiss of data, using Logging Curves, with the method weight of multiple linear regression The indicatrixes such as each mineral constituent content and brittleness index in structure tight rock, and using the indicatrix of reconstruct as constraint item Part predicts rock Mineralss relative amount and brittleness index spatial distribution characteristic with seismic waveform instruction inversion technique And development degree.
Mainly it is considered that brittle mineral relative amount and brittleness index and wave impedance relation unobvious in compact reservoir, Be difficult on seismic profile, although wave impedance inversion and more attribution inversion part well point prediction results are preferable, it is phased not Substantially, no wellblock precision can not meet exploration needs.In this example, with seismic waveform inversion technology, sample for reference space point It is more demanding to the selection of sample number and the selection of radio-frequency component when cloth distance and curve distribution feature establish initial model.Tool Body, it is analyzed based on X diffraction total rock, by log analysis, chooses the well logging to each mineral constituent content sensitivity in tight rock Curve reconstructs indicatrix with the method for multiple linear regression;Each mineral constituent relative amount and brittleness index in rock Prediction needs to carry out sensitivity curves analysis on the basis of well log interpretation, logging explanation and formation testing conclusion, and then selection is suitable Sample number and radio-frequency component establish rational initial model, with each in seismic waveform instruction inversion technique prediction compact reservoir Mineral constituent relative amount and brittleness index spatial feature and development degree.
It is unknown that above-described embodiment solves brittle mineral relative amount and brittleness index and wave impedance relation in compact reservoir Aobvious, the difficult point being difficult on seismic profile is combined by experiment-well logging-seismic inversion, Accurate Prediction brittle mineral Relative amount and brittleness index distribution and development degree, the skill being combined breakthroughly with well log interpretation with seismic inversion Art effectively predicts brittle mineral relative amount and brittleness index distribution characteristics in compact reservoir.The Predicting Technique has obtained well Application effect, play production actual effect well, pass through the brittle mineral relative amount and brittleness index of the technological prediction point Cloth feature and the goodness of fit on well are higher, and poor oil reservoir, oil reservoir feldspar content are high, and brittleness index is high.Especially in commercial oil Jing Chu, This feature becomes apparent, and illustrates that the main reason for brittleness generates in the example is feldspar content.
Based on same inventive concept, a kind of tight rock component relative amount and brittleness are additionally provided in the embodiment of the present invention Index determining device, as described in the following examples.Due to tight rock component relative amount and brittleness index determining device solution Certainly the principle of problem determines that method is similar to tight rock component relative amount and brittleness index, therefore tight rock component is opposite The implementation of content and brittleness index determining device may refer to tight rock component relative amount and brittleness index determines method Implement, overlaps will not be repeated.Used below, term " unit " or " module " can realize the software of predetermined function And/or the combination of hardware.Although following embodiment described device is preferably realized with software, hardware or soft The realization of the combination of part and hardware is also what may and be contemplated.Figure 10, which is that the tight rock component of the embodiment of the present invention is opposite, to be contained A kind of structure diagram of amount and brittleness index determining device, as shown in Figure 10, including:Establish module 1001, reconstructed module 1002 With inverting module 1003, the structure is illustrated below.
Module 1001 is established, for being based on the analysis of X diffraction total rock, it is special to establish each mineral constituent relative amount in tight rock Levy curve;
Reconstructed module 1002, it is special for reconstructing each mineral constituent relative amount indicatrix and brittleness index in tight rock Levy curve;
Inverting module 1003 refers to for each mineral constituent relative amount indicatrix and brittleness in the tight rock with reconstruct Number indicatrix carries out seismic waveform instruction inverting, obtains each mineral constituent relative amount in tight rock as constraints And the inverting data volume of brittleness index.
In one embodiment, reconstructed module 1002 is specific can be used for, by logging sensitivity analy-sis, utilizing Multiple linear regression analysis method reconstructs each mineral constituent relative amount indicatrix in tight rock.
In one embodiment, each mineral constituent relative amount feature is bent in the tight rock that reconstructed module 1002 reconstructs Line can be expressed as:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep side To RS represents shallow lateral, and a, b, c, d, e, f, g are the parameter after reconstruct, is constant.
In one embodiment, reconstructed module 1002 specifically can be used for each mineral group in the tight rock according to reconstruct Split-phase reconstructs the brittleness index indicatrix in tight rock to content characteristics curve, wherein, the brittleness of the tight rock refers to Number indicatrix is expressed as:
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, CARBONATE Represent the relative amount of carbonate rock, CLAY represents the relative amount of clay.
In one embodiment, inverting module 1003 can include:Determination unit, in the tight rock with reconstruct Each mineral constituent relative amount indicatrix and brittleness index indicatrix are as constraints, by analyzing prospect pit, Determine the sample number and radio-frequency component of model foundation;Unit is established, for according to the sample number and radio-frequency component determined, building Vertical initial model;First inverting unit, for based on the initial model, indicating inverting by seismic waveform, being finally inversed by Cross the spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index on well profile and even well profile;Comparison is single Member described crosses well profile and the even spatial distribution of the relative amount of each mineral constituent and crisp on well profile for obtain inverting The spatial distribution of sex index, each mineral constituent relative amount indicatrix and brittleness index indicatrix with the reconstruct, into Row comparison;Second inverting unit, should by the initial model in the case of being less than predetermined error threshold in definite difference Each mineral constituent relative amount and the inverting data volume of brittleness index in rock are obtained for entire seismic data cube inverting.
It can be seen from the above description that the embodiment of the present invention realizes following technique effect:It reconstructs in tight rock Each mineral constituent relative amount indicatrix and brittleness index indicatrix, and with each mineral constituent phase in the tight rock of reconstruct To content characteristics curve and brittleness index indicatrix as constraints, seismic waveform instruction inverting is carried out, obtains fine and close rock Each mineral constituent relative amount and the inverting data volume of brittleness index, solve each mineral in compact reservoir in the prior art in stone Component relative amount and brittleness index relation unobvious and cause be difficult on seismic profile the technical issues of, to have reached letter Single technique effect for accurately determining relation between each mineral constituent relative amount and brittleness index.
Obviously, those skilled in the art should be understood that each module of the above-mentioned embodiment of the present invention or each step can be with It is realized with general computing device, they can concentrate on single computing device or be distributed in multiple computing devices On the network formed, optionally, they can be realized with the program code that computing device can perform, it is thus possible to by it Store and performed in the storage device by computing device, and in some cases, can be held with the order being different from herein They are either fabricated to each integrated circuit modules or will be multiple in them by the shown or described step of row respectively Module or step are fabricated to single integrated circuit module to realize.In this way, the embodiment of the present invention be not restricted to it is any specific hard Part and software combine.
The foregoing is only a preferred embodiment of the present invention, is not intended to limit the invention, for the skill of this field For art personnel, the embodiment of the present invention can have various modifications and variations.Within the spirit and principles of the invention, made Any modification, equivalent substitution, improvement and etc. should all be included in the protection scope of the present invention.

Claims (6)

1. a kind of tight rock component relative amount and brittleness index determine method, which is characterized in that including:
It is analyzed based on X diffraction total rock, establishes each mineral constituent relative amount indicatrix in tight rock;
Reconstruct each mineral constituent relative amount indicatrix and brittleness index indicatrix in tight rock;The densification of the reconstruct Each mineral constituent relative amount indicatrix is expressed as in rock:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep lateral, RS Represent shallow lateral, a, b, c, d, e, f, g are the parameter after reconstruct, are constant;
Each mineral constituent relative amount indicatrix and brittleness index indicatrix are as constraint item using in the tight rock of reconstruct Part carries out seismic waveform instruction inverting, obtains the inverting data of each mineral constituent relative amount and brittleness index in tight rock Body, including:Each mineral constituent relative amount indicatrix and brittleness index indicatrix are as about using in the tight rock of reconstruct Beam condition by analyzing prospect pit, determines the sample number and radio-frequency component of model foundation;According to the sample number determined And radio-frequency component, establish initial model;Based on the initial model, inverting is indicated by seismic waveform, was finally inversed by well The spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index on section and even well profile;Inverting is obtained It is described cross well profile and the even spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index on well profile, With each mineral constituent relative amount indicatrix and brittleness index indicatrix of the reconstruct, compared;If difference is small In predetermined error threshold, then the initial model is obtained into each mineral constituent in rock applied to entire seismic data cube inverting The inverting data volume of relative amount and brittleness index.
2. the method as described in claim 1, which is characterized in that each mineral constituent relative amount feature is bent in reconstruct tight rock Line, including:
By logging sensitivity analy-sis, using multiple linear regression analysis method, it is opposite to reconstruct each mineral constituent in tight rock Content characteristics curve.
3. the method as described in claim 1, which is characterized in that the brittleness index indicatrix in reconstruct tight rock, including:
According to each mineral constituent relative amount indicatrix in the tight rock of reconstruct, the brittleness index reconstructed in tight rock is special Curve is levied, wherein, the brittleness index indicatrix of the tight rock is expressed as:
<mrow> <mi>B</mi> <mi>I</mi> <mo>=</mo> <mfrac> <mrow> <mi>Q</mi> <mo>+</mo> <mi>F</mi> <mo>+</mo> <mi>C</mi> <mi>A</mi> <mi>R</mi> <mi>B</mi> <mi>O</mi> <mi>N</mi> <mi>A</mi> <mi>T</mi> <mi>E</mi> </mrow> <mrow> <mi>Q</mi> <mo>+</mo> <mi>F</mi> <mo>+</mo> <mi>C</mi> <mi>A</mi> <mi>R</mi> <mi>B</mi> <mi>O</mi> <mi>N</mi> <mi>A</mi> <mi>T</mi> <mi>E</mi> <mo>+</mo> <mi>C</mi> <mi>L</mi> <mi>A</mi> <mi>Y</mi> </mrow> </mfrac> <mo>&amp;times;</mo> <mn>100</mn> </mrow>
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, and CARBONATE is represented The relative amount of carbonate rock, CLAY represent the relative amount of clay.
4. a kind of tight rock component relative amount and brittleness index determining device, which is characterized in that including:
Module is established, for being based on the analysis of X diffraction total rock, establishes each mineral constituent relative amount indicatrix in tight rock;
Reconstructed module, for reconstructing each mineral constituent relative amount indicatrix and brittleness index indicatrix in tight rock; Each mineral constituent relative amount indicatrix is expressed as in the tight rock of the reconstructed module reconstruct:
Y=a × AC+b × CN+c × DEN+d × GR+e × RD+f × RS+g
Wherein, AC represents interval transit time, and CN represents neutron, and DEN represents density, and GR represents natural gamma, and RD represents deep lateral, RS Represent shallow lateral, a, b, c, d, e, f, g are the parameter after reconstruct, are constant;
Inverting module, it is bent for each mineral constituent relative amount indicatrix in the tight rock with reconstruct and brittleness index feature Line carries out seismic waveform instruction inverting, obtains each mineral constituent relative amount and brittleness in tight rock and refer to as constraints Several inverting data volumes;Wherein, the inverting module includes:
Determination unit, it is bent for each mineral constituent relative amount indicatrix in the tight rock with reconstruct and brittleness index feature Line, by analyzing prospect pit, determines the sample number and radio-frequency component of model foundation as constraints;
Unit is established, for according to the sample number and radio-frequency component determined, establishing initial model;
First inverting unit, for based on the initial model, indicating inverting by seismic waveform, being finally inversed by well profile And connect the spatial distribution of relative amount of each mineral constituent and the spatial distribution of brittleness index on well profile;
Comparison unit, for the sky for crossing the relative amount of each mineral constituent on well profile and even well profile for obtaining inverting Between the spatial distribution of distribution and brittleness index, it is special with each mineral constituent relative amount indicatrix and brittleness index of the reconstruct Curve is levied, is compared;
Second inverting unit, in the case of being less than predetermined error threshold in definite difference, by the initial model application Each mineral constituent relative amount and the inverting data volume of brittleness index in rock are obtained in entire seismic data cube inverting.
5. device as claimed in claim 4, which is characterized in that the reconstructed module is specifically used for through log sensibility Analysis using multiple linear regression analysis method, reconstructs each mineral constituent relative amount indicatrix in tight rock.
6. device as claimed in claim 4, which is characterized in that the reconstructed module is specifically used for the tight rock according to reconstruct In each mineral constituent relative amount indicatrix, reconstruct tight rock in brittleness index indicatrix, wherein, it is described densification rock The brittleness index indicatrix of stone is expressed as:
<mrow> <mi>B</mi> <mi>I</mi> <mo>=</mo> <mfrac> <mrow> <mi>Q</mi> <mo>+</mo> <mi>F</mi> <mo>+</mo> <mi>C</mi> <mi>A</mi> <mi>R</mi> <mi>B</mi> <mi>O</mi> <mi>N</mi> <mi>A</mi> <mi>T</mi> <mi>E</mi> </mrow> <mrow> <mi>Q</mi> <mo>+</mo> <mi>F</mi> <mo>+</mo> <mi>C</mi> <mi>A</mi> <mi>R</mi> <mi>B</mi> <mi>O</mi> <mi>N</mi> <mi>A</mi> <mi>T</mi> <mi>E</mi> <mo>+</mo> <mi>C</mi> <mi>L</mi> <mi>A</mi> <mi>Y</mi> </mrow> </mfrac> <mo>&amp;times;</mo> <mn>100</mn> </mrow>
Wherein, BI represents brittleness index, and Q represents the relative amount of quartz, and F represents the relative amount of feldspar, and CARBONATE is represented The relative amount of carbonate rock, CLAY represent the relative amount of clay.
CN201511031926.8A 2015-12-31 2015-12-31 Tight rock component relative amount and brittleness index determine method and apparatus Active CN105675635B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511031926.8A CN105675635B (en) 2015-12-31 2015-12-31 Tight rock component relative amount and brittleness index determine method and apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511031926.8A CN105675635B (en) 2015-12-31 2015-12-31 Tight rock component relative amount and brittleness index determine method and apparatus

Publications (2)

Publication Number Publication Date
CN105675635A CN105675635A (en) 2016-06-15
CN105675635B true CN105675635B (en) 2018-06-01

Family

ID=56298489

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511031926.8A Active CN105675635B (en) 2015-12-31 2015-12-31 Tight rock component relative amount and brittleness index determine method and apparatus

Country Status (1)

Country Link
CN (1) CN105675635B (en)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106248494B (en) * 2016-08-29 2019-03-01 中国石油化工股份有限公司江汉油田分公司石油工程技术研究院 A method of for shale gas well reservoir brittleness overall merit
CN106547034B (en) * 2016-11-09 2018-11-02 西南石油大学 A kind of method of determining compact reservoir rock brittleness index
CN108519281B (en) * 2018-02-27 2020-08-11 中国石油天然气股份有限公司 Method, device and system for determining brittleness index of rock
CN110320571B (en) * 2018-03-29 2021-06-15 中国石油化工股份有限公司 Compact sandstone reservoir rock brittleness logging evaluation method
CN108825223B (en) * 2018-06-22 2021-10-22 西南石油大学 Method for extracting bedding characteristics of shale stratum
CN108827774B (en) * 2018-06-23 2019-06-07 东北石油大学 Coal petrography brittleness evaluation method
CN110094198A (en) * 2019-06-05 2019-08-06 中国地质大学(北京) Oil well grey matter content assaying method based on SPSS software
CN110513098B (en) * 2019-07-15 2022-07-01 中国海洋石油集团有限公司 Method for improving longitudinal resolution of element logging curve
CN113125413A (en) * 2019-12-31 2021-07-16 中石化石油工程技术服务有限公司 Shale brittleness online evaluation method based on laser-induced breakdown spectroscopy
CN111429012A (en) * 2020-03-27 2020-07-17 中国石油天然气股份有限公司大港油田分公司 Shale brittle dessert evaluation method
CN112505154B (en) * 2020-11-11 2021-06-29 中国地质大学(北京) Shale reservoir mineral component content analysis and lithofacies identification characterization method
CN114791482B (en) * 2021-01-25 2024-03-08 中国石油化工股份有限公司 Method and device for determining rock mineral content coefficient and storage medium
CN113030440B (en) * 2021-03-22 2022-02-15 西南石油大学 Method for predicting fracture brittleness index of shale in hydraulic fracturing process of shale

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
南方古生界3套富有机质页岩矿物组成与脆性分析;陈吉,肖贤明;《煤炭学报》;20130531;第38卷(第5期);第824页右栏第2-6段 *
基于岩石物理实验的致密油储层脆性指数预测—以柴西南A井区为例;张平 等;《中国石油学会2015年物探技术研讨会论文集》;20150513;第772页右栏第2段-第775页左栏第2段 *
致密地层岩石脆性指数的测井优化建模;范卓颖 等;《石油学报》;20151130;第36卷(第11期);全文 *

Also Published As

Publication number Publication date
CN105675635A (en) 2016-06-15

Similar Documents

Publication Publication Date Title
CN105675635B (en) Tight rock component relative amount and brittleness index determine method and apparatus
CN107817535B (en) The determination method and apparatus of short lap
Grana Joint facies and reservoir properties inversion
CN101930082B (en) Method for distinguishing reservoir fluid type by adopting resistivity data
Verma et al. Estimation of total organic carbon and brittleness volume
CN105510993A (en) Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method
CN102012526A (en) Method for discriminating type of reservoir fluid by using resistivity data
Wu et al. Seismic sedimentology, facies analyses, and high-quality reservoir predictions in fan deltas: a case study of the Triassic Baikouquan Formation on the western slope of the Mahu Sag in China's Junggar Basin
CN105116449B (en) A kind of recognition methods of weak reflection reservoir
CN105301647B (en) The method for distinguishing grey matter mud stone and sandstone
CN106249300B (en) It shakes to combine based on well and determines source rock TOC contents and anisotropism method and apparatus
CN102096098A (en) Seismic attribute extracting method
Qian et al. Prediction and modeling of petrophysical parameters of deep-buried, low permeability glutenite reservoirs in Yubei area, Turpan-Hami Basin, China
Leary et al. Power-law scaling of spatially correlated porosity and log (permeability) sequences from north-central North Sea Brae oilfield well core
Richardson Well correlation and Petrophysical analysis, a case study of “Rickie” field onshore Niger Delta
Ryberg et al. Tomographic V p and V s structure of the California Central Coast Ranges, in the vicinity of SAFOD, from controlled-source seismic data
Khan et al. Three-dimensional Structural Modeling (3D SM) and Joint Geophysical Characterization (JGC) of Hydrocarbon Reservoir: A Case Study of the Kadanwari field in Middle Indus Basin (MIB), Southeastern Pakistan
Gao et al. Seismic structure and texture analyses for fractured reservoir characterization: An integrated workflow
Wu et al. Combination of seismic attributes using clustering and neural networks to identify environments with sandstone-type uranium mineralization
CN112147676A (en) Method for predicting thickness of coal bed and gangue
Hu et al. Seismic interpretation of sandstone-type uranium deposits in the Songliao Basin, Northeast China
CN110174701A (en) A kind of prediction technique and device of Igneous rock invasion range
Yuan et al. Interpretation of fractures and joint inversion using multicomponent seismic data—Marcellus Shale example
Sun et al. 3D spatial characterization of sand body for uranium reservoirs based on geostatistical resistivity inversion
Dong et al. An introduction—Carbonate geophysics

Legal Events

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