CN106707345B - A kind of angle elastic parameter Lithology Identification Methods and device - Google Patents
A kind of angle elastic parameter Lithology Identification Methods and device Download PDFInfo
- Publication number
- CN106707345B CN106707345B CN201611144943.7A CN201611144943A CN106707345B CN 106707345 B CN106707345 B CN 106707345B CN 201611144943 A CN201611144943 A CN 201611144943A CN 106707345 B CN106707345 B CN 106707345B
- Authority
- CN
- China
- Prior art keywords
- angle
- elastic
- impedance
- wave
- shear wave
- 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
Links
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention provides a kind of angle elastic parameter Lithology Identification Methods and devices.This method includes：Obtain the converted shear wave elastic impedance equation with magnitude with longitudinal wave elastic impedance；Obtain converted shear wave elastic impedance puppet log and longitudinal wave elastic impedance puppet log；Obtain angle elastic parameter puppet log；Obtain best angle；The conventional more wave angle degree partial stack seismic channel sets of input；Layer position explanation is carried out to seismic channel set；Obtain converted shear wave elastic impedance data volume, the longitudinal wave elastic impedance data volume with angle information；Obtain angle elastic parameter data body；According to angle elastic parameter data body lithology value range, Lithology Discrimination is carried out, the identification of lithology is completed.The device includes：Conversion module；Log input module；Angle elastic parameter puppet log computing module；Best angle module；Multi-wave seismic data input module；Explanation module；Inverting module；Angle elastic parameter data module；Identification module.
Description
Technical field
The present invention relates to a kind of recognition methods and device more particularly to a kind of angle elastic parameter Lithology Identification Methods and dresses
It sets, belongs to oil earth exploration technical field.
Background technology
The emphasis of China's oil-gas exploration in recent years is gradually shifted to lithological reservoir exploration.Different from structural deposit, this
A little novel reservoirs are influenced by construction and reservoir heterogeneity, and accumulating condition is complicated, and identification difficulty is big, and investment risk is high.Therefore,
Developing a set of high-precision Lithology Discrimination geophysics solution seems particularly urgent.
The conventional technology for carrying out Lithology Discrimination using Prestack seismic data mainly has pre-stack seismic elastic parameter inversion method.
Elastic parameter inversion is based on prestack common offset trace gather or angle part is superimposed trace gather, sparse using maximal possibility estimation, constraint
Pulse isoinversion algorithm is finally inversed by the elastic parameter body comprising lithological information.
Due to the limitation of elastic parameter physical significance itself, there is certain limit, lithology to distinguish the sensibility of lithology
Ability can not meet the requirement of nowadays fining oil-gas exploration and conventional method does not make full use of the angle of seismic data
Converted shear wave earthquake information in information and multi-wave seismic.In fact, the angle information and converted shear wave of seismic data
Shake the lithological information that information includes abundant.Therefore, it is quick to cause current elastic parameter Lithology Identification Methods generally existing
The not strong problem of perception.
Invention content
In order to solve the above-mentioned technical problem, the purpose of the present invention is to provide a kind of angle elastic parameter Lithology Identification Methods
And device, the problem of to overcome the elastic parameter sensibility that conventional Lithology Identification Methods use insufficient, lithology caused to be difficult to differentiate between.
In order to achieve the above technical purposes, the present invention provides angle elastic parameter Lithology Identification Methods, angle elasticity
Parameter Lithology Identification Methods include the following steps：
On the basis of conventional converted shear wave elastic impedance equation, the converted shear wave with magnitude with longitudinal wave elastic impedance is obtained
Elastic impedance equation；
According to the converted shear wave elastic impedance equation and longitudinal wave elastic impedance equation, converted shear wave elastic impedance is obtained
Pseudo- log and longitudinal wave elastic impedance puppet log；
S-wave impedance log, longitudinal wave elastic impedance puppet well curve are substituted using converted shear wave elastic impedance puppet well curve
P-wave impedance log is substituted, brings common elastic parameter calculation formula into, elastic parameter puppet in a certain angular range is obtained and surveys
Well curve；
Sensibility comparison is carried out using the elastic parameter puppet log combination lithologic log curve, chooses sensibility most
The strong corresponding elastic parameter of angle elastic parameter is target angle elastic parameter, and corresponding angle is best angle；
The conventional more wave angle degree partial stack seismic channel sets of input；
Layer position is carried out to seismic channel set to explain；
Converted shear wave elastic impedance inverting is carried out using conventional constraint Sparse Pulse Inversion algorithm, obtains carrying angle information
Converted shear wave elastic impedance data volume, longitudinal wave elastic impedance data volume；
According to log target angle elastic parameter and best angle analysis result, bring into converted shear wave elastic impedance,
In longitudinal wave elastic impedance inverting data volume to corresponding elastic parameter calculation formula, angle elastic parameter data body is obtained；
According to angle elastic parameter data body lithology value range, Lithology Discrimination is carried out, the angle elastic parameter is completed
The identification of lithology.
In the angle elastic parameter Lithology Identification Methods of the present invention, it is preferable that with longitudinal wave elastic impedance turning with magnitude
Shear wave elastic impedance equation is changed to obtain in accordance with the following steps：
It introduces target zone to be averaged shear wave velocity, density, conventional converted shear wave elastic impedance equation is standardized, is obtained
Converted shear wave elastic impedance equation after standardization：
Wherein：
θ indicates that incidence angle, unit are degree；
VpFor velocity of longitudinal wave, unit is meter per second；
VsFor shear wave velocity, unit is meter per second；
ρ is density, and unit is kilograms per cubic meter；
Vs0For purpose interval shear wave velocity, unit is meter per second；
ρ0For the average value of density, unit is kilograms per cubic meter.
In the angle elastic parameter Lithology Identification Methods of the present invention, it is preferable that conventional more wave angle degree partial stack earthquakes
Trace gather includes longitudinal wave angle part superposition seismic channel set, converted shear wave angle part superposition seismic channel set, wherein it is required that superposition road
It is identical as the best angle to collect angle.
In the angle elastic parameter Lithology Identification Methods of the present invention, it is preferable that converted shear wave elastic impedance data volume SEI
(θ) obtains angle elastic parameter data body with conventional longitudinal wave elastic impedance data volume EI (θ) by combination：Pang(θ)=f
(SEI(θ),EI(θ))。
In the angle elastic parameter Lithology Identification Methods of the present invention, it is preferable that Pang(θ)=f (SEI (θ), EI (θ))
In, the structure relationship of common elastic parameter P is expressed as p-wave impedance IpWith S-wave impedance IsFunction P=f (Ip,Is)。
In the angle elastic parameter Lithology Identification Methods of the present invention, it is preferable that best angle determines in accordance with the following steps：
Using longitudinal wave, shear wave, density log curve, according to converted shear wave elastic impedance equation and conventional longitudinal wave elasticity resistance
Anti- equation, is respectively adopted different angle and calculates angle elastic parameter puppet well curve, and sensibility pair is carried out in conjunction with lithologic log curve
Than it is best angle to choose the corresponding angle of the strongest angle elastic parameter of sensibility.
In the angle elastic parameter Lithology Identification Methods of the present invention, function expression Pang(θ)=f (SEI (θ), EI
(θ)) in functional relation established according to the structure relationship of common elastic parameter.The structure relationship of common elastic parameter P can table
It states as p-wave impedance IpWith S-wave impedance IsFunction P=f (Ip,Is).Angle elastic parameter is sought closing using same function
System, p-wave impedance I is replaced with the longitudinal wave elastic impedance EI (θ) with angle informationp, converted shear wave impedance SEI (θ) substitute shear wave
Impedance Is。
For example, Poisson's ratio parameter more sensitive to sandstone and mud stone in common elastic parameter can be expressed asThen angle Poisson's ratio parameter is
In the angle elastic parameter Lithology Identification Methods of the present invention, obtains elastic parameter puppet in a certain angular range and log well
When curve, those skilled in the art can determine specific angular range as needed.
The present invention also provides a kind of angle elastic parameter Lithology Discrimination device, the angle elastic parameter Lithology Discrimination devices
Including：
Conversion module turns with longitudinal wave elastic impedance with magnitude for being converted to conventional converted shear wave elastic impedance equation
Change shear wave elastic impedance equation；
Log input module, for according to converted shear wave elastic impedance equation proposed by the present invention and longitudinal wave elasticity
Impedance equation obtains converted shear wave elastic impedance puppet log and longitudinal wave elastic impedance puppet log；
Angle elastic parameter puppet log computing module, it is horizontal for being substituted using converted shear wave elastic impedance puppet well curve
Wave impedance log, longitudinal wave elastic impedance puppet well curve substitute p-wave impedance log, bring the calculating of common elastic parameter into
Formula obtains elastic parameter puppet log in a certain angular range；
Best angle module chooses the strongest angle of sensibility for combining lithologic log curve to carry out sensibility comparison
The corresponding elastic parameter of elastic parameter is target angle elastic parameter, and corresponding angle is best angle；
Multi-wave seismic data input module, for inputting conventional more wave angle degree partial stack seismic channel sets；
Explanation module is explained for carrying out layer position to seismic channel set；
Inverting module is obtained for carrying out converted shear wave elastic impedance inverting using conventional constraint Sparse Pulse Inversion algorithm
To converted shear wave elastic impedance data volume, longitudinal wave elastic impedance data volume with angle information；
Angle elastic parameter data module, for being analyzed according to log target angle elastic parameter and best angle
As a result, converted shear wave elastic impedance, longitudinal wave elastic impedance inverting data volume are brought into corresponding elastic parameter calculation formula, obtain
To angle elastic parameter data body；
Identification module, for according to angle elastic parameter data body lithology value range, carrying out Lithology Discrimination.
Specific implementation mode according to the present invention, the log of input need to carry out school using conventional logging alignment technique
Just, it is ensured that all curves can accurately reflect formation information.
In the angle elastic parameter Lithology Discrimination device of the present invention, it is preferable that with longitudinal wave elastic impedance turning with magnitude
Shear wave elastic impedance equation is changed to obtain in accordance with the following steps：
It introduces target zone to be averaged shear wave velocity, density, conventional converted shear wave elastic impedance equation is standardized, is obtained
Converted shear wave elastic impedance equation after standardization：
Wherein：
θ indicates that incidence angle, unit are degree；
VpFor velocity of longitudinal wave, unit is meter per second；
VsFor shear wave velocity, unit is meter per second；
ρ is density, and unit is kilograms per cubic meter；
Vs0For purpose interval shear wave velocity, unit is meter per second；
ρ0For the average value of density, unit is kilograms per cubic meter.
In the angle elastic parameter Lithology Discrimination device of the present invention, it is preferable that conventional more wave angle degree partial stack earthquakes
Trace gather includes longitudinal wave angle part superposition seismic channel set, converted shear wave angle part superposition seismic channel set, wherein it is required that superposition road
It is identical as the best angle to collect angle.
In the angle elastic parameter Lithology Discrimination device of the present invention, it is preferable that converted shear wave elastic impedance data volume SEI
(θ) obtains angle elastic parameter data body with conventional longitudinal wave elastic impedance data volume EI (θ) by combination：Pang(θ)=f
(SEI(θ),EI(θ))。
In the angle elastic parameter Lithology Discrimination device of the present invention, it is preferable that Pang(θ)=f (SEI (θ), EI (θ))
In, the structure relationship of common elastic parameter P is expressed as p-wave impedance IpWith S-wave impedance IsFunction P=f (Ip,Is)。
In the angle elastic parameter Lithology Discrimination device of the present invention, it is preferable that best angle carries out in accordance with the following steps：
Using longitudinal wave, shear wave, density log curve, according to converted shear wave elastic impedance equation and conventional longitudinal wave elasticity resistance
Anti- equation, is respectively adopted different angle and calculates angle elastic parameter puppet well curve, and sensibility pair is carried out in conjunction with lithologic log curve
Than it is best angle to choose the corresponding angle of the strongest angle elastic parameter of sensibility.
Angle elastic parameter Lithology Identification Methods provided by the invention and device, for common elastic parameter by physics itself
The limitation of meaning and lithology separating capacity can not meet the problem of nowadays requirement of fining oil-gas exploration, it is contemplated that earthquake
The angle information and converted shear wave earthquake information of data include abundant lithological information.
Angle elastic parameter Lithology Identification Methods provided by the invention and device, it is rich using multi-wave seismic data lithological information
Rich feature, proposes new converted shear wave elastic impedance equation, from multi-wave seismic elastic impedance inverting, using longitudinal wave, turns
It changes shear wave elastic impedance and builds new Lithology Discrimination parameter, there is higher extraction accuracy and ease for use.
Description of the drawings
Fig. 1 is the angle elastic parameter Lithology Discrimination installation drawing in embodiment；
Fig. 2 is the flow chart of the angle elastic parameter Lithology Identification Methods of embodiment；
Fig. 3 be embodiment ash laccolith at, the velocity of longitudinal wave of Mudstone Reservoirs, shear wave velocity, density curve；
P-wave impedance that Fig. 4 is, S-wave impedance curve；
Conventional 10 °, 20 °, the 30 ° of converted shear wave elastic impedance curves that Fig. 5 is；
10 ° of the present invention that Fig. 6 is, 20 °, 30 ° of converted shear wave elastic impedance curves；
Fig. 7 is embodiment converted shear wave elastic impedance sensibility with incidence angle situation of change；
Fig. 8 is common elastic parameter and angle elastic parameter sensibility correlation curve；
Fig. 9 is certain work area practical logging curve and angle elastic parameter sensibility correlation curve；
Figure 10 is certain practical work area angle elastic parameter inversion section.
Specific implementation mode
In order to which technical characteristic, purpose and the advantageous effect to the present invention are more clearly understood, now to the skill of the present invention
Art scheme carry out it is described further below, but should not be understood as to the present invention can practical range restriction.
Embodiment 1
A kind of angle elastic parameter Lithology Discrimination device is present embodiments provided, the structure of the device is as shown in Figure 1, the dress
It sets and specifically includes：
Conversion module turns with longitudinal wave elastic impedance with magnitude for being converted to conventional converted shear wave elastic impedance equation
Change shear wave elastic impedance equation；
Log input module, for according to converted shear wave elastic impedance equation proposed by the present invention and longitudinal wave elasticity
Impedance equation obtains converted shear wave elastic impedance puppet log and longitudinal wave elastic impedance puppet log；
Angle elastic parameter puppet log computing module, it is horizontal for being substituted using converted shear wave elastic impedance puppet well curve
Wave impedance log, longitudinal wave elastic impedance puppet well curve substitute p-wave impedance log, bring the calculating of common elastic parameter into
Formula obtains elastic parameter puppet log in a certain angular range；
Best angle module chooses the strongest angle of sensibility for combining lithologic log curve to carry out sensibility comparison
The corresponding elastic parameter of elastic parameter is target angle elastic parameter, and corresponding angle is best angle；
Multi-wave seismic data input module, for inputting conventional more wave angle degree partial stack seismic channel sets；
Explanation module is explained for carrying out layer position to seismic channel set；
Inverting module is obtained for carrying out converted shear wave elastic impedance inverting using conventional constraint Sparse Pulse Inversion algorithm
To converted shear wave elastic impedance data volume, longitudinal wave elastic impedance data volume with angle information；
Angle elastic parameter data module, for being analyzed according to log target angle elastic parameter and best angle
As a result, converted shear wave elastic impedance, longitudinal wave elastic impedance inverting data volume are brought into corresponding elastic parameter calculation formula, obtain
To angle elastic parameter data body；
Identification module, for according to angle elastic parameter data body lithology value range, carrying out Lithology Discrimination.
Present embodiments provide a kind of angle elastic parameter Lithology Identification Methods, technological process such as Fig. 2 of the recognition methods
It is shown, specifically include following steps：
Step 101, it on the basis of conventional converted shear wave elastic impedance equation, by mathematical derivation, obtains and longitudinal wave bullet
The new converted shear wave elastic impedance equation of property impedance same order；
The general type of conventional converted shear wave elastic impedance equation is：
SEI (θ)=Vs m(θ)ρn(θ)
Wherein：
Wherein, θ indicates incidence angle, VpFor velocity of longitudinal wave, VsIt is density for shear wave velocity, ρ.
The deficiency of the equation is that the order of magnitude of converted shear wave elastic impedance value is very huge with the variation of incident angle,
The comparison being unfavorable between different angle converted shear wave elastic impedance, the combination being especially detrimental between same longitudinal wave elastic impedance make
With.Fig. 3 show velocity of longitudinal wave, shear wave velocity, the density curve of limestone, the three layer model that dolomite forms, the longitudinal wave of calculating
Impedance, S-wave impedance curve are as shown in Figure 4.Model curve shown according to fig. 3 utilizes conventional converted shear wave elastic impedance equation
Converted shear wave elastic impedance curve when incident angle is 10 °, 20 °, 30 ° is calculated separately, as shown in Figure 5.It can be seen by Fig. 5
Arrive, 10 °, 20 °, 30 ° of the order of magnitude be respectively 10 2 powers, 4 powers, 5 powers, the resistance of incidence angles degree converted shear wave elasticity
Order of magnitude difference is larger and different from the S-wave impedance order of magnitude between anti-, is unfavorable for combination and utilizes.
It is averaged shear wave velocity, density present invention introduces target zone, standard is carried out to conventional converted shear wave elastic impedance equation
Change, the converted shear wave elastic impedance equation after being standardized：
Wherein Vs0、ρ0Interval of interest shear wave velocity and the average value of density are indicated respectively.
Converted shear wave elastic impedance is calculated using three layer model curve shown in Fig. 3 using converted shear wave elastic impedance equation
Curve is as shown in Figure 6.By Fig. 6 it can be seen that, incidence angles degree converted shear wave elastic impedance with S-wave impedance in same number
Magnitude is conducive to combination and utilizes.
In the present embodiment, longitudinal wave elastic impedance equation uses equation used by common elastic Impedance Inversion technology.
Step 102, velocity of longitudinal wave, shear wave velocity, density log curve are inputted.The log of input is needed using conventional
Logging correction technology is corrected, it is ensured that all curves can accurately reflect formation information.
Step 103, the survey inputted according to converted shear wave, longitudinal wave elastic impedance equation and the step 102 described in step 101
Well curve calculates converted shear wave, longitudinal wave elastic impedance puppet log.
Step 104, S-wave impedance log, longitudinal wave elastic impedance are substituted using converted shear wave elastic impedance puppet well curve
Pseudo- well curve substitutes p-wave impedance log, brings common elastic parameter calculation formula into, obtains the bullet in a certain angular range
Property parameter puppet well curve.Sensibility comparison is carried out in conjunction with lithologic log explanation results, chooses the strongest angle elasticity ginseng of sensibility
Number is target angle elastic parameter, and corresponding angle is best angle θ.
The present embodiment is using the relative different of elastic parameter value between target lithology and country rock as lithology sensitivity assessment
Parameter, it is best angle to choose the most strong corresponding incidence angle of sensibility.It is illustrated in figure 7 model dolomite and limestone in Fig. 3
Converted shear wave elastic impedance relative different with incidence angle situation of change.It can be seen that converted shear wave elastic impedance is to the quick of lithology
Perception is not the increase with incident angle and monotonic increase, (is in the present embodiment when incidence angle is more than certain angle
35 °), sensibility reduces instead, therefore is needed in actual mechanical process according to the more wave angle degree elasticity ginsengs of practical logging tracing analysis
The case where number sensibility is with incidence angle, chooses best incident angle.
It is best Lithology Discrimination parameter that the embodiment of the present invention, which chooses the most strong corresponding more wave angle degree elastic parameters of sensibility,.
It is illustrated in figure 8 Fig. 3 model common elastic parameters and angle elastic parameter sensibility comparative situation.It can see dolomite interval
Conventional longitudinal wave speed/shear wave velocity parameter, Poisson's ratio parameter and Lame parameter/shearing amount of touching are respectively less than phase with respect to the difference of limestone
The wave angle degree elastic parameter more than 35 ° answered, and the sensibility of wave angle degree Lame parameter/modulus of shearing is most strong wherein more than 35 °.
The present embodiment is using the relative different of elastic parameter value between different lithology as lithology sensitivity assessment parameter, choosing
It is best angle to take the most strong corresponding incidence angle of sensibility.
Step 105, the present embodiment requires the conventional more wave angle degree partial stack seismic channel sets of input, including longitudinal wave angle part
It is superimposed seismic channel set, converted shear wave angle part is superimposed seismic channel set.It is required that superposition trace gather angle chosen with step 104 it is best
Angle is identical.
Step 106, it carries out layer position to seismic channel set to explain, after fine formation information explanation results can be effectively ensured
The precision of renewed treaty beam Sparse Pulse Inversion.
Step 107, in conjunction with above-mentioned steps, elastic impedance inverting is carried out using conventional constraint Sparse Pulse Inversion algorithm.This
In step other than the converted shear wave elastic impedance equation of use is different from conventional method, other all sames.
Step 108, according to step 107 constrained sparse spike inversion inversion result, the converted shear wave bullet with angle information is obtained
Property impedance and longitudinal wave elastic impedance data volume.
Step 109, according to lithology sensitivity analysis in step 104 as a result, bringing converted shear wave elastic impedance SEI (θ) into, indulging
Wave elastic impedance inverting data volume is to corresponding elastic parameter calculation formula PangIn (θ)=f (SEI (θ), EI (θ)), angle is obtained
Elastic parameter data body Pang(θ)。
Step 110, according to the result of calculation of step 109, according to the angle elastic parameter value range on different lithology stratum,
Carry out Lithology Discrimination.
It is bent that Fig. 9 show somewhere practical logging curve, the routine being calculated and angle Lame parameter/modulus of shearing
Line.By Fig. 9 it can be seen that, in dolostone reservoirs section (in figure shown in box) angle Lame parameter/modulus of shearing (shown in dotted line)
More conventional velocity of longitudinal wave, shear wave velocity, density and Lame parameter/modulus of shearing (solid line) amplitude are larger, with limestone country rock area
Clearly demarcated aobvious, angle Lame parameter/modulus of shearing value of dolostone reservoirs is respectively less than 1.Therefore angle Lame parameter/cut can be utilized
Shear modulu carries out Lithology Discrimination.
Figure 10 is somewhere angle Lame parameter/modulus of shearing value inversion result, it can be seen that inverting section intermediate value is less than 1
Interval and lithology interpretation of well logging result (oval frame show dolostone reservoirs) coincide preferable, therefore Lame can be utilized to join
Number/modulus of shearing value identifies advantageous lithology.
Above example illustrates that angle elastic parameter Lithology Identification Methods of the invention and device have higher extraction essence
Degree and ease for use.
System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity,
Or it is realized by the product with certain function.
For convenience of description, it is divided into various units when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each unit is realized can in the same or multiple software and or hardware when application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
It is realized by the mode of software plus required general hardware platform.Based on this understanding, the technical solution essence of the application
On in other words the part that contributes to existing technology can be expressed in the form of software products, in a typical configuration
In, computing device includes one or more processors (CPU), input/output interface, network interface and memory.The computer is soft
Part product may include that some instructions are used so that a computer equipment (can be personal computer, server or network
Equipment etc.) execute method described in certain parts of each embodiment of the application or embodiment.The computer software product can
To store in memory, memory may include the volatile memory in computer-readable medium, random access memory
(RAM) and/or the forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer
The example of readable medium.Computer-readable medium includes that permanent and non-permanent, removable and non-removable media can be by
Any method or technique realizes information storage.Information can be computer-readable instruction, data structure, the module of program or its
His data.The example of the storage medium of computer includes, but are not limited to phase transition internal memory (PRAM), static RAM
(SRAM), dynamic random access memory (DRAM), other kinds of random access memory (RAM), read-only memory
(ROM), electrically erasable programmable read-only memory (EEPROM), fast flash memory bank or other memory techniques, CD-ROM are read-only
Memory (CD-ROM), digital versatile disc (DVD) or other optical storages, magnetic tape cassette, tape magnetic disk storage or
Other magnetic storage apparatus or any other non-transmission medium can be used for storage and can be accessed by a computing device information.According to
Herein defines, and computer-readable medium does not include of short duration computer readable media (transitory media), such as modulation
Data-signal and carrier wave.
Each embodiment in this specification is described in a progressive manner, identical similar portion between each embodiment
Point just to refer each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so description is fairly simple, related place is referring to embodiment of the method
Part explanation.
The application can be used in numerous general or special purpose computing system environments or configuration.Such as：Personal computer, clothes
Business device computer, handheld device or portable device, laptop device, multicomputer system, microprocessor-based system, set
Top box, programmable consumer-elcetronics devices, network PC, minicomputer, mainframe computer including any of the above system or equipment
Distributed computing environment etc..
The application can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Usually, program module includes routines performing specific tasks or implementing specific abstract data types, program, object, group
Part, data structure etc..The application can also be put into practice in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application there are many deformation and
Variation is without departing from spirit herein, it is desirable to which the attached claims include these deformations and change without departing from the application's
Spirit.
Claims (8)
1. a kind of angle elastic parameter Lithology Identification Methods, which is characterized in that the angle elastic parameter Lithology Identification Methods include
Following steps：
On the basis of conventional converted shear wave elastic impedance equation, the converted shear wave elasticity with magnitude with longitudinal wave elastic impedance is obtained
Impedance equation；
According to the converted shear wave elastic impedance equation and longitudinal wave elastic impedance equation, the survey of converted shear wave elastic impedance puppet is obtained
Well curve and longitudinal wave elastic impedance puppet log；
S-wave impedance log is substituted using converted shear wave elastic impedance puppet well curve, longitudinal wave elastic impedance puppet well curve substitutes
P-wave impedance log substitutes into common elastic parameter calculation formula, and it is bent to obtain elastic parameter puppet well logging in a certain angular range
Line；
Sensibility comparison is carried out using the elastic parameter puppet log combination lithologic log curve, it is strongest to choose sensibility
The corresponding elastic parameter of angle elastic parameter is target angle elastic parameter, and corresponding angle is best angle；
The conventional more wave angle degree partial stack seismic channel sets of input；
Layer position is carried out to seismic channel set to explain；
Converted shear wave elastic impedance inverting is carried out using conventional constraint Sparse Pulse Inversion algorithm, obtains turning with angle information
Change shear wave elastic impedance data volume, longitudinal wave elastic impedance data volume；
According to log target angle elastic parameter and best angle analysis result, converted shear wave elastic impedance, longitudinal wave are substituted into
In elastic impedance inverting data volume to corresponding elastic parameter calculation formula, angle elastic parameter data body is obtained；
According to angle elastic parameter data body lithology value range, Lithology Discrimination is carried out, the angle elastic parameter lithology is completed
Identification；
It is obtained in accordance with the following steps with the converted shear wave elastic impedance equation of magnitude with longitudinal wave elastic impedance：
It introduces target zone to be averaged shear wave velocity, density, conventional converted shear wave elastic impedance equation is standardized, standard is obtained
Converted shear wave elastic impedance equation after change：
Wherein：
θ indicates that incidence angle, unit are degree；
VpFor velocity of longitudinal wave, unit is meter per second；
VsFor shear wave velocity, unit is meter per second；
ρ is density, and unit is kilograms per cubic meter；
Vs0For purpose interval shear wave velocity, unit is meter per second；
ρ0For the average value of density, unit is kilograms per cubic meter.
2. angle elastic parameter Lithology Identification Methods according to claim 1, which is characterized in that the more wave angle degree of routine
Partial stack seismic channel set includes that longitudinal wave angle part is superimposed seismic channel set, converted shear wave angle part superposition seismic channel set,
In, it is desirable that superposition trace gather angle is identical as the best angle.
3. angle elastic parameter Lithology Identification Methods according to claim 1, which is characterized in that converted shear wave elastic impedance
Data volume SEI (θ) obtains angle elastic parameter data body with conventional longitudinal wave elastic impedance data volume EI (θ) by combination：Pang
(θ)=f (SEI (θ), EI (θ)), wherein θ indicates that incidence angle, unit are degree.
4. angle elastic parameter Lithology Identification Methods according to claim 3, which is characterized in that Pang(θ)=f (SEI
(θ), EI (θ)) in, the structure relationship of common elastic parameter P is expressed as p-wave impedance IpWith S-wave impedance IsFunction P=f (Ip,
Is)。
5. angle elastic parameter Lithology Identification Methods according to claim 1, which is characterized in that best angle is according to as follows
Step determines：
Using longitudinal wave, shear wave, density log curve, according to converted shear wave elastic impedance equation and conventional longitudinal wave elastic impedance side
Journey is respectively adopted different angle and calculates angle elastic parameter puppet well curve, and sensibility comparison, choosing are carried out in conjunction with lithologic log curve
It is best angle to take the corresponding angle of the strongest angle elastic parameter of sensibility.
6. a kind of angle elastic parameter Lithology Discrimination device, which is characterized in that the angle elastic parameter Lithology Discrimination device includes：
Conversion module, it is horizontal with the conversion of magnitude with longitudinal wave elastic impedance for being converted to conventional converted shear wave elastic impedance equation
Wave elastic impedance equation；
Log input module, for according to converted shear wave elastic impedance equation and longitudinal wave elastic impedance equation, being turned
Change shear wave elastic impedance puppet log and longitudinal wave elastic impedance puppet log；
Angle elastic parameter puppet log computing module is hindered for substituting shear wave using converted shear wave elastic impedance puppet well curve
Anti- log, longitudinal wave elastic impedance puppet well curve substitute p-wave impedance log, substitute into common elastic parameter calculation formula,
Obtain elastic parameter puppet log in a certain angular range；
Best angle module chooses the strongest angle elasticity of sensibility for combining lithologic log curve to carry out sensibility comparison
The corresponding elastic parameter of parameter is target angle elastic parameter, and corresponding angle is best angle；
Multi-wave seismic data input module, for inputting conventional more wave angle degree partial stack seismic channel sets；
Explanation module is explained for carrying out layer position to seismic channel set；
Inverting module obtains band for carrying out converted shear wave elastic impedance inverting using conventional constraint Sparse Pulse Inversion algorithm
Converted shear wave elastic impedance data volume, the longitudinal wave elastic impedance data volume of angled information；
Angle elastic parameter data module, for according to log target angle elastic parameter and best angle analysis knot
Converted shear wave elastic impedance, longitudinal wave elastic impedance inverting data volume are updated in corresponding elastic parameter calculation formula, obtain by fruit
Angle elastic parameter data body；
Identification module, for according to angle elastic parameter data body lithology value range, carrying out Lithology Discrimination；
It is obtained in accordance with the following steps with the converted shear wave elastic impedance equation of magnitude with longitudinal wave elastic impedance：
It introduces target zone to be averaged shear wave velocity, density, conventional converted shear wave elastic impedance equation is standardized, standard is obtained
Converted shear wave elastic impedance equation after change：
Wherein：
θ indicates that incidence angle, unit are degree；
VpFor velocity of longitudinal wave, unit is meter per second；
VsFor shear wave velocity, unit is meter per second；
ρ is density, and unit is kilograms per cubic meter；
Vs0For purpose interval shear wave velocity, unit is meter per second；
ρ0For the average value of density, unit is kilograms per cubic meter.
7. angle elastic parameter Lithology Discrimination device according to claim 6, which is characterized in that the more wave angle degree of routine
Partial stack seismic channel set includes that longitudinal wave angle part is superimposed seismic channel set, converted shear wave angle part superposition seismic channel set,
In, it is desirable that superposition trace gather angle is identical as the best angle.
8. angle elastic parameter Lithology Discrimination device according to claim 6, which is characterized in that converted shear wave elastic impedance
Data volume SEI (θ) obtains angle elastic parameter data body with conventional longitudinal wave elastic impedance data volume EI (θ) by combination：Pang
(θ)=f (SEI (θ), EI (θ)), wherein θ indicates that incidence angle, unit are degree.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611144943.7A CN106707345B (en) | 2016-12-13 | 2016-12-13 | A kind of angle elastic parameter Lithology Identification Methods and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611144943.7A CN106707345B (en) | 2016-12-13 | 2016-12-13 | A kind of angle elastic parameter Lithology Identification Methods and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106707345A CN106707345A (en) | 2017-05-24 |
| CN106707345B true CN106707345B (en) | 2018-10-16 |
Family
ID=58937317
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611144943.7A Active CN106707345B (en) | 2016-12-13 | 2016-12-13 | A kind of angle elastic parameter Lithology Identification Methods and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106707345B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107576985B (en) * | 2017-07-31 | 2019-05-07 | 中国石油天然气集团公司 | A kind of method and apparatus of seismic inversion |
| CN108089228B (en) * | 2017-12-22 | 2019-12-10 | 中国石油天然气股份有限公司 | Interpretation data method and device for determining stratum rock characteristics |
| CN110895348B (en) * | 2019-12-20 | 2022-04-08 | 岭南师范学院 | Method, system and storage medium for extracting low-frequency information of seismic elastic impedance |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140160890A1 (en) * | 2011-08-05 | 2014-06-12 | Xinding Fang | System and method for determining shear wave anisotropy in a vertically transversely isotropic formation |
| CN102749645B (en) * | 2012-03-14 | 2014-11-05 | 中国石油天然气股份有限公司 | Method and device for utilizing angle impedance gradient to carry out reservoir hydrocarbon detection |
| CN103176211B (en) * | 2013-03-07 | 2016-04-20 | 中国石油集团东方地球物理勘探有限责任公司 | Based on gas-bearing reservoir prediction method and the device of many sensibility elasticities parameter |
| CN104375178B (en) * | 2014-11-05 | 2017-10-17 | 中国石油天然气股份有限公司 | A kind of carbonate rock fractured cave method for predicting reservoir and device |
| CN105785446A (en) * | 2016-03-17 | 2016-07-20 | 成都创源油气技术开发有限公司 | Oil shale earthquake identification and evaluation method |
-
2016
- 2016-12-13 CN CN201611144943.7A patent/CN106707345B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106707345A (en) | 2017-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105675635B (en) | Tight rock component relative amount and brittleness index determine method and apparatus | |
| CN104375188B (en) | Seismic wave transmission attenuation compensation method and device | |
| WO2019062655A1 (en) | Method and device for determining thin interlayer | |
| CN101551466B (en) | Method for improving prediction precision of oil and gas reservoir by using seismic attribute related to offset distance | |
| CN105182408B (en) | The preparation method and device of a kind of synthetic seismogram | |
| CN107765298B (en) | A kind of method and device of determining near-surface velocity model | |
| CN105068117B (en) | For the AVO inversion method of Fractured medium, device and equipment | |
| Lu et al. | Prediction of coal seam details and mining safety using multicomponent seismic data: A case history from China | |
| CN103733089B (en) | For including the system and method that the underground of uncertainty estimation characterizes | |
| CN104316965B (en) | Prediction method and system for fissure azimuth and intensity | |
| CN106707345B (en) | A kind of angle elastic parameter Lithology Identification Methods and device | |
| CN104316966B (en) | A kind of Fluid Identification Method and system | |
| CN103775075A (en) | Method for identifying lithology of whole well section | |
| CN105425289A (en) | Method and device of determining low frequency wave impedance | |
| CN107817520A (en) | The pressure coefficient Forecasting Methodology and system of marine facies mud shale stratum | |
| CN106556861A (en) | A kind of azimuthal AVO inversion method based on Omnibearing earthquake auto data | |
| CN104375178B (en) | A kind of carbonate rock fractured cave method for predicting reservoir and device | |
| CN110471113A (en) | Bearing calibration, device and storage medium are moved in inverting based on unstable state seismic data | |
| CN103576197A (en) | Method for extracting converted wave angle channel set | |
| CN102866426A (en) | Method for analyzing oil gas information of rock body by applying amplitude versus offset (AVO) large-angle trace gathers | |
| Liu et al. | Combined petrophysics and 3D seismic attributes to predict shale reservoirs favourable areas | |
| CN111983683B (en) | Prediction method and system for lake-facies limestone reservoir under low-well condition | |
| CN106054252A (en) | Pre-stack time migration method and device | |
| CN106772599B (en) | A kind of method and device calculating formation shear speed | |
| CN105445793A (en) | Method and apparatus for determining bad trace data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |