改进叠前深度偏移的计算效率＆叠前高分辨率成像

An high-resolution prestack imaging technique of seismic data is developed in this thesis. By using this technique, the reflected coefficients of sheet sands can be gained in order to understand and identify thin oil reservoirs. One-way wave equation based migration methods can more accurately model seismic wave propagation effect such as multi-arrivals and obtain almost correct reflected energy in the presence of complex inhomogeneous media, and therefore, achieve more superiorities in imaging complex structure. So it is a good choice to apply the proposed high-resolution imaging to the presatck depth migration gathers. But one of the main shorting of one-way wave equation based migration methods is the low computational efficiency, thus the improvement on computational efficiency is first carried out. The method to improve the computational efficiency of prestack depth migration is first presented in this thesis, that is frequency-dependent varying-step depth exploration scheme plus a table-driven, one-point wavefield interpolation technology for wave equation based migration methods; The frequency-dependent varying-step depth exploration scheme reduces the computational cost of wavefield depth extrapolation, and the a table-driven, one-point wavefield interpolation technology reconstructs the extrapolated wavefield with an equal, desired vertical step with high computational efficiency. The proposed varying-step depth extrapolation plus one-point interpolation scheme results in 2/3 reduction in computational cost when compared to the equal-step depth extrapolation of wavefield, but gives the almost same imaging. The frequency-dependent varying-step depth exploration scheme is presented in theory by using the optimum split-step Fourier. But the proposed scheme can also be used by other wave equation based migration methods of the frequency domain. The proposed method is demonstrated by using impulse response, 2-D Marmousi dataset, 3-D salt dataset and the 3-D field dataset. A method of high-resolution prestack imaging is presented in the 2nd part of this thesis. The seismic interference method to solve the relative reflected coefficients is presented. The high-resolution imaging is obtained by introducing a sparseness- constrained least-square inversion into the reflected coefficient imaging. Gaussian regularization is first imposed and a smoothed solution is obtained by solving equation derived from the least-square inversion. Then the Cauchy regularization is introducing to the least-square inversion , the sparse solution of relative reflected coefficients can be obtained, that is high-resolution solution. The proposed scheme can be used together with other prestack imaging if the higher resolution is needed in a target zone. The seismic interference method in theory and the solution to sparseness-constrained least-square inversion are presented. The proposed method is demonstrated by synthetic examples and filed data.

Fourier有限差分法深度偏移与正演模拟

With the development of seismic exploration, the target becomes more and more complex, which leads to a higher demand for the accuracy and efficiency in 3D exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which keeps the ability of finite-differenc method in dealing with laterally varing media and inherits the predominance of the phase-screen method in stablility and efficiency. In this thesis, the accuracy of the FFD operator is highly improved by using simulated annealing algorithm. This method takes the extrapolation step and band width into account, which is more suitable to various band width and discrete scale than the commonely-used optimized method based on velocity contrast alone. In this thesis, the FFD method is extended to viscoacoustic modeling. Based on one-way wave equation, the presented method is implemented in frequency domain; thus, it is more efficient than two-way methods, and is more convenient than time domain methods in handling attenuation and dispersion effects. The proposed method can handle large velocity contrast and has a high efficiency, which is helpful to further research on earth absorption and seismic resolution. Starting from the frequency dispersion of the acoustic VTI wave equation, this thesis extends the FFD migration method to the acoustic VTI media. Compared with the convetional FFD method, the presented method has a similar computational efficiency, and keeps the abilities of dealing with large velocity contrasts and steep dips. The numerical experiments based on the SEG salt model show that the presented method is a practical migration method for complex acoustical VTI media, because it can handle both large velocity contrasts and large anisotropy variations, and its accuracy is relatively high even in strong anisotropic media. In 3D case, the two-way splitting technique of FFD operator causes artificial azimuthal anisotropy. These artifacts become apparent with increasing dip angles and velocity contrasts, which prevent the application of the FFD method in 3D complex media. The current methods proposed to reduce the azimuthal anisotropy significantly increase the computational cost. In this thesis, the alternating-direction-implicit plus interpolation scheme is incorporated into the 3D FFD method to reduce the azimuthal anisotropy. By subtly utilizing the Fourier based scheme of the FFD method, the improved fast algorithm takes approximately no extra computation time. The resulting operator keeps both the accuracy and the efficiency of the FFD method, which is helpful to the inhancements of both the accuracy and the efficiency for prestack depth migration. The general comparison is presented between the FFD operator and the generalized-screen operator, which is valuable to choose the suitable method in practice. The percentage relative error curves and migration impulse responses show that the generalized-screen operator is much sensiutive to the velocity contrasts than the FFD operator. The FFD operator can handle various velocity contrasts, while the generalized-screen operator can only handle some range of the velocity contrasts. Both in large and weak velocity contrasts, the higher order term of the generalized-screen operator has little effect on improving accuracy. The FFD operator is more suitable to large velocity contrasts, while the generalized-screen operator is more suitable to middle velocity contrasts. Both the one-way implicit finite-difference migration and the two-way explicit finite-differenc modeling have been implemented, and then they are compared with the corresponding FFD methods respectively. This work gives a reference to the choosen of proper method. The FFD migration is illustrated to be more attractive in accuracy, efficiency and frequency dispertion than the widely-used implicit finite-difference migration. The FFD modeling can handle relatively coarse grids than the commonly-used explicit finite-differenc modeling, thus it is much faster in 3D modeling, especially for large-scale complex media.

太阳辐射对热层和电离层变化性的影响

Solar ultraviolet (UV) radiation at wavelengths less than 400 nm is an important source of energy for aeronomic processes throughout the solar system. Solar UV photons are absorbed in planetary atmospheres, as well as throughout the heliosphere, via photodissociation of molecules, photoionization of molecules and atoms, and photoexcitation toexcitation including resonance scattering. In this paper, the solar irradiances data measured by TIMED SEE, as well as the solar proxies such as F10.7 and Mg II, thermosphere neutral density of CHAMP measurements and topside ionospheric plasmas densities from DMSP, are used to analyze solar irradiance effects on the variabilities of the thermosphere and the ionosphere. First, thermosphere densities near 410 km altitude are analyzed for solar irradiance variability effects during the period 2002-2004. Correlations between the densities and the solar irradiances for different spectral lines and wavelength ranges reveal significantly different characteristics. The density correlates remarkably well with all the selected solar irradiances except the lower chromospheric O I (130.4 nm) emission. Among the chosen solar proxies, the Mg II core-to-wing ratio index, EUV (30-120 nm) and F10.7 show the highest correlations with the density for short-term (< ~27 days) variations. For both long- (> ~27 days) and short-term variations, linear correlation coefficients exhibit a decreasing trend from low latitudes towards high latitudes. The density variability can be effectively modeled (capturing 71% of the variance) using multiple solar irradiance indices, including F10.7, SEUV (the EUV 30-120 nm index), and SFUV (the FUV 120-193 nm index), in which a lag time of 1 day was used for both F10.7 and SEUV, and 5 days for SFUV. In our regression formulation SEUV has the largest contribution to the density variation (40%), with the F10.7 having the next largest contribution (32%) and SFUV accounting for the rest (28%). Furthermore, a pronounced period of about 27.2 days (mean period of the Sun's rotation) is present in both density and solar irradiance data of 2003 and 2004, and a pronounced period of about 54.4 days (doubled period of the solar rotation) is also revealed in 2004. However, soft X-ray and FUV irradiances did not present a pronounced 54.4 day period in 2004, in spite of their high correlation with the densities. The Ap index also shows 54-day periodicities in 2004, and magnetic activity, together with solar irradiance, affects the 54-day variation in density significantly. In addition, NRLMSISE00, DTM-2000 and JB2006 model predictions are compared with density measurements from CHAMP to assess their accuracy, and the results show that these models underestimate the response of the thermosphere to variations induced by solar rotation. Next, the equatorial topside ionospheric plasmas densities Ni are analyzed for solar irradiance variability effects during the period 2002-2005. Linear correlations between Ni and the solar irradiances for different wavelength ranges reveal significantly different characteristics. XUV (0-35 nm) and EUV (115-130 nm) show higher correlation with Ni for the long-term variations, whereas EUV (35-115 nm) show higher correlation for the short-term variations. Moreover, partial correlation analysis shows that the long-term variations of Ni are affected by both XUV (0-35 nm) and EUV (35-115 nm), whereas XUV (0-35 nm) play a more important role; the short-term variations of Ni are mostly affected by EUV (35-115 nm). Furthermore, a pronounced period of about 27 days is present in both Ni and solar irradiance data of 2003 and 2004, and a pronounced period of about 54 days is also revealed in 2004. Finally, prompted by previous studies that have suggested solar EUV radiation as a means of driving the semiannual variation, we investigate the intra-annual variation in thermosphere neutral density near 400 km during 2002-2005. The intra-annual variation, commonly referred to as the ‘semiannual variation’, is characterized by significant latitude structure, hemispheric asymmetries, and inter-annual variability. The magnitude of the maximum yearly difference, from the yearly minimum to the yearly maximum, varies by as much as 60% from year to year, and the phases of the minima and maxima also change by 20-40 days from year to year. Each annual harmonic of the intra-annual variation, namely, annual, semiannual, ter-annual and quatra-annual, exhibits a decreasing trend from 2002 through 2005 that is correlated with the decline in solar activity. In addition, some variations in these harmonics are correlated with geomagnetic activity, as represented by the daily mean value of Kp. Recent empirical models of the thermosphere are found to be deficient in capturing most of the latitude dependencies discovered in our data. In addition, the solar flux and geomagnetic activity proxies that we have employed do not capture some latitude and inter-annual variations detected in our data. It is possible that these variations are partly due to other effects, such as seasonal-latitudinal variations in turbopause altitude (and hence O/N2 composition) and ionosphere coupling processes that remain to be discovered in the context of influencing the intra-annual variations depicted here. Our results provide a new dataset to challenge and validate thermosphere-ionosphere general circulation models that seek to delineate the thermosphere intra-annual variation and to understand the various competing mechanisms that may contribute to its existence and variability. We furthermore suggest that the term “intra-annual” variation be adopted to describe the variability in thermosphere and ionosphere parameters that is well-captured through a superposition of annual, semiannual, ter-annual, and quatra-annual harmonic terms, and that “semiannual’ be used strictly in reference to a pure 6-monthly sinusoidal variation. Moreover, we propose the term “intra-seasonal” to refer to those shorter-term variations that arise as residuals from the above Fourier representation.

Estimation of soil organic partition coefficients: from retention factors measured by soil column chromatography with water as eluent

The retention factors (k) of 104 hydrophobic organic chemicals (HOCs) were measured in soil column chromatography (SCC) over columns filled with three naturally occurring reference soils and eluted with Milli-Q water. A novel method for the estimation of soil organic partition coefficient (K-oc) was developed based on correlations with k in soil/water systems. Strong log K-oc versus log k correlations (r>0.96) were found. The estimated K-oc values were in accordance with the literature values with a maximum deviation of less than 0.4 log units. All estimated K-oc values from three soils were consistent with each other. The SCC approach is promising for fast screening of a large number of chemicals in their environmental applications. (C) 2002 Elsevier Science B.V. All rights reserved.

Prediction of soil adsorption coefficients from retention parameters on three reversed-phase liquid chromatographic columns

Reversed-phase high-performance liquid chromatographic (RP-HPLC) retention parameters, which are determined by the intermolecular interactions in retention process, can be considered as the chemical molecular descriptors in linear free energy relationships (LFERs). On the basis of the characterization and comparison of octadecyl-bonded silica gel (ODS), cyano-bonded silica gel (CN), and phenyl-bonded silica gel (Ph) columns with linear solvation energy relationships (LSERs), a new multiple linear regression model using RP-HPLC retention parameters on ODS and CN columns as variables for estimation of soil adsorption coefficients was developed. It was tested on a set of reference substances from various chemical classes. The results showed that the multicolumn method was more promising than a single-column method was for the estimation of soil adsorption coefficients. The accuracy of the suggested model is identical with that of LSERs.

A column method for determination of soil organic partition coefficients of eight pesticides

A column method was developed to conveniently and reliably determine the soil organic partition coefficients (K-oc) of three insecticides (methiocarb, azinphos-methyl, fenthion), four fungicides (triadimenol, fuberidazole, tebuconazole, pencycuron), and one herbicide (atrazine), in which real soil acted as a stationary phase and the water solution of pesticide as an eluent. The processes of sorption equilibrium were directly shown through a breakthrough curve(BTC). The log K-oc values are 1.69, 1.95, 2.25, 2.55, 2.69, 2.67, 3.10, and 3.33 for atrazine, triadimenol, methiocarb, fuberidazole, azinphos-methyl, tebuconazole, fenthion and pencycuron, respectively.

Soil column chromatography for correlation between capacity factors and soil organic partition coefficients for eight pesticides

A soil column chromatographic method was developed to measure the capacity factors (k') of pesticides, in which soil acted as a stationary phase and methanol-water mixture as an eluent. The k' values of eight pesticides, including three insecticides (methiocarb, azinphos-methyl, fenthion), four fungicides (triadimenol, fuberidazole, tebuconazole, pencycuron), and one herbicide (atrazine), were found to be well fitted to a retention equation, ln k'=ln k(w)'-S-phi. Due to similar interactions of solutes with soil and solvent in both sorption determination and retention experiment, log k' has a good linear correlation with log K-oc for the eight pesticides from different classes, in contrast with poor correlation between log k' from C-18 column and log K-oc. So the method provides a tool for rapid estimation of K-oc from experimental k'. (C) 1999 Elsevier Science Ltd. All rights reserved.

A new method for determination of adsorption and desorption coefficients of pesticides with soil column liquid chromatography

The adsorption and desorption coefficients of atrazine, methiocarb and simazine on a sandy loam soil were measured in this study with soil column liquid chromatographic (SCLC) technique. The adsorption and desorption data of all the three pesticides followed Freundlich isotherms revealing the existence of hysteresis. In comparing with other methods, SCLC method showed some characteristics such as rapidity, online and accuracy.

Prediction of soil organic partition coefficients by a soil leaching column chromatographic method

The soil organic partition coefficient (K-oc) is one of the most important parameters to depict the transfer and fate of a chemical in the soil-water system. Predicting K-oc by using a chromatographic technique has been developing into a convenient and low-cost method. In this paper, a soil leaching column chromatograpy (SLCC) method employing the soil column packed with reference soil GSE 17201 (obtained from Bayer Landwirtschaftszentrum, Monheim, Germany) and methanol-water eluents was developed to predict the K-oc of hydrophobic organic chemicals (HOCs), over a log K-oc range of 4.8 orders of magnitude, from their capacity factors. The capacity factor with water as an eluent (k(w)') could be obtained by linearly extrapolating capacity factors in methanol-water eluents (k') with various volume fractions of methanol (phi). The important effects of solute activity coefficients in water on k(w)' and K-oc were illustrated. Hence, the correlation between log K-oc and log k(w)' (and log k') exists in the soil. The correlation coefficient (r) of the log K-oc vs. log k(w)' correlation for 58 apolar and polar compounds could reach 0.987, while the correlation coefficients of the log K-oc-log k' correlations were no less than 0.968, with phi ranging from 0 to 0.50. The smaller the phi, the higher the r. Therefore, it is recommended that the eluent of smaller phi, such as water, be used for accurately estimating K-oc. Correspondingly, the r value of the log K-oc-log k(w)' correlation on a reversed-phase Hypersil ODS (Thermo Hypersil, Kleinostheim, Germany) column was less than 0.940 for the same solutes. The SLCC method could provide a more reliable route to predict K-oc indirectly from a correlation with k(w)' than the reversed-phase liquid chromatographic (RPLC) one.