COORDINATION CHEMISTRY AND SUPERCONDUCTIVITY

Analysing the coordination state of copper ions in cuprate superconductors, it is found that the larger the energy splitting between d(x2-y2) and d(z2) orbitals of Cu or the higher the energy of the d(x2-y2) orbital, the higher the Tc. Thus, appropriate coordination structures and strong-field ligands must be chosen for expanding the energy splitting and increasing the energy of the d(x2-y2) orbital when searching for new high-Tc superconductors. Summarizing the experimental results of ESR and XPS, it is considered that the [Cu2+ - O open-square-box 2- - Cu3+] resonance exists in cuprate superconductors and the electron field breathing mode is present. Analysing the mechanism and the relationship between the coordination state of Cu and Tc, we consider that the two dimensional Cu-O planes are responsible for the superconductivity of YBa2Cu3O7-y.

OBSERVATION OF MOSSBAUER-EFFECT IN THE HIGH-TEMPERATURE SUPERCONDUCTOR EUBA2CU3O7-X WITH DC ELECTRIC-CURRENT

The difference between the Mossbauer parameters for EuBa2Cu3O7-x with dc electric current and those without dc electric current at 83 K has been observed. The change in isomer shift, electric quadrupole splitting and the asymmetry parameter of the electric field gradient at the Eu-151 nucleus may be caused by the movement of a mass of conduction electrons along a certain direction in the EuBa2Cu3O7-x crystal with a layered structure.

复杂条件地层滤波预测的理论与关键技术研究

The theory researches of prediction about stratigraphic filtering in complex condition are carried out, and three key techniques are put forward in this dissertation. Theoretical aspects: The prediction equations for both slant incidence in horizontally layered medium and that in laterally variant velocity medium are expressed appropriately. Solving the equations, the linear prediction operator of overlaid layers, then corresponding reflection/transmission operators, can be obtained. The properties of linear prediction operator are elucidated followed by putting forward the event model for generalized Goupillaud layers. Key technique 1: Spectral factorization is introduced to solve the prediction equations in complex condition and numerical results are illustrated. Key technique 2: So-called large-step wavefield extrapolation of one-way wave under laterally variant velocity circumstance is studied. Based on Lie algebraic integral and structure preserving algorithm, large-step wavefield depth extrapolation scheme is set forth. In this method, the complex phase of wavefield extrapolation operator’s symbol is expressed as a linear combination of wavenumbers with the coefficients of this linear combination in the form of the integral of interval velocity and its derivatives over depth. The exponential transform of the complex phase is implemented through phase shifting, BCH splitting and orthogonal polynomial expansion. The results of numerical test show that large-step scheme takes on a great number of advantages as low accumulating error, cheapness, well adaptability to laterally variant velocity, small dispersive, etc. Key technique 3: Utilizing large-step wavefield extrapolation scheme and based on the idea of local harmonic decomposition, the technique generating angle gathers for 2D case is generalized to 3D case so as to solve the problems generating and storing 3D prestack angle gathers. Shot domain parallel scheme is adopted by which main duty for servant-nodes is to compute trigonometric expansion coefficients, while that for host-node is to reclaim them with which object-oriented angle gathers yield. In theoretical research, many efforts have been made in probing into the traits of uncertainties within macro-dynamic procedures.

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.

煤系地层三维三分量地震数据解释方法研究

Conventional 3D seismic exploration cannot meet the demand of high yield and high efficiency safe production in coal mine any more. Now it is urgent to improve the discovery degree of coal mine geological structures for coal production in China. Based on 3D3C seismic exploration data, multi-component seismic information is fully excavated. First systematic research on 3D3C seismic data interpretation of coal measure strata is carried out. Firstly, by analyzing the coal measure strata, the seismic-geologic model of coal measure strata is built. Shear wave logging is built by using regression analysis. Horizon calibration methods of PP-wave and PS-wave are studied and the multi-wave data are used together to interpret small faults. Using main amplitude analysis technology, small faults which cannot be found from PP-wave sections can be interpreted from the low frequency PS-wave sections. Thus, the purpose to applying PS-wave data to fine structure assistant interpretation is achieved. Secondly, PP- and PS-wave post-stack well constrained inversion methods of coal measure strata are studied. Joint PP- and PS-wave post-stack inversion flow is established. More attribute parameters, which are applied in fine lithology interpretation of coal measure strata, are obtained from combinations of the inversion results. Exploring the relation between rock with negative Poisson’s ratio and anisotropy, fracture development in coal seam are predicted. Petrophysical features of coal measure strata are studied, and the relations between elastic parameters and lithology, fluid and physical properties are established. Inversions of the physical parameters such as porosity, permeability and water saturation, which reflect lithology and fluid property, are obtained. Finally, the approaches of shear wave splitting and Thomsen parameters inversion, which provide new ideas for seismic anisotropy interpretation of coal measure strata, are studied to predict fracture development. The results of practical application indicate that the methods in this paper have good feasibility and applicability. They have positive significance for high yield and high efficiency safe production in coal mine.

三维非均匀各向异性介质中的纵波走时反演研究

Many observations show that seismic anisotropy is very common in the crust and upper mantle of the Earth. Seismic anisotropy can provide some clue about the changing and transporting process inside the earth. in recent years, abundant earthquake travel time data are accumulated, computers become more powerful, and these make the inversion of earthquake travel time data practical. In this thesis we studied the theory of elastic wave in anisotropic media, some formule for travel time inversion were derived. We present an iterative procedure to determine 21 elastic parameters from qP wave travel times. No a priori assumptions about heterogeneity and anisotropy of the model are made. The procedure is suitable for the case when we know nothing about the symmetry of anisotropy of the media, as well as for the case of earthquake travel time inversion which may contain various symmetry of anisotropy. The procedure is tested with a synthetic multiple-source offset VSP experiment. The results proved that the formulae are correct, and the procedure is practical. The results and the related theory indicate that the anisotropic inversion needs more rays than isotropic case. For a 2-D weak anisotropic (WA) medium, we need at least 5 rays in different directions to retrieve the elastic parameters on one grid point, and for a 3-D WA medium we need at least 15 rays in different directions to retrieve the elastic parameters on one grid point. The results also indicate that the starting background velocity has no influence on the final results, at least for the model we specified. Our results also show that insufficient illumination coverage will slow down the convergence rate, and make the results more sensitive to noise. We apply the procedure to a set of field travel time data. The data is from an artificial seismic observation. This observation is for locating micro-seismic events around a tunnel, its purpose is to find out if the digging process and the stress condition around the tunnel can generate micro-cracks. The size of this area is around 100m. The anisotropy derived from qP travel times is the same as the anisotropy showed by apparent velocities, and is also consistent with the anisotropy derived from S-wave splitting phenomena.

各向异性介质中的地震面波和S波分裂

Surface wave propagation in the anisotropic media and S-wave splitting in China mainland are focused in this M.S. dissertation. We firstly introduced Anderson parameters in the research of surface wave propagation in the anisotropic media were deduced, respectively. By applying the given initial model to the forward calculation of Love wave, we compared dispersion curves of Love wave in the anisotropic media with the one in the isotropic media. the results show that, although the two kind of results are similar with each other, the effect of anisotropy can not be neglected. Furthermore, the variation of anisotropy factors will result in the variation of dispersion curves, especially for high-mode one. The method of grid dispersion inversion was then described for further tectonic inversion. We also deduced inversion equation on the condition that the layered media is anisotropic, and calculated the phase-velocity partial derivatives with respect to the model parameters, P- and S-wave velocities, density, anisotropic parameters for Rayleigh wave and Love wave. Having analyzed the results of phase-velocity partial derivatives, we concluded that the derivatives within each period decreased with the depth increasing, the phase-velocity of surface wave is sensitive to the S-wave velocities and anisotropic factors and is not sensitive to the densities of layers. Dispersion data of Love wave from the events occurred during the period from 1991 to 1998 around the Qinghai and Tibet Plateau, which magnitudes are more than 5.5, have been used in the grid dispersion inversion. Those data have been preprocessed and analyzed in the F-T domain. Then the results of 1°*1° grid dispersion inversion, the pure path dispersion data, in the area of Qianghai and Tibet Plateau were obtained. As an example, dispersion data have been input for the tectonic inversion in the anisotropic media, and the results of anisotropic factors under the region of Qianghai and Tibet Plateau were initially discussed. As for the other part of this dissertation. We first introduced the phenomena of S-wave splitting and the methods for calculation the splitting parameters. Then, We applied Butterworth band-pass filter to S-wave data recorded at 8 stations in China mainland, and analyzed S-wave splitting at different frequency bands. The results show the delay time and the fast polarization directions of S-wave splitting depend upon the frequency bands. There is an absence of S-wave splitting at the station of Wulumuqi (WMQ) for the band of 0.1-0.2Hz. With the frequency band broaden, the delay time of S-wave splitting decreases at the stations of Beijing (BJI), Enshi (ENH), Kunming (KMI) and Mudanjiang (MDJ); the fast polarization direction at Enshi (ENH) changes from westward to eastward, and eastward to westward at Hailaer (HIA). The variations of delay time with bands at Lanzhou (LZH) and qiongzhong (QIZ) are similar, and there is a coherent trend of fast polarization directions at BJI, KMI and MDJ respectively. Initial interpretations to the results of frequency band-dependence of S-wave splitting were also presented.

指数矩阵近似计算在地震波场延拓中的应用

Seismic wave field numerical modeling and seismic migration imaging based on wave equation have become useful and absolutely necessarily tools for imaging of complex geological objects. An important task for numerical modeling is to deal with the matrix exponential approximation in wave field extrapolation. For small value size matrix exponential, we can approximate the square root operator in exponential using different splitting algorithms. Splitting algorithms are usually used on the order or the dimension of one-way wave equation to reduce the complexity of the question. In this paper, we achieve approximate equation of 2-D Helmholtz operator inversion using multi-way splitting operation. Analysis on Gauss integral and coefficient of optimized partial fraction show that dispersion may accumulate by splitting algorithms for steep dipping imaging. High-order symplectic Pade approximation may deal with this problem, However, approximation of square root operator in exponential using splitting algorithm cannot solve dispersion problem during one-way wave field migration imaging. We try to implement exact approximation through eigenfunction expansion in matrix. Fast Fourier Transformation (FFT) method is selected because of its lowest computation. An 8-order Laplace matrix splitting is performed to achieve a assemblage of small matrixes using FFT method. Along with the introduction of Lie group and symplectic method into seismic wave-field extrapolation, accurate approximation of matrix exponential based on Lie group and symplectic method becomes the hot research field. To solve matrix exponential approximation problem, the Second-kind Coordinates (SKC) method and Generalized Polar Decompositions (GPD) method of Lie group are of choice. SKC method utilizes generalized Strang-splitting algorithm. While GPD method utilizes polar-type splitting and symmetric polar-type splitting algorithm. Comparing to Pade approximation, these two methods are less in computation, but they can both assure the Lie group structure. We think SKC and GPD methods are prospective and attractive in research and practice.

多波分量地震资料处理解释方法

To deal with some key problems in multi-component seismic exploration, some methods are introduced in this thesis based on reading amounts of papers about multi-component seismic theories and methods. First, to find a solution for the detection of the fracture density and orientation in igneous, carbonate and shale reservoirs, a large amount of which exist in domestic oil fields with low exploration and development degree, a new fast and slow shear waves separation method called Ratio Method based on S-wave splitting theory is discussed in this thesis, through which the anisotropy coefficient as well as fracture parameters such as density and azimuthal angle can be acquired. Another main point in this thesis involves the application of seismic velocity ratio (Vp/Vs) to predict the Hthological parameters of subsurface medium. To deal with the unfeasibility of velocity ratio calculation method based on time ratio due to the usually low single-noise ratio of S-wave seismic data acquired on land, a new method based on detailed velocity analysis is introduced. Third, pre-stack Kirchhoff integral migration is a new method developed in recent years, through which both S and P component seismic data as well as amplitude ratio of P/S waves can be acquired. In this thesis, the research on untilizing the P and S wave sections as well as amplitude ratio sections to interpret low-amplitude structures and lithological traps is carried out. The fast and slow shear wave separation method is then be applied respectively to detect the density and azimuthal angle of fractures in an igneous rock gas reservoir and the coal formation in a coal field. Two velocity ratio-calculating methods are applied respectively in the lithological prediction at the gas and coal field after summarizing a large amount of experimental results draw domestically and abroad. P and S wave sections as well as amplitude ratio sections are used to identify low-amplitude structures and lithological traps in the slope area of a oil-bearing sedimentary basin. The calculated data concerning fracture density and azimuthal angle through the introduced method matches well with the regional stress and actual drilling data. The predicted lithological data reflects the actual drilling data. Some of the low-amplitude and lithological traps determined by Kirchhoff migration method are verified by the actual drilling data. These results indicate that these methods are very meaningful when dealing with complex oil and gas reservoir, and can be applied in other areas.

Formation of the V center on H-type zeolites with thermal treatment

The effect of thermal treatment on H-MCM-22 and H-ZSM-5 zeolites was investigated using the electron spin resonance technique. A six-line signal (denoted as A, g = 2.048, A = 22. 15 G) was detected on H-MCM-22 after He purging at high temperatures, whose intensities increased with the treating temperature. The same signal was also found on H-ZSM-5 zeolites with different crystal sizes. The paramagnetic center was identified as a V center, namely, a hole of an electron trapped on an oxygen atom bonding to a nearby aluminum atom. These signals appeared only on a dealuminated sample or a sample concomitantly with dealumination. The formation of the hole might involve an electron transferring from the lattice oxygen to a nonframework aluminum species, and the hyperfine splitting is caused by the interaction between the electron hole locating on the p orbit of oxygen and the framework aluminum bonding with the oxygen. The signal disappeared after the sample was exposed to air or oxygen at room temperature. However, the process was reversible. A new set of signals (denoted as B, g(1) = 2.008, g(2) = 2.003, g(3) = 1.9985) was observed after oxygen adsorption on the H-MCM-22 pretreated with He at 973 K or He purging at 973 K on the H-MCM-22 pretreated with oxygen at 813 K, which was attributed to the O- species.