519 results about "Wave equation" patented technology

The instant invention is a method for numerically propagating waves or solving wave equations on a digital computer. This invention can be used to compute highly accurate solutions to the wave equation, in fact in some cases it computes the analytical solution, something previously considered impossible. The instant method can also propagate waves that are not described by differential equations, such as anisotropic scalar waves. The invention has the advantages that it is computationally efficient, accurate, and flexible. Of importance is the ability to propagate waves that simulate the P-wave arrivals in both isotropic and anisotropic media with a scalar as opposed to a vector equation.

Method for interpolating aliased seismic data to generate shot records at regular shot points. The method involves sorting the data into common offsets to obtain the finer sampling in the common-offset domain, then using DMO to transform the data to zero offset to eliminate azimuth and offset variations, from which densely sampled shot records are reconstructed by an efficient Log-Stretch Inverse Shot Record DMO. The method has ready application to techniques that benefit from input data that are regularly and densely sampled, such as 3-D surface related multiple elimination and shot record wave equation migration.

A method for estimating near-surface material properties in the vicinity of a locally dense group of seismic receivers is disclosed. The method includes receiving seismic data that has been measured by a locally dense group of seismic receivers. Local derivatives of the wavefield are estimated such that the derivatives are centered at a single location preferably using the Lax-Wendroff correction. Physical relationships between the estimated derivatives including the free surface condition and wave equations are used to estimate near-surface material properties in the vicinity of the receiver group. Another embodiment of the invention is disclosed wherein the physical relationships used to estimate material properties are derived from the physics of plane waves arriving at the receiver group, and the group of receivers does not include any buried receivers.

The invention discloses a method and a device for full-wave-shape inversion based on a time domain-step velocity-stress elastic wave equation, wherein the method includes: adopting an adjoint method based on the perturbation theory, and confirming the adjoint equation of the time domain-step velocity-stress elastic wave equation and a corresponding gradient expression, in regard to model parameters, of an objective functional; conforming a forward-propagating wave field according to the time domain-step velocity-stress elastic wave equation, and conforming an inverse time extrapolation adjoint wave field according to the adjoint equation; confirming the gradient, in regard to the model parameters, of the objective functional according to the forward-propagating wave field, the inverse time extrapolation adjoint wave field and the gradient expression; and using a gradient iterative algorithm to carry out multi-scale full-wave-shape inversion. According to the method for the full-wave-shape inversion, before transmitted in an inverse-time mode, a wave field residual error does not need pretreatment, the conversion problems of particle velocity and particle displacement do not need to be considered, and inversion calculation efficiency can be improved.

The invention discloses an elastic reverse time migration imaging method by combining seismic multi-component and belongs to the field of exploration geophysics. The method comprises the following steps: directly inputting seismic multi-component data; without performing the wave field separation on the input data, extrapolating forwards at inverse time based on an elastic wave equation; building an underground elastic vector seismic wave field by combining with multi-component; acquiring four elastic wave imaging results having specific physical significance, such as longitudinal wave, transverse wave, transition longitudinal wave and transition transverse wave, by applying an elastic wave cross-correlation imaging condition of lighting compensation; and suppressing the low wave number noise in the elastic reverse time migration by performing low-pass angle filtering, thereby acquiring an imaging result of final elastic wave migration. By directly inputting the seismic multi-component data and building the underground elastic vector seismic wave field by combining with multi-component, the method can be used for accurately imaging for complex earth medium.

Method for determining reservoir permeability from Stoneley wave attenuation extracted from conventional sonic logs by inversion of the full Biot wave equations for a porous medium. Frequency-dependent Stoneley-wave attenuation is extracted by analyzing array sonic measurements. Then, based on Biot's full theory applied to a borehole model and the standard logs (gamma ray, caliper, density, neutron, resistivity, sonic, etc.), a simulation model with the same parameters as the Stoneley-wave measurements is built. Next, a theoretical Stoneley-wave attenuation is computed for a given permeability. Finally, reservoir permeability is determined by comparing the modeled Stoneley-wave attenuation with the measured Stoneley-wave attenuation by an iterative inversion process.

The invention relates to a method and a system for oil pumping control of an oil pumping unit with a sucker rod pump. In the method, a ground indicating diagram is generated by collecting suspension center load and displacement data of the oil pumping unit; an underground pump power diagram is obtained by calculation according to a wave equation; according to the graphic characteristic of the pump power diagram, fuzzy logic is applied to carry out analysis and control to a traveling valve leakage loss and other pump working conditions; if serious fault occurs, the machine is stopped for detection; otherwise, operating parameters are adjusted according to the full degree of the pump obtained by analysis; if an indirect pumping condition is satisfied, the pumping stop time is calculated or adjusted for carrying out indirect pumping control; otherwise, a new stroke parameter is calculated according to the full degree of the pump. An ARM embedded type oil pumping controller adopting the method can combine the pump power diagram diagnosis with the real-time oil pumping control, realize fault protection, pumping stop, indirect pumping or stroke parameter adjustment and remarkably improve working efficiency.

A system and method for shallow water simulation may provide a framework for solving General Shallow Wave Equations (GSWE) to efficiently simulate 3D fluid effects on arbitrary surfaces using a height field representation. The height field representation may include height columns constructed along surface normals, which may be dependent on a condition of boundary cells adjacent to fluid cells and / or artificial viscosity effects. The framework may provide implicit schemes for solving for the effects of external forces applied to the fluid, including gravity and surface tension, and explicit schemes for solving for advection effects. The system and method may be implemented on general-purpose CPU(s) and / or GPU(s) and may be capable of simulating a variety of fluid effects including: waves, rivulets and streams, drops, and capillary events. In some embodiments, the system and method may achieve real-time fluid control and fluid shape design through user-interaction (e.g., in a graphical painting application).

An instrumentation system for assessing operation of a reciprocating pump system which produces hydrocarbons from a non-vertical or a vertical wellbore. The instrumentation system periodically produces a downhole pump card as a function of a directly or indirectly measured surface card and a friction law function from a wave equation which describes the linear vibrations in a long slender rod. A control signal or command signal is generated based on characteristics of the downhole pump card for controlling the pumping system. It also generates a pump and well analysis report that is useful for a pump operation and determination of its condition.

The present invention relates to a method of processing data representing energy propagating through a medium (e.g., acoustic, elastic or electromagnetic energy) and describes an efficient and flexible e approach to forward modeling and inversion of such energy for a given medium. The representation theorem for the wave-equation is used, in combination with time-reversal invariance and reciprocity, to express the Green's function between two points in the interior of the model as an integral over the response in those points due to sources regularly distributed on a surface surrounding the medium and the points.

A staggered grid 3D seismic prestack reverse time migration imaging method based on GPU small storage capacity is a very important technology for 3D prestack depth migration imaging of complex structures. The reverse time migration method based on the two-way wave equation has been realized, which can accurately image the rotary wave and the large-angle reflected wave, and at the same time solve the problem of the large amount of calculation of this method. Utilizing parallel computing technology based on GPU hardware, staggered grid difference format, checkpoint storage and other technologies, and relying on Laplace imaging conditions, the pre-stack reverse time migration imaging of 3D seismic data is realized, which greatly shortens the computing time. The efficiency has been increased hundreds of times, and the imaging accuracy has been significantly improved.

Solving wave equation with damp calculates rate of change from swabbing parameters including stroke, number of strokes, diameter of pump, depth of pump, diameter and length of sucker rod etc. to system effectiveness. Sorting rate of change is carried out according to magnitude of absolute value. Under limiting conditions of well yield, load, torque moment etc, analyzing degree of sensitivity of swabbing parameters, the invention optimizes a group of optimal swabbing parameters for oil pumping well system. The method can calculate out polished rod horsepower, efficiencies of sucker rod, and oil pipe difficult to calculate out before. Carrying out quantitative analysis for each factor affecting system effectiveness, the invention solves issue of designing complex parameters so as to raise system effectiveness. The invention is applicable to oil pumping well system composed of inclined well or straight well, different typed pumping units and electric motor.

A migration method using hybrid one-way and full-way (HOF) wave equation propagation. The HOF method extrapolates seismic wavefields in less complex media with a one-way wave equation propagator and extrapolates seismic wavefields in extremely complex media with a full-way wave equation propagator. For prestack depth migration, the HOF extrapolates source-side and receiver-side wavefields independently. Frequency-space domain and time-space domain imaging conditions are applied to the one-way and the full-way extrapolated wavefields, respectively. A suitable amplitude matching factor is introduced to combine the one-way and full-way images. The HOF method is a cost-effective migration that produces superior image quality with less noises and less computational resources.

The invention discloses ultrasound dynamic elastic imaging probe and method. The probe internally comprises an imaging transducer with higher frequency and exciting transducers with lower frequency, wherein the imaging transducer is arranged in the middle, and the exciting transducers are arranged at both sides of the imaging transducer; the imaging transducer and the exciting transducers synchronously work in a coupling mode; ultrasonic radiometric force generated by the exciting transducers generates shear waves in a tissue, and the transverse propagation of the shear waves causes the longitudinal displacement of the tissue; and the imaging transducer transmits ultrasonic waves to detect the longitudinal displacement of the tissue and then reconstructs an elastic physical quantity according to a wave equation so as to acquire an elastic image of the tissue. In the invention, the transducers with lower frequency are used for exciting the tissue to generate the shear waves, and the generated shear waves have strong intensity and wider effective exciting range; and meanwhile, the traditional transducer with higher frequency is used for imaging the shear waves, and the elastic image has high signal to noise ratio and high image quality and conforms to the standard of national acoustic power.

According to a preferred aspect of the instant invention, there is provided herein a system and method for extending zero-offset or stacked wave-equation illumination analysis into the angle-gather domain, where it becomes an appropriate tool for assessing the effects of complex overburden on AVA response. A preferred method for doing this involves first creating an angle gather that has a perfect AVA response (i.e. a constant amplitude as a function of angle). This gather is then preferably used as a reflectivity map that is fed into a demigration process which creates modeled data that by construction carries with it a completely flat reflectivity signature. Remigration of such a data set then results in a gather on which any amplitude variation is more likely to be a measure of illumination effects alone. The resulting AVA signature on the gather can then be used to assess the validity of the AVA response on modeled or actual data, resulting in a useful AVA risk analysis.

A method of three-dimensional preserved-amplitude pre-stack time migration which is applied to reflected seismic data treatment in seismic prospecting is a pre-stack migration imaging method aiming at the three-dimensional collected seismic data. The technical scheme comprises three aspects: 1) performing the migration amplitude calculation by using a weight coefficient acquired according to a one-way wave equation and the reverse convolution imaging condition of the depth migration, 2) confirming the time varying three-dimensional migration aperture according to the biggest incidence formed underground, and 3) automatically suppressing the migration voice by performing strip attenuation method on the edge of the migration aperture. The method of the invention realizes high-efficient, high SNR (Signal to Noise Ratio) and preserved-amplitude migration imaging by improving the migration aperture, providing a weight coefficient of the preserved-amplitude and suppressing the migration voice in the imaging process. The preserved-amplitude reflective spot collection generated by using the method of the invention can better serve for the oil gas and fluid detection technique for pre-stack inversion. The method of the invention has high application value to oil gas and mine resource exploration.

The invention discloses a method for inverting a near-surface velocity model by utilizing preliminary waveforms. The method comprises acoustic wave equation-based wave field forward modeling and steepest descent-based waveform inversion technologies, and comprises the following steps of 1, extracting time-domain preliminary waveform records and an initial model; 2, calculating a simulated wave field and a wave field residual by utilizing acoustic wave equation staggered grid finite-difference forward modeling simulation; 3, reversely propagating the wave field residual to obtain a retransmission wave field; 4, calculating a gradient of a target function by utilizing the retransmission wave field and a forward propagation wave field, and calculating an updating step length; 5, updating a speed model; 6, inspecting whether the speed model is consistent with an iteration stopping condition, outputting the speed model if the speed model is consistent with the iteration stopping condition, otherwise returning to the step 2, and continuing iterative updating. According to the method, a wave equation theory-based full-waveform inversion technology is used as reference, and preliminary waves with higher energy and more stable waveforms are used for inversion, so that the multiplicity of solutions of full-waveform inversion is reduced, and the inversion stability and the calculation efficiency are improved; the accuracy of static correction and shallow depth imaging is improved.

The optimized design method of geological geophysical model based three-dimensional seismological observation system includes the following steps: A). creating geological geophysical model; B). designing two kinds of several kinds of observation systems with three-dimentional geological data acquisition; C). adopting imaginary spectrum method and combinatino of three-dimensional wave equation numerical analog and three-dimensional geophysical model analog two methods of implement data acquisitino work; D). using common-reflection point as basis, making depth deviation treatment of acquired data before superposition to obtain several kinds of imaging results of seismic treatment; and E). making comprehensive comparison and evaluation of imaging quality of said imaging results and defining one optimum seismic imaging result. The invention is extensive in application range, can obtain good effect for making seismic exploration of petroleum and natural gas.

The invention provides a rock longitudinal wave speed prediction method based on an ellipsoid double porosity model. The method comprises the following steps: rock physical parameters are obtained by acquiring well-logging data and experimental observation data, and a dry skeleton model and a porosity fluid model of the rock are generated; a double porosity medium ellipsoid plaque saturation model is established, the potential energy / kinetic energy of the porosity medium and kinetic energy and dissipative equation of the fluid in an inner embedded body are calculation, an Lagrange equation system is derived, and the longitudinal wave and transverse wave speeds are solved; according to the plane wave analysis method, the dispersion relationship of the wave equations is obtained, and the longitudinal wave speed dispersion and attenuation computational formula is obtained. According to the invention, the rock longitudinal wave speed prediction accuracy can be improved.

A method for generating an illumination map. The method includes receiving an earth model that represents one or more properties of the earth. The method then includes receiving a target area in the earth model and propagating one or more wavefields from the target area to a potential seismic data acquisition region in the earth model. After propagating the wavefields, the method includes decomposing the propagated wavefields into one or more subsets of the propagated wavefields and creating the illumination map based on the subsets of the propagated wavefields. The illumination map may indicate one or more contributions of one or more synthetic source and receiver pairs for illuminating the target area in the earth model. The method may then include creating a seismic survey design based on the contributions.

A one-return wave equation migration is used to extrapolate both down-going and up-going waves. Followed by a properly designed imaging condition, the partial image contributed form turning waves is correctly reconstructed. Numerical examples show that this method can significantly enhance definition of an overhanging salt boundary and a geological structure with vertical features.

The invention is a prestack depth migration method for seismic data processing of a relief land surface, determining an imaging reference plane according to the highest point of a collection region, and projecting the relief land surface to a regular grid to fill a constant speed approaching to the near land surface in the relief land surface; using Fourier finite extrapolation operators of one-way wave equation to make wave field extrapolation on shot seismic source wave fields and detector received wave fields along the depth axis, and keeping the wave field between input and output datum equal to zero; imaging the extrapolated shot and detector wave field record by conjugate related imaging conditions, and stacking the related results of each frequency and processing the imaging result into a section image by routine display software, and the invention has an obvious imaging effect of complex construction on the condition of relief land surface.

The invention relates to an anisotropic reverse time migration method for a quasi-P wave equation in a transverse isotropy with a vertical axis of symmetry (VTI) medium. The method comprises the following steps of: (1) performing discrete differentiation on a first-order quasi-P wave equation in a two-dimensional VTI medium and a perfect matched layer (PML) absorption layer boundary equation by adopting staggered meshes to obtain high-order difference formats of forward continuation and reverse continuation of the two equations; (2) performing numeric calculation to obtain a forwards-continued wave field of a shot point and a reverse time continued wave field of a receiver point, and performing normalized correlation operation on the two wave fields to obtain a migration imaging result of each imaging point in a model; and (3) extracting common imaging point gathers from migration results to obtain a final migration profile. By the method, a problem about the migration imaging of intensive transverse speed variation and a high dip angle stratum can be solved; and the influence of anisotropy of the medium is also taken into account, and the good imaging effect of longitudinal wave data acquired from an anisotropic region can be achieved by an anisotropy-theory-based migration method.

The invention provides an interpolation and denoising method and system for seismic data. The method comprises the steps that a preset sampling matrix is acquired; observation data are collected according to the sampling matrix; an objective functional is constructed; according to the observation data and the objective functional, seismic data volume, after interpolation and denoising are performed, corresponding to the observation data is determined. By studying an interpolation and denoising integration scheme of low-signal-to-noise-ratio seismic data, denoising is performed through a threshold value after each time of interpolation, and finally the complete and regular seismic data volume high in signal to noise ratio is obtained. Accordingly, the interpolation and denoising method and system lay a data foundation for multiple processes such as multiple attenuation based on wave equation surface correlation, migration imaging and wave equation inversion.

The invention provides a method for improving the imaging effect of wave equation prestack depth migration, wherein the method is used for processing seismic reflection data in seismic exploration so as to improve the application effect of the wave equation prestack depth migration. The method can be used for improving the imaging effect of the wave equation prestack depth migration by constructing the relationship of residual normal moveout between a depth domain and an angle domain and an interval velocity inversion method to directly update a migration velocity model of the wave equation prestack depth migration based on the residual normal moveout of an angle gather. Based on the relationship of residual normal moveout between the depth domain and the angle domain, the method can be used for removing the residual normal moveout, noise and stretching of a migration gather of the wave equation prestack depth migration, thereby improving the signal-to-noise ratio and resolution of a migration stack profile and improving the quality of the angle gather applied to the prestack inversion. Accordingly, the method can be used for well achieving the direct recognition of underground oil, gas or water, can be used for 2D (two-dimensional) and 3D (three-dimensional) wave equation prestack depth migration of the seismic reflection data, and has important application value for exploration of oil-gas and mineral resources.

The invention provides a time domain analysis method for a transient response of lossy nonuinform multi-conductor transmission lines and aims to solve the problem that the coupling among the lossy nonuniform multi-conductor transmission lines interferes with the computation in the field of integrity of electronic circuit signals. The method comprises the following steps of: performing differential dispersion on space differential operators in an electric wave equation; and integrating time differential operators in the electrical wave equation by a trapezoidal integration method to ensure that the oscillation, which is caused by a central difference method, of a computed result can be eliminated effectively, voltage and current transient response waveforms of any point on the transmissionlines can be acquired; and the transient analysis of the transmission lines in a longer time period is effective, the waveforms can be analyzed in a longer duration without large errors, and the coupled transmission lines are not needed to be decoupled when a model of the transmission lines is built. By the method, the computational load is reduced, and the simulation efficiency is improved. The method can be used for analyzing the transient response of the lossy nonuinform multi-conductor transmission lines.

The present invention provides an efficient time domain full waveform inversion method, belonging to the field of efficient, low storage, high precision and high speed modeling of petroleum seismography data. The efficient time domain full waveform inversion method comprises the following steps: (1) employing original seismic observation records to randomly form supper-shot gathers with different frequency bands through adoption of a source polarity randomizing scheme, a source position randomizing scheme and a source number randomizing scheme; (2) accelerating the forward modeling of high order staggered mesh finite differences of a sound wave equation system and an elastic wave system of the one-form basis of the time domain through adoption of a multi-GPU parallel computing technique on the basis of CUDA, and obtaining a supper-shot gather with the ifreq frequency band simulated by the GPU; and (3) simulating propagations with different seismic sources through multi-GPU nucleus and through adoption of the collaborative technology of the CPU and the GPU at the forward modeling, and storing boundary data of a seismic source wave field and all the wave fields at the last moment to a computer memory.

The invention relates to the prestack depth migration technique during the data processing of petroleum seismic prospecting, including the following concrete steps: establishing depth rate pattern of prestack seismic data; well designing observation grids, and solving wave equation in sphere coordinate system to simulate the communication process of seismic waves; changing and converting the calculated seismic wave field of a frequency domain into a time domain, and fitting the Green function energy spectrum of the seismic wave field in the time domain to detect the arriving time and the vibration amplitude when the maximum energy arrives; transferring the travel time and the vibration amplitude field calculated in the sphere coordinate system to a rectangular coordinate system; accomplishing the maximum energy integral prestack depth migration according to the well calculated travel time and the vibration amplitude; analyzing the operation of speed to the imaging gather of the prestack depth migration to modify the depth rate pattern; and modifying the depth rate pattern by iteration, and outputting the ultimate result until the migration result meets the precision requirement. The invention provides the self-adapting variation difference grid calculation technique so as to achieve the purpose that the difference grids are automatically and gradually thinned along with radius increase, thereby guaranteeing the finite difference calculation precision.