CN106707345A - Method and device for identifying lithology of angle elastic parameter - Google Patents
Method and device for identifying lithology of angle elastic parameter Download PDFInfo
- Publication number
- CN106707345A CN106707345A CN201611144943.7A CN201611144943A CN106707345A CN 106707345 A CN106707345 A CN 106707345A CN 201611144943 A CN201611144943 A CN 201611144943A CN 106707345 A CN106707345 A CN 106707345A
- Authority
- CN
- China
- Prior art keywords
- angle
- elastic
- wave
- impedance
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 238000006243 chemical reaction Methods 0.000 claims abstract 2
- 238000001615 p wave Methods 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000010008 shearing Methods 0.000 description 8
- 238000003860 storage Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 5
- 235000019738 Limestone Nutrition 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 239000006028 limestone Substances 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 239000010459 dolomite Substances 0.000 description 3
- 229910000514 dolomite Inorganic materials 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 241000287196 Asthenes Species 0.000 description 1
- 241001274660 Modulus Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000010297 mechanical methods and process Methods 0.000 description 1
- 230000005226 mechanical processes and functions Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 239000004575 stone Substances 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/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 invention provides a method and a device for identifying lithology of an angle elastic parameter. The method comprises the following steps: obtaining a converted transverse wave elastic impedance equation with the same magnitude as the longitudinal wave elastic impedance; obtaining a converted transverse wave elastic impedance pseudo logging curve and a longitudinal wave elastic impedance pseudo logging curve; obtaining angle elastic parameter pseudo logging curves; obtaining an optimal angle; inputting a conventional multi-wave angle partial stack seismic gather; performing horizon interpretation on the seismic channel set; obtaining a converted transverse wave elastic impedance data body and a longitudinal wave elastic impedance data body with angle information; obtaining an angle elastic parameter data body; and according to the lithological dereferencing range of the angle elastic parameter data, performing lithological identification to finish the identification of the lithological. The device includes: a conversion module; a logging curve input module; an angle elastic parameter pseudo logging curve calculation module; an optimal angle module; a multi-wave seismic data input module; an interpretation module; an inversion module; an angle elastic parameter data volume module; and identifying the 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 dress
Put, belong 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 influenceed by construction and reservoir heterogeneity, and accumulating condition is complicated, and identification difficulty is big, and investment risk is high.Therefore,
Develop a set of high-precision Lithology Discrimination geophysics solution and seem 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 road collection or angle part superposition road collection, using maximal possibility estimation, constrains sparse
Pulse isoinversion algorithm, is finally inversed by the elastic parameter body comprising lithological information.
Due to the limitation of elastic parameter physical significance in itself, its sensitiveness to lithology has certain limit, and lithology is distinguished
Ability can not meet the requirement of the oil-gas exploration that nowadays becomes more meticulous and conventional method does not make full use of the angle of geological data
Converted shear wave earthquake information in information and multi-wave seismic.In fact, the angle information of geological data and converted shear wave ground
Shake information includes abundant lithological information.Therefore, current elastic parameter Lithology Identification Methods generally existing is caused quick
The not strong problem of perception.
The content of the invention
In order to solve the above-mentioned technical problem, it is an object of the invention to provide a kind of angle elastic parameter Lithology Identification Methods
And device, it is not enough with the elastic parameter sensitiveness for overcoming conventional Lithology Identification Methods to use, cause the problem that lithology is difficult to differentiate between.
In order to realize above-mentioned technical purpose, the invention provides angle elastic parameter Lithology Identification Methods, angle elasticity
Parameter Lithology Identification Methods are comprised the following steps：
On the basis of conventional converted shear wave elastic impedance equation, the converted shear wave with magnitude with compressional wave elastic impedance is obtained
Elastic impedance equation；
According to the converted shear wave elastic impedance equation and compressional wave elastic impedance equation, converted shear wave elastic impedance is obtained
Pseudo- log and compressional wave elastic impedance puppet log；
S-wave impedance log, compressional wave elastic impedance puppet well curve are substituted using converted shear wave elastic impedance puppet well curve
P-wave impedance log is substituted, common elastic parameter calculation formula is brought into, the pseudo- survey of elastic parameter in a certain angular range is obtained
Well curve；
Sensitiveness contrast is carried out using elastic parameter puppet log combination lithologic log curve, sensitiveness is chosen most
The strong corresponding elastic parameter of angle elastic parameter is angle on target elastic parameter, and corresponding angle is best angle；
The conventional many 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, is obtained with angle information
Converted shear wave elastic impedance data volume, compressional wave elastic impedance data volume；
According to log angle on target elastic parameter and best angle analysis result, bring into converted shear wave elastic impedance,
Compressional wave elastic impedance inverting data volume obtains angle elastic parameter data body in corresponding elastic parameter computing formula；
According to angle elastic parameter data body lithology span, Lithology Discrimination is carried out, complete the angle elastic parameter
The identification of lithology.
In angle elastic parameter Lithology Identification Methods of the invention, it is preferable that with compressional wave elastic impedance turning with magnitude
Shear wave elastic impedance equation is changed to obtain in accordance with the following steps：
The average shear wave velocity of target zone, density are introduced, conventional converted shear wave elastic impedance equation is standardized, obtained
Converted shear wave elastic impedance equation after standardization：
Wherein：
θ represents incidence angle, and unit is degree；
VpIt is velocity of longitudinal wave, unit is meter per second；
VsIt is shear wave velocity, unit is meter per second；
ρ is density, unit for kilograms per cubic meter；
Vs0It is purpose interval shear wave velocity, unit is meter per second；
ρ0Be the average value of density, unit for kilograms per cubic meter.
In angle elastic parameter Lithology Identification Methods of the invention, it is preferable that conventional many wave angle degree partial stack earthquakes
Road collection includes compressional wave angle part superposition seismic channel set, converted shear wave angle part superposition seismic channel set, wherein, it is desirable to it is superimposed
Collection angle is identical with the best angle.
In angle elastic parameter Lithology Identification Methods of the invention, it is preferable that converted shear wave elastic impedance data volume SEI
(θ), with conventional compressional wave elastic impedance data volume EI (θ), angle elastic parameter data body is obtained by combination：Pang(θ)=f
(SEI(θ),EI(θ))。
In angle elastic parameter Lithology Identification Methods of the invention, it is preferable that Pang(θ)=f (SEI (θ), EI (θ))
In, the structure relation of common elastic parameter P is expressed as p-wave impedance IpWith S-wave impedance IsFunction P=f (Ip,Is)。
In angle elastic parameter Lithology Identification Methods of the invention, it is preferable that best angle determines in accordance with the following steps：
Using compressional wave, shear wave, density log curve, according to converted shear wave elastic impedance equation and conventional compressional wave elasticity resistance
Anti- equation, is respectively adopted different angle calculation angle elastic parameter puppet well curves, and sensitiveness pair is carried out with reference to lithologic log curve
Than it is best angle to choose the most strong corresponding angle of angle elastic parameter of sensitiveness.
In angle elastic parameter Lithology Identification Methods of the invention, function expression Pang(θ)=f (SEI (θ), EI
(θ)) in functional relation set up according to the structure relation of common elastic parameter.The structure relation of common elastic parameter P can table
It is p-wave impedance I to statepWith S-wave impedance IsFunction P=f (Ip,Is).The asking for of angle elastic parameter is closed using same function
System, p-wave impedance I is replaced with compressional 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 angle elastic parameter Lithology Identification Methods of the invention, the pseudo- well logging of elastic parameter in a certain angular range is obtained
During curve, those skilled in the art can as needed determine specific angular range.
Present invention also offers a kind of angle elastic parameter Lithology Discrimination device, the angle elastic parameter Lithology Discrimination device
Including：
Modular converter, for conventional converted shear wave elastic impedance equation to be converted into the turning with magnitude with compressional wave elastic impedance
Change shear wave elastic impedance equation；
Log input module, for according to converted shear wave elastic impedance equation proposed by the present invention and compressional wave elasticity
Impedance equation, obtains converted shear wave elastic impedance puppet log and compressional 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, compressional 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, for carrying out sensitiveness contrast with reference to lithologic log curve, chooses the most strong angle of sensitiveness
The corresponding elastic parameter of elastic parameter is angle on target elastic parameter, and corresponding angle is best angle；
Multi-wave seismic data input module, for being input into conventional many wave angle degree partial stack seismic channel sets；
Explanation module, explains for carrying out layer position to seismic channel set；
Inverting module, for carrying out converted shear wave elastic impedance inverting using conventional constraint Sparse Pulse Inversion algorithm, obtains
To converted shear wave elastic impedance data volume, compressional wave elastic impedance data volume with angle information；
Angle elastic parameter data module, for being analyzed according to log angle on target elastic parameter and best angle
As a result, converted shear wave elastic impedance, compressional wave elastic impedance inverting data volume are brought into corresponding elastic parameter computing formula, are obtained
To angle elastic parameter data body；
Identification module, for according to angle elastic parameter data body lithology span, carrying out Lithology Discrimination.
Specific embodiment of the invention, the log of input needs to carry out school using conventional logging alignment technique
Just, it is ensured that all curves can exactly reflect formation information.
In angle elastic parameter Lithology Discrimination device of the invention, it is preferable that with compressional wave elastic impedance turning with magnitude
Shear wave elastic impedance equation is changed to obtain in accordance with the following steps：
The average shear wave velocity of target zone, density are introduced, conventional converted shear wave elastic impedance equation is standardized, obtained
Converted shear wave elastic impedance equation after standardization：
Wherein：
θ represents incidence angle, and unit is degree；
VpIt is velocity of longitudinal wave, unit is meter per second；
VsIt is shear wave velocity, unit is meter per second；
ρ is density, unit for kilograms per cubic meter；
Vs0It is purpose interval shear wave velocity, unit is meter per second；
ρ0Be the average value of density, unit for kilograms per cubic meter.
In angle elastic parameter Lithology Discrimination device of the invention, it is preferable that conventional many wave angle degree partial stack earthquakes
Road collection includes compressional wave angle part superposition seismic channel set, converted shear wave angle part superposition seismic channel set, wherein, it is desirable to it is superimposed
Collection angle is identical with the best angle.
In angle elastic parameter Lithology Discrimination device of the invention, it is preferable that converted shear wave elastic impedance data volume SEI
(θ), with conventional compressional wave elastic impedance data volume EI (θ), angle elastic parameter data body is obtained by combination：Pang(θ)=f
(SEI(θ),EI(θ))。
In angle elastic parameter Lithology Discrimination device of the invention, it is preferable that Pang(θ)=f (SEI (θ), EI (θ))
In, the structure relation of common elastic parameter P is expressed as p-wave impedance IpWith S-wave impedance IsFunction P=f (Ip,Is)。
In angle elastic parameter Lithology Discrimination device of the invention, it is preferable that best angle is carried out in accordance with the following steps：
Using compressional wave, shear wave, density log curve, according to converted shear wave elastic impedance equation and conventional compressional wave elasticity resistance
Anti- equation, is respectively adopted different angle calculation angle elastic parameter puppet well curves, and sensitiveness pair is carried out with reference to lithologic log curve
Than it is best angle to choose the most strong corresponding angle of angle elastic parameter of sensitiveness.
Angle elastic parameter Lithology Identification Methods and device that the present invention is provided, physics itself is received for common elastic parameter
The limitation of meaning and lithology separating capacity can not meet the problem of the requirement of the oil-gas exploration that nowadays becomes more meticulous, 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 and device that the present invention is provided, it is rich using multi-wave seismic data lithological information
Rich the characteristics of, new converted shear wave elastic impedance equation is proposed, from multi-wave seismic elastic impedance inverting, using compressional wave, turned
Change shear wave elastic impedance and build new Lithology Discrimination parameter, with extraction accuracy and ease for use higher.
Brief 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 into, the velocity of longitudinal wave of Mudstone Reservoirs, shear wave velocity, density curve；
Fig. 4 is p-wave impedance, the S-wave impedance curve for obtaining；
Fig. 5 is conventional 10 °, 20 °, the 30 ° of converted shear wave elastic impedance curves for obtaining；
Fig. 6 is 10 °, 20 °, the 30 ° of converted shear wave elastic impedance curves of the invention for obtaining；
Fig. 7 is embodiment converted shear wave elastic impedance sensitiveness with incidence angle situation of change；
Fig. 8 is common elastic parameter and angle elastic parameter sensitiveness correlation curve；
Fig. 9 is certain work area practical logging curve and angle elastic parameter sensitiveness correlation curve；
Figure 10 is certain actual work area angle elastic parameter inversion section.
Specific embodiment
In order to be more clearly understood to technical characteristic of the invention, purpose and beneficial effect, now to skill of the invention
Art scheme carry out it is described further below, but it is not intended that to it is of the 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
Put and specifically include：
Modular converter, for conventional converted shear wave elastic impedance equation to be converted into the turning with magnitude with compressional wave elastic impedance
Change shear wave elastic impedance equation；
Log input module, for according to converted shear wave elastic impedance equation proposed by the present invention and compressional wave elasticity
Impedance equation, obtains converted shear wave elastic impedance puppet log and compressional 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, compressional 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, for carrying out sensitiveness contrast with reference to lithologic log curve, chooses the most strong angle of sensitiveness
The corresponding elastic parameter of elastic parameter is angle on target elastic parameter, and corresponding angle is best angle；
Multi-wave seismic data input module, for being input into conventional many wave angle degree partial stack seismic channel sets；
Explanation module, explains for carrying out layer position to seismic channel set；
Inverting module, for carrying out converted shear wave elastic impedance inverting using conventional constraint Sparse Pulse Inversion algorithm, obtains
To converted shear wave elastic impedance data volume, compressional wave elastic impedance data volume with angle information；
Angle elastic parameter data module, for being analyzed according to log angle on target elastic parameter and best angle
As a result, converted shear wave elastic impedance, compressional wave elastic impedance inverting data volume are brought into corresponding elastic parameter computing formula, are obtained
To angle elastic parameter data body；
Identification module, for according to angle elastic parameter data body lithology span, 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, on the basis of conventional converted shear wave elastic impedance equation, by mathematical derivation, obtains and compressional 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, θ represents incidence angle, VpIt is velocity of longitudinal wave, VsFor shear wave velocity, ρ are density.
The deficiency of the equation be the order of magnitude of converted shear wave elastic impedance value as the change of incident angle is very huge,
It is unfavorable for the comparing between different angular transition shear wave elastic impedances, the combination being especially detrimental between same compressional wave elastic impedance makes
With.Fig. 3 show limestone, velocity of longitudinal wave, shear wave velocity, the density curve of the three layer model of dolomite composition, the compressional wave of calculating
Impedance, S-wave impedance curve are as shown in Figure 4.The model curve according to Fig. 3, using conventional converted shear wave elastic impedance equation
Calculate respectively incident angle be 10 °, 20 °, 30 ° when converted shear wave elastic impedance curve, as shown in Figure 5.Can be seen by Fig. 5
Arrive, 10 °, 20 °, 30 ° of the order of magnitude be respectively 10 2 powers, 4 powers, 5 powers, incidence angles degree converted shear wave elasticity resistance
Order of magnitude difference is larger and different from the S-wave impedance order of magnitude between anti-, is unfavorable for that combination is utilized.
Present invention introduces the average shear wave velocity of target zone, density, 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、ρ0The average value of interval of interest shear wave velocity and density is represented respectively.
Using converted shear wave elastic impedance equation, using three layer model curve shown in Fig. 3, converted shear wave elastic impedance is calculated
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, utilizes beneficial to combination.
In the present embodiment, the equation that compressional wave elastic impedance equation is used using common elastic Impedance Inversion technology.
Step 102, input velocity of longitudinal wave, shear wave velocity, density log curve.The log of input is needed using conventional
Logging correction technology is corrected, it is ensured that all curves can exactly reflect formation information.
Step 103, the survey of converted shear wave, compressional wave elastic impedance equation and step 102 input according to step 101
Well curve calculates converted shear wave, compressional wave elastic impedance puppet log.
Step 104, S-wave impedance log, compressional 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.Sensitiveness contrast is carried out with reference to lithologic log explanation results, the most strong angle elasticity ginseng of sensitiveness is chosen
Number is angle on target elastic parameter, and corresponding angle is best angle θ.
The present embodiment is by the use of 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 sensitiveness.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 monotonic increase with the increase of incident angle, (is in the present embodiment when incidence angle is more than certain angle
35 °), sensitiveness is reduced on the contrary, therefore is needed in actual mechanical process according to many wave angle degree elasticity ginsengs of practical logging tracing analysis
Number sensitiveness chooses optimal incident angle with the situation of incidence angle.
It is optimal Lithology Discrimination parameter that the embodiment of the present invention chooses the most strong corresponding many wave angle degree elastic parameters of sensitiveness.
It is illustrated in figure 8 Fig. 3 model common elastic parameters and contrasts situation with angle elastic parameter sensitiveness.It can be seen that dolomite interval
The difference respectively less than phase of conventional compressional wave speed/shear wave velocity parameter, Poisson's ratio parameter and Lame parameter/shearing amount of touching with respect to limestone
The wave angle degree elastic parameter more than 35 ° answered, and the sensitiveness of wave angle degree Lame parameter/modulus of shearing is most strong wherein more than 35 °.
The present embodiment is selected by the use of the relative different of elastic parameter value between different lithology as lithology sensitivity assessment parameter
It is best angle to take the most strong corresponding incidence angle of sensitiveness.
Step 105, the conventional many wave angle degree partial stack seismic channel sets of the present embodiment requirement input, including compressional wave angle part
Superposition seismic channel set, converted shear wave angle part superposition seismic channel set.It is required that superposition road collection angle is optimal with what step 104 was chosen
Angle is identical.
Step 106, layer position is carried out to seismic channel set and is explained, after fine formation information explanation results can be effectively ensured
The precision of renewed treaty beam Sparse Pulse Inversion.
Step 107, with reference to above-mentioned steps, elastic impedance inverting is carried out using conventional constraint Sparse Pulse Inversion algorithm.This
In step in addition to the converted shear wave elastic impedance equation for using is different from conventional method, other all sames.
Step 108, according to step 107 constrained sparse spike inversion inversion result, obtains the converted shear wave bullet with angle information
Property impedance and compressional wave elastic impedance data volume.
Step 109, according to lithology sensitivity analysis result in step 104, brings converted shear wave elastic impedance SEI (θ) into, indulges
Ripple elastic impedance inverting data volume is to corresponding elastic parameter computing 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 span on different lithology stratum,
Carry out Lithology Discrimination.
It is bent that Fig. 9 show somewhere practical logging curve, the routine that is 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 square frame) 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, the angle Lame parameter/modulus of shearing value of dolostone reservoirs is respectively less than 1.Therefore can be using angle Lame parameter/cut
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 can be joined using Lame
Number/modulus of shearing value recognizes favourable lithology.
Above example illustrates that angle elastic parameter Lithology Identification Methods of the invention and device have extraction essence higher
Degree and ease for use.
System, device, module or unit that above-described embodiment is illustrated, can specifically be realized by computer chip or entity,
Or realized by the product with certain function.
For convenience of description, it is divided into various units with function during description apparatus above to describe respectively.Certainly, this is being implemented
The function of each unit can be realized in same or multiple softwares and/or hardware during application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
Realized by the mode of software plus required general hardware platform.Based on such understanding, the technical scheme essence of the application
On the part that is contributed to prior art in other words can be embodied in the form of software product, in a typical configuration
In, computing device includes one or more processors (CPU), input/output interface, network interface and internal memory.The computer is soft
Part product can be used to including some instructions so that a computer equipment (can be personal computer, server, or network
Equipment etc.) perform method described in some parts of each embodiment of the application or embodiment.The computer software product can
To store in internal memory, internal memory potentially includes the volatile memory in computer-readable medium, random access memory
And/or the form, such as read-only storage (ROM) or flash memory (flash RAM) such as Nonvolatile memory (RAM).Internal memory is computer
The example of computer-readable recording medium.Computer-readable medium includes that permanent and non-permanent, removable and non-removable media can be by
Any method or technique realizes information Store.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 storage
(ROM), Electrically Erasable Read Only Memory (EEPROM), fast flash memory bank or other memory techniques, read-only optical disc are read-only
Memory (CD-ROM), digital versatile disc (DVD) or other optical storages, magnetic cassette tape, tape magnetic rigid disk storage or
Other magnetic storage apparatus or any other non-transmission medium, can be used to store the information that can be accessed by a computing device.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 by the way of progressive, identical similar portion between each embodiment
Divide mutually referring to what each embodiment was stressed is the difference with other embodiment.Especially for system reality
Apply for example, because it is substantially similar to embodiment of the method, so description is fairly simple, related part 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.For example：Personal computer, clothes
Business device computer, handheld device or portable set, laptop device, multicomputer system, the system based on microprocessor, put
Top box, programmable consumer-elcetronics devices, network PC, minicom, mainframe computer, including any of the above system or equipment
DCE etc..
The application can be described in the general context of computer executable instructions, such as program
Module.Usually, program module includes performing particular task or realizes routine, program, object, the group of particular abstract data type
Part, data structure etc..The application can also be in a distributed computing environment put into practice, in these DCEs, by
Remote processing devices connected by communication network perform task.In a distributed computing environment, program module can be with
In local and remote computer-readable storage medium including including storage device.
Although depicting the application by embodiment, it will be appreciated by the skilled addressee that the application have it is many deformation and
Change is without deviating from spirit herein, it is desirable to which appended claim includes these deformations and changes without deviating from the application's
Spirit.
Claims (10)
1. a kind of angle elastic parameter Lithology Identification Methods, it is characterised 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 compressional wave elastic impedance is obtained
Impedance equation；
According to the converted shear wave elastic impedance equation and compressional wave elastic impedance equation, the pseudo- survey of converted shear wave elastic impedance is obtained
Well curve and compressional wave elastic impedance puppet log；
S-wave impedance log, compressional wave elastic impedance puppet well curve are substituted using converted shear wave elastic impedance puppet well curve to substitute
P-wave impedance log, brings common elastic parameter calculation formula into, obtains the pseudo- well logging of elastic parameter in a certain angular range bent
Line；
Sensitiveness contrast is carried out using elastic parameter puppet log combination lithologic log curve, sensitiveness is chosen most strong
The corresponding elastic parameter of angle elastic parameter is angle on target elastic parameter, and corresponding angle is best angle；
The conventional many 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, turning with angle information is obtained
Change shear wave elastic impedance data volume, compressional wave elastic impedance data volume；
According to log angle on target elastic parameter and best angle analysis result, converted shear wave elastic impedance, compressional wave are brought into
Elastic impedance inverting data volume obtains angle elastic parameter data body in corresponding elastic parameter computing formula；
According to angle elastic parameter data body lithology span, Lithology Discrimination is carried out, complete the angle elastic parameter lithology
Identification.
2. angle elastic parameter Lithology Identification Methods according to claim 1, it is characterised in that same with compressional wave elastic impedance
The converted shear wave elastic impedance equation of magnitude is obtained in accordance with the following steps：
The average shear wave velocity of target zone, density are introduced, conventional converted shear wave elastic impedance equation is standardized, obtain standard
Converted shear wave elastic impedance equation after change：
Wherein：
θ represents incidence angle, and unit is degree；
VpIt is velocity of longitudinal wave, unit is meter per second；
VsIt is shear wave velocity, unit is meter per second；
ρ is density, unit for kilograms per cubic meter；
Vs0It is purpose interval shear wave velocity, unit is meter per second；
ρ0Be the average value of density, unit for kilograms per cubic meter.
3. angle elastic parameter Lithology Identification Methods according to claim 1, it is characterised in that many wave angle degree of routine
Partial stack seismic channel set includes compressional wave angle part superposition seismic channel set, converted shear wave angle part superposition seismic channel set, its
In, it is desirable to it is superimposed road collection angle identical with the best angle.
4. angle elastic parameter Lithology Identification Methods according to claim 1, it is characterised in that converted shear wave elastic impedance
Data volume SEI (θ), with conventional compressional wave elastic impedance data volume EI (θ), angle elastic parameter data body is obtained by combination：Pang
(θ)=f (SEI (θ), EI (θ)).
5. angle elastic parameter Lithology Identification Methods according to claim 4, it is characterised in that Pang(θ)=f (SEI
(θ), EI (θ)) in, the structure relation of common elastic parameter P is expressed as p-wave impedance IpWith S-wave impedance IsFunction P=f (Ip,
Is)。
6. angle elastic parameter Lithology Identification Methods according to claim 1, it is characterised in that best angle is according to as follows
Step determines：
Using compressional wave, shear wave, density log curve, according to converted shear wave elastic impedance equation and conventional compressional wave elastic impedance side
Journey, is respectively adopted different angle calculation angle elastic parameter puppet well curves, and sensitiveness contrast is carried out with reference to lithologic log curve, selects
The most strong corresponding angle of angle elastic parameter of sensitiveness is taken for best angle.
7. a kind of angle elastic parameter Lithology Discrimination device, it is characterised in that the angle elastic parameter Lithology Discrimination device includes：
Modular converter, it is horizontal with the conversion of magnitude with compressional wave elastic impedance for conventional converted shear wave elastic impedance equation to be converted to
Ripple elastic impedance equation；
Log input module, for according to converted shear wave elastic impedance equation and compressional wave elastic impedance equation, being turned
Change shear wave elastic impedance puppet log and compressional wave elastic impedance puppet log；
Angle elastic parameter puppet log computing module, for substituting shear wave resistance using converted shear wave elastic impedance puppet well curve
Anti- log, compressional wave elastic impedance puppet well curve substitute p-wave impedance log, bring common elastic parameter calculation formula into,
Obtain elastic parameter puppet log in a certain angular range；
Best angle module, for carrying out sensitiveness contrast with reference to lithologic log curve, chooses the most strong angle elasticity of sensitiveness
The corresponding elastic parameter of parameter is angle on target elastic parameter, and corresponding angle is best angle；
Multi-wave seismic data input module, for being input into conventional many wave angle degree partial stack seismic channel sets；
Explanation module, explains for carrying out layer position to seismic channel set；
Inverting module, for carrying out converted shear wave elastic impedance inverting using conventional constraint Sparse Pulse Inversion algorithm, obtains band
The converted shear wave elastic impedance data volume of angled information, compressional wave elastic impedance data volume；
Angle elastic parameter data module, for according to log angle on target elastic parameter and best angle analysis knot
Really, converted shear wave elastic impedance, compressional wave elastic impedance inverting data volume are brought into corresponding elastic parameter computing formula, are obtained
Angle elastic parameter data body；
Identification module, for according to angle elastic parameter data body lithology span, carrying out Lithology Discrimination.
8. angle elastic parameter Lithology Discrimination device according to claim 7, it is characterised in that same with compressional wave elastic impedance
The converted shear wave elastic impedance equation of magnitude is obtained in accordance with the following steps：
The average shear wave velocity of target zone, density are introduced, conventional converted shear wave elastic impedance equation is standardized, obtain standard
Converted shear wave elastic impedance equation after change：
Wherein：
θ represents incidence angle, and unit is degree；
VpIt is velocity of longitudinal wave, unit is meter per second；
VsIt is shear wave velocity, unit is meter per second；
ρ is density, unit for kilograms per cubic meter；
Vs0It is purpose interval shear wave velocity, unit is meter per second；
ρ0Be the average value of density, unit for kilograms per cubic meter.
9. angle elastic parameter Lithology Discrimination device according to claim 7, it is characterised in that many wave angle degree of routine
Partial stack seismic channel set includes compressional wave angle part superposition seismic channel set, converted shear wave angle part superposition seismic channel set, its
In, it is desirable to it is superimposed road collection angle identical with the best angle.
10. angle elastic parameter Lithology Discrimination device according to claim 7, it is characterised in that converted shear wave elasticity resistance
Anti- data volume SEI (θ), with conventional compressional wave elastic impedance data volume EI (θ), angle elastic parameter data body is obtained by combination：
Pang(θ)=f (SEI (θ), EI (θ)).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611144943.7A CN106707345B (en) | 2016-12-13 | 2016-12-13 | Method and device for identifying lithology of angle elastic parameter |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611144943.7A CN106707345B (en) | 2016-12-13 | 2016-12-13 | Method and device for identifying lithology of angle elastic parameter |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106707345A true CN106707345A (en) | 2017-05-24 |
| CN106707345B 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 | Method and device for identifying lithology of angle elastic parameter |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106707345B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107576985A (en) * | 2017-07-31 | 2018-01-12 | 中国石油天然气集团公司 | A kind of method and apparatus of seismic inversion |
| CN108089228A (en) * | 2017-12-22 | 2018-05-29 | 中国石油天然气股份有限公司 | Interpretation data method and device for determining stratum rock characteristics |
| CN110895348A (en) * | 2019-12-20 | 2020-03-20 | 岭南师范学院 | A method, system and storage medium for extracting low frequency information of seismic elastic impedance |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102749645A (en) * | 2012-03-14 | 2012-10-24 | 中国石油天然气股份有限公司 | Method and device for detecting reservoir hydrocarbons by using angle impedance gradient |
| WO2013022421A1 (en) * | 2011-08-05 | 2013-02-14 | Halliburton Energy Services, Inc | System and method for determining shear wave anisotropy in a vertically transversely isotropic formation |
| CN103176211A (en) * | 2013-03-07 | 2013-06-26 | 中国石油集团东方地球物理勘探有限责任公司 | Method and device, for predicting gas reservoir, based on multi-sensitive elastic parameters |
| CN104375178A (en) * | 2014-11-05 | 2015-02-25 | 中国石油天然气股份有限公司 | Carbonate rock fracture-cave reservoir prediction method 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
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013022421A1 (en) * | 2011-08-05 | 2013-02-14 | Halliburton Energy Services, Inc | System and method for determining shear wave anisotropy in a vertically transversely isotropic formation |
| CN102749645A (en) * | 2012-03-14 | 2012-10-24 | 中国石油天然气股份有限公司 | Method and device for detecting reservoir hydrocarbons by using angle impedance gradient |
| CN103176211A (en) * | 2013-03-07 | 2013-06-26 | 中国石油集团东方地球物理勘探有限责任公司 | Method and device, for predicting gas reservoir, based on multi-sensitive elastic parameters |
| CN104375178A (en) * | 2014-11-05 | 2015-02-25 | 中国石油天然气股份有限公司 | Carbonate rock fracture-cave reservoir prediction method and device |
| CN105785446A (en) * | 2016-03-17 | 2016-07-20 | 成都创源油气技术开发有限公司 | Oil shale earthquake identification and evaluation method |
Non-Patent Citations (1)
| Title |
|---|
| 张树林: "基于纵波地震和四分量地震的弹性波阻抗反演", 《海洋地质与第四纪地质》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107576985A (en) * | 2017-07-31 | 2018-01-12 | 中国石油天然气集团公司 | A kind of method and apparatus of seismic inversion |
| CN108089228A (en) * | 2017-12-22 | 2018-05-29 | 中国石油天然气股份有限公司 | Interpretation data method and device for determining stratum rock characteristics |
| CN110895348A (en) * | 2019-12-20 | 2020-03-20 | 岭南师范学院 | A method, system and storage medium for extracting low frequency information of seismic elastic impedance |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106707345B (en) | 2018-10-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| 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 | |
| CN104375188B (en) | Seismic wave transmission attenuation compensation method and device | |
| CN105675635B (en) | Method and device for determining relative content of components of compact rock and brittleness index | |
| 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 | |
| CN110501744A (en) | Hydrocarbon source rock organic carbon geophysics quantitative forecasting technique, device, equipment and storage medium | |
| US11221425B1 (en) | Generating a model for seismic velocities in a subsurface region using inversion with lateral variations | |
| CN104316965A (en) | Prediction method and system for fissure azimuth and intensity | |
| CN104965224A (en) | PP wave and PS wave joint AVO inversion method using average angle of incidence gather | |
| CN102866426B (en) | A kind of method utilizing AVO wide-angle road set analysis rock mass hydrocarbon information | |
| Pech et al. | Quartic moveout coefficient: 3D description and application to tilted TI media | |
| CN106707345B (en) | Method and device for identifying lithology of angle elastic parameter | |
| CN104375178B (en) | Carbonate rock fracture-cave reservoir prediction method and device | |
| CN105445793B (en) | A kind of method and device for determining bad track data | |
| CN112684498A (en) | Reservoir fracture prediction method and system based on wide-azimuth seismic data | |
| CN111812716B (en) | Pre-stack quantitative prediction method, device and equipment for total organic carbon content of shale gas reservoir | |
| CN105527648B (en) | Computational methods and system for the sensitivity matrix of Anisotropic parameters inversion | |
| US11531129B2 (en) | Picking seismic stacking velocity based on structures in a subterranean formation | |
| CN110007349B (en) | Elastic parameter inversion method | |
| CN107991705B (en) | Well logging curve correction method and device based on two-dimensional statistical characteristics | |
| CN106845086A (en) | formation pressure calculation method and device | |
| Copley | Seismic characterization of Niobrara fluid and rock properties: A 4D study and multicomponent (3C) analysis | |
| CN107329169A (en) | Method and device for extracting angle gather | |
| CN107422375A (en) | The determination method and apparatus of the CRP gather of subsea node |
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 |