CN108333626A - A kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy - Google Patents
A kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy Download PDFInfo
- Publication number
- CN108333626A CN108333626A CN201810021231.9A CN201810021231A CN108333626A CN 108333626 A CN108333626 A CN 108333626A CN 201810021231 A CN201810021231 A CN 201810021231A CN 108333626 A CN108333626 A CN 108333626A
- Authority
- CN
- China
- Prior art keywords
- population
- individual
- wave impedance
- impedance inversion
- genetic algorithm
- 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
- 230000002068 genetic effect Effects 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000014759 maintenance of location Effects 0.000 title claims abstract description 16
- 238000010353 genetic engineering Methods 0.000 claims abstract description 8
- 230000003044 adaptive effect Effects 0.000 claims abstract description 6
- 230000035772 mutation Effects 0.000 claims abstract description 6
- 230000008569 process Effects 0.000 claims abstract description 4
- 108090000623 proteins and genes Proteins 0.000 claims description 4
- 230000003595 spectral effect Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 3
- 230000006978 adaptation Effects 0.000 claims description 2
- 230000010429 evolutionary process Effects 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims 2
- 230000002028 premature Effects 0.000 abstract description 4
- 239000000203 mixture Substances 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6226—Impedance
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
A kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy, belong to Geophysics Inversion field, especially with regard to the wave impedance inversion technique in oil geophysical exploration, it is intended to provide a kind of improved adaptive GA-IAGA Optimum Impedance Inversion Method, for solving premature problem of the standard genetic algorithm in wave impedance inversion, so that inversion result is relatively reliable, mainly include the following steps：1. building the object function of wave impedance inversion according to convolution model；2. the estimation of seismic wavelet by homomorphic theory；3. being encoded to wave impedance using binary coding mode；4. calculating the fitness value of each individual with object function, and quantitatively evaluating is carried out to individual according to this；5. generating population of new generation according to the selection mode of best retention strategy；6. carrying out genetic manipulation according to the intersection of design, mutation operator；7. idiotype is converted to phenotype according to corresponding decoding process and realizes that algorithm recycles；8. seeking wave impedance using recurrence method.
Description
Technical field
The present invention provides a kind of genetic algorithm Optimum Impedance Inversion Methods based on best retention strategy, belong to geophysics
It is calculated especially with regard to the wave impedance inversion technique in oil geophysical exploration for solving standard genetic in inverting field
Premature problem of the method in wave impedance inversion.
Background technology
There are the weakness that Premature Convergence and convergence rate are slow, current improved methods in practical application for standard genetic algorithm
Mostly it is directed to macrotechnique, genetic manipulation, genetic operator improvement and the parallelization operation etc. of population.In coding mode,
Dynamic coding mode, floating-point encoding mode, Gray code mode etc., in terms of the macrotechnique of population, Yao Wei powder are developed
(2015) propose dynamic microhabitat coevolution model, and in terms of the improvement of selection opertor, Zhang Jing (2015) is proposed
Adaptive sequencing selection mode, in terms of crossover operator improvement, Davis (1991) proposes serial number crossover operator and uniformly row
Sequence crossover operator, in terms of the improvement of mutation operator, Liu Li (2015) proposes gene position TSP question genetic operator, on
The search performance that innovatory algorithm improves genetic algorithm to varying degrees is stated, but is directed to the particular problem of wave impedance inversion or more
Compared to standard genetic algorithm, implementation process is cumbersome for corrective measure, poor for different problem adaptability, in terms of convergence not
Clear superiority can be embodied.The present invention sets about from selection opertor, proposes a kind of genetic algorithm based on optimal selection strategy, can have
Effect overcomes the problems, such as the premature problem in wave impedance inversion, and adaptability is good, easy to implement.
Invention content
The present invention is intended to provide a kind of overcome standard genetic algorithm to be lost in the improvement of wave impedance inversion mid-early maturity convergence problem
Propagation algorithm, it is to ensure that other genetic operators are constant on the basis of standard genetic algorithm, roulette wheel selection mode is abandoned, using one
The best retention strategy selection mode of kind.
The present invention specific steps include：
(1) initialization operation, setting control parameter and generation initial population, and calculate the fitness value of population.
(2) individual in parent population is ranked up by fitness value size, seeks the average fitness value of population, will fits
The individual that should be worth more than average fitness value is genetic directly in next-generation population.
(3) using highest fitness value as template, fitness value makees cross-correlation judgement with the individual with highest fitness, will
The individual that fitness value is high and difference in correlation is larger forms new population.
(4) according to principle in (3), gradually with the high individual of fitness value for template, the individual composition of different templates is selected
New population.
(5) judge whether to reach population scale, if it is, carrying out the genetic manipulations such as next step intersection, variation, otherwise will
The individual of removal sequentially supplies the scarce quantity of population institute by fitness value size, until reaching population scale.
The present invention is a kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy, has following features：
(1) genetic algorithm based on best retention strategy takes a kind of be based on by the sequence of individual adaptation degree size, mutually
The selection mode of the operations such as pass, selection ensure that in each evolutionary process, filial generation can retain the optimized individual in parent,
The diversity and intersection, the stability of mutation operation for ensureing population gene are avoided and are lost based on the standard under roulette wheel selection mode
The precocious phenomenon that propagation algorithm generates, makes algorithm that may finally search globally optimal solution.
(2) compared to other Revised genetic algorithums, realize that simply algorithm stability is good for wave impedance inversion problem,
Calculating speed is fast.
Description of the drawings
Fig. 1 and Fig. 2 is respectively the single-channel seismic record of standard genetic algorithm and improved adaptive GA-IAGA inverting, from inversion result
On see, Fig. 1 is that standard genetic algorithm iteration 11 times restrains obtained inverting record, is 94.4% with original record related coefficient, partially
Difference is larger；Fig. 2 is 32 inverting records for restraining to obtain of improved adaptive GA-IAGA iteration, is with original seismic data related coefficient
99.6%, the wave impedance that recurrence method acquires is coincide with actual well drilled.It theoretically analyzes, standard genetic algorithm is due to using roulette wheel
Selection mode generates higher fitness individual x at evolution initial stage, and other individuals are eliminated rapidly, most of individual and x phases
Together, intersected, mutation operation individual be in a disadvantageous position in competition, be easy to be eliminated, in initial stage population of evolving owning
Individual is absorbed in same extreme value and stops evolving.And improved adaptive GA-IAGA is by pressing the sequence of fitness size, cross-correlation judges, template
The operation such as selection improves the diversity of population gene under the premise of ensureing population's fitness, ensure that follow-up intersection, variation
The stability of operation may finally search globally optimal solution, inversion accuracy higher.
Specific implementation mode
A kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy, specific implementation step are as follows：
(1) object function is built
The object function of wave impedance inversion is built according to convolution model：
In formula, D is real seismic record, and W (t) is seismic wavelet, and Z is wave impedance, and R is reflectance factor, and N is reflectance factor
Sequence length, Δ t are the sampling interval.
Build wave impedance recursive function：
R is reflectance factor in formula, and Z is wave impedance.
(2) the estimation of seismic wavelet by homomorphic theory, by reflection coefficient sequence intermediary heat spectrum and earthquake in intermediary heat spectral domain
The intermediary heat spectrum of wave separates, and then obtains the intermediary heat spectral sequence of seismic wavelet, is then transformed to time-domain to get then
Between domain seismic wavelet.
(3) wave impedance is encoded using binary coding mode, determines the genotype X of individual.
(4) object function that formula (1) is built is used to calculate the fitness value of each individual and the average fitness of population
Value, and quantitatively evaluating is carried out to individual according to this.
(5) population of new generation is generated according to the selection mode of best retention strategy.
(6) genetic manipulation is carried out according to the intersection of design, mutation operator.
(7) idiotype is converted to phenotype and substitutes into end condition according to corresponding decoding process and judged,
Follow-up genetic manipulation is terminated if meeting end condition, otherwise returns to (4) step.
(8) wave impedance is asked using recurrence method.
Claims (5)
1. a kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy, it is characterised in that foring a set of completely newly has
The Optimum Impedance Inversion Method of effect, such as following steps：
(1) object function of wave impedance inversion is built according to convolution model：
In formula, D is real seismic record, and W (t) is seismic wavelet, and Z is wave impedance, and R is reflectance factor, and N is reflection coefficient sequence
Length, Δ t are the sampling interval.
Build wave impedance recursive function：
R is reflectance factor in formula, and Z is wave impedance.
(2) the estimation of seismic wavelet by homomorphic theory, by reflection coefficient sequence intermediary heat spectrum and seismic wavelet in intermediary heat spectral domain
Intermediary heat spectrum separates, and then obtains the intermediary heat spectral sequence of seismic wavelet, and wavelet is separated from seismic data, then will
It transforms to time-domain to get to time-domain seismic wavelet.
(3) initialization operation, setting control parameter and generation initial population, and calculate the fitness value of population.
(4) individual in parent population is ranked up by fitness value size, the average fitness value of population is sought, by adaptive value
Individual more than average fitness value is genetic directly in next-generation population.
(5) using highest fitness value as template, fitness value makees cross-correlation judgement with the individual with highest fitness, will adapt to
The individual that angle value is high and difference in correlation is larger forms new population.
(6) according to principle in (5), gradually the individual of different templates is selected to form newly for template with the high individual of fitness value
Population.
(7) judge whether to reach population scale, if it is, the genetic manipulations such as next step intersection, variation are carried out, it otherwise will removal
Individual sequentially supply the scarce quantity of population institute by fitness value size, until reaching population scale.
(8) genetic manipulation is carried out according to the intersection of design, mutation operator.
(9) idiotype is converted to phenotype and substitutes into end condition according to corresponding decoding process and judged, if
Meet end condition and then terminate follow-up genetic manipulation, otherwise returns to (4) step.
(10) wave impedance is asked using recurrence method.
2. a kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy according to claim 1, feature
It is：The object function established based on convolution model is compared to the object function that the wave equation based on wave theory is established
It is strong with noiseproof feature, the characteristics of algorithmic stability.
3. a kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy according to claim 1, feature
It is：Wavelet extraction method is not by borehole restraint, and the wavelet precision of extraction is higher, and fitting effect is good.
4. a kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy according to claim 1, feature
It is：The recurrence method inverting wave impedance of use has well without that can be used under the conditions of well, is suitable for exploration initial stage Wu Jing and well is few
The case where.
5. a kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy according to claim 1, feature
It is：Based on the best retention strategy selection mode by operations such as the sequence of individual adaptation degree size, cross-correlation, selections, each
In secondary evolutionary process, filial generation always remains individual best in parent, ensure that diversity and intersection, the variation of population gene
The stability of operation may finally search globally optimal solution.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810021231.9A CN108333626B (en) | 2018-01-10 | 2018-01-10 | A kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810021231.9A CN108333626B (en) | 2018-01-10 | 2018-01-10 | A kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108333626A true CN108333626A (en) | 2018-07-27 |
| CN108333626B CN108333626B (en) | 2019-06-18 |
Family
ID=62923504
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810021231.9A Active CN108333626B (en) | 2018-01-10 | 2018-01-10 | A kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108333626B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109918659A (en) * | 2019-02-28 | 2019-06-21 | 华南理工大学 | A method of based on not retaining optimum individual genetic algorithm optimization term vector |
| CN112115642A (en) * | 2020-09-14 | 2020-12-22 | 四川航天燎原科技有限公司 | High maneuvering platform SAR imaging parameter optimization design method |
| CN113504568A (en) * | 2021-07-09 | 2021-10-15 | 吉林大学 | Median filtering method based on niche differential evolution algorithm |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1710446A (en) * | 2005-06-21 | 2005-12-21 | 中国石油大学(北京) | Method for inversion constituting virtual well data using before-folded seismic wave form |
| US20070090836A1 (en) * | 2001-10-25 | 2007-04-26 | Intematix Corporation | Detection with evanescent wave probe |
| CN104977609A (en) * | 2014-04-11 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | Prestack longitudinal wave and transverse wave combined inversion method based on rapid simulated annealing |
| CN105445791A (en) * | 2015-11-25 | 2016-03-30 | 成都理工大学 | Stratum aperture pressure prediction method based on variety earthquake attributes |
| CN107462924A (en) * | 2017-07-27 | 2017-12-12 | 西安交通大学 | A kind of absolute wave impedance inversion method independent of well-log information |
-
2018
- 2018-01-10 CN CN201810021231.9A patent/CN108333626B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070090836A1 (en) * | 2001-10-25 | 2007-04-26 | Intematix Corporation | Detection with evanescent wave probe |
| CN1710446A (en) * | 2005-06-21 | 2005-12-21 | 中国石油大学(北京) | Method for inversion constituting virtual well data using before-folded seismic wave form |
| CN104977609A (en) * | 2014-04-11 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | Prestack longitudinal wave and transverse wave combined inversion method based on rapid simulated annealing |
| CN105445791A (en) * | 2015-11-25 | 2016-03-30 | 成都理工大学 | Stratum aperture pressure prediction method based on variety earthquake attributes |
| CN107462924A (en) * | 2017-07-27 | 2017-12-12 | 西安交通大学 | A kind of absolute wave impedance inversion method independent of well-log information |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109918659A (en) * | 2019-02-28 | 2019-06-21 | 华南理工大学 | A method of based on not retaining optimum individual genetic algorithm optimization term vector |
| CN109918659B (en) * | 2019-02-28 | 2023-06-20 | 华南理工大学 | Method for optimizing word vector based on unreserved optimal individual genetic algorithm |
| CN112115642A (en) * | 2020-09-14 | 2020-12-22 | 四川航天燎原科技有限公司 | High maneuvering platform SAR imaging parameter optimization design method |
| CN113504568A (en) * | 2021-07-09 | 2021-10-15 | 吉林大学 | Median filtering method based on niche differential evolution algorithm |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108333626B (en) | 2019-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108333626B (en) | A kind of genetic algorithm Optimum Impedance Inversion Method based on best retention strategy | |
| CN110346831B (en) | Intelligent seismic fluid identification method based on random forest algorithm | |
| CN109360604B (en) | Ovarian cancer molecular typing prediction system | |
| US9146903B2 (en) | Method of using spatially independent subsets of data to calculate vertical trend curve uncertainty of spatially correlated reservoir data | |
| WO2005098647A2 (en) | Genetic algorithm based selection of neural network ensemble for processing well logging data | |
| Mahmoudabadi et al. | A hybrid method for grade estimation using genetic algorithm and neural networks | |
| CN114723095A (en) | Missing well logging curve prediction method and device | |
| CN113553780A (en) | Stratum pore pressure prediction method based on machine learning | |
| CN115308799B (en) | Method and system for identifying free gas structure of seismic imaging | |
| CN111767674A (en) | Well logging lithology identification method based on active domain adaptation | |
| CN111914478A (en) | Comprehensive geological drilling well logging lithology identification method | |
| CN108663711B (en) | A kind of Bayes's seismic inversion method based on τ distribution | |
| CN111058840A (en) | Organic carbon content (TOC) evaluation method based on high-order neural network | |
| CN103278849B (en) | The method and system of higher-order spectra are carried out based on seismic data and well-log information | |
| WO2020096572A1 (en) | Dictionary generation for downhole signal compression | |
| CN108828680B (en) | Method and device for determining formation resistivity and polarizability | |
| CN111751878A (en) | Method and device for predicting transverse wave velocity | |
| CN115620831A (en) | Method for generating sequence mutation fitness through loop iteration optimization and related device | |
| CN112034515B (en) | Volcanic channel identification method based on unsupervised neural network | |
| CN116335648A (en) | Method for quantitatively evaluating dessert quality of conglomerate horizontal well based on logging while drilling data | |
| CN112649869A (en) | Reservoir characteristic parameter prediction method and system based on GA-WNN | |
| CN115759792A (en) | Horizontal well profile control and flooding effect analysis method based on whitening weighted gray model | |
| CN114075973A (en) | Stratum element logging curve reconstruction method and device | |
| RU2808168C1 (en) | Method for determining reservoir pressure in volume of field based on artificial neural networks | |
| CN116187206B (en) | COD spectrum data migration method based on generation countermeasure network |
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 |