排序方式： 共有23条查询结果，搜索用时 15 毫秒
1.
A nutrient dynamic model coupled with a 3D physical model has been developed to study the annual cycle of phytoplankton production in the Yellow Sea. The biological model involves interactions between inorganic nitrogen (nitrate and ammonium), phosphate and phytoplankton biomass. The model successfully reproduces the main features of phytoplankton-nutrient variation and dynamics of production. 1. The well-mixed coastal water is characterized by high primary production, as well as high new production. 2. In summer, the convergence of tidal front is an important hydrodynamic process, which contributes to high biomass at frontal areas. 3. The evolution of phytoplankton blooms and thermocline in the central region demonstrate that mixing is a dominant factor to the production in the Yellow Sea. In this simulation, nitrate- and ammonium-based productions are estimated regionally and temporally. The northern Yellow Sea is one of the highly ranked regions in the Yellow Sea for the capability of fixing carbon and nitrogen. The annual averaged f-ratio of 0.37 indicates that regenerated production prevails over the Yellow Sea. The result also shows that phosphate is the major nutrient, limiting phytoplankton growth throughout the year and it can be an indicator to predict the bloom magnitude. Finally, the relative roles of external nutrient sources have been evaluated, and benthic fluxes might play a significant role in compensating 54.6% of new nitrogen for new production consumption. 相似文献
2.
This is the first in a series of three papers focused on using variants of a logarithmic objective function approach to full waveform inversion. In this article, we investigate waveform inversion using full logarithmic principles and compare the results with the conventional least squares approach. We demonstrate theoretically that logarithmic inversion is computational similar to the conventional method in the sense that it uses exactly the same back‐propagation technology as used in least‐squares inversion. In the sense that it produces better results for each of three numerical examples, we conclude that logarithmic inversion is also more robust. We argue that a major reason for the inherent robustness is the fact that the logarithmic approach produces a natural scaling of the amplitude of the residual wavefield by the amplitude of the modelled wavefield that tends to stabilize the computations and consequently improve the final result. We claim that any superiority of the logarithmic inversion is based on the fact that it tends to be tomographic in the early stage of the inversion and more dependent on amplitude differences in the latter stages. 相似文献
3.
In the second paper of this three part series, we studied the case of conventional and logarithmic phase‐only approaches to full‐waveform inversion. Here, we concentrate on deriving amplitude‐only approaches for both conventional‐ and logarithmic‐based methods. We define two amplitude‐only objective functions by simply assuming that the phase of the modelled wavefield is equal to that of the observed wavefield. We do this for both the conventional least‐squares approach and the logarithmic approach of Shin and Min. We show that these functions can be optimized using the same reverse‐time propagation algorithm of the full conventional methodology. Although the residuals in this case are not really residual wavefields, they can both be considered and utilized in that sense. In contrast to the case for our phase‐only algorithms, we show through numerical tests that the conventional amplitude‐only inversion is better than the logarithmic method. 相似文献
4.
Improved amplitude preservation for prestack depth migration by inverse scattering theory 总被引：6，自引：1，他引：6
A prestack reverse time-migration image is not properly scaled with increasing depth. The main reason for the image being unscaled is the geometric spreading of the wavefield arising during the back-propagation of the measured data and the generation of the forward-modelled wavefields. This unscaled image can be enhanced by multiplying the inverse of the approximate Hessian appearing in the Gauss–Newton optimization technique. However, since the approximate Hessian is usually too expensive to compute for the general geological model, it can be used only for the simple background velocity model.We show that the pseudo-Hessian matrix can be used as a substitute for the approximate Hessian to enhance the faint images appearing at a later time in the 2D prestack reverse time-migration sections. We can construct the pseudo-Hessian matrix using the forward-modelled wavefields (which are used as virtual sources in the reverse time migration), by exploiting the uncorrelated structure of the forward-modelled wavefields and the impulse response function for the estimated diagonal of the approximate Hessian. Although it is also impossible to calculate directly the inverse of the pseudo-Hessian, when using the reciprocal of the pseudo-Hessian we can easily obtain the inverse of the pseudo-Hessian. As examples supporting our assertion, we present the results obtained by applying our method to 2D synthetic and real data collected on the Korean continental shelf. 相似文献
5.
6.
Ho Seuk Bae Sukjoon Pyun Wookeen Chung Seung‐Goo Kang Changsoo Shin 《Geophysical Prospecting》2012,60(3):413-432
We developed a frequency‐domain acoustic‐elastic coupled waveform inversion based on the Gauss‐Newton conjugate gradient method. Despite the use of a high‐performance computer system and a state‐of‐the‐art parallel computation algorithm, it remained computationally prohibitive to calculate the approximate Hessian explicitly for a large‐scale inverse problem. Therefore, we adopted the conjugate gradient least‐squares algorithm, which is frequently used for geophysical inverse problems, to implement the Gauss‐Newton method so that the approximate Hessian is calculated implicitly. Thus, there was no need to store the Hessian matrix. By simultaneously back‐propagating multi‐components consisting of the pressure and displacements, we could efficiently extract information on the subsurface structures. To verify our algorithm, we applied it to synthetic data sets generated from the Marmousi‐2 model and the modified SEG/EAGE salt model. We also extended our algorithm to the ocean‐bottom cable environment and verified it using ocean‐bottom cable data generated from the Marmousi‐2 model. With the assumption of a hard seafloor, we recovered both the P‐wave velocity of complicated subsurface structures as well as the S‐wave velocity. Although the inversion of the S‐wave velocity is not feasible for the high Poisson's ratios used to simulate a soft seafloor, several strategies exist to treat this problem. Our example using multi‐component data showed some promise in mitigating the soft seafloor effect. However, this issue still remains open. 相似文献
7.
Waveform inversion in the Laplace domain 总被引：12，自引：0，他引：12
8.
Elastic waves, such as Rayleigh and mode‐converted waves, together with amplitude versus offset variations, serve as noise in full waveform inversion using the acoustic approximation. Heavy preprocessing must be applied to remove elastic effects to invert land or marine data using the acoustic inversion method in the time or frequency domains. Full waveform inversion using the elastic wave equation should be one alternative; however, multi‐parameter inversion is expensive and sensitive to the starting velocity model. We implement full acoustic waveform inversion of synthetic land and marine data in the Laplace domain with minimum preprocessing (i.e., muting) to remove elastic effects. The damping in the Laplace transform can be thought of as an automatic time windowing. Numerical examples show that Laplace‐domain acoustic inversion can yield correct smooth velocity models even with the noise originating from elastic waves. This offers the opportunity to develop an accurate smooth starting model for subsequent inversion in the frequency domain. 相似文献
9.
Eunjin Park Wansoo Ha Wookeen Chung Changsoo Shin Dong-Joo Min 《Pure and Applied Geophysics》2013,170(12):2075-2085
The wavefield in the Laplace domain has a very small amplitude except only near the source point. In order to deal with this characteristic, the logarithmic objective function has been used in many Laplace domain inversion studies. The Laplace-domain waveform inversion using the logarithmic objective function has fewer local minima than the time- or frequency domain inversion. Recently, the power objective function was suggested as an alternative to the logarithmic objective function in the Laplace domain. Since amplitudes of wavefields are very small generally, a power <1 amplifies the wavefields especially at large offset. Therefore, the power objective function can enhance the Laplace-domain inversion results. In previous studies about synthetic datasets, it is confirmed that the inversion using a power objective function shows a similar result when compared with the inversion using a logarithmic objective function. In this paper, we apply an inversion algorithm using a power objective function to field datasets. We perform the waveform inversion using the power objective function and compare the result obtained by the logarithmic objective function. The Gulf of Mexico dataset is used for the comparison. When we use a power objective function in the inversion algorithm, it is important to choose the appropriate exponent. By testing the various exponents, we can select the range of the exponent from 5 × 10?3 to 5 × 10?8 in the Gulf of Mexico dataset. The results obtained from the power objective function with appropriate exponent are very similar to the results of the logarithmic objective function. Even though we do not get better results than the conventional method, we can confirm the possibility of applying the power objective function for field data. In addition, the power objective function shows good results in spite of little difference in the amplitude of the wavefield. Based on these results, we can expect that the power objective function will produce good results from the data with a small amplitude difference. Also, it can partially be utilized at the sections where the amplitude difference is very small. 相似文献
10.
Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix 
下载免费PDF全文
下载免费PDF全文 To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems. 相似文献