3-D multi-parameter type traveltime tomography in a spherical coordinate frame: comparison of double and triple class simultaneous inversions

It is now common practice to perform simultaneous traveltime inversion for the velocity field and the reflector geometry in reflection/refraction tomography, or the velocity field and the hypocenter locations in regional earthquake tomography, but seldom are all three classes of model parameters updated simultaneously. This is mainly due to the trade-off between the different types of model parameters and the lack of different seismic phases to constrain the model parameters. Using a spherical-coordinate ray tracing algorithm for first and later（primary reflected） arrival tracing algorithm in combination with a popular linearized inversion solver, it is possible to simultaneously recover the three classes of model parameters in regional or global tomographic studies. In this paper we incorporate the multistage irregular shortest-path ray tracing algorithm（in a spherical coordinate system） with a subspace inversion solver to formulate a simultaneous inversion algorithm for triple model parameters updating using direct and later arrival time information.Comparison tests for two sets of data（noise free and added noise） indicate that the new triple-class parameter inversion algorithm is capable of obtaining nearly the same results as the double-class parameter inversion scheme. Furthermore,the proposed multi-parameter type inversion method is not sensitive to a modest level of picking error in the traveltime data, and also performs well with a relatively large uncertainty in earthquake hypocentral locations. This shows it to be a feasible and promising approach in regional or global tomographic applications.

Simultaneous Waveform Inverse Modelling for Litho-Fluid Prediction in an Old Marginal, “Agbbo”Field, Onshore Niger Delta, Nigeria

Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with uncertainties, especially in marginal fields. An approach is employed in this study that integrated rock physics and waveform inverse modelling for lithology and fluid-type characterization to appropriately identify potential hydrocarbon saturated zones and their corresponding lithology. Seismic and well-log data were analyzed using Hampson Russel software. The method adopted includes lithofacies and fluid content analysis using rock physics parameters and seismic simultaneous inverse modelling. Rock physics analysis identified 2 broad reservoirs namely: HDZ1 and HDZ2 reservoirs. Results from the inverse modelling showed that low values of acoustic impedance from 19,743 to 20,487 (ft/s)(g/cc) reflect hydrocarbon-bearing reservoirs while medium to high values shows brine and shale respectively, with brine zone ranging from 20,487 to 22,531 (ft/s)(g/cc) and shale above 22,531 (ft/s)(g/cc). Two lithofacies were identified from inversion analysis of Vp/Vs and Mu-Rho, namely: sand and shale with VpVs 1.95 values respectively. Mu-Rho > 12.29 (GPa)(g/cc) and <12.29 (GPa) (g/cc) represent sand and shale respectively. From 3D volume, it was observed that a high accumulation of hydrocarbon was observed to be saturated at the north to the eastern part of the field forming a meandering channel. Sands were mainly distributed around the northeastern to the southwestern part of the field, that tends to be away from Well 029. This was also validated by the volume of rigidity modulus (Mu-Rho) showing high values indicating sands fall within the northeastern part of the field.

AVA Simultaneous Inversion of Prestack Seismic Data Using Particle Swarm Optimization

A new prestack AVA simultaneous inversion using particle swarm optimization algorithm is proposed, which can obtain the elastic parameters such as P-wave and S-wave impedance from P-wave reflection data simultaneously. Compared with the conventional AVA inversion based on generalized linear technique, this method does not depend on the initial model and can reach the global minimum. In order to increase the stability of the inversion, low-frequency trends of P-wave and S-wave impedances are built into the inversion. This method has been successfully applied to synthetic and field data. The estimated P-wave and S-wave impedances can be combined to derive other elastic parameters, which are sensitive for lithology identification and fluid prediction.

A high resolution inversion method for fluid factor with dynamic dryrock V_(P)/V_(S) ratio squared

As an important indicator parameter of fluid identification,fluid factor has always been a concern for scholars.However,when predicting Russell fluid factor or effective pore-fluid bulk modulus,it is necessary to introduce a new rock skeleton parameter which is the dry-rock VP/VS ratio squared(DVRS).In the process of fluid factor calculation or inversion,the existing methods take this parameter as a static constant,which has been estimated in advance,and then apply it to the fluid factor calculation and inversion.The fluid identification analysis based on a portion of the Marmousi 2 model and numerical forward modeling test show that,taking the DVRS as a static constant will limit the identification ability of fluid factor and reduce the inversion accuracy.To solve the above problems,we proposed a new method to regard the DVRS as a dynamic variable varying with depth and lithology for the first time,then apply it to fluid factor calculation and inversion.Firstly,the exact Zoeppritz equations are rewritten into a new form containing the fluid factor and DVRS of upper and lower layers.Next,the new equations are applied to the four parameters simultaneous inversion based on the generalized nonlinear inversion(GNI)method.The testing results on a portion of the Marmousi 2 model and field data show that dynamic DVRS can significantly improve the fluid factor identification ability,effectively suppress illusion.Both synthetic and filed data tests also demonstrate that the GNI method based on Bayesian deterministic inversion(BDI)theory can successfully solve the above four parameter simultaneous inversion problem,and taking the dynamic DVRS as a target inversion parameter can effectively improve the inversion accuracy of fluid factor.All these results completely verified the feasibility and effectiveness of the proposed method.

Simultaneous inversion of configurations of C_(22),C_(23) in 24S-ethyl-3,5-cyclo-cholestane-22S,23S-diol-6-one

In this paper,obtainment of 3 by simultaneous inversion of configurations of C_(22),C_(23) in 2 via its ditosylate or dimesylate was unsuccessful,leading to elimination products instead.However, monotosylate of 2 can be readily obtained with exclusive regioselectivity,which in turn gave a single epoxide 7b Hydrolysis of 7b gave 3 and 2 is a ratio of 2:1,with a combined yield of 80%.The recovered 2 can of course be recycled again.

Joint PP and PS AVO inversion based on Zoeppritz equations

Considering Zoeppritz equations, reflections of PP and PS are only the function of ratios of density and velocity. So the inversion results will be the same if the ratios are the same but values of density, velocities of P- wave and S-wave are different without strict constraint. This paper makes efforts to explore nonlinear simultaneous PP and PS inversion with expectation to reduce the ambiguity of AVO analysis by utilizing the redundancy of multi-component AVO measurements. Accurate estimation of ratio parameters depends on independence of input data. There are only two independent AVO attributes for PP reflectivity （i.e. intercept and gradient） and two for PS reflectivity （i.e. pseudo-intercept and pseudo-gradient or extreme amplitude）, respectively. For individual PP and PS inversion, the values of least-squares objective function do not converge around a large neighborhood of chosen true model parameters. Fortunately for joint PP and PS inversion the values of the least-squares objective function show closed contours with single minima. Finally the power function fitting is used to provide a higher precision AVO attributes than traditional polynomial fitting. By using the four independent fitting attributes （two independent attributes for PP and PS respectively）, the inversion of four ratio parameters （velocities and densities） would be estimated with less errors than that in traditional method.

ON SIMULTANEOUS APPROXIMATION TO A DIFFERENTIABLE FUNCTION AND ITS DERIVATIVE BY INVERSE PAL-TYPE INTERPOLATION POLYNOMIALS

Let ξn-1<ξn-2 <ξn-2 <… < ξ1 be the zeros of the the (n -1)-th Legendre polynomial Pn-1(x) and - 1 = xn < xn-1 <… < x1 = 1 the zeros of the polynomial W n(x) =- n(n - 1) Pn-1(t)dt = (1 -x2)P'n-1(x). By the theory of the inverse Pal-Type interpolation, for a function f(x) ∈ C[-1 1], there exists a unique polynomial Rn(x) of degree 2n - 2 (if n is even) satisfying conditions Rn(f,ξk) = f(∈ek)(1≤ k≤ n - 1) ;R'n(f,xk) = f'(xk)(1≤ k≤ n). This paper discusses the simultaneous approximation to a differentiable function f by inverse Pal-Type interpolation polynomial {Rn(f,x)} (n is even) and the main result of this paper is that if f ∈ C'[1,1], r≥2, n≥ + 2> and n is even thenholds uniformly for all x ∈ [- 1,1], where h(x) = 1 +

Earthquake relocation and 3-dimensional crustal structure of P-wave velocity in central-western China

A simultaneous inversion of earthquake relocation and three-dimensional crustal structure of P-wave velocity in central-western China (21癗~36癗, 98癊~112癊) were performed in this paper. The crustal P-wave velocity model and earthquake relocation for this region are obtained using Pg and Sg phase readings of 9 988 earthquakes from 1992 to 1999 recorded at 193 seismic stations within central-western China by SPHYPIT90 and SPHREL3D90 programs. A lateral inhomogeneous structure of P-wave velocity in this region was obtained. Ob-vious contrast of P-wave velocities was revealed on both sides of active fault zones. Relocated epicenters of 6 459 events show clear lineation along active faults, which indicated a close correlation between seismicity and the active faults in this region. Focal depths of 82% relocated events ranged from 0 to 20 km, which is in good agreement with that from double-difference earthquake location algorithm.

Preliminary estimation for the focal mechanisms of 1996 Lijiang main shock in Yunnan Province and studies of related problems

The mechanisms of the February 3, 1996 Lijiang main shock, Yunnan Province, are estimated by using the principle to inverse the mechanisms of two point sources simultaneously. The results are that the main shock of Lijiang consists of two large ruptures, the time difference and the distance between the two ruptures are about 12 s (by the inversion) and about 26 km respectively. An extensional normal with strike-slip fault in about the north-south direction was formed by the first rupture, the mechanism of the second rupture is to be further studied. The method to inverse mechanisms of two point sources at the same time and the results obtained by directly analyzing P waveform records of the main shock are introduced, some related problems are also discussed. The Wuding earthquakes of October, 1995 and the Lijiang earthquake are considered to be the manifestation of the same dynamic process at different temporal and spatial points and the occurrence order of the two earthquakes is related to the direction of dynamics transmission.