Three‐term inversion of prestack seismic data using a weighted l2, 1 mixed norm          下载免费PDF全文

Daniel O. Pérez  Danilo R. Velis  Mauricio D. Sacchi 《Geophysical Prospecting》2017,65(6):1477-1495

We present a new inversion method to estimate, from prestack seismic data, blocky P‐ and S‐wave velocity and density images and the associated sparse reflectivity levels. The method uses the three‐term Aki and Richards approximation to linearise the seismic inversion problem. To this end, we adopt a weighted mixed l_{2, 1}‐norm that promotes structured forms of sparsity, thus leading to blocky solutions in time. In addition, our algorithm incorporates a covariance or scale matrix to simultaneously constrain P‐ and S‐wave velocities and density. This a priori information is obtained by nearby well‐log data. We also include a term containing a low‐frequency background model. The l_{2, 1} mixed norm leads to a convex objective function that can be minimised using proximal algorithms. In particular, we use the fast iterative shrinkage‐thresholding algorithm. A key advantage of this algorithm is that it only requires matrix–vector multiplications and no direct matrix inversion. The latter makes our algorithm numerically stable, easy to apply, and economical in terms of computational cost. Tests on synthetic and field data show that the proposed method, contrarily to conventional l₂‐ or l₁‐norm regularised solutions, is able to provide consistent blocky and/or sparse estimators of P‐ and S‐wave velocities and density from a noisy and limited number of observations.  相似文献   

Some remarks on Q-compensated sparse deconvolution without knowing the quality factor Q     

Xintao Chai  Ronghua Peng  Genyang Tang  Wei Chen  Jingnan Li 《Geophysical Prospecting》2019,67(8):2003-2021

The subsurface media are not perfectly elastic, thus anelastic absorption, attenuation and dispersion (aka Q filtering) effects occur during wave propagation, diminishing seismic resolution. Compensating for anelastic effects is imperative for resolution enhancement. Q values are required for most of conventional Q-compensation methods, and the source wavelet is additionally required for some of them. Based on the previous work of non-stationary sparse reflectivity inversion, we evaluate a series of methods for Q-compensation with/without knowing Q and with/without knowing wavelet. We demonstrate that if Q-compensation takes the wavelet into account, it generates better results for the severely attenuated components, benefiting from the sparsity promotion. We then evaluate a two-phase Q-compensation method in the frequency domain to eliminate Q requirement. In phase 1, the observed seismogram is disintegrated into the least number of Q-filtered wavelets chosen from a dictionary by optimizing a basis pursuit denoising problem, where the dictionary is composed of the known wavelet with different propagation times, each filtered with a range of possible values. The elements of the dictionary are weighted by the infinity norm of the corresponding column and further preconditioned to provide wavelets of different values and different propagation times equal probability to entry into the solution space. In phase 2, we derive analytic solutions for estimates of reflectivity and Q and solve an over-determined equation to obtain the final reflectivity series and Q values, where both the amplitude and phase information are utilized to estimate the Q values. The evaluated inversion-based Q estimation method handles the wave-interference effects better than conventional spectral-ratio-based methods. For Q-compensation, we investigate why sparsity promoting does matter. Numerical and field data experiments indicate the feasibility of the evaluated method of Q-compensation without knowing Q but with wavelet given.  相似文献   

Effect of Polianthes tuberosa (Amaryllidaceae) on the Reproduction and Biochemical Parameters in the Ovotestis of Snail Lymnaea acuminata (Mollusca,Pulmonata)     

K. Singh  V. K. Singh  D. K. Singh 《洁净——土壤、空气、水》1999,27(1):32-37

The present study elucidates the effect of powder of Polianthes tuberosa bulb and its sapogenins tigogenin and hecogenin on the reproduction of the snail Lymnaea acuminata. It has been reported that the two sapogenins tigogenin and hecogenin caused a significant reduction in fecundity, hatchability, and survival of young snails. Exposure to sublethal concentration caused a small but significant reduction in protein, amino acids, nucleic acids, and phosphatase activities in gonadal tissue of treated snails. Withdrawal after 7 days indicates that the effect of powder of P. tuberosa bulb, tigogenin, and hecogenin on the reproduction and different biochemical parameters is reversible.  相似文献   

Spatial constrained inverse rock physics modelling   总被引：1，自引：0，他引：1       下载免费PDF全文

Erling H. Jensen  Tor Arne Johansen 《Geophysical Prospecting》2015,63(1):183-191

Predicting reservoir parameters, such as porosity, lithology, and saturations, from geophysical parameters is a problem with non‐unique solutions. The variance in solutions can be extensive, especially for saturation and lithology. However, the reservoir parameters will typically vary smoothly within certain zones—in vertical and horizontal directions. In this work, we integrate spatial correlations in the predicted parameters to constrain the range of predicted solutions from a particular type of inverse rock physics modelling method. Our analysis is based on well‐log data from the Glitne field, where vertical correlations with depth are expected. It was found that the reservoir parameters with the shortest depth correlation (lithology and saturation) provided the strongest constraint to the set of solutions. In addition, due to the interdependence between the reservoir parameters, constraining the predictions by the spatial correlation of one parameter also reduced the number of predictions of the other two parameters. Moreover, the use of additional constraints such as measured log data at specific depth locations can further narrow the range of solutions.  相似文献   

Joint parameter estimation from magnetic resonance and vertical electric soundings using a multi‐objective genetic algorithm     

İrfan Akca  Thomas Günther  Mike Müller‐Petke  Ahmet T. Başokur  Ugur Yaramanci 《Geophysical Prospecting》2014,62(2):364-376

Magnetic resonance sounding (MRS) has increasingly become an important method in hydrogeophysics because it allows for estimations of essential hydraulic properties such as porosity and hydraulic conductivity. A resistivity model is required for magnetic resonance sounding modelling and inversion. Therefore, joint interpretation or inversion is favourable to reduce the ambiguities that arise in separate magnetic resonance sounding and vertical electrical sounding (VES) inversions. A new method is suggested for the joint inversion of magnetic resonance sounding and vertical electrical sounding data. A one‐dimensional blocky model with varying layer thicknesses is used for the subsurface discretization. Instead of conventional derivative‐based inversion schemes that are strongly dependent on initial models, a global multi‐objective optimization scheme (a genetic algorithm [GA] in this case) is preferred to examine a set of possible solutions in a predefined search space. Multi‐objective joint optimization avoids the domination of one objective over the other without applying a weighting scheme. The outcome is a group of non‐dominated optimal solutions referred to as the Pareto‐optimal set. Tests conducted using synthetic data show that the multi‐objective joint optimization approximates the joint model parameters within the experimental error level and illustrates the range of trade‐off solutions, which is useful for understanding the consistency and conflicts between two models and objectives. Overall, the Levenberg‐Marquardt inversion of field data measured during a survey on a North Sea island presents similar solutions. However, the multi‐objective genetic algorithm method presents an efficient method for exploring the search space by producing a set of non‐dominated solutions. Borehole data were used to provide a verification of the inversion outcomes and indicate that the suggested genetic algorithm method is complementary for derivative‐based inversions.  相似文献   

A framework for the seismic assessment of existing masonry buildings accounting for different sources of uncertainty          下载免费PDF全文

M. Rota  A. Penna  G. Magenes 《地震工程与结构动力学》2014,43(7):1045-1066

The current formulation of Eurocode 8 Part 3 and the Italian building code for the seismic assessment of existing buildings accounts for epistemic (knowledge‐based) uncertainties by means of the identification of knowledge levels with associated values of the so‐called confidence factors, applied only as a reduction of material strengths. This formulation does not always produce consistent results and it does not explicitly account for other sources of uncertainty. The paper proposes a probabilistic methodology for the quantification of appropriately defined factors, allowing consideration of the different sources of uncertainty involved in the seismic assessment of masonry buildings by means of nonlinear static analyses. This simple approach, also including an alternative formulation of the confidence factors related with material properties, allows to obtain results which are consistent with the acquired level of knowledge and correctly account for the different sources of uncertainty without requiring to carry out any stochastic nonlinear analysis. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

草原生态系统模式中参数不确定性导致的模拟结果不确定性研究   总被引：1，自引：1，他引：0  

谢东东  孙国栋  邵爱梅  穆穆 《气候与环境研究》2013,18(3):375-386

基于五变量草原生态系统理论模式,应用与参数有关的条件非线性最优扰动(CNOP-P)方法,探讨了由参数不确定性导致的草原生态系统模式模拟结果的不确定性问题。参数的不确定性可能来源于观测和（或）对物理过程描述等的不确定性。选取了五变量草原生态系统模式中具有物理意义的32个模式参数进行数值试验。试验结果表明,对所考察的32个模式参数,在一定的不确定性和给定的优化时刻范围内,单独优化每个参数所得CNOP-Ps的联合模态与同时优化32个参数所得CNOP-P的模态并不相同。比较了上述两类参数误差以及随机参数误差对草原生态系统模拟的差异。随机参数误差与上述优化方法所得参数误差的不确定性范围大小相同。数值结果表明,同时优化32个参数所得 CNOP-P 类型参数误差使得草原生态系统模拟的不确定性程度最大。这种影响表现在使得草原生态系统转变为沙漠生态系统,或者使得草原生态系统转变为具有更多生草量的草原生态系统。上述数值结果不依赖于优化时间和参数不确定性程度的大小。这些数值结果建议我们应当考虑多参数的非线性相互作用来研究草原生态系统模式模拟的不确定性问题,并且揭示出CNOP-P方法是讨论上述问题的一个有用的工具。  相似文献