Genome-wide survey of putative serine/threonine protein kinases in cyanobacteria

Background: Serine/threonine kinases (STKs) have been found in an increasing number of prokaryotes, showing important roles in signal transduction that supplement the well known role of two-component system. Cyanobacteria are photoautotrophic prokaryotes able to grow in a wide range of ecological environments, and their signal transduction systems are important in adaptation to the environment. Sequence information from several cyanobacterial genomes offers a unique opportunity to conduct a comprehensive comparative analysis of this kinase family. In this study, we extracted information regarding Ser/Thr kinases from 21 species of sequenced cyanobacteria and investigated their diversity, conservation, domain structure, and evolution. Results: 286 putative STK homologues were identified. STKs are absent in four Prochlorococcus strains and one marine Synechococcus strain and abundant in filamentous nitrogen-fixing cyanobacteria. Motifs and invariant amino acids typical in eukaryotic STKs were conserved well in these proteins, and six more cyanobacteria- or bacteria-specific conserved residues were found. These STK proteins were classified into three major families according to their domain structures. Fourteen types and a total of 131 additional domains were identified, some of which are reported to participate in the recognition of signals or substrates. Cyanobacterial STKs show rather complicated phylogenetic relationships that correspond poorly with phylogenies based on 16S rRNA and those based on additional domains. Conclusion: The number of STK genes in different cyanobacteria is the result of the genome size, ecophysiology, and physiological properties of the organism. Similar conserved motifs and amino acids indicate that cyanobacterial STKs make use of a similar catalytic mechanism as eukaryotic STKs. Gene gain-and-loss is significant during STK evolution, along with domain shuffling and insertion. This study has established an overall framework of sequence-structure-function interactions for the STK gene family, which may facilitate further studies of the role of STKs in various organisms.

Genome-wide survey of putative serine/threonine protein kinases in cyanobacteria

Background: Serine/threonine kinases (STKs) have been found in an increasing number of prokaryotes, showing important roles in signal transduction that supplement the well known role of two-component system. Cyanobacteria are photoautotrophic prokaryotes able to grow in a wide range of ecological environments, and their signal transduction systems are important in adaptation to the environment. Sequence information from several cyanobacterial genomes offers a unique opportunity to conduct a comprehensive comparative analysis of this kinase family. In this study, we extracted information regarding Ser/Thr kinases from 21 species of sequenced cyanobacteria and investigated their diversity, conservation, domain structure, and evolution. Results: 286 putative STK homologues were identified. STKs are absent in four Prochlorococcus strains and one marine Synechococcus strain and abundant in filamentous nitrogen-fixing cyanobacteria. Motifs and invariant amino acids typical in eukaryotic STKs were conserved well in these proteins, and six more cyanobacteria- or bacteria-specific conserved residues were found. These STK proteins were classified into three major families according to their domain structures. Fourteen types and a total of 131 additional domains were identified, some of which are reported to participate in the recognition of signals or substrates. Cyanobacterial STKs show rather complicated phylogenetic relationships that correspond poorly with phylogenies based on 16S rRNA and those based on additional domains. Conclusion: The number of STK genes in different cyanobacteria is the result of the genome size, ecophysiology, and physiological properties of the organism. Similar conserved motifs and amino acids indicate that cyanobacterial STKs make use of a similar catalytic mechanism as eukaryotic STKs. Gene gain-and-loss is significant during STK evolution, along with domain shuffling and insertion. This study has established an overall framework of sequence-structure-function interactions for the STK gene family, which may facilitate further studies of the role of STKs in various organisms.

Geophysical signatures associated with fluid flow and gas hydrate occurrence in a tectonically quiescent sequence, Qiongdongnan Basin, South China Sea

Interpretation of high-resolution two-dimensional (2D) and three-dimensional (3D) seismic data collected in the Qiongdongnan Basin, South China Sea reveals the presence of polygonal faults, pockmarks, gas chimneys and slope failure in strata of Pliocene and younger age. The gas chimneys are characterized by low-amplitude reflections, acoustic turbidity and low P-wave velocity indicating fluid expulsion pathways. Coherence time slices show that the polygonal faults are restricted to sediments with moderate-amplitude, continuous reflections. Gas hydrates are identified in seismic data by the presence of bottom simulating reflectors (BSRs), which have high amplitude, reverse polarity and are subparallel to seafloor. Mud diapirism and mounded structures have variable geometry and a great diversity regarding the origin of the fluid and the parent beds. The gas chimneys, mud diapirism, polygonal faults and a seismic facies-change facilitate the upward migration of thermogenic fluids from underlying sediments. Fluids can be temporarily trapped below the gas hydrate stability zone, but fluid advection may cause gas hydrate dissociation and affect the thickness of gas hydrate zone. The fluid accumulation leads to the generation of excess pore fluids that release along faults, forming pockmarks and mud volcanoes on the seafloor. These features are indicators of fluid flow in a tectonically-quiescent sequence, Qiongdongnan Basin. Geofluids (2010) 10, 351-368.

The characteristics of fluid potential in mud diapirs associated with gas hydrates in the Okinawa Trough

Many mud diapirs have been identified in the southern Okinawa Trough from a seismic survey using R/V KEXUE I in 2001. The movement and accumulation of free gas related to mud diapirs are discussed in detail by an analysis of fluid potential which is based upon velocity data. It can be found that free gas moves from the higher fluid potential strata to the lower ones and the gas hydrate comes into being during free gas movement meeting the proper criteria of temperature and pressure. In fact, gas hydrates have been found in the upper layers above the mud diapirs and in host rocks exhibiting other geophysical characteristics. As the result of the formation of the gas hydrate, the free gas bearing strata are enclosed by the gas hydrate bearing strata. Due to the high pressure anomalies of the free gas bearing strata the fluid potential increases noticeably. It can then be concluded that the high fluid potential anomaly on the low fluid potential background may be caused by the presence of the free gas below the gas hydrate bearing strata.

Geophysical evaluation methods for buried hill reservoirs in the Jiyang superdepression of the Bohai Bay basin, eastern China

The Jiyang superdepression is one of the richest hydrocarbon accumulations in the Bohai Bay basin, eastern China. Comprehensive seismic methods have been used in buried hill exploration in Jiyang to describe these fractured reservoirs better. Accurate seismic stratigraphic demarcation and variable-velocity mapping were applied to reveal the inner structure of the buried hills and determine the nature of the structural traps more precisely. Based on the analysis of rock properties and the characteristics of well-developed buried hill reservoirs, we have successfully linked the geology and seismic response by applying seismic forward technology. Log-constrained inversion, absorption coefficient analysis and tectonic forward-inversion with FMI loggings were applied to analyse and evaluate the buried hill reservoirs and gave satisfying results. The reservoir prediction was successful, which confirmed that the comprehensive utilization of these methods can be helpful in the exploration of buried hill reservoirs.

On SAR wind speed ambiguities and related geophysical model functions

Recent investigations show that normalized radar cross sections for C-band microwave sensors decrease under high wind conditions with certain incident angles instead of increase, as is the case for low to moderate wind speeds. This creates the problem of ambiguities in high wind speed retrievals from synthetic aperture radar (SAR). In the present work, four geophysical model functions (GMFs) are studied, namely the high wind C-band model 4 (CMOD4HW), C-band model 5 (CMOD5), the high wind vertical polarized GMF (HWGMF_VV), and the high wind horizontal polarized GMF (HWGMF_HH). Our focus is on model behaviours relative to wind speed ambiguities. We show that, except for CMOD4HW, the other GMFs exhibit the wind speed ambiguity problem. To consider this problem in high wind speed retrievals from SAR, we focus on hurricanes and propose a method to remove the speed ambiguity using the dominant hurricane wind structure.

A survey for paralytic shellfish poisoning (PSP) in Vancouver Harbour

Shellfish samples were collected from seven inter-tidal and two sub-tidal sites between 23 May and 8 June 1999 in Vancouver Harbour and were analysed for paralytic shellfish poisoning (PSP) using a mouse bioassay. PSP was detected in mussels collected at four sampling sites in English Bay and Burrard Inlet, at a concentration below 20 mug saxitoxin equivalents (STXeq)/100 g wet weight. (C) 2003 Elsevier Ltd. All rights reserved.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

内蒙古中部道郎和都格铜多金属矿的地质特征矿床类型及找矿模式研究

Daolangheduge copper polymetallic deposit is located on east edge of Ondor Sum-Bainaimiao metallogenic belt, which is a prospective area of porphyry copper deposit, in Xianghuangqi of central Inner Mongolia. Geotectonically, it occurred in the continental margin accretion belt along the north margin of North China Plate, south of the suture zone between North China Plate and Siberian Plate. The intrusive rocks in this area mainly consist of intermediate-acid magmatic rocks, and the quartz veins, tourmaline veins and the transitional phase are comparatively developed. According to our research, the ore-bearing rock body is mainly quartz diorite while the surrounding rock is mainly biotite granite. Besides, the wall rock alteration are mainly propylitization, pyritization and silicification, which consist of epidotization, actinolitization, chloritzation and so on. The metallic minerals are mainly chalcopyrite and pyrite. In addition, the primary ore is mainly of quartz-chalcopyrite-pyrite type. Above all, Daolangheduge copper polymetallic deposit is suggested to be categorized in the porphyry copper type. With isotopic dating and geochemical research on quartz diorite of ore-bearing rock body, the zircon LA-ICP-MS U-Pb dating of two samples yields an age of 266±2 Ma, falling into the range of late Permian Epoch. It is the first accurate age data in Xianghuangqi area, so it should play a key role in the research of deposit and magmatic rocks in this area. With the major elements and trace elements analysis of 14 samples, the quartz diorite should be among the calc-alkaline series, the geochemical characteristics show higher large-ion lithophile elements of Rb, Sr and LREE, low high-field strength elements of Nb, Ta and high transition elements of Cu, Cr . Also, the REE patterns have negative Eu anomalies. With the same analysis of 4 sample for the biotite granite, the geochemical characteristics show higher Rb, Th,, Zr, Hf and LREE, low Nb, Sm and HREE and Eu has no anomaly. It should be among the calc-alkaline series, over aluminum quality and has characteristics of Adakites. According to isotopic dating and geochemical characteristics of ore-bearing rock body, it is suggested that its materials mainly derived from upper mantle that had fractional crystallization and its magma source region may be affected by fluid metasomatism of paleo-asian ocean. It should be an extensional process of post-orogeny according to regional tectonic evolution. Consequently, because of the decrease of temperature and pressure, the ore forming fluid was raised to surface and mineralized accompanied by magmatic activity which might occur in south of the suture zone. By geological survey, further geophysical and geochemical work is needed. In this area, we have accomplished high precision magnetic prospecting, high density electrical survey, gravity prospecting, soil geochemical prospecting, X-ray fluorescence analyzer prospecting and so on. According to geophysical and geochemical abnormal and surface occurrence, 11 drills are arranged to verification. The type of ores are mainly quartz-chalcopyrite-pyrite ores within 3 drills by drill core logging. Although the grade as well as the scale of already-found Cu deposits are insufficient for industrial exploitation, the mineralization prospect in this region is supposed to be great and the potential in mineral exploration at depth is excellent.

活动构造区区域地壳稳定性研究——以金沙江水电梯级工程规划区为例

The study of regional crustal stability of active tectonic region basically includes analysis of recent activity of Earth's crust, single factor assessment, study of complexity, and comprehensive assessment of crust stability. In this thesis, some work are made as follows: · Based on abundant data from gravity field, aeromagnetic survey, magnetism, magnetotelluric deep sounding, remote sensing and geotectonic as well as earthquakes observed in recent years around this region and adjacent zones, we can get a through understanding about the structural features and activity of the earth's crust in Chuan-Dian region. The results from explosion earthquake and telluric electromagnetic sounding are consistent with the structural features of the crust manifested by the geophysical field. The data of deep geologic structures are important for us to work out a vivid three-dimensional model of the earth's crustal structure of the Jinsha River region. According to a synthesis, the author of this thesis proposes some indicators for dividing the faulted blocks. It can also be inferred that the movement of the Chuan-Dian faulted block, which is the relatively active part of southwestern China, is controlled by the boundary faults, and the intensive activity and deformation are concentrated along the boundaries of the block. · Mainly discussing respectively the mechanism and laws of active faults, earthquakes, and geological hazards activity, and their influences on the stability and security of engineering, also trying to probe into the way to assess the risk of single factor in this section. Especially with the method of fractal geometry, the thesis has discussed how to study the complexity of each factor. These geologic hazards which are distributed at the uppermost part of the crust in this region form a typical mountainous set of the active tectonic areas. The results of survey show that some slopes are liable- to -sliding with a weak layer of low shear strength. Occurrences of landslides are to a great extent related to local geological structures, in particular active faults. This is why numerous landslides have occurred simultaneously around the epicenter of a strong earthquake or the center of a strong rainfall, which are related to active faults. · The analysis of the crustal stability is based on a regional grid division, and a fuzzy comprehensive analysis method is used to determine the grade of the quality in each grid. The evaluation factors and their weights are taken from the results of the hierarchical analysis. The evaluation indexes consist of qualitative and quantitative ones. The qualitative ones can be quantified through the experts weighing system, while the quantitative ones can be obtained from statistical analysis. For quality grades, four levels are used: stable, essentially stable, sub-stable, and unstable. The results of the evaluation on Jinshajiang region demonstrate that the crustal stability become distinctly worse in the areas controlled by active deep faults. Therefore, detailed investigations on the active faulting and geologic hazards, include earthquake activity are especially necessary for those areas adjacent to the deep fault belts. On the bases of the data available and the survey results, we have made a preliminary assessment for the construction conditions and adaptability of every planned site in the middle or lower reaches of Jinsha River. Finally, the thesis prospected the vista of the study of crustal stability.

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

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.

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.

基于信息重建技术的地震信号处理方法研究

The dissertation addressed the problems of signals reconstruction and data restoration in seismic data processing, which takes the representation methods of signal as the main clue, and take the seismic information reconstruction (signals separation and trace interpolation) as the core. On the natural bases signal representation, I present the ICA fundamentals, algorithms and its original applications to nature earth quake signals separation and survey seismic signals separation. On determinative bases signal representation, the paper proposed seismic dada reconstruction least square inversion regularization methods, sparseness constraints, pre-conditioned conjugate gradient methods, and their applications to seismic de-convolution, Radon transformation, et. al. The core contents are about de-alias uneven seismic data reconstruction algorithm and its application to seismic interpolation. Although the dissertation discussed two cases of signal representation, they can be integrated into one frame, because they both deal with the signals or information restoration, the former reconstructing original signals from mixed signals, the later reconstructing whole data from sparse or irregular data. The goal of them is same to provide pre-processing methods and post-processing method for seismic pre-stack depth migration. ICA can separate the original signals from mixed signals by them, or abstract the basic structure from analyzed data. I surveyed the fundamental, algorithms and applications of ICA. Compared with KL transformation, I proposed the independent components transformation concept (ICT). On basis of the ne-entropy measurement of independence, I implemented the FastICA and improved it by covariance matrix. By analyzing the characteristics of the seismic signals, I introduced ICA into seismic signal processing firstly in Geophysical community, and implemented the noise separation from seismic signal. Synthetic and real data examples show the usability of ICA to seismic signal processing and initial effects are achieved. The application of ICA to separation quake conversion wave from multiple in sedimentary area is made, which demonstrates good effects, so more reasonable interpretation of underground un-continuity is got. The results show the perspective of application of ICA to Geophysical signal processing. By virtue of the relationship between ICA and Blind Deconvolution , I surveyed the seismic blind deconvolution, and discussed the perspective of applying ICA to seismic blind deconvolution with two possible solutions. The relationship of PC A, ICA and wavelet transform is claimed. It is proved that reconstruction of wavelet prototype functions is Lie group representation. By the way, over-sampled wavelet transform is proposed to enhance the seismic data resolution, which is validated by numerical examples. The key of pre-stack depth migration is the regularization of pre-stack seismic data. As a main procedure, seismic interpolation and missing data reconstruction are necessary. Firstly, I review the seismic imaging methods in order to argue the critical effect of regularization. By review of the seismic interpolation algorithms, I acclaim that de-alias uneven data reconstruction is still a challenge. The fundamental of seismic reconstruction is discussed firstly. Then sparseness constraint on least square inversion and preconditioned conjugate gradient solver are studied and implemented. Choosing constraint item with Cauchy distribution, I programmed PCG algorithm and implement sparse seismic deconvolution, high resolution Radon Transformation by PCG, which is prepared for seismic data reconstruction. About seismic interpolation, dealias even data interpolation and uneven data reconstruction are very good respectively, however they can not be combined each other. In this paper, a novel Fourier transform based method and a algorithm have been proposed, which could reconstruct both uneven and alias seismic data. I formulated band-limited data reconstruction as minimum norm least squares inversion problem where an adaptive DFT-weighted norm regularization term is used. The inverse problem is solved by pre-conditional conjugate gradient method, which makes the solutions stable and convergent quickly. Based on the assumption that seismic data are consisted of finite linear events, from sampling theorem, alias events can be attenuated via LS weight predicted linearly from low frequency. Three application issues are discussed on even gap trace interpolation, uneven gap filling, high frequency trace reconstruction from low frequency data trace constrained by few high frequency traces. Both synthetic and real data numerical examples show the proposed method is valid, efficient and applicable. The research is valuable to seismic data regularization and cross well seismic. To meet 3D shot profile depth migration request for data, schemes must be taken to make the data even and fitting the velocity dataset. The methods of this paper are used to interpolate and extrapolate the shot gathers instead of simply embedding zero traces. So, the aperture of migration is enlarged and the migration effect is improved. The results show the effectiveness and the practicability.

共反射面元叠加

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. As one important seismic data processing method, Stacking is an effective means to suppress the records noise. Broadening area of surface stacked is more important to enhance genuine reflection signals and suppressing unwanted energy in the form of coherent and random ambient noise. Common reflection surface stack is a macro-model independent seismic imaging method. Based on the similarity of CRP trace gathers in one coherent zone, CRS stack effectively improves S/N ratio by using more CMP trace gathers to stack. It is regarded as one important method of seismic data processing. Performing CRS stack depends on three attributes. However, the equation of CRS is invalid under condition of great offset. In this thesis, one method based on velocity model in depth domain is put forward. Ray tracing is used to determine the traveltime of CRP in one common reflection surface by the least squares method to regress the equation of CRS. Then we stack in the coherent seismic data set according to the traveltime, and get the zero offset section. In the end of flowchart of implementing CRS stack, one method using the dip angle to enhance the ratio of S/N is used. Application of the method on synthetic examples and field seismic records, the results of this method show an excellent performance of the algorithm both in accuracy and efficiency.