超参数控制下的波阻抗约束反演

The seismic survey is the most effective geophysical method during exploration and development of oil/gas. As a main means in processing and interpreting seismic data, impedance inversion takes up a special position in seismic survey. This is because the impedance parameter is a ligament which connects seismic data with well-logging and geological information, while it is also essential in predicting reservoir properties and sand-body. In fact, the result of traditional impedance inversion is not ideal. This is because the mathematical inverse problem of impedance is poor-pose so that the inverse result has instability and multi-result, so it is necessary to introduce regularization. Most simple regularizations are presented in existent literature, there is a premise that the image(or model) is globally smooth. In fact, as an actual geological model, it not only has made of smooth region but also be separated by the obvious edge, the edge is very important attribute of geological model. It's difficult to preserve these characteristics of the model and to avoid an edge too smooth to clear. Thereby, in this paper, we propose a impedance inverse method controlled by hyperparameters with edge-preserving regularization, the inverse convergence speed and result would be improved. In order to preserve the edge, the potential function of regularization should satisfy nine conditions such as basic assumptions edge preservation and convergence assumptions etc. Eventually, a model with clear background and edge-abnormity can be acquired. The several potential functions and the corresponding weight functions are presented in this paper. The potential functionφLφHL andφGM can meet the need of inverse precision by calculating the models. For the local constant planar and quadric models, we respectively present the neighborhood system of Markov random field corresponding to the regularization term. We linearity nonlinear regularization by using half-quadratic regularization, it not only preserve the edge, and but also simplify the inversion, and can use some linear methods. We introduced two regularization parameters (or hyperparameters) λ2 and δ in the regularization term. λ2 is used to balance the influence between the data term and the transcendental term; δ is a calibrating parameter used to adjust the gradient value at the discontinuous position(or formation interface). Meanwhile, in the inverse procedure, it is important to select the initial value of hyperparameters and to change hyperparameters, these will then have influence on convergence speed and inverse effect. In this paper, we roughly give the initial value of hyperparameters by using a trend- curve of φ-(λ2, δ) and by a method of calculating the upper limit value of hyperparameters. At one time, we change hyperparameters by using a certain coefficient or Maximum Likelihood method, this can be simultaneously fulfilled with the inverse procedure. Actually, we used the Fast Simulated Annealing algorithm in the inverse procedure. This method overcame restrictions from the local extremum without depending on the initial value, and got a global optimal result. Meanwhile, we expound in detail the convergence condition of FSA, the metropolis receiving probability form Metropolis-Hasting, the thermal procession based on the Gibbs sample and other methods integrated with FSA. These content can help us to understand and improve FSA. Through calculating in the theoretic model and applying it to the field data, it is proved that the impedance inverse method in this paper has the advantage of high precision practicability and obvious effect.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

地震反射界面形态与介质各向异性属性联立成像：地震反射波旅行时与偏振信息综合反演

It has been a difficult problem faced by seismologists for long time that how exactly to reconstruct the earth's geometric structure and distribution of physical attributes according to seismic wave's kinematical and dynamic characteristics, obtained in seismological observation. The jointing imaging of seismic reflector and anisotropy attributes in the earth interior is becoming the research hot spot. The limitation of shoot and observation system makes that the obtained seismic data are too scarce to exactly reconstruct the geological objects. It is popular that utilizing only seismic reflection traveltimes or polarizations information make inversion of the earth's velocity distribution by fixing seismic reflector configuration (vice versa), these will lead to the serious non-uniqueness reconstruction due to short of effective data, the non-uniqueness problem of reconstructing anisotropy attributes will be more serious than in isotropy media. Obviously it is not enough to restrict the media structure only by information of seismic reflection traveltimes or polarizations, which even sometimes will lead to distorted images and misinterpretation of subsurface structure. So we try to rebuild seismic reflection structure (geometry) and media anisotropic structure (physics) in the earth interior by jointing data of seismic wave kinematics and dynamics characteristics, we carry out the new experiment step by step, and the research mainly comprises of two parts: one is the reconstruction of P-wave vertical velocity and anisotropic structure(Thomsen parameter s and 8) in the transversely isotropic media with vertical symmetrical axis(VTI) by fixing geometrical structure, and the other is the simultaneous inversion of the reflector surface conformation and seismic anisotropic structure by jointing seismic reflection traveltimes and polarizations data. Simulated annealing method is used to the first research part, linear inversion based on BG theory and Simulated annealing are applied to the second one. All the research methods are checked by model experiments, then applied to the real data of the wide-angle seismic profile from Tunxi, Anhui Province, to Wenzhou, Zhejiang Province. The results are as following The inversion results based on jointing seismic PP-wave or PSV-wavereflection traveltimes and polarizations data are more close to real model than themodels based simply on one of the two data respectively. It is shown that the methodwe present here can effectively reconstruct the anisotropy attributes in the earth'sinterior when seismic reflector structure is fixed.The layer thickness, P-wave vertical velocity and Thomsen anisotropicparameters {s and 8) could be resolved simultaneously by jointing inversion ofseismic reflection traveltimes and polarizations with the linear inversion methodbased on BG theory.The image of the reflector structure, P-wave vertical velocity and theanisotropy parameters in the crust could be obtained from the wide-angle seismicprofile from Tunxi (in Anhui Province), to Wenzhou (in Zhejiang Province). Theresults reveal the difference of the reflector geometrical structure and physicalattributes in the crust between Yangtze block and Cathaysia block, and attempt tounderstand the characteristics of the crustal stress field in the areas.

Effect of Cold Rolling Ratio and Heating Rate of Batch Annealing on Textures of Micro-Carbon Steel Sheet

在四辊冷轧试验机和Gleeble-1500试验机上进行了热轧微碳钢板的冷轧和退火试验。用D/max-RC衍射仪测量了试样的,/”层织构,并用Roe软件进行了ODF分析。研究表明,所研究的热轧微碳深冲板压下率约为75%,退火升温速度为20-40℃/h时,试样为{111}织构特征；压下率较大(80%)时,退火织构为较弱的{111}组分。无论{111}织构还是非{111}织构都是在形核阶段开始形成,在晶粒长大优先长大,受到定向形核和选择生长双重机制的作用。

Stability of ZrTiCuNiBe Bulk Metallic Glass Upon Isothermal Annealing Near the Glass Transition Temperature

The stability of Zr41Ti14Cu12.5Ni10Be22.5 bulk metallic glass (BMG) upon isothermal annealing near the glass transition temperature has been investigated by using x-ray diffraction, differential scanning calorimetry, and the pulse echo overlap method. The density, elastic constants, and thermodynamic parameters as well as their annealing time dependence have been determined. The microstructural and properties changes of the annealed BMG were checked by acoustic measurement. Obvious structural and property changes were observed with prolonged annealing of the BMG near the glass transition temperature.

Work Softening And Annealing Hardening Of Deformed Nanocrystalline Nickel

We reported that work softening takes place during room-temperature rolling of nanocrystalline Ni at an equivalent strain of around 0.30. The work softening corresponds to a strain-induced phase transformation from a face-centered cubic (fcc) to a body-centered cubic (bcc) lattice. The hardness decreases with increasing volume fraction of the bcc phase. When the deformed samples are annealed at 423 K, a hardening of the samples takes place. This hardening by annealing can be attributed to a variety of factors including the recovery transformation from the bcc to the fcc phase, grain boundary relaxation, and retardation of dislocation gliding by microtwins.

Simulated pathogenic conformational switch regions matched well with the biochemical findings

Pathogenic conformational conversion is a general causation of many disease, such as transmissible spon- giform encephalopathy (TSE) caused by misfolding of prion, sickle cell anemia, and etc. In such structural changes, misfolding occurs in regions important for the stability of native structure firstly. This destabi- lizes the normal conformation and leads to subsequent errors in folding pathway. Sites involved in the first stage can be deemed switch regions of the protein, and are vital for conformational conversion. Namely it could be a switch of disease at residue level. Here we report an algorithm that can identify such sites computationally with an accuracy of 93%, by calculating the probability of the native structure of a short segment jumping to a mistake one. Knowledge of such switch sites could be used to target clinical therapy, study physiological and pathologic mechanism of protein, and etc.

Influence Of Gravity On Structure Of Colloidal Crystal Using Simulated Microgravity

Liquid mixtures of water and deuterium oxide as the liquid phase, were used to match the density of charged colloidal particles. Kossel diffraction method was used to detect the crystal structures. The experiments under the density-matched (g=0) and unmatched (g=1) conditions are compared to examine the influence of gravity on the crystal structures formed by self-assembly of 110 nm (in diameter) polystyrene microspheres. The result shows that die gravity tends to make the lattice constants of colloidal crystals smaller at lower positions, which indicates that the effect of gravity should be taken into account in the study of the colloidal crystals.

Influence of thermal annealing on optical properties and surface morphology of NiOx thin films

NiOx thin films were deposited by reactive DC-magnetron sputtering from a nickel metal target in Ar + O-2 with the relative O-2 content of 5%. Thermal annealing effects on optical properties and surface morphology of NiOx, films were investigated by X-ray photoelectron spectroscopy, thermogravimetric analysis, scanning electron microscope and optical measurement. The results showed that the changes in optical properties and surface morphology depended on the temperature. The surface morphology of the films changed obviously as the annealing temperature increased due to the reaction NiOx -> NiO + O-2 releasing O-2. The surface morphology change was responsible for the variation of the optical properties of the films. The optical contrast between the as-deposited films and 400 degrees C annealed films was about 52%. In addition, the relationship of the optical energy band gap with the variation of annealing temperature was studied. (c) 2006 Elsevier B.V. All rights reserved.

Effects of annealing on the color, absorption spectra, and light yield of Ce : YAlO3 single crystal grown by the temperature gradient technique

YAlO3 single crystal doped with Ce3+ at concentration 1% was grown by the temperature gradient technique. The as-grown crystal was pink. After H-2 annealing or air annealing at 1400degreesC for 20 h, the crystal was turned into colorless. We concluded there were two kinds of color centers in the as-grown crystal. One is F+ center attributed to absorption band peaking at about 530 nm, the other is O- center attributed to absorption band peaking at about 390 nm. This color centers model can be applied in explaining the experiment phenomena including the color changes, the absorption spectra changes, and the light yield changes of Ce:YAP crystals before and after annealing. (C) 2004 American Institute of Physics.

Effects of sapphire substrate annealing on ZnO epitaxial films grown by MOCVD

The annealing effects of sapphire substrate on the quality of epitaxial ZnO films grown by metalorganic chemical vapor deposition (MOCVD) were studied. The atomic steps formed on (0 0 0 1) sapphire (alpha-Al2O3) substrate surface by annealing at high temperature was analyzed by atomic force microscopy (AFM). The annealing effects of sapphire substrate on the ZnO films were examined by X-ray diffraction (XRD), AFM and photoluminescence (PL) measurements. Experimental results indicate that the film quality is strongly affected by annealing treatment of the sapphire substrate surface., The optimum annealing temperature of sapphire substrates is given. (c) 2006 Elsevier B.V. All rights reserved.

Annealing effect on properties of ZnO thin films grown on LiNbO3 substrates by MOCVD

ZnO films were grown on (0 0 0 1) LiNbO3 substrates by metal organic chemical vapor deposition (MOCVD). Annealing of ZnO films was performed in air for I h at 800 degrees C. The effects of annealing on the structural and optical properties of ZnO thin films on LiNbO3 substrates were investigated using X-ray diffraction (XRD), atomic force microscopy (AFM) and photoluminescence (PL) measurements. XRD patterns and AFM showed that the as-grown and the annealed ZnO films grown on LiNbO3 substrates had c-axis preferential orientation, the crystallinity of the ZnO films grown on LiNbO3 Substrates was improved, and the grain size increased by annealing. The PL spectra showed that the intensity of the UV near-band-edge peak was increased after annealing, while the intensity of visible peak (deep-level emission) decreased. (c) 2005 Elsevier B.V. All rights reserved.

Effects of gamma irradiation and air annealing on Yb-doped YAlO3 single crystals

The effects of gamma irradiation on as-grown 5 at% Yb:YAlO3 (YAP) and air annealing on gamma-irradiated 5 at% Yb: YAP have been studied by the difference in the absorption spectra before and after treatment. The gamma irradiation and air annealing led to opposite changes of the absorption properties of the Yb: YAP crystal. After air annealing, the gamma-irradiation effects were totally removed over the wavelength range 390-800 nm and the concentrations of Fe3+ and Yb3+ were slightly increased. For the first time, the gamma-irradiation-induced valence changes between Yb3+ and Yb2+ ions in Yb: YAP crystals have been observed. (c) 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Effect of air annealing on the spectral properties of Ce : Y3Al5O12 single crystals grown by the temperature gradient technique

Five absorption hands, at 227, 300 340, 370 and 457nm, were observed in the optical absorption spectrum of Ce:Y3Al5O12 (Ce:YAG) crystals grown by the temperature gradient technique (TGT). The absorption bands at 227, 340, and 457 nm were identified Lis belonging to the Ce3+ -ion in the YAG crystal. A near UV optical emission band at 398nm was observed. with an excitation spectrum containing two bands, at 235 and 370nm. No fluorescence was detected under 300 nm excitation. The pair of absorption bands at 235 and 370 nm and the absorption band at 300 nm were attributed to the F- and F+-type color centers, respectively. The color centers model was also applied to explain the spectral changes in the Ce:YAG (TGT) crystal, including the reduction in the Ce 31 -ion absorption intensity, after annealing in an oxidizing atmosphere (air). (C) 2004 Elsevier B.V. All rights reserved.