CN106484989B - A method of coal rank type is quickly divided automatically using well-log information computer - Google Patents
A method of coal rank type is quickly divided automatically using well-log information computer Download PDFInfo
- Publication number
- CN106484989B CN106484989B CN201610871853.1A CN201610871853A CN106484989B CN 106484989 B CN106484989 B CN 106484989B CN 201610871853 A CN201610871853 A CN 201610871853A CN 106484989 B CN106484989 B CN 106484989B
- Authority
- CN
- China
- Prior art keywords
- coal
- rank type
- coal rank
- resistivity
- compensated neutron
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
Abstract
A method of coal rank type is quickly divided automatically using well-log information computer, first select the sensibility log parameter of different rank type: recycling makes full use of interval transit time, compensation density, compensated neutron and resistivity to construct view four kinds of wave impedance, compensated neutron enhancement coefficient and resistivity contrasts index combination logging parameters: and then establish three-dimensional coal rank type identification plate, coal rank type identification equation is resettled, finally identifies coal rank type；The present invention is based on the well-log informations such as the natural gamma of coal petrography, natural potential, resistivity, interval transit time, compensation density and compensated neutron, on the basis of analyzing the sensibility log parameter of different rank type, establish the combination loggings parameters such as coal petrography view wave impedance, compensated neutron enhancement coefficient, resistivity contrasts index, and three-dimensional coal rank type identification plate is constructed with this combination logging parameter, the sequencing plate and then division coal rank type will exploit for safety of coal mines and cbm development provide borehole logging technical support.
Description
Technical field
The invention belongs to the borehole logging assessment techniques during coal mine, coal-bed gas exploitation, in particular to a kind of to utilize survey
The method that well material computation machine quickly divides coal rank type automatically.
Background technique
In order to accurately divide coal rank type, it is automatic using well-log information computer that the purpose of the present invention is to provide a kind of
Quickly divide the method for coal rank type.Existing coal rank Type division method, mostly according to the vitrinite of coal petrography heart laboratory test
Reflectance value calculates the reflectance of vitrinite of coal petrography to divide, or using well-log information to divide coal rank type.Due to coal petrography
Heart data is very limited, and well-log information is commonly used in actual production to divide coal rank type.However, due to the vitrinite reflectance of coal petrography
Rate is influenced by many geologic(al) factors, and well-log information is difficult to accurately calculate the reflectance of vitrinite of coal petrography.In consideration of it, actually answering
It is relatively limited in.For this purpose, the scientific worker in the field turns at visual angle on crossplot chart.Vitrinite is carried out to the coal petrography heart
Albedo measurement and playback, and then two log parameter cross plots of different rank were established, but the method has highlighted two parameter and has been difficult to
The problems such as precisely dividing coal rank type, furthermore the method is not tried weak without the sensibility log parameter for different rank of investigating yet
Change or reduce the influence factor of log parameter, this makes troubles to safety of coal mines exploitation and coal bed gas Efficient Development detail design.
Summary of the invention
In order to overcome the shortcomings of above-mentioned existing method, well-log information computer is utilized the purpose of the present invention is to provide a kind of
The automatic method for quickly dividing coal rank type, natural gamma, natural potential, resistivity, interval transit time, compensation based on coal petrography are close
The well-log informations such as degree and compensated neutron establish coal petrography view on the basis of analyzing the sensibility log parameter of different rank type
The combination loggings parameters such as wave impedance, compensated neutron enhancement coefficient, resistivity contrasts index, and constructed with this combination logging parameter
Three-dimensional coal rank type identification plate, the sequencing plate and then division coal rank type will be opened for safety of coal mines exploitation and coal bed gas
Hair provides borehole logging technical support.
A method of coal rank type is quickly divided automatically using well-log information computer, comprising the following steps:
Step 1: the sensibility log parameter of different rank type selects: utilizing the coal rank of lab analysis chemical examination identification
Type, on the basis of coal core Location, natural gamma corresponding to the different rank type of extraction and analysis chemical examination, naturally electricity
Position, resistivity, compensated neutron, interval transit time and compensated density log data intersect its above-mentioned six kinds of log parameter two-by-two, system
To make two log parameter cross plot of coal rank type, is obtained by cross plot, natural potential, natural gamma identification coal rank type effect are poor,
Interval transit time, compensation density, compensated neutron and resistivity can identify coal rank type to a certain extent；
Step 2: establishing combination logging parameter: interval transit time, compensation density, compensated neutron and resistivity being made full use of
Building view four kinds of wave impedance, compensated neutron enhancement coefficient and resistivity contrasts index combination logging parameters:
RC=LOG (RLLD-RLLDS) (3)
In formula: Z is the view wave impedance of coal petrography, dimensionless；DEN is the density on stratum, g/cm3；When Δ t is the sound wave on stratum
Difference, μ s/m；NCFor the compensated neutron enhancement coefficient of coal petrography, dimensionless；CNL,CNLSRespectively the compensated neutron value of coal petrography, top bottom
The compensated neutron value of plate, %；RCFor resistivity contrasts index, dimensionless；RLLD,RLLDSThe respectively resistivity of coal petrography and top bottom
The resistivity of plate, Ω m；
Step 3: establishing three-dimensional coal rank type identification plate: regarding wave impedance, compensation relatively using what above-mentioned steps two constructed
Neutron enhancement coefficient and resistivity contrasts index construct three-dimensional coal rank type identification plate；
Step 4: coal rank type identification establishing equation: based on three-dimensional coal rank type identification plate constructed by step 3, with
Midpoint between different rank type is interface, and each coal rank type interface corresponds to an equation, then, above-mentioned three-dimensional coal rank
Type identification plate can determine that seven discriminant equations, as follows respectively:
Z=0.34629-0.00097432 × RC+0.01463×NC (4)
Z=0.37225-0.01829 × RC+0.042×NC (5)
Z=0.33984+0.002 × RC+0.02493×NC (6)
Z=0.08079+0.07686 × RC+0.03956×NC (7)
Z=0.394-0.00824 × RC-0.08756×NC (8)
Z=0.24223+0.04762 × RC-0.09047×NC (9)
Z=0.23031+0.02505 × RC+0.03576×NC (10)
Step 5: coal rank type identification: according to referring to depending on wave impedance, compensated neutron enhancement coefficient and resistivity contrasts for coal petrography
The determining data point of number falls in which delimited area of above-mentioned 7 equations, then can determine that coal rank type；If shown in above-mentioned steps two
Combination logging parameter computation model, the coal rank type identification equation computer sequencing in step 4, be articulated in logging data processing
On interpretation platform, it can realize that the computer of coal rank type automatically processes division.
The present invention is in order to find out the sensibility log parameter of different rank, and using intersection, drawing method is full and accurate analyzes different coals
Inherent pass between the parameters such as rank and natural potential, natural gamma, compensation density, resistivity, interval transit time and compensated neutron log
System；Based on log parameters such as interval transit time, compensation density, compensated neutron and the resistivity more sensitive to coal rank, energy is constructed
Three enough view wave impedance for effectively distinguishing coal rank type, resistivity contrasts coefficient and compensated neutron enhancement coefficient well logging combination ginsengs
Number, and organize well logging combination parameter using this and construct coal rank type three dimensional identification plate, and then the sequencing plate is next automatic fast
Speed identification coal rank type.Sensibility well logging of the method from acoustics, electricity and the preferred different rank type of radioactivity many aspects
Parameter, therefore, the invention the method more can accurately divide coal rank type, and then can be safety of coal mines exploitation and coal bed gas
Efficient Development provides logging technique and supports.
Detailed description of the invention
Fig. 1 is the method flow diagram for quickly dividing coal rank type automatically using well-log information computer in the present invention.
Fig. 2 is different rank type natural gamma and interval transit time relational graph in the present invention.
Fig. 3 is different rank type compensated neutron and compensation density relational graph in the present invention.
Fig. 4 is different rank type natural potential and resistivity relation figure in the present invention.
Fig. 5, which is that the different rank type in the present invention is opposite, regards wave impedance, compensated neutron enhancement coefficient and resistivity contrasts
Index three-dimensional identification figure.
Fig. 6 is that the coal rank type computer in the present invention divides result map automatically.
Specific embodiment
Technical solution of the present invention is described in detail with reference to the accompanying drawing.
Referring to Fig.1, a method of quickly dividing coal rank type, including following step automatically using well-log information computer
It is rapid:
Step 1: the sensibility log parameter of different rank type selects: utilizing the coal rank of lab analysis chemical examination identification
Type, on the basis of coal core Location, natural gamma corresponding to the different rank type of extraction and analysis chemical examination, naturally electricity
Position, resistivity, compensated neutron, interval transit time and compensated density log data intersect its above-mentioned six kinds of log parameter two-by-two, system
Make Fig. 1~two log parameter cross plot of coal shown in Fig. 3 rank type, thus organizes cross plot and obtain, remove natural potential, natural gamma
Identify that coal rank type effect is poor outer, interval transit time, compensation density, compensated neutron and resistivity can identify to a certain extent
Coal rank type, but the mixed class point between different rank type is more.It follows that identifying coal rank class using conventional two log parameter methods
Type precision is lower, it is difficult to meet the needs of actual production.
Step 2: establishing combination logging parameter: the precision for fully taking into account conventional two log parameter methods identification coal rank is lower
This objective factor, and interval transit time, compensation density, compensated neutron and resistivity can identify coal rank class to a certain extent
Type makes full use of interval transit time, compensation density, in compensation from the purpose for reducing logging error, promoting coal rank accuracy of identification
Sub and resistivity regards four kinds of wave impedance, compensated neutron enhancement coefficient and resistivity contrasts index combination logging parameters to construct.
RC=LOG (RLLD-RLLDS) (3)
In formula: Z is the view wave impedance of coal petrography, dimensionless；DEN is the density on stratum, g/cm3；When Δ t is the sound wave on stratum
Difference, μ s/m；NCFor the compensated neutron enhancement coefficient of coal petrography, dimensionless；CNL,CNLSRespectively the compensated neutron value of coal petrography, top bottom
The compensated neutron value of plate, %；RCFor resistivity contrasts index, dimensionless；RLLD,RLLDSThe respectively resistivity of coal petrography and top bottom
The resistivity of plate, Ω m.
Step 3: establishing three-dimensional coal rank type identification plate: referring to Fig. 5, based on the thought of geology clustering, identical coal
The logging response character of rank type coal petrography is almost the same, then using in the opposite view wave impedance of the building of above-mentioned steps two, compensation
Sub- enhancement coefficient and resistivity contrasts index construct three-dimensional coal rank type identification plate.
Step 4: coal rank type identification establishing equation: based on three-dimensional coal rank type identification plate constructed by step 3, with
Midpoint between different rank type is interface, and each coal rank type interface corresponds to an equation, then, above-mentioned three-dimensional coal rank
Type identification plate can determine that seven discriminant equations, as follows respectively:
Z=0.34629-0.00097432 × RC+0.01463×NC (4)
Z=0.37225-0.01829 × RC+0.042×NC (5)
Z=0.33984+0.002 × RC+0.02493×NC (6)
Z=0.08079+0.07686 × RC+0.03956×NC (7)
Z=0.394-0.00824 × RC-0.08756×NC (8)
Z=0.24223+0.04762 × RC-0.09047×NC (9)
Z=0.23031+0.02505 × RC+0.03576×NC (10)
Step 5: coal rank type identification: referring to Fig. 5, if the view wave impedance of coal petrography, compensated neutron enhancement coefficient and resistance
The data point that rate differential index (di) determines falls in the region that equation (4) and equation (5) define, then can determine that as lignite；Successively similarly,
It can determine other kinds of coal petrography.By the coal in combination logging parameter computation model, step 4 shown in above-mentioned steps two
Rank type identification equation computer sequencing, is articulated on logging data processing interpretation platform, can realize the computer of coal rank type
Automatically process division.
The present invention is tried out in practical coalfield.In the coal rank Type division application of X well, referring to Fig. 6, the well main force coal
Interval 573.5-577.4m, with a thickness of 3.9m, without obvious dirt band in layer.From the figure the 6th it is found that coalbed methane reservoir 573.5-
575.8m well section is identified as meager coal, and 575.8-577.4m well section is identified as lean coal.Compare the 6th coal rank type identification result of figure
With the 7th the coal petrography heart determine coal rank type it is found that the method for the present invention identification coal rank type completely with coal core analysis
The coal rank type tested is completely the same.This method had not only analyzed preferably the sensibility log parameter of different rank type, but also constructed
Combination logging parameter and three-dimensional coal rank type identification plate, the coal rank type divided and practical coal rank are more coincide.In the method
Each combination logging parameter can be sought from coalfield borehole logging data, and almost all of coalfield all has largely
Borehole logging data.Therefore, coal rank Type division method of the present invention has good popularization and application foreground and value.
It will be understood by those of skill in the art that the coal seam buried depth and degree of metamorphism due to different regions differ greatly, it is
Coal rank type is more accurately divided, is very necessary for the different corresponding three-dimensional coal lithotypes identification plates of area building
, and well-log information has to pass through the correction method of surroundings effecting.
The method that the present invention divides coal rank type using well-log information makes full use of the well logging more sensitive with coal rank type
Parameter organically blends its log parameter to construct combination logging parameter, and the coal rank type divided and lab analysis are chemically examined
Coal rank type it is almost the same, coal rank Type division precision is greatly improved compared with conventional method.
Claims (1)
1. a kind of method for quickly dividing coal rank type automatically using well-log information computer, which is characterized in that including following step
It is rapid:
Step 1: the sensibility log parameter of different rank type selects: the coal rank type identified using lab analysis chemical examination,
On the basis of coal core Location, natural gamma, natural potential, resistance corresponding to the different rank type of extraction and analysis chemical examination
Rate, compensated neutron, interval transit time and compensated density log data intersect its above-mentioned six kinds of log parameter two-by-two, make coal rank
Thus two log parameter cross plot of type is organized cross plot and is obtained, natural potential, natural gamma identification coal rank type effect are poor, sound
The wave time difference, compensation density, compensated neutron and resistivity can identify coal rank type；
Step 2: establishing combination logging parameter: interval transit time, compensation density, compensated neutron and resistivity being made full use of to construct
Depending on three kinds of wave impedance, compensated neutron enhancement coefficient and resistivity contrasts index combination logging parameters:
RC=LOG (RLLD-RLLDS) (3)
In formula: Z is the view wave impedance of coal petrography, dimensionless；DEN is the density on stratum, g/cm3；Δ t is the interval transit time on stratum, μ
s/m；NCFor the compensated neutron enhancement coefficient of coal petrography, dimensionless；CNL,CNLSRespectively the compensated neutron value of coal petrography, roof and floor
Compensated neutron value, %；RCFor resistivity contrasts index, dimensionless；RLLD,RLLDSThe respectively resistivity of coal petrography and roof and floor
Resistivity, Ω m；
Step 3: establishing three-dimensional coal rank type identification plate: opposite the view wave impedance, compensated neutron constructed using above-mentioned steps two
Enhancement coefficient and resistivity contrasts index construct three-dimensional coal rank type identification plate；
Step 4: coal rank type identification establishing equation: based on three-dimensional coal rank type identification plate constructed by step 3, with difference
Midpoint between coal rank type is interface, and each coal rank type interface corresponds to an equation, then, above-mentioned three-dimensional coal rank type
Identification plate can determine that seven discriminant equations, as follows respectively:
Z=0.34629-0.00097432 × RC+0.01463×NC (4)
Z=0.37225-0.01829 × RC+0.042×NC (5)
Z=0.33984+0.002 × RC+0.02493×NC (6)
Z=0.08079+0.07686 × RC+0.03956×NC (7)
Z=0.394-0.00824 × RC-0.08756×NC (8)
Z=0.24223+0.04762 × RC-0.09047×NC (9)
Z=0.23031+0.02505 × RC+0.03576×NC (10)
Step 5: coal rank type identification: true according to the view wave impedance of coal petrography, compensated neutron enhancement coefficient and resistivity contrasts index
Fixed data point falls in which delimited area of above-mentioned seven discriminant equations, then can determine that coal rank type；If two institute of above-mentioned steps
Coal rank type identification equation computer sequencing in the combination logging parameter computation model that shows, step 4, is articulated at well logging
It manages on interpretation platform, can realize that the computer of coal rank type automatically processes division.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610871853.1A CN106484989B (en) | 2016-09-30 | 2016-09-30 | A method of coal rank type is quickly divided automatically using well-log information computer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610871853.1A CN106484989B (en) | 2016-09-30 | 2016-09-30 | A method of coal rank type is quickly divided automatically using well-log information computer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106484989A CN106484989A (en) | 2017-03-08 |
| CN106484989B true CN106484989B (en) | 2019-06-07 |
Family
ID=58269036
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610871853.1A Expired - Fee Related CN106484989B (en) | 2016-09-30 | 2016-09-30 | A method of coal rank type is quickly divided automatically using well-log information computer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106484989B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109031461A (en) * | 2018-09-21 | 2018-12-18 | 西南石油大学 | Carbonate porosity type oil and gas reservoir quantitative identification method and region quantitative method |
| CN111141866B (en) * | 2020-01-14 | 2022-11-11 | 山东科技大学 | Test device and method for simulating spontaneous combustion process of coal seam and detection of high-temperature abnormal area |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1412575A (en) * | 2001-10-11 | 2003-04-23 | 中国石油天然气集团公司西北地质研究所 | Oil-gas detection method based on multiphase medium theory |
| CN103592687A (en) * | 2013-11-22 | 2014-02-19 | 中国石油化工集团公司 | Quantitative calculation method for coal and rock adsorbed gas content |
| CN104704089A (en) * | 2012-10-01 | 2015-06-10 | 格雷特波因特能源公司 | Agglomerated particulate low-rank coal feedstock and uses thereof |
| CN105761159A (en) * | 2016-03-17 | 2016-07-13 | 成都创源油气技术开发有限公司 | Stratum fracture and cavern development single well section modeling method |
-
2016
- 2016-09-30 CN CN201610871853.1A patent/CN106484989B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1412575A (en) * | 2001-10-11 | 2003-04-23 | 中国石油天然气集团公司西北地质研究所 | Oil-gas detection method based on multiphase medium theory |
| CN104704089A (en) * | 2012-10-01 | 2015-06-10 | 格雷特波因特能源公司 | Agglomerated particulate low-rank coal feedstock and uses thereof |
| CN103592687A (en) * | 2013-11-22 | 2014-02-19 | 中国石油化工集团公司 | Quantitative calculation method for coal and rock adsorbed gas content |
| CN105761159A (en) * | 2016-03-17 | 2016-07-13 | 成都创源油气技术开发有限公司 | Stratum fracture and cavern development single well section modeling method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106484989A (en) | 2017-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104612675B (en) | A kind of carbonate formation Lithology while drilling method for quickly identifying | |
| Verma et al. | Estimation of total organic carbon and brittleness volume | |
| CN101738639B (en) | Method for improving computing precision of rock fracture parameters | |
| CN103993871B (en) | Method and device for processing well logging information of thin interbed stratums in standardization mode | |
| CN109085663A (en) | A kind of tight sandstone reservoir stratification seam recognition methods | |
| CN110318745B (en) | Particle size lithology logging evaluation method under deposition microphase constraint | |
| CN105697002A (en) | Method for recognizing coal measure strata lithology | |
| Chen et al. | Identification of thin-layer coal texture using geophysical logging data: Investigation by wavelet transform and linear discrimination analysis | |
| Euzen et al. | Well log cluster analysis: an innovative tool for unconventional exploration | |
| CN104991286A (en) | Sedimentary facies characterization method based on sedimentary modes | |
| CN106484989B (en) | A method of coal rank type is quickly divided automatically using well-log information computer | |
| CN112946751B (en) | Method and system for predicting plane distribution of different reservoirs | |
| CN105350959B (en) | The method that shale gas reservoir gas-bearing saturation degree is determined by WELL LITHOLOGY density | |
| CN110529106B (en) | Method for determining content of coal seam micro-components by using logging information | |
| CN104234706B (en) | A kind of coal bearing property well logging sensitive parameter evaluation method | |
| CN112505754A (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
| CN109826623A (en) | Knowledge method is sentenced in a kind of geophysical log of tight sandstone reservoir stratification seam | |
| CN112946782B (en) | Earthquake fine depicting method for dense oil-gas storage seepage body | |
| Stadtműller et al. | Characterization of the carbonate formation fracture system based on well logging data and results of laboratory measurements | |
| CN113484907B (en) | Method for predicting distribution on different types of reservoir planes | |
| CN114861515A (en) | Method, device, equipment and medium for calculating layer speed data volume | |
| Kato et al. | Sweet spot mapping in the Montney tight gas reservoir | |
| CN112528106A (en) | Volcanic lithology identification method | |
| Zhang et al. | Identifying minor faults on top of coalfield Ordovician limestone stratum using seismic attributes derived from azimuthally stacked data | |
| CN110454155A (en) | A kind of method determining the Quaternary Stratigraphic age using magnetic susceptibility log method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190607 Termination date: 20190930 |