Generating velocity models using full waveform inversion related to subsurface azimuth and reflection angles
阅读说明:本技术 使用与地下方位角和反射角相关的全波形反演产生速度模型 (Generating velocity models using full waveform inversion related to subsurface azimuth and reflection angles ) 是由 康斯坦丁诺斯·钦戈斯 郑宇敦 金永成 于 2018-05-08 设计创作,主要内容包括:本公开描述了用于产生针对地下区域的速度模型的方法和系统,包括计算机实现的方法、计算机程序产品和计算机系统。一种计算机实现的方法包括:在数据处理装置处获得源波场和残余波场,其中,基于与地下区域相关联的地震数据来计算源波场和残余波场;由数据处理装置将源波场和残余波场分解为地下方位角和反射角域;由数据处理装置基于分解后的源波场和残余波场计算多个梯度方向分量;由数据处理装置基于与多个梯度方向分量中的每个梯度方向分量相关联的方位角和角度来缩放多个梯度方向分量;基于缩放后的梯度方向分量计算梯度方向;以及使用梯度方向产生速度模型。(The present disclosure describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for generating a velocity model for a subsurface region. A computer-implemented method comprising: obtaining, at a data processing device, a source wavefield and a residual wavefield, wherein the source wavefield and the residual wavefield are calculated based on seismic data associated with a subsurface region; decomposing, by the data processing device, the source wavefield and the residual wavefield into subsurface azimuth and reflection angle domains; calculating, by the data processing device, a plurality of gradient direction components based on the decomposed source and residual wavefields; scaling, by the data processing apparatus, the plurality of gradient direction components based on the azimuth angle and the angle associated with each of the plurality of gradient direction components; calculating a gradient direction based on the scaled gradient direction component; and generating a velocity model using the gradient direction.)
1. A computer-implemented method for generating a velocity model for a subsurface region, comprising:
obtaining, at a data processing device, a source wavefield and a residual wavefield, wherein the source wavefield and the residual wavefield are calculated based on seismic data associated with the subsurface region;
decomposing, by the data processing device, the source wavefield and the residual wavefield into subsurface azimuth and reflection angle domains;
calculating, by the data processing device, a plurality of gradient direction components based on the decomposed source and residual wavefields;
scaling, by the data processing apparatus, the plurality of gradient direction components based on the azimuth angle and the angle associated with each of the plurality of gradient direction components;
calculating, by the data processing apparatus, a gradient direction based on the scaled gradient direction component; and
generating the velocity model using the gradient direction.
2. The method of claim 1, further comprising: updating the velocity model using an iterative process, wherein a different scaling factor is used in each iteration of the iterative process.
3. The method of claim 2, wherein the iterative process includes a first iteration in which a first scaling factor corresponding to a first angle is used and a second iteration in which a second scaling factor corresponding to a second angle is used, and the first angle is greater than the second angle.
4. The method of claim 1, wherein the decomposing is performed using optical flow processing.
5. The method of claim 1, wherein the scaling comprises: applying a filter to the plurality of gradient direction components based on the range of azimuth angles and angles.
6. The method of claim 1, wherein the scaling comprises: applying a weighting function that weights each of the plurality of gradient direction components with a respective weighting factor.
7. The method of claim 1, wherein the decomposing comprises calculating a reflection angle using optical flow processing.
8. A non-transitory computer-readable medium storing instructions that, when executed, cause a computing device to perform operations comprising:
obtaining, at a data processing device, a source wavefield and a residual wavefield, wherein the source wavefield and the residual wavefield are calculated based on seismic data associated with a subsurface region;
decomposing, by the data processing device, the source wavefield and the residual wavefield into subsurface azimuth and reflection angle domains;
calculating, by the data processing device, a plurality of gradient direction components based on the decomposed source and residual wavefields;
scaling, by the data processing apparatus, the plurality of gradient direction components based on the azimuth angle and the angle associated with each of the plurality of gradient direction components;
calculating, by the data processing apparatus, a gradient direction based on the scaled gradient direction component; and
a velocity model is generated using the gradient direction.
9. The non-transitory computer-readable medium of claim 8, the operations further comprising: updating the velocity model using an iterative process, wherein a different scaling factor is used in each iteration of the iterative process.
10. The non-transitory computer-readable medium of claim 9, wherein the iterative process includes a first iteration in which a first scaling factor corresponding to a first angle is used and a second iteration in which a second scaling factor corresponding to a second angle is used, and the first angle is greater than the second angle.
11. The non-transitory computer-readable medium of claim 8, wherein the decomposing is performed using optical flow processing.
12. The non-transitory computer-readable medium of claim 8, wherein the scaling comprises: applying a filter to the plurality of gradient direction components based on the range of azimuth angles and angles.
13. The non-transitory computer-readable medium of claim 8, wherein the scaling comprises: applying a weighting function that weights each of the plurality of gradient direction components with a respective weighting factor.
14. The non-transitory computer-readable medium of claim 8, wherein the decomposing comprises calculating a reflection angle using optical flow processing.
15. An apparatus, comprising:
at least one hardware processor; and
a non-transitory computer-readable storage medium coupled to the at least one hardware processor and storing programming instructions for execution by the at least one hardware processor, wherein the programming instructions, when executed, cause the at least one hardware processor to perform operations comprising:
obtaining a source wavefield and a residual wavefield, wherein the source wavefield and the residual wavefield are calculated based on seismic data associated with a subsurface region;
decomposing the source wavefield and the residual wavefield into subsurface azimuth and reflection angle domains;
calculating a plurality of gradient direction components based on the decomposed source and residual wavefields;
scaling the plurality of gradient direction components based on an azimuth angle and an angle associated with each of the plurality of gradient direction components;
calculating a gradient direction based on the scaled gradient direction component; and
a velocity model is generated using the gradient direction.
16. The device of claim 15, the operations further comprising: updating the velocity model using an iterative process, wherein a different scaling factor is used in each iteration of the iterative process.
17. The device of claim 16, wherein the iterative process includes a first iteration in which a first scaling factor corresponding to a first angle is used and a second iteration in which a second scaling factor corresponding to a second angle is used, and the first angle is greater than the second angle.
18. The apparatus of claim 15, wherein the decomposition is performed using optical flow processing.
19. The apparatus of claim 15, wherein the scaling comprises: applying a filter to the plurality of gradient direction components based on the range of azimuth angles and angles.
20. The apparatus of claim 15, wherein the scaling comprises: applying a weighting function that weights each of the plurality of gradient direction components with a respective weighting factor.
Technical Field
The present disclosure relates to generating a velocity model of a subsurface structure.
Background
In geophysical analysis, velocity models are used to generate images of subsurface structures. Velocity models are important factors that affect the quality of images of subsurface structures. A good velocity model may provide a better understanding of the subsurface structure and improve the efficiency of oil and gas exploration or production operations within a geographic region.
Disclosure of Invention
The present disclosure describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for generating a velocity model. A computer-implemented method for generating a velocity model comprising: obtaining, at a data processing device, a source wavefield and a residual wavefield, wherein the source wavefield and the residual wavefield are calculated based on seismic data associated with a subsurface region; decomposing, by the data processing device, the source wavefield and the residual wavefield into subsurface azimuth and reflection angle domains; calculating, by the data processing device, a plurality of gradient direction components based on the decomposed source and residual wavefields; scaling, by the data processing apparatus, the plurality of gradient direction components based on the azimuth angle and the angle associated with each of the plurality of gradient direction components; calculating, by the data processing apparatus, a gradient direction based on the scaled gradient direction component; and generating a velocity model using the gradient direction.
Other implementations of this aspect include corresponding computer systems, apparatus, and computer programs recorded on one or more computer storage devices, each configured to perform the actions of the methods. A system of one or more computers may be configured to perform particular operations or actions by software, firmware, hardware, or a combination of software, firmware, or hardware installed on the system that in operation cause the system to perform the actions. One or more computer programs may be configured to perform particular operations or actions by including instructions that, when executed by data processing apparatus, cause the apparatus to perform the actions.
The details of one or more implementations of the subject matter of this specification are set forth in the accompanying drawings and the description below. Other features, aspects, and advantages of the subject matter will become apparent from the description, the drawings, and the claims.
Drawings
This patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the patent and trademark office upon request and payment of the necessary fee.
FIG. 1 (including FIGS. 1A and 1B) illustrates a subsurface image according to an implementation.
FIG. 2 illustrates an example speed model generation process using the AADFWI algorithm according to an implementation.
Fig. 3 is a schematic diagram illustrating the calculation of the reflection angle according to an implementation.
FIG. 4 is a schematic diagram illustrating application of subsurface azimuth and reflection angle filtering during an inversion process according to an implementation.
Fig. 5 (including fig. 5A, 5B, and 5C) illustrates an example sensitivity kernel for gradient directions according to an implementation.
FIG. 6 illustrates an example weighting function according to an implementation.
Fig. 7 (including fig. 7A, 7B, 7C, and 7D) illustrates example numerical tests of the AADFWI according to implementations.
FIG. 8 illustrates an example velocity model generation process according to an implementation.
Figure 9 is a high level architectural block diagram of a geophysical imaging system according to an implementation.
Like reference numbers and designations in the various drawings indicate like elements.
Detailed Description
The following description is presented to enable any person skilled in the art to make and use the disclosed subject matter, and is provided in the context of one or more particular implementations. Various modifications to the disclosed implementations will be readily apparent to those skilled in the art, and the general principles defined herein may be applied to other implementations and applications without departing from the scope of the present disclosure. Thus, the present disclosure is not intended to be limited to the implementations described and/or shown, but is to be accorded the widest scope consistent with the principles and features disclosed herein.
The present disclosure generally describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for generating a velocity model. In some cases, seismic signals may be transmitted into the earth's subsurface by a source device at a source location. Examples of seismic signals include acoustic signals. Seismic signals propagate through the subsurface and may be received by receiver devices placed at receiver locations. In some cases, the source device, the receiver device, or a combination thereof may be placed at the surface. The signal may travel downward until reaching the reflective structure and then reflect upward toward the surface. Because the signal is refracted and reflected by the subsurface structures, the characteristics of the received signal contain information about the subsurface structures. The received signals may be analyzed to generate an image of the subsurface structure.
In marine surveys, air guns and hydrophones may be used as source and receiver devices, respectively. During acquisition, a seismic source is fired from an air gun array. The reflected and refracted signals are acquired by the streamers of the hydrophones. In land acquisition, explosives may be used as the source of the detonation and geophones employed as the receiver equipment. In another example, a vibrating truck may be used as the source device. Other devices that generate and receive seismic signals may also be used.
In some cases, seismic data including sampling of received signals may be collected to establish Common Imaging Gathers (CIGs) or Common Reflection Points (CRPs). In geophysical analysis, a CIG refers to a subset of an entire image with fixed surface locations for measuring variations between local images at fixed image points. If the local image is a function of the angle of reflection, the corresponding CIG is referred to as a co-angle gather (CAG) or an angle domain co-imaging gather (ADCIG).
In some cases, the CAG may be used to determine amplitude variation using angle analysis, offset velocity quality control, and offset velocity analysis. When the shift speed is correct, the images from various angles of incidence are concentrated at the same depth, resulting in a flat event with respect to the CAG. Conversely, when there is an error in velocity, the event in the CAG will become non-flat. A measurement of the residual moveout from a non-flat event can be used to estimate the offset velocity. Such velocity estimation may be performed by using a computed tomography algorithm. The estimated offset velocity can be used to refine the velocity model and determine the correct velocity of the signal propagating through the subsurface structure. With the correct velocity model, an image of the subsurface structure may be generated based on the collected seismic data.
In some implementations, Reverse Time Migration (RTM) may be used in analyzing the migration. The RTM algorithm includes calculations for: forward propagation of the source wavefield, backward propagation of the receiver wavefield, and associated imaging conditions between the two calculated wavefields. In some cases, the source and receiver wavefields may be referred to as forward and backward wavefields, respectively.
Using the two-way wave equation, RTM can produce images of complex subsurface structures (sub-salt reservoirs), focus on steep dip events (salt flank regions), and provide more accurate image amplitudes for seismic interpretation. Accordingly, RTM CAG may be used in the offset velocity analysis rather than using one-way wave equation migration, kirchhoff, or beam migration methods. A cost-effective method of generating CAG based on advanced RTM technology would help to improve the efficiency of subsurface analysis.
In some cases, a seismic Full Waveform Inversion (FWI) may be used to build a velocity model for RTM. FWI is an iterative optimization process that minimizes the error between observed data and modeled data. Conventional FWI typically uses a gradient-based or steepest descent-based approach. In some cases, local minima and highly nonlinear effects can be a problem with this approach. Thus, the successful application of FWI may depend on the following factors: low frequency information, good start model and long offset data (offset data). Wherein the low frequency information and the long offset data may be associated with a data acquisition phase, while the start model may be associated with a data processing phase. If there are no low frequency components in the observed data and the offset range is limited, a good starting velocity model that well delineates the background velocity can be used to solve these problems. Therefore, the reconstructed background velocity (long wavelength velocity) is one of the important factors for FWI application.
The method of reconstructing the long wavelength velocity model may include: the method comprises the steps of using a large damping coefficient to generate extremely low frequency information, obtaining envelope information from a signal by using a Hilbert transform, and adopting a wave field decomposition method to suppress backscattered noise. Separating the gradient direction into a tomographic term and an offset term and enhancing the intrinsic tomographic term of the gradient direction can help to restore the long wavelength structure. In some cases, wavenumber filtering based on scattering angle analysis may be used. Velocity models of different scales are selectively updated by angle domain wavenumber filters from large to small scales during the inversion process. Gradually relaxing the scatter angle allows the velocity model to be approximated in a more stable and accurate manner and avoids falling into local minima.
In some cases, azimuth and angle dependent fwi (aadfwi) may be used to decompose the forward and backward wavefields into the angular domain and the spatial and temporal domains by using optical flow. By avoiding performing transform computations (e.g., hilbert transforms), the method may reduce computational overhead and be more suitable for 3-dimensional (3D) applications. Optical flow is used to determine each vector of the wavefield relative to each direction. Subsurface azimuth and reflection angles between the two wavefields may then be calculated. The new gradient direction compensates for the evanescent long wave information when reflection dominated data is given without low frequency components. The use of large reflection angles helps to restore the long wavelength structure and mitigates the periodic jumps and offset footprints (migrationfotprints) along the gradient direction of the FWI. When the long wavelength structure is restored, the angular range can be widened so that the FWI can be concentrated on the short wavelength structure. Thus, a high resolution velocity model can be obtained.
FIG. 1 (including FIGS. 1A and 1B) illustrates a subsurface image according to an implementation. Fig. 1A and 1B include an
Fig. 2 illustrates an example speed
At 202, forward time domain propagation modeling is performed on the source signal data. The source signal data includes sample data of a source signal transmitted by a source device into the subsurface of the geographic region. In some cases, the source signal data is a record of a source wavelet generated by a source device. Alternatively or in combination, the source signal data may be a computer-generated source excitation using a mathematical function (e.g., a Ricker wavelet function) or other simulation and approximation of the source signal. During the forward time domain propagation modeling process, a plurality of source wavefields are calculated from a source location to a receiver location using the source signal data and an initial velocity model. In one example of an initial velocity model, the velocity at a location point may increase linearly with the depth of the location point. Other examples of initial velocity models may be obtained based on other assumptions or derivation processes. Each source wavefield may be calculated corresponding to a particular time step and a particular position between the source position and the receiver position. Source wavefields corresponding to receiver position are referred to as modeled receiver data
At 204, a residual vector is calculated. The residual vector may be calculated by comparing the modeled receiver data and the observed data. The observed data may be time domain samples of receiver signal data received by the receiver device at the receiver location.
At 206, backward time-domain propagation modeling is performed on the residual vector. In FWI, the residual vector is used as a source for generating the backward wavefield. For example, during a backward time-domain propagation modeling process, a plurality of residual wavefields from a receiver position to a source position are calculated using the residual vectors and the initial velocity model. Each residual wavefield corresponds to a particular time step and a particular position between the receiver position and the source position.
At 210, the source wavefield and the residual wavefield are decomposed using optical flow. The source and residual wavefields may be decomposed into subsurface azimuth and reflection angle domains and time and space domains.
By applying the optical flow to the source and residual wavefields, a motion vector for each wavefield at each time step may be obtained. These motion vectors represent the instantaneous wave propagation at the corresponding time step. The angle between these vectors can be calculated by calculating the dot product of the vectors. For 2-dimensional (2D) analysis, the dot products of the 2 vectors on the x-axis and z-axis are calculated, respectively. For 3-D analysis, the dot products of the 3 vectors are calculated on the x, y and z axes, respectively.
The following is an example optical flow equation for calculating motion vector M:
Dxux+Dyuy+Dzuz+Dt=0 (1)
wherein
Here, DiIndicating the partial derivative of the wavefield P with respect to each propagation direction and time, u represents a vector for each direction, which isProblem M (u)x,uy,uz) The solution of (1). By solving equation (1) using horns and schunck (hs), the solution becomes an optimization process that minimizes an objective function C, which is as follows:
wherein
And α is a weighting factor, then, an iterative solution is given as follows:
where the upper dash symbols represent the weighted average obtained from neighboring points according to the laplacian operator.
Finally, the subsurface reflection angle θ between the two wavefields at a single point is calculated by the following equationrefAnd underground azimuth
And
where p, represented by subscripts s and r, indicates the source and residual wavefields, θ, respectivelyoIs the opening angle between two vectors, nxIs the unit vector on the x-axis and is (1, 0, 0).
Fig. 3 is a schematic diagram illustrating the calculation of the reflection angle according to an implementation. As shown in fig. 3, a source wavefield is computed along a
As shown in the right enlarged view, a
Returning to FIG. 2, the gradient direction of the FWI may be obtained by correlation between the forward (source) and backward (residual) propagating wavefields, as follows:
where C is the objective function of FWI and v represents velocity. Equation (7) can also be expressed in an angular domain approach:
at 212, the source and residual wavefields are filtered. The conventional gradient direction in equation 7 and its decomposed form in equation 8 can be calculated by using the correlation of the forward and backward wavefields at each single point at each time step. In order to extract long wavelength information according to those equations, a filter for angle filtering may be used. The filter for the AADFWI can be expressed as:
wherein
wherein
At 214a, a Fourier transform is performed on the filtered source wavefield to obtain a monochromatic virtual source wavefield relative to a set of frequency components. Similarly, at 214b, a Fourier transform is performed on the filtered residual wavefield to obtain a monochromatic residual wavefield with respect to the set of frequency components. The set of frequency components may be selected to target a range of frequencies to be analyzed in the FWI process.
At 216, gradient directions and Hessian matrices are computed based on the filtered wavefield. Fig. 4 is a schematic diagram 400 illustrating the application of azimuth and angle filtering during an inversion process according to an implementation. The red dots indicate the values obtained by equation 12. When superimposed for all axes, the cube decreases at a single point of the conventional gradient direction cube. If the filter is applied over a range of angles (subsurface reflection angles) (e.g. at θ)1And theta2In between) and azimuth (subsurface azimuth) (e.g., in
Andin between), only gradients within that range of the filter are superimposed for each time step to produce a gradient direction, values outside that range of the filter being padded with zeros. The Hessian matrix may be calculated by automatically correlating the virtual source wavefields.At 218, the velocity model is updated with the filtered gradient directions.
Fig. 5 (including fig. 5A, 5B, and 5C) illustrates an example sensitivity kernel for gradient directions according to an implementation. In FWI, the sensitivity kernel describes how model parameters are perturbed in the earth model with respect to data residual changes. Because the sensitivity kernel includes both low and high wavenumbers (long and short wavelengths) in the conventional gradient direction, FWI processing tends to invert both simultaneously (invert). The results are typically affected by the RTM footprint of the high wavenumber components and when the initial guess is insufficient, the solution may be trapped at local extreme minima. As previously mentioned, the angular filtering can be considered as a sensitivity kernel decomposition. In the conventional gradient direction, it may be difficult to update the background velocity because the high wavenumber structure is concentrated on the offset term in the gradient direction (which is typically much stronger than the tomographic term in the reflection-dominated data). Fig. 5A-5C include
In some cases, a weighting function may be used to attenuate values with respect to azimuth and angle, rather than using the filtering discussed above. In this case, the
Applying a variable weighting to equation 11, the equation can be rewritten as:
wherein
For convenience, the range of the azimuth is omitted and equation 13 can be rewritten as:
wherein
Equations 13 to 16 show that the gradient direction for each angular range can be individually weighted by a weighting factor w1、w2...、wnThe weighting is performed. When the weighting factor is set to 1, the result degenerates to the normal gradient direction. FIG. 6 illustrates example weighting functions 602, 604, and 606 according to an implementation. The weighting function can be expressed as:
y=1-a cos(x) (17)
where a denotes the coefficient of the cosine function. For the
In some cases, a weighting function may be used in addition to the filtering process. For example, instead of padding wavefields outside the filtered range with zeros, weighting factors may be used. For wavefields outside the range of filtering, the value of the weighting factor may be gradually attenuated to approach zero.
Fig. 7 (including fig. 7A, 7B, 7C, and 7D) illustrates example numerical tests of the AADFWI according to implementations. Here, the synthetic Marmousi model was used and the water layer was added from the original model. Time-domain finite difference modeling (Time domain finite difference modeling) has been used to generate synthetic shot gathers for observed data, and 192 sources are used with Ricker wavelets with a dominant frequency of 10 Hz. The maximum offset is 6 km. The time-frequency domain-mixing algorithm uses 9 discrete frequencies from 5Hz to 8Hz for the inversion process. The source wavelet is estimated for each shot and each frequency. Fig. 7A-7D include
FIG. 8 illustrates an example velocity
At 802, a source wavefield and a residual wavefield are obtained at a data processing device. A source wavefield and a residual wavefield are calculated based on seismic data associated with the subsurface region. The seismic data includes receiver signal data and source signal data.
At 804, the source and residual wavefields are decomposed into subsurface azimuth and reflection angle domains using optical flow processing. In some cases, the decomposition is performed based on equations (5) and (6).
At 806, a plurality of gradient direction components are calculated based on the source wavefield and the residual wavefield. In some cases, a plurality of gradient direction components are calculated based on equation (12).
At 808, each of the plurality of gradient direction components is scaled based on the azimuth angle and the angle associated with the respective gradient direction component. In some cases, scaling may be performed using a filtering process based on equation (10). Alternatively, the scaling may be performed using a weighting function based on equation (15).
At 810, a gradient direction is calculated based on the scaled gradient direction component. At 812, a velocity model of the subsurface region is generated based on the gradient directions. In some cases, the velocity model may be updated in an iterative process. The final velocity model may be used to generate a subsurface image. Subsurface images can be used to provide a better understanding of subsurface structures and to improve the efficiency of oil and gas exploration or production operations within a geographic region.
FIG. 9 is a high-level architectural block diagram of a
The described illustration is only one possible implementation of the described subject matter and is not intended to limit the disclosure to a single described implementation. Those of ordinary skill in the art will appreciate the fact that the described components may be connected, combined, or used in an alternative manner in light of the present disclosure.
The
The
The
At a high level, the
The
Each component of the
The
The
The
There may be any number of
Implementations of the described subject matter may include one or more features, either alone or in combination.
For example, in a first implementation, a computer-implemented method for generating a velocity model for a subsurface region, comprises: obtaining, at a data processing device, a source wavefield and a residual wavefield, wherein the source wavefield and the residual wavefield are calculated based on seismic data associated with a subsurface region; decomposing, by the data processing device, the source wavefield and the residual wavefield into subsurface azimuth and reflection angle domains; calculating, by the data processing device, a plurality of gradient direction components based on the decomposed source and residual wavefields; scaling, by the data processing apparatus, the plurality of gradient direction components based on the azimuth angle and the angle associated with each of the plurality of gradient direction components; calculating, by the data processing apparatus, a gradient direction based on the scaled gradient direction component; and generating a velocity model using the gradient direction.
The foregoing and other described implementations may each optionally include one or more of the following features:
a first feature combinable with any of the following features, the method further comprising: the velocity model is updated using an iterative process, where a different scaling factor is used in each iteration of the iterative process.
A second feature combinable with any of the preceding or following features, wherein the iterative process comprises a first iteration in which a first scaling factor corresponding to the first angle is used and a second iteration in which a second scaling factor corresponding to the second angle is used, and the first angle is larger than the second angle.
A third feature combinable with any of the preceding or following features wherein the decomposition is performed using optical flow processing.
Fourth features combinable with any of the preceding or following features, wherein the scaling comprises: a filter is applied to the plurality of gradient direction components based on the range of azimuth angles and angles.
A fifth feature which can be combined with any of the preceding or following features, wherein the scaling comprises: a weighting function is applied that weights each of the plurality of gradient direction components with a respective weighting factor.
A sixth feature combinable with any of the preceding or following features wherein the decomposing comprises calculating the reflection angle using optical flow processing.
In a second implementation, a non-transitory computer-readable medium storing instructions that, when executed, cause a computing device to: obtaining, at a data processing device, a source wavefield and a residual wavefield, wherein the source wavefield and the residual wavefield are calculated based on seismic data associated with a subsurface region; decomposing, by the data processing device, the source wavefield and the residual wavefield into azimuth and angle domains; calculating, by the data processing device, a plurality of gradient direction components based on the decomposed source and residual wavefields; scaling, by the data processing apparatus, the plurality of gradient direction components based on the azimuth angle and the angle associated with each of the plurality of gradient direction components; calculating, by the data processing apparatus, a gradient direction based on the scaled gradient direction component; and generating a velocity model using the gradient direction.
The foregoing and other described implementations may each optionally include one or more of the following features:
a first feature combinable with any of the following features, the operations further comprising: the velocity model is updated using an iterative process, where a different scaling factor is used in each iteration of the iterative process.
A second feature combinable with any of the preceding or following features, wherein the iterative process comprises a first iteration in which a first scaling factor corresponding to the first angle is used and a second iteration in which a second scaling factor corresponding to the second angle is used, and the first angle is larger than the second angle.
A third feature combinable with any of the preceding or following features wherein the decomposition is performed using optical flow processing.
Fourth features combinable with any of the preceding or following features, wherein the scaling comprises: a filter is applied to the plurality of gradient direction components based on the range of azimuth angles and angles.
A fifth feature which can be combined with any of the preceding or following features, wherein the scaling comprises: a weighting function is applied that weights each of the plurality of gradient direction components with a respective weighting factor.
A sixth feature combinable with any of the preceding or following features wherein the decomposing comprises calculating the reflection angle using optical flow processing.
In a third implementation, an apparatus includes: at least one hardware processor; and a non-transitory computer-readable storage medium coupled to the at least one hardware processor and storing programming instructions for execution by the at least one hardware processor, wherein the programming instructions, when executed, cause the at least one hardware processor to perform operations comprising: obtaining a source wavefield and a residual wavefield, wherein the source wavefield and the residual wavefield are calculated based on seismic data associated with the subsurface region; decomposing the source wavefield and the residual wavefield into subsurface azimuth and reflection angle domains; calculating a plurality of gradient direction components based on the decomposed source and residual wavefields; scaling the plurality of gradient direction components based on the azimuth angle and the angle associated with each of the plurality of gradient direction components; calculating a gradient direction based on the scaled gradient direction component; and generating a velocity model using the gradient direction.
The foregoing and other described implementations may each optionally include one or more of the following features:
a first feature combinable with any of the following features, the operations further comprising: the velocity model is updated using an iterative process, where a different scaling factor is used in each iteration of the iterative process.
A second feature combinable with any of the preceding or following features, wherein the iterative process comprises a first iteration in which a first scaling factor corresponding to the first angle is used and a second iteration in which a second scaling factor corresponding to the second angle is used, and the first angle is larger than the second angle.
A third feature combinable with any of the preceding or following features wherein the decomposition is performed using optical flow processing.
Fourth features combinable with any of the preceding or following features, wherein the scaling comprises: a filter is applied to the plurality of gradient direction components based on the range of azimuth angles and angles.
A fifth feature which can be combined with any of the preceding or following features, wherein the scaling comprises: a weighting function is applied that weights each of the plurality of gradient direction components with a respective weighting factor.
A sixth feature combinable with any of the preceding or following features wherein the decomposing comprises calculating the reflection angle using optical flow processing.
The implementation of the subject matter and the functional operations described in this specification can be realized in the following forms: digital electronic circuitry, tangibly embodied computer software or firmware, computer hardware, including the structures disclosed in this specification and their structural equivalents, or combinations of one or more of them. Implementations of the subject matter described in this specification can be implemented as one or more computer programs, i.e., one or more modules of computer program instructions, encoded on a tangible, non-transitory computer storage medium for execution by, or to control the operation of, data processing apparatus. Alternatively or additionally, the program instructions may be encoded on an artificially generated propagated signal (e.g., a machine-generated electrical, optical, or electromagnetic signal) that is generated to encode information for transmission to suitable receiver apparatus for execution by a data processing apparatus. The computer storage medium may be a machine-readable storage device, a machine-readable storage substrate, a random or serial access memory device, or a combination of one or more of them.
The terms "data processing apparatus," "computer," or "electronic computer device" (or equivalents thereof as understood by those of ordinary skill in the art) refer to data processing hardware and include various apparatus, devices, and machines for processing data, including by way of example a programmable processor, a computer, or multiple processors or computers. The apparatus can also be or include special purpose logic circuitry, e.g., a Central Processing Unit (CPU), an FPGA (field programmable gate array), or an ASIC (application-specific integrated circuit). In some implementations, the data processing apparatus and/or dedicated logic circuits may be hardware-based and/or software-based. Alternatively, the apparatus can include code that creates an execution environment for the computer program, e.g., code that constitutes processor firmware, a protocol stack, a database management system, an operating system, or a combination of one or more of them. The present disclosure contemplates the use of data processing devices with or without a conventional operating system (e.g., LINUX, UNIX, WINDOWS, MAC OS, ANDROID, IOS, or any other suitable conventional operating system).
A computer program (also can be called or described as a program, software application, module, software module, script, or code) can be written in any form of programming language, including: a compiled or interpreted language, or a declarative or procedural language, and the computer program may be deployed in any form, including as a stand-alone program or as a module, component, subroutine, or other unit suitable for use in a computing environment. The computer program may, but need not, correspond to a file in a file system. A program can be stored in a portion of a file that holds other programs or data (e.g., one or more scripts stored in a markup language document), in a single file dedicated to the program in question, or in multiple coordinated files (e.g., files that store one or more modules, sub programs, or portions of code). A computer program can be deployed to be executed on one computer or on multiple computers that are located at one site or distributed across multiple sites and interconnected by a communication network. While the portions of the program shown in the figures are illustrated as individual modules implementing various features and functions through various objects, methods, or other processes, the program may alternatively include multiple sub-modules, third party services, components, libraries, and the like, as appropriate. Rather, the features and functionality of the various components may be combined into a single component as appropriate.
The processes and logic flows described in this specification can be performed by one or more programmable computers executing one or more computer programs to perform functions by operating on input data and generating output. The processes and logic flows can also be performed by, and apparatus can also be implemented as, special purpose logic circuitry, e.g., a CPU, FPGA or ASIC.
A computer suitable for executing a computer program may be based on a general purpose or special purpose microprocessor, both, or any other type of CPU. Generally, a CPU will receive instructions and data from a read-only memory (ROM) or a Random Access Memory (RAM) or both. The essential elements of a computer are a CPU for executing instructions and one or more memory devices for storing instructions and data. Generally, a computer will also include, or be operatively coupled to receive data from or transfer data to, or both, one or more mass storage devices for storing data, e.g., magnetic, magneto-optical disks, or optical disks. However, the computer need not have these devices. Moreover, a computer may be embedded in another device, e.g., a mobile telephone, a Personal Digital Assistant (PDA), a mobile audio or video player, a game player, a Global Positioning System (GPS) receiver, or a portable storage device (e.g., a Universal Serial Bus (USB) flash drive), to name a few.
Computer-readable media (transitory or non-transitory) suitable for storing computer program instructions and data include all forms of non-volatile memory, media and memory devices, including by way of example semiconductor memory devices (e.g., erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), and flash memory devices); magnetic disks (e.g., internal hard disks or removable disks); magneto-optical disks; and CD ROM, DVD +/-R, DVD-RAM and DVD-ROM disks. The memory may store various objects or data, including: caches, classes (classes), frames, applications, backup data, jobs, web pages, web page templates, database tables, knowledge bases storing business and/or dynamic information, and any other suitable information including any parameters, variables, algorithms, instructions, rules, constraints, references thereto. Further, the memory may include any other suitable data, such as well logs, policies, security or access data, reporting files, and so forth. The processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry.
To provide for interaction with a user, the subject matter described in this specification can be implemented on a computer having a display device (e.g., a CRT (cathode ray tube), LCD (liquid crystal display), LED (light emitting diode), or plasma monitor) for displaying information to the user and a keyboard and a pointing device (e.g., a mouse, a trackball, or a trackpad) by which the user can provide input to the computer. Touch screens (e.g., a tablet computer surface with pressure sensitivity, a multi-touch screen using capacitive or electrical sensing, or other types of touch screens) may also be used to provide input to the computer. Other types of devices may also be used to provide for interaction with a user; for example, feedback provided to the user can be any form of sensory feedback, e.g., visual feedback, auditory feedback, or tactile feedback; and input from the user may be received in any form, including acoustic, speech, or tactile input. Further, the computer may interact with the user by sending and receiving documents to and from the device used by the user; for example, a user is interacted with by sending a web page to a web browser on a user client device in response to a request received from the web browser.
The term "graphical user interface" or GUI may be used in the singular or plural to describe one or more graphical user interfaces and each display of a particular graphical user interface. Thus, the GUI may represent any graphical user interface, including but not limited to a web browser, touch screen, or Command Line Interface (CLI) that processes information and effectively presents the results of the information to the user. In general, the GUI may include a number of User Interface (UI) elements, some or all of which are associated with a web browser, such as interactive fields, drop-down lists, and buttons that are operable by a business suite user. These and other UI elements may be related to or represent functionality of a web browser.
Implementations of the subject matter described in this specification can be implemented in a computing system that includes a back-end component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front-end component (e.g., a client computer having a graphical user interface or a web browser through which a user can interact with an implementation of the subject matter described in this specification), or any combination of one or more such back-end, middleware, or front-end components. The components of the system can be interconnected by any form or medium of wired and/or wireless digital data communication (e.g., a communication network). Examples of communication networks include a Local Area Network (LAN), a Radio Access Network (RAN), a Metropolitan Area Network (MAN), a Wide Area Network (WAN), Worldwide Interoperability for Microwave Access (WIMAX), a Wireless Local Area Network (WLAN) using, for example, 802.11a/b/g/n and/or 802.20, all or a portion of the Internet, and/or any other communication system at one or more locations. The network may communicate, for example, Internet Protocol (IP) packets, frame relay frames, Asynchronous Transfer Mode (ATM) cells, voice, video, data, and/or other suitable information between network addresses.
The computing system may include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other.
In some implementations, any or all of the components (hardware and/or software) of the computing system may interface with each other and/or with interfaces using an Application Programming Interface (API) and/or a service layer. The API may include specifications for routines, data structures, and object classes. An API may be independent or dependent on the computer language and refers to a complete interface, a single function, or even a collection of APIs. The service layer provides software services to the computing system. The functionality of the various components of the computing system may be accessible to all service consumers via the service layer. Software services provide reusable, defined business functions through defined interfaces. For example, the interface may be software written in JAVA, C + +, or other suitable language that provides data in an extensible markup language (XML) format or other suitable format. The API and/or service layer may be an integrated component and/or a stand-alone component related to other components of the computing system. Further, any or all portions of the service layer may be implemented as sub-modules or sub-modules of another software module, enterprise application, or hardware module without departing from the scope of the present disclosure.
While this specification contains many specific implementation details, these should not be construed as limitations on the scope of any inventions or of what may be claimed, but rather as descriptions of features specific to particular implementations of particular inventions. Certain features that are described in this specification in the context of separate implementations can also be implemented in combination in a single implementation. Conversely, various features that are described in the context of a single implementation can also be implemented in multiple implementations separately or in any suitable subcombination. Furthermore, although features may be described above as acting in certain combinations and even initially claimed as such, one or more features from a claimed combination can in some cases be excised from the combination, and the claimed combination may be directed to a subcombination or variation of a subcombination.
Particular implementations of the present subject matter have been described. Other implementations, changes, and substitutions to the described implementations are apparent to those of skill in the art and are within the scope of the appended claims. Although operations are depicted in the drawings or claims in a particular order, this should not be understood as: it may be desirable to perform the operations in the particular order shown, or in sequential order, or to perform all of the operations shown (some of which may be viewed as optional) in order to achieve the desired results. In some scenarios, multitasking and parallel processing may be advantageous.
Moreover, the separation and/or integration of various system modules and components in the implementations described above should not be understood as requiring such separation and/or integration in all implementations, and it should be understood that the described program components and systems can generally be integrated together in a single software product or packaged into multiple software products.
Accordingly, the above description of example implementations does not define or constrain this disclosure. Other changes, substitutions, and alterations are also possible without departing from the spirit and scope of this disclosure.
- 上一篇:一种医用注射器针头装配设备
- 下一篇:抑制地震数据中的噪声