SYNTHESES AND CHARACTERIZATION OF [(MU-CF3CO2)(2)LN(MU-CF3CHO2)ALR(2)CENTER-DOT-2THF](2) AND THEIR CATALYTIC ACTIVITIES FOR POLYMERIZATION OF SOME POLAR MONOMERS

Three new bimetallic complexes were synthesized and crystalized by reactions of (CF3CO2)(3)Ln With R(1) AlR(2)(Ln=Nd and Y, R(1)=H, R=i-C4H9; Ln=Eu, R=R(1)=C2H5) in tetrahydrofuran solution, and their crystal structures were determined using a X-ray diffraction method. The structures and the questions on valence state and noncoplanarity in the structures were confirmed and cracked by means of H-1 NMR and C-13 NMR spectra, especially by C-13-H-1 COSY 2D NMR technique. A general formula of molecules of the three rare earth complexes was defined as follows: [(mu-CF3CO2)(2)Ln(mu-CF3CHO2)AlR(2) . 2THF](2) A mechanism on the formation of the new complexes was also proposed through the following five steps: alkylating, beta-elimination (or hydrogenation), hydrogen transfer, linkage and association. Both Y-Al and Eu-Al complexes function as a catalyst in polymerization of MMA and ECH. The polymer obtained from the first monomer is mainly syndiotactic chain structure and the polymerization of the last monomer shows higher catalytic activity. The Y-Al complex also capable of ring-opening polymerization of THF in case of adding-vary small amount of ECH and a oxonium ion mechanism of THF polymerization was suggested from the analysis of THF polymer terminal.

APPLICATION OF FACTOR-ANALYSIS TO CORRECTION OF SPECTRAL OVERLAP INTERFERENCE IN INDUCTIVELY COUPLED PLASMA ATOMIC EMISSION-SPECTROMETRY

Correction of spectral overlap interference in inductively coupled plasma atomic emission spectrometry by factor analysis is attempted. For the spectral overlap of two known lines, a data matrix can be composed from one or two pure spectra and a spectrum of the mixture. The data matrix is decomposed into a spectra matrix and a concentration matrix by target transformation factor analysis. The component concentration of interest in a binary mixture is obtained from the concentration matrix and interference from the other component is eliminated. This method is applied to correcting spectral interference of yttrium on the determination of copper and aluminium: satisfactory results are obtained. This method may also be applied to correcting spectral overlap interference for more than two lines. Like other methods of correcting spectral interferences, factor analysis can only be used for additive spectral overlap. Results obtained from measurements on copper/yttrium mixtures with different white noise added show that random errors in measurement data do not significantly affect the results of the correction method.

SPECTROPHOTOMETRIC DETERMINATION OF LANTHANUM AND SUM OF OTHER LIGHT RARE-EARTHS BY FLEXIBLE TOLERANCE SIMPLEX-METHOD

The absorptivities of color elements in a mixture can be obtained by using Gauas' elimination with selection of principal element in matrix to the standards. These values can be applied to flexible tolerance simplex method to give the composition of samples. In the exprimental design and data treatment, an effort was made to minimize the errors of results according the principal of optimization. When the difference of absorptivities of color material is significant to the exprimental error, the pr...

Comparative expression analysis of transcription factor genes in the endostyle of invertebrate chordates

The endostyle of invertebrate chordates is a pharyngeal organ that is thought to be homologous with the follicular thyroid of vertebrates. Although thyroid-like features such as iodine-concentrating and peroxidase activities are located in the dorsolateral part of both ascidian and amphioxus endostyles, the structural organization and numbers of functional units are different. To estimate phylogenetic relationships of each functional zone with special reference to the evolution of the thyroid, we have investigated, in ascidian and amphioxus, the expression patterns of thyroid-related transcription factors such as TTF-2/MoxE4 and Pax2/5/8, as well as the forkhead transcription factors FoxQ1 and FoxA. Comparative gene expression analyses depicted an overall similarity between ascidians and amphioxus endostyles, while differences in expression patterns of these genes might be specifically related to the addition or elimination of a pair of glandular zones. Expressions of Ci-FoxE and BbFoxE4 suggest that the ancestral FoxE class might have been recruited for the formation of thyroid-like region in a possible common ancestor of chordates. Furthermore, coexpression of FoxE4, Pax2/5/8, and TPO in the dorsolateral part of both ascidian and amphioxus endostyles suggests that genetic basis of the thyroid function was already in place before the vertebrate lineage. (c) 2005 Wiley-Liss, Inc.

Genetic characterization of asymmetric reciprocal hybridization between the flatfishes Paralichthys olivaceus and Paralichthys dentatus

Interspecific reciprocal crosses between the two flatfishes Paralichthys olivaceus and P. dentatus yielded hybrids with different viabilities. Specifically, the hybrids of P. olivaceus female and P. dentatus male (HI) were found to be viable, while the reciprocal hybrids from P. dentatus female and P. olivaceus male (HII) were completely inviable. All the HII individuals showed morphological deformities and died before first feeding. The chromosome analysis showed that HI individuals had the same chromosome number as parents. However, two chromosomes were missing in HII offspring indicating that the latter were aneuploids. Genomic inheritance from the parents to F-1 progeny was also examined by amplified fragment length polymorphism (AFLP) analyses, and the results showed differences between reciprocal hybrids. Almost all AFLP bands (97.71%) observed in parents were passed on to HI individuals. In contrast, only 86.64% of the AFLP bands from parents were scored in HII individuals. Frequency of lost parental bands was thus significantly higher in HII than that in HI and intraspecific crosses, which was probably associated with chromosomal elimination. In addition, higher segregation distortions were found in hybrids than in controls, although these differences were not significant. The present study indicates that chromosomal elimination and loss of AFLP loci occurred in inviable HII individuals, while such genomic changes were not found in viable HI individuals. Possible implications of such difference on genomic changes for asymmetric viability in reciprocal hybrids are discussed.

中国明对虾免疫系统中抗氧化相关基因的克隆与表达分析

对虾病害在世界范围内的广泛传播，给水产养殖和沿海农村经济造成了重大损失。深入开展对虾免疫机制研究并在此基础上寻找对虾疾病防治的有效方法已成为当务之急。研究表明，当对虾等甲壳动物受到外界病原刺激时，其体内的吞噬细胞在吞噬活动中会激活磷酸己糖支路的代谢，引起呼吸爆发，产生多种活性氧分子。另外，受到病原侵染的对虾还会产生其他多种免疫反应，这些免疫反应将消耗大量的能量（ATP），产能的呼吸链会加速运转，由此也会引发大量活性氧的产生。这些活性氧分子可以杀灭入侵的病原微生物，但同时由于活性氧分子反应的非特异性，它们也会对宿主的细胞、组织和器官造成严重伤害，进而导致对虾生理机能的损伤和免疫系统的破坏。所以，消除对虾体内因过度免疫反应产生的过量氧自由基将能够增强其抵御病原侵染的能力，提高免疫力。本论文从中国明对虾体内克隆了线粒体型超氧化物歧化酶（mMnSOD）、胞质型超氧化物歧化酶（cMnSOD）、过氧化氢酶（Catalase）和过氧化物还原酶（Peroxiredoxin）等四种与免疫系统相关的抗氧化酶基因，分析了它们的分子结构特征，组织分布及应答不同病原刺激的表达变化模式，并对其中的mMnSOD基因和Peroxiredoxin基因进行了体外重组表达、分离纯化和酶活性分析。 采用RACE技术从中国明对虾血细胞中克隆了两个超氧化物歧化酶（SOD）基因，通过序列比对分析发现，其中一个为mMnSOD基因，另一个为cMnSOD基因。mMnSOD基因的cDNA全长为1185个碱基，其中开放阅读框为660个碱基，编码220个氨基酸，其中推测的信号肽为20个氨基酸。多序列比对结果显示中国明对虾mMnSOD基因的推导氨基酸序列与罗氏沼虾、蓝蟹的推导氨基酸序列同源性分别为88%和82%。Northern blot结果表明，该基因在对虾的肝胰脏、血细胞、淋巴器官、肠、卵巢、肌肉和鳃等组织中均有表达。半定量RT-PCR结果显示，对虾感染病毒3 h时，该基因在血细胞和肝胰脏中的转录水平显著升高。此外，通过构建原核表达载体，本研究对该基因进行了体外重组表达，并对纯化的重组蛋白进行了质谱鉴定和酶活分析。cMnSOD基因的cDNA全长为1284个碱基，其中开放阅读框为861个碱基，编码287个氨基酸。多序列比对结果显示中国明对虾cMnSOD基因的推导氨基酸序列与斑节对虾和凡纳滨对虾的同源性高达98%和94%。组织半定量结果显示，cMnSOD基因在对虾被检测的各个组织中均有表达。 另外，半定量RT-PCR结果表明，对虾感染病毒23h时，该基因在肝胰脏中的转录上升到正常水平的3.5倍；而感染后59 h时，该基因在血细胞中的转录上升到正常水平的2.5倍。 利用根据其他生物过氧化氢酶保守氨基酸序列设计的简并引物，结合RACE技术，从中国明对虾肝胰脏中克隆到了过氧化氢酶基因的部分片段，片段长1725个碱基。多序列比对结果发现目前所得中国明对虾Catalase基因部分片段的推导氨基酸序列与罗氏沼虾和皱纹盘鲍Catalase氨基酸序列的同源性分别达到95%和73%。通过实时荧光定量PCR技术对中国明对虾Catalase基因在各个组织中的分布情况及病毒感染后该基因在血细胞和肝胰脏中的转录变化进行了研究。结果发现，该基因在肝胰脏、鳃、肠和血细胞中表达水平较高，在卵巢、淋巴器官和肌肉中的表达水平相对较弱；感染病毒23 h和37 h时，对虾血细胞和肝胰脏中该基因mRNA的表达量分别出现显著性上升。 依据中国明对虾头胸部cDNA文库提供的部分片段信息，结合SMART-RACE技术，从中国明对虾肝胰脏中克隆到了过氧化物还原酶基因（Peroxiredoxin）， 该基因的cDNA全长为942个碱基，其中开放阅读框为594个碱基，编码198个氨基酸。中国明对虾Peroxiredoxin基因的推断氨基酸序列与伊蚊、文昌鱼和果蝇等Peroxiredoxin基因的推断氨基酸序列同源性分别为77%、76%和73%。其蛋白理论分子量为22041.17 Da，pI为5.17。Northern blot结果表明，Peroxiredoxin基因在对虾的肝胰脏、血细胞、淋巴器官、肠、卵巢、肌肉和鳃等组织中均有表达。实时荧光定量PCR结果显示，弧菌感染后，该基因在对虾血细胞和肝胰脏中的转录水平都有明显变化并且表达模式不同。另外，对该基因进行了体外重组表达，并对纯化的重组蛋白进行了质谱鉴定和酶活性分析。酶活性分析表明，复性后的重组蛋白能在DTT存在的条件下还原H2O2。

非均质储层叠前地震反演技术研究及应用

(1) I research on the relationship between elastic parameters, lithology and liquid. It is a physical base for pre-stack seismic inversion. I research all kinds of approximate expressions of Zoeppritz function. Then the relation of all kinds of approximate expressions can be confirmed. The geological model of water sand and gas sand in different depth was designed. Moreover I research on precision of all kinds of approximate expressions. (2) In process of seismic data which aim at amplitude recovery and apply in pre-stack seismic inversion, I advance to adopt double flow chart for different aim. Pre-stack noise elimination, real amplitude recovery and NMO correction of long offset are the key taches. (3) I made a systemic expatiate for the thinking and applicability about all kinds of expressions of elastic impedance. And mathematical model was applied to compare the precision with all kinds of expressions of elastic impedance. I propose a new pre-stack simultaneous inversion which is based on the Zoppritz function and simulated annealing algorithm. This method can ensure calculation precision of reflection coefficient from different incident angle and get a global optimum solution. Therefore this method improves the precision of pre-stack seismic inversion. (4) The object function of P-S wave pre-stack simultaneous inversion was established. I compared the precision and convergence between simultaneous inversion and P-wave inversion. And the results show that simultaneous inversion is superior to P-wave inversion. Through the study of AVO event of transformed wave, AVO characters of different kinds of gas sand were analyzed. (5) I carried out the study work of pre-stack seismic inversion for carbonate reservoir in middle of Tarim basin and sand shale reservoir in Sulige Area of Erdos Basin. The method and technology in this paper was applied to practical work. And I made a prediction for heterogeneous reservoir. Moreover it acquires a good application effect. Key Word: reflection coefficient， amplitude recovery， pre-stack seismic inversion， Heterogeneous reservoir，prediction.

反射地震转换波成像方法研究

Multi-waves and multi-component get more and more attentions from oil industry. On the basis of existent research results, My research focuses on some key steps of OBC 4C datum processing. OBC datum must be preprocessed quite well for getting a good image. We show a flow chart of preprocess including attenuation of noise on multi-component datum、elimination ghost by summing P and Z and rotation of horizontal components. This is a good foundation for the coming steps about OBC processing. How to get exact converted point location and to analyze velocity are key points in processing reflection seismic converted wave data. This paper includes computing converted point location, analyzing velocity and nonhyperbolic moveout about converted waves. Anisotropic affects deeply the location of converted wave and the nonhyperbolic moveout. Supposed VTI, we research anisotropic effect on converted wave location and the moveout. Since Vp/Vs is important, we research the compute method of Vp/Vs from post-stack data and pre-stack data. It is a part of the paper that inversing anisotropic parameter by traveltime. Pre-stack time migration of converted wave is an focus, using common-offset Kirchhoff migration, we research the velocity model updating in anisotropic media. I have achieved the following results: 1) using continued Fractions, we proposed a new converted point approximate equation, when the offset is long enough ,the thomsen’s 2 order equation can’t approximate to the exact location of converted point, our equation is a good approximate for the exact location. 2) our new methods about scanning nonhyperbolic velocity and Vp/Vs can get a high quality energy spectrum. And the new moveout can fit the middle and long offset events. Processing the field data get a good result. 3) a new moveout equation, which have the same form as Alkhalifah’s long offset P wave moveout equation, have the same degree preciseness as thomsen’s moveout equation by testing model data. 4) using c as a function of the ratio offset to depth, we can uniform the Li’s and thomsen’s moveout equation in a same equation, the model test tell us choice the reasonable function C can improve the exact degree of Li’s and thomsen’s equation. 5) using traveltime inversion ,we can get anisotropic parameter, which can help to flat the large offset event and propose a model of anisotropic parameter which will useful for converted wave pre-stack time migration in anisotropic media. 6）using our pre-stack time migration method and flow, we can update the velocity model and anisotropic parameter model then get good image. Key words: OBC, Common converted Point (CCP), Nonhyperbolic moveout equation, Normal moveout correction, Velocity analysis, Anisotropic parameters inversion, Kirchhoff anisotropic pre-stack time migration, migration velocity model updating

表面相关多次波衰减方法与应用研究

The function of seismic data in prospecting and exploring oil and gas has exceeded ascertaining structural configuration early. In order to determine the advantageous target area more exactly, we need exactly image the subsurface media. So prestack migration imaging especially prestack depth migration has been used increasingly widely. Currently, seismic migration imaging methods are mainly based on primary energy and most of migration methods use one-way wave equation. Multiple will mask primary and sometimes will be regarded as primary and interferes with the imaging of primary, so multiple elimination is still a very important research subject. At present there are three different wavefield prediction and subtraction methods: wavefield extrapolation; feedback loop; and inverse-scattering series. I mainly do research on feedback loop method in this paper. Feedback loop method includs prediction and subtraction.Currently this method has some problems as follows. Firstly, feedback loop method requires the seismic data used to predict multiple is full wavefield data, but usually the original seismic data don’t meet this assumption, so seismic data must be regularized. Secondly, Multiple predicted through feedback loop method usually can’t match the real multiple in seismic data and they are different in amplitude, phase and arrrival time. So we need match the predicted multiple and that in seismic data through estimating filtering factors and subtract multiple from seismic data. It is the key for multiple elimination how to select a correct matching filtering method. There are many matching filtering methods and I put emphasis on Least-square adaptive matching filtering and L1-norm minimizing adaptive matching filtering methods. Least-square adaptive matching filtering method is computationally very fast, but it has two assumptions: the signal has minimum energy and is orthogonal to the noise. When seismic data don’t meet the two assumptions, this method can’t get good matching results and then can’t attenuate multiple correctly. L1-norm adaptive matching filtering methods can avoid these two assumptions and then get good matching results, but this method is computationally a little slow. The results of my research are as follows: 1. Proposed a method that interpolates seismic traces based on F-K migration and demigration. The main advantage of this method is that it can interpolate seismic traces in any offsets. It shows this method is valid through a simple model. 2. Comparing different Least-square adaptive matching filtering methods. The results show that equipose multi-channel adaptive matching filtering methods can get better results of multiple elimination than other matcing methods through three model data and two field data. 3. Proposed equipose multi-channel L1-norm adaptive matching filtering method. Because L1-norm is robust to large amplitude differences, there are no assumption on the signal has minimum energy and orthogonality, this method can get better results of multiple elimination. 4. Research on multiple elimination in inverse data space. The method is a new multiple elimination method and it is different from those methods mentioned above.The advantages of this method is that it is simple in theory and no need for the adaptive subtraction and computationally very fast. The disadvantage of this method is that it is not stabilized in its solution. The results show that equipose multi-channel and equipose pesudo-multi-channel least-square matching filtering and equipose multi-channel and equipose pesudo-multi-channel L1-norm matching filtering methods can get better results of multiple elimination than other matcing methods through three model data and many field data.

地震波传播模拟的边界元―体积元积分方程法及其近似求解

To improve the efficiency of boundary-volume integral equation technique, this paper is involved in the approximate solutions of boundary-volume integral equation technique. Firstly, based on different interpretations of the self-interaction and extrapolation operators of the resulting boundary integral equation matrix, two different hybrid BEM+Born series modeling schemes are formulated and validated through comparisons with the full-waveform BE numerical solutions for wave propagation simulation in a semicircular alluvial valley and a complex fault model respectively. Numerical experiments indicate that both the BEM+Born series modeling schemes are suitable for complex geological structures and significantly improve computational efficiency especially for the cases of high frequencies and multisource seismic survey. Then boundary-volume integral equation technique is illuminated in detail and verified by modeling wave propagation in complex media. Furthermore, the first-order and second-order Born approximate solutions for the volume-scattering waves are studied and quantified by numerical simulation in different random medium models. Finally, preconditioning generalized minimal residual method is applied to solve boundary-volume integral equation and compared with Gaussian elimination method. Numerical experiments indicate this method makes the calculations more efficient.

隐伏金矿床预测中的地震勘探险方法

Based on brief introduction of seismic exploration and it's general developing situation, the seismic exploration method in field work implementation and some problems frequently encountered in field, which should be pay attention to, are analyzed in detail. The most economic field work techniques are emphasized. Then the seismic data processing flow and it's interpretation technique about the processing results are presented. At last four examples of seismic prospecting in gold deposits are showed. The main conclusions of our research are: 1. Seismic prospecting technique is a very efficient method in the prediction of concealed gold deposits. With appropriate application, it can absolutely reflect the detail underground geological structure in the condition of rugged area and complicated geological environment. 2. The field geometry should be designed and changed according to different kinds of objective exploration depth and ground situation. The best field implementing parameters which include offset, the distance between two adjacent traces, the quantity of dynamite and the depth of hole for explosion, should be determined with examination. Only this way, the high quality original seismic data can be gotten. 3. In seismic data processing, the edition of invalid trace and source gather, signal enhancement, velocity analysis and migration are the key steps. It has some different points with conventional processing and needs a new processing flow and methods which is suitable to the data acquired in rugged area and complicated geological environment. 4. The new common reflection area stacking method in crooked line data processing is an efficient method to improve the signal to noise ratio of seismic data The innovations of our research work are: 1. In the areas which were considered to be forbidden zone, we implement the seismic exploration in several gold deposits in China through our application. All acquire distinguished effects. This show the seismic exploration method is a new effective method in the prediction of concealed gold deposits. 2. We developed a set of seismic field work techniques and data processing which is suitable to complex environment, especially find a effective method in stacking and noise elimination in crooked line data processing. 3. In the field of seismic profile interpretation, through our research work, we are convinced of that: in different kinds of geological condition, the seismic reflection character are not same. For example the lava, the intrusion rock and sediment layers are different in the character of reflection structure and strength. So we accumulate some experience about seismic data interpretation in the area of gold deposits.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.

f-x-y域去除随机噪音方法研究

The seismic survey is the most effective prospecting geophysical method during exploration and development of oil/gas. The structure and the lithology of the geological body become increasingly complex now. So it must assure that the seismic section own upper resolution if we need accurately describe the targets. High signal/noise ratio is the precondition of high-resolution. For the sake of improving signal/noise ratio, we put forward four methods for eliminating random noise on the basis of detailed analysis of the technique for noise elimination using prediction filtering in f-x-y domain. The four methods are put forward for settling different problems, which are in the technique for noise elimination using prediction filtering in f-x-y domain. For weak noise and large filters, the response of the noise to the filter is little. For strong noise and short filters, the response of the noise to the filter is important. For the response of the noise, the predicting operators are inaccurate. The inaccurate operators result in incorrect results. So we put forward the method using prediction filtering by inversion in f-x-y domain. The method makes the assumption that the seismic signal comprises predictable proportion and unpredictable proportion. The transcendental information about predicting operator is introduced in the function. The method eliminates the response of the noise to filtering operator, and assures that the filtering operators are accurate. The filtering results are effectively improved by the method. When the dip of the stratum is very complex, we generally divide the data into rectangular patches in order to obtain the predicting operators using prediction filtering in f-x-y domain. These patches usually need to have significant overlap in order to get a good result. The overlap causes that the data is repeatedly used. It effectively increases the size of the data. The computational cost increases with the size of the data. The computational efficiency is depressed. The predicting operators, which are obtained by general prediction filtering in f-x-y domain, can not describe the change of the dip when the dip of the stratum is very complex. It causes that the filtering results are aliased. And each patch is an independent problem. In order to settle these problems, we put forward the method for eliminating noise using space varying prediction filtering in f-x-y domain. The predicting operators accordingly change with space varying in this method. Therefore it eliminates the false event in the result. The transcendental information about predicting operator is introduced into the function. To obtain the predicting operators of each patch is no longer independent problem, but related problem. Thus it avoids that the data is repeatedly used, and improves computational efficiency. The random noise that is eliminated by prediction filtering in f-x-y domain is Gaussian noise. The general method can't effectively eliminate non-Gaussian noise. The prediction filtering method using lp norm (especially p=l) can effectively eliminate non-Gaussian noise in f-x-y domain. The method is described in this paper. Considering the dip of stratum can be accurately obtained, we put forward the method for eliminating noise using prediction filtering under the restriction of the dip in f-x-y domain. The method can effectively increase computational efficiency and improve the result. Through calculating in the theoretic model and applying it to the field data, it is proved that the four methods in this paper can effectively solve these different problems in the general method. Their practicability is very better. And the effect is very obvious.

利用分频与反演提高地震资料信噪比与分辨率方法研究

This thesis mainly talks about the wavelet transfrom and the frequency division method. It describes the frequency division processing on prestack or post-stack seismic data and application of inversion noise attenuation, frequency division residual static correction and high resolution data in reservoir inversion. This thesis not only describes the frequency division and inversion in theory, but also proves it by model calculation. All the methods are integrated together. The actual data processing demonstrates the applying results. This thesis analyzes the differences and limitation between t-x prediction filter and f-x prediction filter noise attenuation from wavelet transform theory. It considers that we can do the frequency division attenuation process of noise and signal by wavelet frequency division theory according to the differences of noise and signal in phase, amplitude and frequency. By comparison with the f-x coherence noise, removal method, it approves the effects and practicability of frequency division in coherence and random noise isolation. In order to solve the side effects in non-noise area, we: take the area constraint method and only apply the frequency division processing in the noise area. So it can solve the problem of low frequency loss in non-noise area. The residual moveout differences in seismic data processing have a great effect on stack image and resolutions. Different frequency components have different residual moveout differences. The frequency division residual static correction realizes the frequency division and the calculation of residual correction magnitude. It also solves the problems of different residual correction magnitude in different frequency and protects the high frequency information in data. By actual data processing, we can get good results in phase residual moveout differences elimination of pre-stack data, stack image quality and improvement of data resolution. This thesis analyses the characters of the random noises and its descriptions in time domain and frequency domain. Furthermore it gives the inversion prediction solution methods and realizes the frequency division inversion attenuation of the random noise. By the analysis of results of the actual data processing, we show that the noise removed by inversion has its own advantages. By analyzing parameter's about resolution and technology of high resolution data processing, this thesis describes the relations between frequency domain and resolution, parameters about resolution and methods to increase resolution. It also gives the processing flows of the high resolution data; the effect and influence of reservoir inversion caused by high resolution data. Finally it proves the accuracy and precision of the reservoir inversion results. The research results of this thesis reveal that frequency division noise attenuation, frequency residual correction and inversion noise attenuation are effective methods to increase the SNR and resolution of seismic data.

基于结构分析地震资料去噪方法研究

With the development of oil and gas exploration, the exploration of the continental oil and gas turns into the exploration of the subtle oil and gas reservoirs from the structural oil and gas reservoirs in China. The reserves of the found subtle oil and gas reservoirs account for more than 60 percent of the in the discovered oil and gas reserves. Exploration of the subtle oil and gas reservoirs is becoming more and more important and can be taken as the main orientation for the increase of the oil and gas reserves. The characteristics of the continental sedimentary facies determine the complexities of the lithological exploration. Most of the continental rift basins in East China have entered exploration stages of medium and high maturity. Although the quality of the seismic data is relatively good, this areas have the characteristics of the thin sand thickness, small faults, small range of the stratum. It requests that the seismic data have high resolution. It is a important task how to improve the signal/noise ratio of the high frequency of seismic data. In West China, there are the complex landforms, the deep embedding the targets of the prospecting, the complex geological constructs, many ruptures, small range of the traps, the low rock properties, many high pressure stratums and difficulties of boring well. Those represent low signal/noise ratio and complex kinds of noise in the seismic records. This needs to develop the method and technique of the noise attenuation in the data acquisition and processing. So that, oil and gas explorations need the high resolution technique of the geophysics in order to solve the implementation of the oil resources strategy for keep oil production and reserves stable in Ease China and developing the crude production and reserves in West China. High signal/noise ratio of seismic data is the basis. It is impossible to realize for the high resolution and high fidelity without the high signal/noise ratio. We play emphasis on many researches based on the structure analysis for improving signal/noise ratio of the complex areas. Several methods are put forward for noise attenuation to truly reflect the geological features. Those can reflect the geological structures, keep the edges of geological construction and improve the identifications of the oil and gas traps. The ideas of emphasize the foundation, give prominence to innovate, and pay attention to application runs through the paper. The dip-scanning method as the center of the scanned point inevitably blurs the edges of geological features, such as fault and fractures. We develop the new dip scanning method in the shap of end with two sides scanning to solve this problem. We bring forward the methods of signal estimation with the coherence, seismic wave characteristc with coherence, the most homogeneous dip-sanning for the noise attenuation using the new dip-scanning method. They can keep the geological characters, suppress the random noise and improve the s/n ratio and resolution. The rutine dip-scanning is in the time-space domain. Anew method of dip-scanning in the frequency-wavenumber domain for the noise attenuation is put forward. It use the quality of distinguishing between different dip events of the reflection in f-k domain. It can reduce the noise and gain the dip information. We describe a methodology for studying and developing filtering methods based on differential equations. It transforms the filtering equations in the frequency domain or the f-k domain into time or time-space domains, and uses a finite-difference algorithm to solve these equations. This method does not require that seismic data be stationary, so their parameters can vary at every temporal and spatial point. That enhances the adaptability of the filter. It is computationally efficient. We put forward a method of matching pursuits for the noise suppression. This method decomposes any signal into a linear expansion of waveforms that are selected from a redundant dictionary of functions. These waveforms are chosen in order to best match the signal structures. It can extract the effective signal from the noisy signal and reduce the noise. We introduce the beamforming filtering method for the noise elimination. Real seismic data processing shows that it is effective in attenuating multiples and internal multiples. The s/n ratio and resolution are improved. The effective signals have the high fidelity. Through calculating in the theoretic model and applying it to the real seismic data processing, it is proved that the methods in this paper can effectively suppress the random noise, eliminate the cohence noise, and improve the resolution of the seismic data. Their practicability is very better. And the effect is very obvious.