CN109459790A - For coal measure strata seismic velocity field method for building up and system - Google Patents

For coal measure strata seismic velocity field method for building up and system Download PDF

Info

Publication number
CN109459790A
CN109459790A CN201710797077.XA CN201710797077A CN109459790A CN 109459790 A CN109459790 A CN 109459790A CN 201710797077 A CN201710797077 A CN 201710797077A CN 109459790 A CN109459790 A CN 109459790A
Authority
CN
China
Prior art keywords
velocity
coal seam
coal
seismic
interval
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
Application number
CN201710797077.XA
Other languages
Chinese (zh)
Other versions
CN109459790B (en
Inventor
时磊
刘俊州
夏红敏
温立峰
王震宇
许云书
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201710797077.XA priority Critical patent/CN109459790B/en
Publication of CN109459790A publication Critical patent/CN109459790A/en
Application granted granted Critical
Publication of CN109459790B publication Critical patent/CN109459790B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters

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

It discloses a kind of for coal measure strata seismic velocity field method for building up and system.This method may include: stack velocity to be converted to interval velocity, and carry out correction process with except abnormality processing, obtain initial optimization interval velocity；Wavelet is extracted to seismic data, and carries out seismic inversion, on the Wave Impedance Data Volume of seismic inversion, the layer position information on prediction coal seam top and coal seam bottom；Velocity of longitudinal wave is filtered, obtains velocity of longitudinal wave and velocity of longitudinal wave under coal seam on coal seam, and then be based on initial optimization interval velocity, obtains interval velocity and interval velocity under coal seam on coal seam；According to coal seam top and the layer position information at coal seam bottom, coal inter-layer information is obtained, and then practical coal seam speed is set as speed between coal seam；By speed merges between interval velocity and coal seam under interval velocity, coal seam on coal seam, final earthquake velocity field is obtained.The present invention is calculated by seismic inversion and layering section, and the available true seismic velocity of coal measure strata ensure that the computational accuracy of angle and the reasonability of AVO inversion result.

Description

For coal measure strata seismic velocity field method for building up and system
Technical field
The present invention relates to oil gas technical field of physical geography, are directed to coal measure strata seismic velocity more particularly, to one kind Field method for building up and system.
Background technique
From Ostrander since nineteen eighty-two proposes AVO technology, numerous scholars conduct in-depth research AVO theory With analysis, it is applied to oil-gas exploration, have passed through development in more than 20 years, AVO technology, which has become, predicts oil using pre-stack data The effective ways of gas, in terms of FRACTURE PREDICTION, pressure prediction, deposit dynamic monitoring, petroleum-gas prediction, reservoir heterogeneity It is widely applied.AVO is to describe amplitude variation with Offset and a technology changing, can be with by AVO inversion technique Various information of description underground medium variation are obtained, the information of lithology and oil gas is disclosed, are had to oil-gas exploration very heavy Want meaning.And seismic data is all to describe as the function of geophone offset, and elastic wave is anti-from processing, reflected amplitude is collected It is then using incidence angle as variable that drill problem reduction, which be approximately the Zoeppritz equation more easily solved and its approximate expression, It is described, by the elastic parameter information of actual seismic data estimation rock.In order to preferably utilize AVO, AVA technological prediction Oil gas needs that offset distance information is switched to angle information using interval velocity, therefore the speed of underground medium is a very crucial Factor, accurate speed are sought, and are always the key problem of seismic prospecting, speed seek it is accurate whether directly affect it is final Predict achievement.
It is a variety of to the method for the interval velocity optimization of earthquake velocity field at present:
Conventional interval velocity acquiring method is sought by Dix formula (nineteen fifty-five), and this method assumed condition is in speed Field is horizontally unchanged, and calculation amount is small, convenient and efficient, is easy to implement, and is suitable for less than 15 ° situations of stratigraphic dip, complex area Error is big, and stack velocity precision itself is also limited；
For tilted stratum, interval velocity precision is improved frequently with modeling tomographic and modeling inversion, wherein model chromatographs Method is asked on the basis of establishing the space geological model of To and velocity field curved surface X, Y-direction using ray communication theory Lead, obtain refraction point the angle of emergence and reflection point deviate incidence point spatial offset, thus seek each layer interval velocity and really Determine reflecting interface；Dix formula is abandoned in the calculating of its interval velocity, solves the problems such as tilted stratum and speed reversing.And model inversion Method is to utilize computer simulation field actual arrangement, it is assumed that underground is two dimension inclination layered medium, by horizontal stack velocity spectrum And the slope of corresponding lineups is finally inversed by underground structure, including interval velocity, interface depth and inclination angle.Basic ideas are: given Initial model synthesizes stack velocity spectrum, just with the difference of practical stack velocity spectrum, modifies model, layer-by-layer iteration is superimposed synthesis Speed and actual speed spectrum error are minimum；
For high steep complex area, the method for inversion that can be concerned with using stack velocity or the relevant method of inversion (Landa of CMP trace gather Et al. 1991), input CMP trace gather, influenced speed is sought on by pre-processing factor it is small, it is theoretical based on ray tracing, it is applicable Range is wide, but needs the road prestack CMP set information, and operand is big.The stack velocity method of inversion is according to the self-excitation at the CMP trace gather Interval velocity is calculated from paper T0 figure and stack velocity is received, each puts given interval velocity range to certain target zone, passes through ray The actual time curve of tracing computation.The time curve calculated by ray tracing time curve corresponding with stack velocity is carried out Compare, then matching the best interval velocity that the corresponding interval velocity of best (closest) time curve is exactly stratum.By above-mentioned side Method node-by-node algorithm speed then can be obtained given time layer along Interval velocity spectrum.And CMP trace gather be concerned with the method for inversion, Method And Principle with Stack velocity is concerned with, and the method for inversion is essentially identical, and implementation method is almost the same with step, but needs to input prestack CMP trace gather, fortune Calculation amount is big, and the period is long.Therefore, it is necessary to develop a kind of for coal measure strata seismic velocity field method for building up and system.
The information for being disclosed in background of invention part is merely intended to deepen the reason to general background technique of the invention Solution, and it is known to those skilled in the art existing to be not construed as recognizing or imply that the information is constituted in any form Technology.
Summary of the invention
The invention proposes one kind to be directed to coal measure strata seismic velocity field method for building up and system, can be by earthquake about Beam Sparse Pulse Inversion and layering section calculate, the available true seismic velocity of coal measure strata, ensure that the calculating essence of angle The reasonability of degree and AVO inversion result.
According to an aspect of the invention, it is proposed that a kind of be directed to coal measure strata seismic velocity field method for building up.The method It may include: stack velocity to be converted into interval velocity, and carry out correction process with except abnormality processing, obtain initial optimization layer speed Degree；Wavelet is extracted to seismic data, and carries out seismic inversion, on the Wave Impedance Data Volume of the seismic inversion, predicts coal seam The layer position information on top and coal seam bottom；Velocity of longitudinal wave is filtered, velocity of longitudinal wave and longitudinal wave speed under coal seam on coal seam are obtained Degree, and then it is based on the initial optimization interval velocity, obtain interval velocity and interval velocity under coal seam on coal seam；According to the coal seam top with The layer position information at coal seam bottom obtains coal inter-layer information, and then practical coal seam speed is set as speed between coal seam；By the coal seam Speed merges between interval velocity and the coal seam under upper interval velocity, the coal seam, obtains final earthquake velocity field.
Preferably, stack velocity is converted to by interval velocity by Dix formula.
Preferably, the Dix formula are as follows:
Wherein, VnIndicate the interval velocity of n-th layer, VR,nIndicate the stack velocity of n-th layer, t0,nIndicate that n-th layer is arrived on ground Time.
Preferably, be filtered to the velocity of longitudinal wave includes: to divide the velocity of longitudinal wave on coal seam, under coal seam It is not filtering to the frequency band of earthquake overlap speed, and then obtains velocity of longitudinal wave and velocity of longitudinal wave under the coal seam on the coal seam.
Preferably, obtaining interval velocity and interval velocity under the coal seam on the coal seam includes: building coal seam overlying strata screen work With Stratigraphic framework under coal seam, well speed interpolation is carried out under corresponding Stratigraphic framework, and then based on the initial optimization interval velocity Processing obtains interval velocity and interval velocity under the coal seam on the coal seam.
Preferably, the correction process includes: singular value correction, low-pass filtering, space interpolation and smooth correction.
According to another aspect of the invention, it is proposed that a kind of establish system for coal measure strata seismic velocity field, can wrap Include: memory is stored with computer executable instructions；Processor, the computer that the processor is run in the memory can It executes instruction, executes following steps: stack velocity is converted into interval velocity, and carry out correction process with except abnormality processing, obtain Initial optimization interval velocity；Wavelet is extracted to seismic data, and carries out seismic inversion, in the Wave Impedance Data Volume of the seismic inversion On, the layer position information on prediction coal seam top and coal seam bottom；Velocity of longitudinal wave is filtered, velocity of longitudinal wave and coal on coal seam are obtained The lower velocity of longitudinal wave of layer, and then it is based on the initial optimization interval velocity, obtain interval velocity and interval velocity under coal seam on coal seam；According to institute Coal seam top and the layer position information at coal seam bottom are stated, obtains coal inter-layer information, and then practical coal seam speed is set as speed between coal seam； By speed merges between interval velocity and the coal seam under interval velocity, the coal seam on the coal seam, final earthquake velocity field is obtained.
Preferably, stack velocity is converted to by interval velocity, the Dix formula by Dix formula are as follows:
Wherein, VnIndicate the interval velocity of n-th layer, VR,nIndicate the stack velocity of n-th layer, t0,nIndicate that n-th layer is arrived on ground Time.
Preferably, be filtered to the velocity of longitudinal wave includes: to divide the velocity of longitudinal wave on coal seam, under coal seam It is not filtering to the frequency band of earthquake overlap speed, and then obtains velocity of longitudinal wave and velocity of longitudinal wave under the coal seam on the coal seam.
Preferably, obtaining interval velocity and interval velocity under the coal seam on the coal seam includes: building coal seam overlying strata screen work With Stratigraphic framework under coal seam, well speed interpolation is carried out under corresponding Stratigraphic framework, and then based on the initial optimization interval velocity Processing obtains interval velocity and interval velocity under the coal seam on the coal seam.
Methods and apparatus of the present invention has other characteristics and advantages, these characteristics and advantages are attached from what is be incorporated herein It will be apparent in figure and subsequent specific embodiment, or will be in the attached drawing being incorporated herein and subsequent specific reality It applies in mode and is stated in detail, the drawings and the detailed description together serve to explain specific principles of the invention.
Detailed description of the invention
Exemplary embodiment of the present is described in more detail in conjunction with the accompanying drawings, of the invention is above-mentioned and other Purpose, feature and advantage will be apparent, wherein in exemplary embodiments of the present invention, identical reference label is usual Represent same parts.
Fig. 1 shows the flow chart of the step of being directed to coal measure strata seismic velocity field method for building up according to the present invention.
Fig. 2 a and Fig. 2 b respectively illustrate the seismic velocity field according to an embodiment of the invention being not optimised and final The schematic diagram of seismic velocity field.
Fig. 3 a and Fig. 3 b respectively illustrate the seismic velocity field conversion according to an embodiment of the invention being not optimised The schematic diagram of incident angle gathers and the incident angle gathers of final earthquake velocity field conversion.
Specific embodiment
The present invention will be described in more detail below with reference to accompanying drawings.Although showing the preferred embodiment of the present invention in attached drawing, However, it is to be appreciated that may be realized in various forms the present invention and should not be limited by the embodiments set forth herein.On the contrary, providing These embodiments are of the invention more thorough and complete in order to make, and can will fully convey the scope of the invention to ability The technical staff in domain.
Fig. 1 shows the flow chart of the step of being directed to coal measure strata seismic velocity field method for building up according to the present invention.
In this embodiment, according to the present invention to may include: for coal measure strata seismic velocity field method for building up
Step 101, stack velocity is converted into interval velocity, and carries out correction process with except abnormality processing, obtained initial excellent Change interval velocity；In one example, stack velocity is converted to by interval velocity by Dix formula.In one example, correction process It include: singular value correction, low-pass filtering, space interpolation and smooth correction.
In one example, Dix formula are as follows:
Wherein, VnIndicate the interval velocity of n-th layer, VR,nIndicate the stack velocity of n-th layer, t0,nIndicate that n-th layer is arrived on ground Time.
Step 102, wavelet is extracted to seismic data, and carries out seismic inversion, on the Wave Impedance Data Volume of seismic inversion, Predict the layer position information on coal seam top and coal seam bottom.
Step 103, velocity of longitudinal wave is filtered, obtains velocity of longitudinal wave and velocity of longitudinal wave under coal seam on coal seam, into And it is based on initial optimization interval velocity, obtain interval velocity and interval velocity under coal seam on coal seam；In one example, to velocity of longitudinal wave into Row filtering processing includes: velocity of longitudinal wave to be filtering to the frequency band of earthquake overlap speed respectively on coal seam, under coal seam, and then obtain Obtain velocity of longitudinal wave and velocity of longitudinal wave under coal seam on coal seam.In one example, interval velocity and interval velocity under coal seam on coal seam are obtained Include: Stratigraphic framework under building coal seam overlying strata screen work and coal seam, under corresponding Stratigraphic framework, and then is based on initial optimization layer Speed carries out the processing of well speed interpolation, obtains interval velocity and interval velocity under coal seam on coal seam.
Step 104, according to the layer position information on coal seam top and coal seam bottom, coal inter-layer information is obtained, and then practical coal seam is fast Degree is set as speed between coal seam.
Step 105, by speed merges between interval velocity and coal seam under interval velocity, coal seam on coal seam, final seismic velocity is obtained ?.
Specifically, since the velocity anomaly in coal seam is low, the velocity field precision that seismic processing is established is inadequate, can not be apparent The interval velocity of coal output layer is portrayed, therefore the angle changed by offset distance is caused false information occur, inversion result is directly affected, is This, can depict the top in coal seam and bottom to come in conjunction with constrained sparse spike inversion inverting, then be layered section (on coal seam under coal seam and Coal seam) Lai Jinhang interval velocity body constraint well interpolation optimal speed, fully consider coal measure strata interval velocity precision problem.
Stack velocity is converted into interval velocity using Dix formula according to formula (1), and carries out singular value correction, low pass filtered The correction process such as wave, space interpolation and smooth correction remove outlier and velocity anomaly, obtain initial optimization interval velocity；In conjunction with survey Well information extracts wavelet to seismic data, seismic constraint Sparse Pulse Inversion is completed, in the Wave Impedance Data Volume of seismic inversion On, the coal seam of low wave impedance part is carved out, coal output layer spatial distribution is portrayed, and predicts the layer position on coal seam top and coal seam bottom Information；Velocity of longitudinal wave is filtering to the frequency band of earthquake overlap speed respectively on coal seam, under coal seam, constructs coal seam overlying strata lattice Stratigraphic framework under frame and coal seam carries out Jing Suduchazhichu under corresponding Stratigraphic framework, and then based on initial optimization interval velocity Reason, and then obtain velocity of longitudinal wave and velocity of longitudinal wave under coal seam on coal seam；According to the layer position information on coal seam top and coal seam bottom, coal is obtained Inter-layer information, and then practical coal seam speed is set as speed between coal seam；By interval velocity and coal seam under interval velocity, coal seam on coal seam Between speed merge, obtain final earthquake velocity field.
This method is calculated by seismic constraint Sparse Pulse Inversion and layering section, the available true earthquake of coal measure strata Speed ensure that the computational accuracy of angle and the reasonability of AVO inversion result.
Using example
A concrete application example is given below in the scheme and its effect of the embodiment of the present invention for ease of understanding.This field It should be understood to the one skilled in the art that the example is only for the purposes of understanding the present invention, any detail is not intended to be limited in any way The system present invention.
Stack velocity is converted into interval velocity using Dix formula according to formula (1), and carries out singular value correction, low pass filtered The correction process such as wave, space interpolation and smooth correction remove outlier and velocity anomaly, obtain initial optimization interval velocity；In conjunction with survey Well information extracts wavelet to seismic data, seismic constraint Sparse Pulse Inversion is completed, in the Wave Impedance Data Volume of seismic inversion On, the coal seam of low wave impedance part is carved out, coal output layer spatial distribution is portrayed, and predicts the layer position on coal seam top and coal seam bottom Information；Velocity of longitudinal wave is filtering to the frequency band of earthquake overlap speed respectively on coal seam, under coal seam, as 0-2Hz constructs coal Stratigraphic framework under layer overlying strata screen work and coal seam carries out well under corresponding Stratigraphic framework, and then based on initial optimization interval velocity Speed interpolation processing, and then obtain velocity of longitudinal wave and velocity of longitudinal wave under coal seam on coal seam；According to the layer position on coal seam top and coal seam bottom Information obtains coal inter-layer information, and then practical coal seam speed 2420m/s is set as speed between coal seam；By interval velocity on coal seam, Speed merges between interval velocity and coal seam under coal seam, obtains final earthquake velocity field.
Fig. 2 a and Fig. 2 b respectively illustrate the seismic velocity field according to an embodiment of the invention being not optimised and final The schematic diagram of seismic velocity field.
Fig. 3 a and Fig. 3 b respectively illustrate the seismic velocity field conversion according to an embodiment of the invention being not optimised The schematic diagram of incident angle gathers and the incident angle gathers of final earthquake velocity field conversion.
Applied to the DaNiuDi gas field of Ordos Basin, the gas field stratum is more gentle, generally develops coal seam, and there is speed Spend Inversion Problem, from well velocity of longitudinal wave and seismic interval velocity compare it is both visible there are certain trend differences, especially In the place of the coal seam 1.5s or so development, logging speed and seismic velocity field difference are big, as shown in Figure 2 a, the quarter that can not be apparent Draw the speed in coal seam.After the method for the present invention, it can be very good solve speed Inversion Problem, while what is be apparent depicts The interval velocity in coal seam guarantees to be consistent with well under coal seam with the trend under coal seam, as shown in Figure 2 b, and ensure that survey Well speed and the good consistency in seismic velocity field.
The trace gather conversion of AVO is carried out using velocity field after optimization of the invention, the seismic velocity field conversion being not optimised enters Firing angle trace gather is at the position of development position H1 layers of coal seam and under H1, because seismic interval velocity is much larger than coal seam true velocity, after conversion Incident angle it is small, as shown in Figure 3a；The incident angle gathers of final earthquake velocity field conversion H1 layers of position and its under, reduce The incidence angle of coal seam true velocity, conversion is more, and angle is big, as shown in Figure 3b, conducive to the AVO attributive analysis in later period and prestack inversion Research.
In conclusion the present invention is calculated by seismic constraint Sparse Pulse Inversion and layering section, available coal measure strata True seismic velocity ensure that the computational accuracy of angle and the reasonability of AVO inversion result.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
System is established for coal measure strata seismic velocity field according to an embodiment of the invention, providing one kind, may include: Memory is stored with computer executable instructions；Processor, the computer executable instructions in processor run memory, holds Row following steps: being converted to interval velocity for stack velocity, and carries out correction process with except abnormality processing, obtains initial optimization layer speed Degree；Wavelet is extracted to seismic data, and carries out seismic inversion, on the Wave Impedance Data Volume of seismic inversion, prediction coal seam top with The layer position information at coal seam bottom；Velocity of longitudinal wave is filtered, velocity of longitudinal wave and velocity of longitudinal wave under coal seam on coal seam are obtained, into And it is based on initial optimization interval velocity, obtain interval velocity and interval velocity under coal seam on coal seam；According to the layer position on coal seam top and coal seam bottom Information obtains coal inter-layer information, and then practical coal seam speed is set as speed between coal seam；It will be under interval velocity, coal seam on coal seam Speed merges between interval velocity and coal seam, obtains final earthquake velocity field.
In one example, stack velocity is converted to by interval velocity, Dix formula by Dix formula are as follows:
Wherein, VnIndicate the interval velocity of n-th layer, VR,nIndicate the stack velocity of n-th layer, t0,nIndicate that n-th layer is arrived on ground Time.
In one example, be filtered to velocity of longitudinal wave includes: to distinguish velocity of longitudinal wave on coal seam, under coal seam It is filtering to the frequency band of earthquake overlap speed, and then obtains velocity of longitudinal wave and velocity of longitudinal wave under coal seam on coal seam.
In one example, obtain on coal seam interval velocity under interval velocity and coal seam include: building coal seam overlying strata screen work with Stratigraphic framework under coal seam carries out the processing of well speed interpolation under corresponding Stratigraphic framework, and then based on initial optimization interval velocity, obtains Obtain interval velocity and interval velocity under coal seam on coal seam.
The present invention is calculated by seismic constraint Sparse Pulse Inversion and layering section, the available true earthquake of coal measure strata Speed ensure that the computational accuracy of angle and the reasonability of AVO inversion result.
It will be understood by those skilled in the art that above to the purpose of the description of the embodiment of the present invention only for illustratively saying The beneficial effect of bright the embodiment of the present invention is not intended to limit embodiments of the invention to given any example.
Various embodiments of the present invention are described above, above description is exemplary, and non-exclusive, and It is not limited to disclosed each embodiment.Without departing from the scope and spirit of illustrated each embodiment, for this skill Many modifications and changes are obvious for the those of ordinary skill in art field.

Claims (10)

1. one kind is directed to coal measure strata seismic velocity field method for building up, comprising:
Stack velocity is converted into interval velocity, and carries out correction process with except abnormality processing, obtains initial optimization interval velocity；
Wavelet is extracted to seismic data, and carries out seismic inversion, on the Wave Impedance Data Volume of the seismic inversion, predicts coal seam The layer position information on top and coal seam bottom；
Velocity of longitudinal wave is filtered, obtains velocity of longitudinal wave and velocity of longitudinal wave under coal seam on coal seam, and then based on described first Begin optimization interval velocity, obtains interval velocity and interval velocity under coal seam on coal seam；
According to the layer position information on coal seam top and coal seam bottom, coal inter-layer information is obtained, and then practical coal seam speed is set as Speed between coal seam；
By speed merges between interval velocity and the coal seam under interval velocity, the coal seam on the coal seam, final seismic velocity is obtained ?.
2. according to claim 1 be directed to coal measure strata seismic velocity field method for building up, wherein will be folded by Dix formula Acceleration is converted to interval velocity.
3. according to claim 2 be directed to coal measure strata seismic velocity field method for building up, wherein the Dix formula are as follows:
Wherein, VnIndicate the interval velocity of n-th layer, VR,nIndicate the stack velocity of n-th layer, t0,nTime of the expression ground to n-th layer.
4. it is according to claim 1 be directed to coal measure strata seismic velocity field method for building up, wherein to the velocity of longitudinal wave into Row is filtered
The velocity of longitudinal wave is filtering to the frequency band of earthquake overlap speed respectively on coal seam, under coal seam, and then described in acquisition Velocity of longitudinal wave and velocity of longitudinal wave under the coal seam on coal seam.
5. according to claim 1 be directed to coal measure strata seismic velocity field method for building up, wherein obtain the coal seam upper layer Speed includes: with interval velocity under the coal seam
Stratigraphic framework under coal seam overlying strata screen work and coal seam is constructed, under corresponding Stratigraphic framework, and then based on described initial excellent Change interval velocity and carry out the processing of well speed interpolation, obtains interval velocity and interval velocity under the coal seam on the coal seam.
6. according to claim 1 be directed to coal measure strata seismic velocity field method for building up, wherein the correction process packet It includes: singular value correction, low-pass filtering, space interpolation and smooth correction.
7. one kind establishes system for coal measure strata seismic velocity field, which is characterized in that the system includes:
Memory is stored with computer executable instructions；
Processor, the processor run the computer executable instructions in the memory, execute following steps:
Stack velocity is converted into interval velocity, and carries out correction process with except abnormality processing, obtains initial optimization interval velocity；
Wavelet is extracted to seismic data, and carries out seismic inversion, on the Wave Impedance Data Volume of the seismic inversion, predicts coal seam The layer position information on top and coal seam bottom；
Velocity of longitudinal wave is filtered, obtains velocity of longitudinal wave and velocity of longitudinal wave under coal seam on coal seam, and then based on described first Begin optimization interval velocity, obtains interval velocity and interval velocity under coal seam on coal seam；
According to the layer position information on coal seam top and coal seam bottom, coal inter-layer information is obtained, and then practical coal seam speed is set as Speed between coal seam；
By speed merges between interval velocity and the coal seam under interval velocity, the coal seam on the coal seam, final seismic velocity is obtained ?.
8. according to claim 7 establish system for coal measure strata seismic velocity field, wherein will be folded by Dix formula Acceleration is converted to interval velocity, the Dix formula are as follows:
Wherein, VnIndicate the interval velocity of n-th layer, VR,nIndicate the stack velocity of n-th layer, t0,nTime of the expression ground to n-th layer.
9. according to claim 7 establish system for coal measure strata seismic velocity field, wherein to the velocity of longitudinal wave into Row is filtered
The velocity of longitudinal wave is filtering to the frequency band of earthquake overlap speed respectively on coal seam, under coal seam, and then described in acquisition Velocity of longitudinal wave and velocity of longitudinal wave under the coal seam on coal seam.
10. according to claim 7 establish system for coal measure strata seismic velocity field, wherein obtain on the coal seam Interval velocity includes: with interval velocity under the coal seam
Stratigraphic framework under coal seam overlying strata screen work and coal seam is constructed, under corresponding Stratigraphic framework, and then based on described initial excellent Change interval velocity and carry out the processing of well speed interpolation, obtains interval velocity and interval velocity under the coal seam on the coal seam.
CN201710797077.XA 2017-09-06 2017-09-06 Method and system for establishing seismic velocity field of coal measure strata Active CN109459790B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710797077.XA CN109459790B (en) 2017-09-06 2017-09-06 Method and system for establishing seismic velocity field of coal measure strata

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710797077.XA CN109459790B (en) 2017-09-06 2017-09-06 Method and system for establishing seismic velocity field of coal measure strata

Publications (2)

Publication Number Publication Date
CN109459790A true CN109459790A (en) 2019-03-12
CN109459790B CN109459790B (en) 2021-02-09

Family

ID=65606069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710797077.XA Active CN109459790B (en) 2017-09-06 2017-09-06 Method and system for establishing seismic velocity field of coal measure strata

Country Status (1)

Country Link
CN (1) CN109459790B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110471112A (en) * 2019-09-20 2019-11-19 西南石油大学 Dipping bed interval velocity exception inversion method based on stack velocity variation
CN112147700A (en) * 2019-06-28 2020-12-29 中国石油化工股份有限公司 Low-frequency model construction method and system for speed abnormal area

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130028049A1 (en) * 2011-07-25 2013-01-31 Naide Pan Method for handling rough sea and irregular recording conditions in multi-sensor towed streamer data
CN102918422A (en) * 2010-03-31 2013-02-06 斯塔特伊石油公司 Estimating interval velocities
CN104459782A (en) * 2014-12-10 2015-03-25 中国石油天然气股份有限公司 Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit
US20160061973A1 (en) * 2014-08-29 2016-03-03 Pgs Geophysical As Methods and systems that determine a velocity wavefield from a measured pressure wavefield
CN106405643A (en) * 2016-08-30 2017-02-15 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Well-to-seismic integration velocity field construction method and device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102918422A (en) * 2010-03-31 2013-02-06 斯塔特伊石油公司 Estimating interval velocities
US20130028049A1 (en) * 2011-07-25 2013-01-31 Naide Pan Method for handling rough sea and irregular recording conditions in multi-sensor towed streamer data
US20160061973A1 (en) * 2014-08-29 2016-03-03 Pgs Geophysical As Methods and systems that determine a velocity wavefield from a measured pressure wavefield
CN104459782A (en) * 2014-12-10 2015-03-25 中国石油天然气股份有限公司 Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit
CN106405643A (en) * 2016-08-30 2017-02-15 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Well-to-seismic integration velocity field construction method and device

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
孙铁军，等: "DIVA技术在Tenere地区的应用", 《石油地球物理勘探》 *
张水山，等: "基于综合信息的三维层速度场建立及其在涪陵页岩气田应用研究", 《石油物探》 *
易远元，等: "三维地震速度场建立技术-以饶阳凹陷河间南地区为例", 《石油学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112147700A (en) * 2019-06-28 2020-12-29 中国石油化工股份有限公司 Low-frequency model construction method and system for speed abnormal area
CN110471112A (en) * 2019-09-20 2019-11-19 西南石油大学 Dipping bed interval velocity exception inversion method based on stack velocity variation

Also Published As

Publication number Publication date
CN109459790B (en) 2021-02-09

Similar Documents

Publication Publication Date Title
CN107817535B (en) The determination method and apparatus of short lap
WO2018010628A1 (en) Seismic rock physics inversion method based on a large area tight reservoir
US20150066460A1 (en) Stratigraphic function
CN109738945B (en) Method for directly generating construction diagram by using prestack depth migration result
EA020635B1 (en) Method for seismic survey for hydrocarbons using average velocity of wave velocity field construction
WO2017035104A1 (en) Velocity model seismic static correction
CN103439740B (en) The method predicted based on the relative impedances of dipole seismic wavelet multiple integral and device
US20140336940A1 (en) Estimation of q-factor in time domain
CN107045143A (en) A kind of method and device of Prediction of fracture
WO2016135505A1 (en) A method of producing and utilising high resolution impedance logs derived from vsp data for use in assessing an oilfield subterranean formation
CN111007567A (en) Sand shale thin interbed prediction method and system based on seismic waveform inversion
CN105549082A (en) Establishing method and system of three-dimensional geomechanical field of extra-deep carbonate reservoir
CN107942379A (en) A kind of method for improving complex fault block rate pattern precision
CN111751885A (en) Shale gas volume fracturing micro-seismic monitoring method
CN104975851B (en) For amplitude with the reservoir model optimization method of geophone offset variation road set analysis
CN109459790A (en) For coal measure strata seismic velocity field method for building up and system
CN111983683B (en) Prediction method and system for lake-facies limestone reservoir under low-well condition
CN111722284B (en) Method for establishing speed depth model based on gather data
US20120099396A1 (en) System and method for characterization with non-unique solutions of anisotropic velocities
CN106908838A (en) The method for building target area stratum inelastic attenuation quality factor three-dimensional model
CN109143335B (en) Method, system, medium and apparatus for producing synthetic seismic record
CN106597547A (en) Method for accurately describing earthquake in thin reservoir
US20220236435A1 (en) Low-Frequency Seismic Survey Design
CN109839659B (en) Method for carrying out iterative optimization on prestack depth migration profile
CN103424774B (en) A kind of earthquake hydrocarbon-bearing pool formation method and device

Legal Events

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