复杂介质边界元法数值模拟研究

A major impetus to study the rough surface and complex structure in near surface model is because accuracy of seismic observation and geophysical prospecting can be improved. Wave theory study about fluid-satuated porous media has important significance for some scientific problems, such as explore underground resources, study of earth's internal structure, and structure response of multi-phase porous soil under dynamic and seismic effect. Seismic wave numerical modeling is one of the effective methods which understand seismic propagation rules in complex media. As a numerical simulation method, boundary element methods had been widely used in seismic wave field study. This paper mainly studies randomly rough surface scattering which used some approximation solutions based on boundary element method. In addition, I developed a boundary element solution for fluid saturated porous media. In this paper, we used boundary element methods which based on integral expression of wave equation to study the free rough surface scattering effects of Kirchhoff approximation method, Perturbation approximation method, Rytov approximation method and Born series approximation method. Gaussian spectrum model of randomly rough surfaces was chosen as the benchmark model. The approximation methods result were compared with exact results which obtained by boundary element methods, we study that the above approximation methods were applicable how rough surfaces and it is founded that this depends on and ( here is the wavenumber of the incident field, is the RMS height and is the surface correlation length ). In general, Kirchhoff approximation which ignores multiple scatterings between any two surface points has been considered valid for the large-scale roughness components. Perturbation theory based on Taylor series expansion is valid for the small-scale roughness components, as and are .Tests with the Gaussian topographies show that the Rytov approximation methods improves the Kirchhoff approximation in both amplitude and phase but at the cost of an extra treatment of transformation for the wave fields. The realistic methods for the multiscale surfaces come with the Born series approximation and the second-order Born series approximation might be sufficient to guarantee the accuracy of randomly rough surfaces. It could be an appropriate choice that a complex rough surface can be divided into large-, medium-, and small-scale roughness components with their scattering features be studied by the Kirchhoff or Rytov phase approximations, the Born series approximation, and the perturbation theory, respectively. For this purpose, it is important to select appropriate parameters that separate these different scale roughness components to guarantee the divided surfaces satisfy the physical assumptions of the used approximations, respectively. In addition, in this paper, the boundary element methods are used for solving the porous elastic wave propagation and carry out the numerical simulation. Based on the fluid-saturated porous model, this paper analyses and presents the dynamic equation of elastic wave propagation and boundary integral equation formulation of fluid saturated porous media in frequency domain. The fundamental solutions of the elastic wave equations are obtained according to the similarity between thermoelasticity and poroelasticity. At last, the numerical simulation of the elastic wave propagation in the two-phase isotropic media is carried out by using the boundary element method. The results show that a slow quasi P-wave can be seen in both solid and fluid wave-field synthetic seismograms. The boundary element method is effective and feasible.

波动方程辛几何算法三维叠前深度偏移及其并行实现

Oil and scientific groups have been focusing on the 3D wave equation prestack depth migration since it can solve the complex problems of the geologic structure accurately and maintain the wave information, which is propitious to lithology imaging. The symplectic method was brought up by Feng Kang firstly in 1984 and became the hotspot of numerical computation study. It will be widely applied in many scientific field of necessity because of its great virtue in scientific sense. This paper combines the Symplectic method and the 3-D wave equation prestack depth migration to bring up an effectual numerical computation method of wave field extrapolatation technique under the scientific background mentioned above. At the base of deep analysis of computation method and the performance of PC cluster, a seismic prestack depth migration flow considering the virtue of both seismic migration method and Pc cluster has formatted. The software, named 3D Wave Equation Prestack Depth Migration of Symplectic Method, which is based on the flow, has been enrolled in the National Bureau of Copyright (No. 0013767). Dagang and Daqing Oil Field have now put it into use in the field data processing. In this paper, the one way wave equation operator is decompounded into a phase shift operator and a time shift operator and the correct item with high rank Symplectic method when approaching E exponent. After reviewing eliminating alias frequency of operator, computing the maximum angle of migration and the imaging condition, we present the test result of impulse response of the Symplectic method. Taking the imaging results of the SEG/EAGE salt and overthrust models for example and seeing about the imaging ability with complex geologic structure of our software system, the paper has discussed the effect of the selection of imaging parameters and the effectuation on the migration result of the seismic wavelet and compared the 2-D and 3-D prestack depth migration result of the salt mode. We also present the test result of impulse response with the overthrust model. The imaging result of the two international models indicates that the Symplectic method of 3-D prestack depth migration accommodates great transversal velocity variation and complex geologic structure. The huge computing cost is the key obstruction that 3-D prestack depth migration wave equation cannot be adopted by oil industry. After deep analysis of prestack depth migration flow and the character of PC cluster ,the paper put forward :i)parallel algorithms in shot and frequency domain of the common shot gather 3-D wave equation prestack migration; ii)the optimized setting scheme of breakpoint in field data processing; iii)dynamic and static load balance among the nodes of the PC cluster in the 3-D prestack depth migration. It has been proven that computation periods of the 3-D prestack depth migration imaging are greatly shortened given that adopting the computing method mentioned in the paper. In addition,considering the 3-D wave equation prestack depth migration flow in complex medium and examples of the field data processing, the paper put the emphasis on: i)seismic data relative preprocessing, ii) 2.5D prestack depth migration velocity analysis, iii)3D prestack depth migration. The result of field data processing shows satisfied application ability of the flow put forward in the paper.

双聚焦偏移速度估计与可分表示叠前弹性反演

In exploration seismology, the geologic target of oil and gas reservoir in complex medium request the high accuracy image of the structure and lithology of the medium. So the study of the prestack image and the elastic inversion of seismic wave in the complex medium come to the leading edge. The seismic response measured at the surface carries two fundamental pieces of information: the propagation effects of the medium and the reflections from the different layer boundaries in the medium. The propagation represent the low-wavenumber component of the medium, it is so-called the trend or macro layering, whereas the reflections represent the high-wavenumber component of the medium, it is called the detailed or fine layering. The result of migration velocity analysis is the resolution of the low-wavenumber component of the medium, but the prestack elastic inversion provided the resolution of the high-wavvenumber component the medium. In the dissertation, the two aspects about the migration velocity estimation and the elastic inversion have been studied.Firstly, any migration velocity analysis methods must include two basic elements: the criterion that tell us how to know whether the model parameters are correct and the updating that tell us how to update the model parameters when they are incorrect, which are effected on the properties and efficiency of the velocity estimation method. In the dissertation, a migration velocity analysis method based on the CFP technology has been presented in which the strategy of the top-down layer stripping approach are adapted to avoid the difficult of the selecting reduce .The proposed method has a advantage that the travel time errors obtained from the DTS panel are defined directly in time which is the difference with the method based on common image gather in which the residual curvature measured in depth should be converted to travel time errors.In the proposed migration velocity analysis method, the four aspects have been improved as follow:? The new parameterization of velocity model is provided in which the boundaries of layers are interpolated with the cubic spline of the control location and the velocity with a layer may change along with lateral position but the value is calculated as a segmented linear function of the velocity of the lateral control points. The proposed parameterization is suitable to updating procedure.? The analytical formulas to represent the travel time errors and the model parameters updates in the t-p domain are derived under local lateral homogeneous. The velocity estimations are iteratively computed as parametric inversion. The zero differential time shift in the DTS panel for each layer show the convergence of the velocity estimation.? The method of building initial model using the priori information is provided to improve the efficiency of velocity analysis. In the proposed method, Picking interesting events in the stacked section to define the boundaries of the layers and the results of conventional velocity analysis are used to define the velocity value of the layers? An interactive integrate software environment with the migration velocity analysis and prestack migration is built.The proposed method is firstly used to the synthetic data. The results of velocity estimation show both properties and efficiency of the velocity estimation are very good.The proposed method is also used to the field data which is the marine data set. In this example, the prestack and poststack depth migration of the data are completed using the different velocity models built with different method. The comparison between them shows that the model from the proposed method is better and improves obviously the quality of migration.In terms of the theoretical method of expressing a multi-variable function by products of single-variable functions which is suggested by Song Jian (2001), the separable expression of one-way wave operator has been studied. A optimization approximation with separable expression of the one-way wave operator is presented which easily deal with the lateral change of velocity in space and wave number domain respectively and has good approach accuracy. A new prestack depth migration algorithm based on the optimization approximation separable expression is developed and used to testing the results of velocity estimation.Secondly, according to the theory of the seismic wave reflection and transmission, the change of the amplitude via the incident angle is related to the elasticity of medium in the subsurface two-side. In the conventional inversion with poststack datum, only the information of the reflection operator at the zero incident angles can be used. If the more robust resolutions are requested, the amplitudes of all incident angles should be used.A natural separable expression of the reflection/transmission operator is represented, which is the sum of the products of two group functions. One group function vary with phase space whereas other group function is related to elastic parameters of the medium and geological structure.By employing the natural separable expression of the reflection/transmission operator, the method of seismic wave modeling with the one-way wave equation is developed to model the primary reflected waves, it is adapt to a certain extent heterogeneous media and confirms the accuracy of AVA of the reflections when the incident angle is less than 45'. The computational efficiency of the scheme is greatly high.The natural separable expression of the reflection/transmission operator is also used to construct prestack elastic inversion algorithm. Being different from the AVO analysis and inversion in which the angle gathers formed during the prstack migration are used, the proposed algorithm construct a linear equations during the prestack migration by the separable expression of the reflection/transmission operator. The unknowns of the linear equations are related to the elasticity of the medium, so the resolutions of them provided the elastic information of the medium.The proposed method of inversion is the same as AVO inversion in , the difference between them is only the method processing the amplitude via the incident angle and computational domain.

Optical recording properties of metal-azo dye as recording medium with super-resolution near-field structure

A novel metallized azo dye has been synthesized. The absorption spectra of the thin film and thermal characteristic are measured. Static optical recording properties with and without the Bi mask layer super-resolution near-field structure (Super-RENS) of the metal-azo dye are investigated. The results show that the metal-azo dye film has a broad absorbance band in the region of 450-650 nm and the maximum absorbance wavelength is located at 603 nm. It is also found that the new metallized azo dye occupies excellent thermal stability, initiatory decomposition temperature is at 270 degrees C and the mass loss is about 48% in a narrow temperature region (15 degrees C). The complex refractive index N (N = n + ik) is measured. High refractive index (n = 2.45) and low extinction coefficient (k = 0.2) at the recording wavelength 650nm are attained. Static optical recording tests with and without Super-RENS are carried out using a 650nm semiconductor diode laser with recording power of 7mW and laser pulse duration of 200ns. The AFM images show that the diameter of recording mark on the dye film with the Bi mask layer is reduced about 42%, compared to that of recorded mark on the dye film without Super-RENS. It is indicated that Bi can well performed as a mask layer of the dye recording layer and the metallized azo dye can be a promising candidate for recording media with the super-resolution near-field structure.

The diphasic nc-Si/a-Si : H thin film with improved medium-range order

A series of silicon film samples were prepared by plasma enhanced chemical vapor deposition (PECVD) near the threshold from amorphous to nanocrystalline state by adjusting the plasma parameters and properly increasing the reactions between the hydrogen plasma and the growing surface. The microstucture of the films was studied by micro-Raman and Fourier transform infrared (FTIR) spectroscopy. The influences of the hydrogen dilution ratio of silane (R-H = [H-2]/[SiH4]) and the substrate temperature (T-s) on the microstructural and photoelectronic properties of silicon films were investigated in detail. With the increase of RH from 10 to 100, a notable improvement in the medium-range order (MRO) of the films was observed, and then the phase transition from amorphous to nanocrystalline phase occurred, which lead to the formation of diatomic hydrogen complex, H-2* and their congeries. With the increase of T-s from 150 to 275 degreesC, both the short-range order and the medium range order of the silicon films are obviously improved. The photoconductivity spectra and the light induced changes of the films show that the diphasic nc-Si/a-Si:H films with fine medium-range order present a broader light spectral response range in the longer wavelength and a lower degradation upon illumination than conventional a-Si:H films. (C) 2004 Elsevier B.V. All rights reserved.

The diphasic nc-Si/a-Si : H thin film with improved medium-range order

A series of silicon film samples were prepared by plasma enhanced chemical vapor deposition (PECVD) near the threshold from amorphous to nanocrystalline state by adjusting the plasma parameters and properly increasing the reactions between the hydrogen plasma and the growing surface. The microstucture of the films was studied by micro-Raman and Fourier transform infrared (FTIR) spectroscopy. The influences of the hydrogen dilution ratio of silane (R-H = [H-2]/[SiH4]) and the substrate temperature (T-s) on the microstructural and photoelectronic properties of silicon films were investigated in detail. With the increase of RH from 10 to 100, a notable improvement in the medium-range order (MRO) of the films was observed, and then the phase transition from amorphous to nanocrystalline phase occurred, which lead to the formation of diatomic hydrogen complex, H-2* and their congeries. With the increase of T-s from 150 to 275 degreesC, both the short-range order and the medium range order of the silicon films are obviously improved. The photoconductivity spectra and the light induced changes of the films show that the diphasic nc-Si/a-Si:H films with fine medium-range order present a broader light spectral response range in the longer wavelength and a lower degradation upon illumination than conventional a-Si:H films. (C) 2004 Elsevier B.V. All rights reserved.

Solvent extraction study of rare earth elements from chloride medium by mixtures of sec-nonylphenoxy acetic acid with Cyanex301 or Cyanex302

The extraction of rare earth elements from chloride medium by mixtures of sec-nonylphenoxy acetic acid (CA100) with bis(2,4,4-trimethylpentyl) dithiophosphinic acid (Cyanex301) or bis(2,4,4-trimethylpentyl) monothiophosphinic acid (Cyanex302) in n-heptane has been studied. The synergistic enhancement of the extraction of lanthanum (III) by mixtures of CA100 with Cyanex301 has been investigated using the methods of slope analysis and constant mole. The extracted complex of lanthanum (III) is determined. The logarithm of the equilibrium constant is calculated as - 1.41. The formation constants and the thermodynamic functions, Delta H, Delta G, and Delta S have also been determined.

Solvent extraction of rare earths from chloride medium with mixtures of 1-phenyl-3-methyl-4-benzoyl-pyrazalone-5 and sec-octylphenoxyacetic acid

The extraction of rare earth elements from chloride medium by mixtures of sec-nonylphenoxy acetic acid (CA100) with bis(2,4,4-trimethylpentyl) dithiophosphinic acid (Cyanex301) or bis(2,4,4-trimethylpentyl) monothiophosphinic acid (Cyanex302) in n-heptane has been studied. The synergistic enhancement of the extraction of lanthanum (III) by mixtures of CA100 with Cyanex301 has been investigated using the methods of slope analysis and constant mole. The extracted complex of lanthanum (III) is determined. The logarithm of the equilibrium constant is calculated as - 1.41. The formation constants and the thermodynamic functions, Delta H, Delta G, and Delta S have also been determined.

Solvent extraction of lanthanides and yttrium from nitrate medium with CYANEX 925 in heptane

The extraction behavior of lanthanides and yttrium usinsg CYANEX 925 (mixture of branched chain alkylated phosphine oxides) in n-heptane from nitrate medium has been studied. The effects of aqueous phase ionic strength, CYANEX 925 concentration in the organic phase, and temperature on Sm3+, Nd3+ and Y3+ extraction have been investigated. The extractability of the lanthanides and yttrium increases with increasing nitrate concentration, as well as with increasing CYANEX 925 concentration. An extraction mechanism is proposed based on slope analysis. Furthermore, the infra-red spectra of CYANEX 925 saturated with lanthanides are employed to provide evidence of the composition of the complex. The relationship between the logarithm of the distribution ratio and lanthanide atomic number is also discussed which indicates that yttrium can be separated from fight lanthanides. In addition separation of the light and heavy lanthanide groups is also possible using CYANEX 925. From the temperature dependence data, the thermodynamic parameters values (Delta H, Delta S and Delta G) are calculated.

Extraction and stripping of ytterbium (III) from H2SO4 medium by Cyanex 923

The extraction and stripping of ytterbium (III) from sulfuric acid medium using Cyanex 923 in heptane solution was investigated. The effects of extractant concentration, pH and sulfate ion as well as stripping agents, acidity and temperature on the extraction and stripping were studied. The equilibrium constants and thermodynamic parameters, such as Delta H (10.76 kJ(.)mol(-1)), Delta G (-79.26 kJ(.)mol(-1)) and Delta S (292.41 J(.)K(-1.)mol(-1)), were calculated. The extraction mechanism and the complex species extracted were determined by slope analysis and FrIR spectra. Furthermore, it was found that the extraction of Yb (III) from sulfuric acid medium by Cyanex 923 increased with pH, concentration of SO42-, HSO4-, and extractant concentration, and approximately a quantitative extraction of Yb (III) was achieved at an equilibrium pH near 3.0, and the extracted complex was YbSO4(HSO4)(.)2Cyanex923((o)).

KINETIC-STUDY OF ELECTROCATALYTIC REDUCTION OF HYDROGEN-PEROXIDE AT IRON(III) PORPHYRIN MODIFIED GLASSY-CARBON ELECTRODE IN ALKALINE-MEDIUM USING THE ROTATION-SCAN METHOD

Reduction of hydrogen peroxide at a glassy carbon (GC) electrode modified with sigma-bonded pyrrole iron(III) octaethylporphyrin complex, (OEP)Fe(Pyr), was studied by cyclic voltammetry and a rotating disk electrode. In 0.1N NaOH solution, it is shown that such an (OEP)Fe(Pyr)/GC electrode has a significant catalytic activity towards hydrogen peroxide reduction (E(D) = -0.80 V, k = 0.066 cm s(-1)); however, the electrode stability is low. The deactivation is observed when the reaction charge (Q) is passing through the (OEP)Fe(Pyr)/GC disk electrode. A linear rotation scan method is applied to study the kinetic process by determining the disk electrochemical response (i(D)) to rotation rate (omega) at a definite disk potential (E(D)). Considering that the number of adsorbed electroreduced catalyst molecules (Red) varies according to the disk potential, a factor theta(= Gamma(Red)/(Gamma(Red) + Gamma(Ox))) is introduced to describe the electrode surface area fraction for electroreduced species. The obtained Koutecky-Levich equation is applicable whatever the potential is.