小波域地震资料的处理方法（WSSP）研究

The modeling formula based on seismic wavelet can well simulate zero - phase wavelet and hybrid-phase wavelet, and approximate maximal - phase and minimal - phase wavelet in a certain sense. The modeling wavelet can be used as wavelet function after suitable modification item added to meet some conditions. On the basis of the modified Morlet wavelet, the derivative wavelet function has been derived. As a basic wavelet, it can be sued for high resolution frequency - division processing and instantaneous feature extraction, in acoordance with the signal expanding characters in time and scale domains by each wavelet structured. Finally, an application example proves the effectiveness and reasonability of the method. Based on the analysis of SVD (Singular Value Decomposition) filter, by taking wavelet as basic wavelet and combining SVD filter and wavelet transform, a new de - noising method, which is Based on multi - dimension and multi-space de - noising method, is proposed. The implementation of this method is discussed the detail. Theoretical analysis and modeling show that the method has strong capacity of de - noising and keeping attributes of effective wave. It is a good tool for de - noising when the S/N ratio is poor. To give prominence to high frequency information of reflection event of important layer and to take account of other frequency information under processing seismic data, it is difficult for deconvolution filter to realize this goal. A filter from Fourier Transform has some problems for realizing the goal. In this paper, a new method is put forward, that is a method of processing seismic data in frequency division from wavelet transform and reconstruction. In ordinary seismic processing methods for resolution improvement, deconvolution operator has poor part characteristics, thus influencing the operator frequency. In wavelet transform, wavelet function has very good part characteristics. Frequency - division data processing in wavelet transform also brings quite good high resolution data, but it needs more time than deconvolution method does. On the basis of frequency - division processing method in wavelet domain, a new technique is put forward, which involves 1) designing filter operators equivalent to deconvolution operator in time and frequency domains in wavelet transform, 2) obtaining derivative wavelet function that is suitable to high - resolution seismic data processing, and 3) processing high resolution seismic data by deconvolution method in time domain. In the method of producing some instantaneous characteristic signals by using Hilbert transform, Hilbert transform is very sensitive to high - frequency random noise. As a result, even though there exist weak high - frequency noises in seismic signals, the obtained instantaneous characteristics of seismic signals may be still submerged by the noises. One method for having instantaneous characteristics of seismic signals in wavelet domain is put forward, which obtains directly the instantaneous characteristics of seismic signals by taking the characteristics of both the real part (real signals, namely seismic signals) and the imaginary part (the Hilbert transfom of real signals) of wavelet transform. The method has the functions of frequency division and noise removal. What is more, the weak wave whose frequency is lower than that of high - frequency random noise is retained in the obtained instantaneous characteristics of seismic signals, and the weak wave may be seen in instantaneous characteristic sections (such as instantaneous frequency, instantaneous phase and instantaneous amplitude). Impedance inversion is one of tools in the description of oil reservoir. one of methods in impedance inversion is Generalized Linear Inversion. This method has higher precision of inversion. But, this method is sensitive to noise of seismic data, so that error results are got. The description of oil reservoir in researching important geological layer, in order to give prominence to geological characteristics of the important layer, not only high frequency impedance to research thin sand layer, but other frequency impedance are needed. It is difficult for some impedance inversion method to realize the goal. Wavelet transform is very good in denoising and processing in frequency division. Therefore, in the paper, a method of impedance inversion is put forward based on wavelet transform, that is impedance inversion in frequency division from wavelet transform and reconstruction. in this paper, based on wavelet transform, methods of time - frequency analysis is given. Fanally, methods above are in application on real oil field - Sansan oil field.

中国北方草本植物及表土有机质碳同位素组成

Recently, more and more attention has been paid to stable isotope ratios in terrestrial depositional systems. Among them, δ~(13)C value is mainly determined by the surface vegetation, while vegetation is directly related to climate, therefore, carbon isotope ratio in soil organic matter and pedogenic carbonate has been employed as an important paleoecological indicator. In order to test the paleoecological information extracted from stable isotope ratios in terrestrial depositional systems, it is necessary to study the relationships between δ~(13)C value in standing terrestrial plants and today climate, as well as between δ~(13)C value in modern surface soil organic matter and standing vegetation. Thus, these relationships were studied in this paper by means of analysing δ~(13)C in standing plants and modem surface soil organic matter in North China. The main results and conclusions are presented as following: 1. According to their δ~(13)C values, 40 C-4 species represent about 16% of the 257 plant species sarnpled from the North China. C-4 photosynthesis mainly occurs in Poaceae, Cyperaceae and Chenopidaceae families, and percentage representation of C-4 photosynthesis is up to 56% in Poaceae family. 2. The δ~(13)C values of C-3 plant species in North China vary from -21.7‰ to -32.0‰ with an average of -27.1‰, and 93% focus on the range of -24.0‰ ～ -30.0‰; δ~(13)C values of C-4 plant species in North China are between -10.0‰ ～ -15.5‰ with an average of -12.9‰, and 90% concentrate on the range of -11.0‰ ～ -15.0‰. 3. The δ~(13)C composition of C-3 plant species collected from Beijing, a semi-moist district, mainly vary between -27.0‰ ～ -30.0‰, and the average is -28.7‰; the δ ~(13)C values of plants in the semi-arid district, east and west to the Liu Pan Moutain, focus on the range of-26.0‰ ～ -29.0‰ and -25.0‰ ～ -28.0‰, respectively, with the mean value of -27.6‰ and -26.6‰, respectively; the δ~(13)C composition in the arid district dominantly vary from -24.0‰ to -29.0‰, with the average of -26.2‰, and among them, the δ~(13)C values of C-3 plant species in deserts are often between - 22‰ ～ -24‰; the δ~(13)C values in the cold mountain district concentrate on the range of -24.0‰ to -29.0‰, with the average of -26.3‰. 4. The main range of δ~(13)C composition of C-4 plant species, derived from Beijing, a semi-moist district, are -13.0‰ ～ -15.5‰; the semi-arid district, -11.0‰ ～ -14.0‰; the arid district, -11.0‰ ～ -14.0‰. The mean values of them are -14.0‰, -12.4‰,-12.7‰, respectively. 5. From east to west in North China, δ~(13)C values of C-3 plant species increase with longitude. The correlation between δ~(13)C ratios of C-3 plant species and longitude is linear. Changing temperate and precipitation and changing atmosphere pressure are spossible explanations. 6. Almost all C-3 plant species have the trends that their δ~(13)C values gradually increase with decreasing precipitation, decreasing temperature and increasing altitude. Our results show the increases of the δ~(13)C value by 0.30 ～ 0.45‰, 0.19 ～ 0.27‰ and 1.1 ～ 1.2‰ per 100 mm, I℃ and 1000 m, respectively, for all C-3 plant species together. 7. The δ~(13)C values of all C-3 plant species together and a part of C-3 species show highly significant linear correlation with the mean annual temperature, the mean annual precipitation and the altitude, and the results suggest that they can be used as proxies of these environmental variables, while, those without highly significant correlation, may be not suitable as the proxies. 8. The extent, which of responses of δ~(13)C composition to environmental variables, is different for each C-3 plant specie. 9. The δ~(13)C variations along altitude and longitude may be non-linear for C-4 p1ant species in North China. The mean annual temperature may be not important influential factor, thus, it suggests that the δ~(13)C composition of C-4 plant species may be not suitable as the proxy of the mean annual temperature. The influences of summer temperature on δ~(13)C values are much bigger than that of annual temperature, among them, the influence of September temperature is biggest. The mean annual precipitation may be one of the dominant influential factors, and it shows a highly significant non-linear relationship with δ~(13)C values, and the result indicates that δ~(13) C composition of C-4 plant species can be employed as the proxy of the mean annual precipitation. 10. The variations of δ~(13)C ratios do not show systematic trends along longitude, latitude and altitude for modern surface soil organic in Northwest China. ll. The δ~(13)C ratios of modern surface soil organic do not exhibit systematic patterns with temperature and precipitation in Northwest China, it suggests that, unless soil organic is transferred from pure C-3 or C-4 vegetation, the δ~(13)C composition of soil organic may be not used as proxies of climatic variables. 12. The δ~(13)C values of modem surface soil organic are heavier than that of standing vegetation, and the difference ofrnean δ~(13)C between them is -2.18‰. 13. Without considering the δ~(13)C difference between vegetation and soil organic, as well as the δ~(13)C drift in various enviromnent, we may not obtain the valuable information of C-3, C-4 relative biomass in vegetation. 14. The C-4 biomass contribution in vegetation increase with decreasing latitude, increasing longitude and decreasing altitude in Northwest China. The C-4 biomass almost are zero in those regions north to 38 ° N, or west to 100°E, or above 2400 m. 15. The C-4 relative biomass in vegetation increase with growing temperature and precipitation. and, C-4 plants are rare at those regions where the mean annual temperature is less 4 ℃, or the mean annual precipitation is less 200 mm, and their biomass contribution in vegetation are almost zero. Both the mean annual temperature and the mean annual average precipitation may be the important influential factors of C-4 distribution, but the dominant factors.

利用重力与地形数据探讨我国大陆若干地区岩石圈有效弹性厚度

The lithosphere behaves as an elastic thin plate overlying a weak fluid. When loaded by topography, it deflects to compensate the topography and isostasy is achieved. The response of the plate is often characterized by the flexural rigidity, or equivalently, by the effective elastic thickness(EET). The relationships between gravity and topography have commonly been employed to investigate quantitively the isostatic compensation of the lithosphere. Cross-spectral techniques have been developed to estimate the admittance and coherence between gravity and topography. The observed admittance and coherence functions, when compared with theoretical admittance and coherence computed for various models and parameters, can provide an estimate of the effective elastic thickness. As for coherence, the wavelength at which the coherence drops from approaching unity at long wavelengths to value approaching zero at short wavelengths is a measure of the rigidity of the lithosphere. This research takes advantage of the high-resolution gravity and topography data in Erdos. Using the coherence technique, we have estimated the effective elastic thickness of this region. Subsurface loads are also considered in our calculation. When the admittance between topography and gravity is obtained, we can filter the topography to give an estimate of gravity anomaly. This would be very important when the gravity data is scarce. Several other regions of China have been selected to investigate the effective elastic thickness of lithosphere. We compare those result with lithosphere thickness obtained through seismological technique and heat flow of the region. We find the effective elastic thickness is always smaller than the seismogenic thickness. This is due to what the effective elastic thickness reflects is only the up elastic part of the lithosphere. And we also find there is some degree of correlation between effective elastic thickness and heat flow. This suggests EET is probably controlled by the thermal state of the lithosphere and correlates with their tectonic age. Thus the estimations of the effective elastic thickness of lithosphere can help to investigate geophysical features of active tectonics of continental lithosphere.

地震波场延拓+菲涅尔体叠加方法及成像效果

In the increasingly enlarged exploration target, deep target layer(especially for the reservoir of lava) is a potential exploration area. As well known, the reflective energy becomes weak because the seismic signals of reflection in deep layer are absorbed and attenuate by upper layer. Caustics and multi-values traveltime in wavefield are aroused by the complexity of stratum. The ratio of signal to noise is not high and the fold numbers are finite(no more than 30). All the factors above affect the validity of conventional processing methods. So the high S/N section of stack can't always be got with the conventional stack methods even if the prestack depth migration is used. So it is inevitable to develop another kind of stack method instead. In the last a few years, the differential solution of wave equation was hold up by the condition of computation. Kirchhoff integral method rose in the initial stages of the ninetieth decade of last century. But there exist severe problems in it, which is are too difficult to resolve, so new method of stack is required for the oil and gas exploration. It is natural to think about upgrading the traditionally physic base of seismic exploration methods and improving those widely used techniques of stack. On the other hand, great progress is depended on the improvement in the wave differential equation prestack depth migration. The algorithm of wavefield continuation in it is utilized. In combination with the wavefield extrapolation and the Fresnel zone stack, new stack method is carried out It is well known that the seismic wavefield observed on surface comes from Fresnel zone physically, and doesn't comes from the same reflection points only. As to the more complex reflection in deep layer, it is difficult to describe the relationship between the reflective interface and the travel time. Extrapolation is used to eliminate caustic and simplify the expression of travel time. So the image quality is enhanced by Fresnel zone stack in target. Based on wave equation, high-frequency ray solution and its character are given to clarify theoretical foundation of the method. The hyperbolic and parabolic travel time of the reflection in layer media are presented in expression of matrix with paraxial ray theory. Because the reflective wave field mainly comes from the Fresnel Zone, thereby the conception of Fresnel Zone is explained. The matrix expression of Fresnel zone and projected Fresnel zone are given in sequence. With geometrical optics, the relationship between object point in model and image point in image space is built for the complex subsurface. The travel time formula of reflective point in the nonuniform media is deduced. Also the formula of reflective segment of zero-offset and nonzero offset section is provided. For convenient application, the interface model of subsurface and curve surface derived from conventional stacks DMO stack and prestack depth migration are analyzed, and the problem of these methods was pointed out in aspects of using data. Arc was put forward to describe the subsurface, thereby the amount of data to stack enlarged in Fresnel Zone. Based on the formula of hyperbolic travel time, the steps of implementation and the flow of Fresnel Zone stack were provided. The computation of three model data shows that the method of Fresnel Zone stack can enhance the signal energy and the ratio of signal to noise effectively. Practical data in Xui Jia Wei Zhi, a area in Daqing oilfield, was processed with this method. The processing results showed that the ability in increasing S/N ratio and enhancing the continuity of weak events as well as confirming the deep configuration of volcanic reservoir is better than others. In deeper target layer, there exists caustic caused by the complex media overburden and the great variation of velocity. Travel time of reflection can't be exactly described by the formula of travel time. Extrapolation is bring forward to resolve the questions above. With the combination of the phase operator and differential operator, extrapolating operator adaptable to the variation of lateral velocity is provided. With this method, seismic records were extrapolated from surface to any different deptlis below. Wave aberration and caustic caused by the inhomogenous layer overburden were eliminated and multi-value curve was transformed into the curve.of single value. The computation of Marmousi shows that it is feasible. Wave field continuation extends the Fresnel Zone stack's application.

内蒙古草地植被群落生物量的梯度特征研究

Relationship between biology and environment is always the theme of ecology. Transect is becoming one of the important methods in studies on relationship between global change and terrestrial ecosystems, especially for analysis of its driving factors. Inner Mongolia Grassland is the most important in China Grassland Transect brought forward by Yu GR. In this study, changes in grassland community biomass along gradients of weather conditions in Inner Mongolia was researched by the method of transect. Methods of regression about biomass were also compared. The transect was set from Eerguna county to Alashan county (38° 07' 35" ~50° 12' 20" N, 101° 55' 25" -120° 20' 46" E) in Inner Mongolia, China. The sample sites were mainly chosen along the gradient of grassland type, meadow steppe-* typical steppe-*desert steppe-*steppification desert-^desert. The study was carried out when grassland community biomass got the peak in August or September, 2003 and 2004. And data of 49 sample sites was gotten, which included biomass, mean annual temperature, annual precipitation, accumulated temperature above zero, annual hours of sunshine and other statistical and descriptive data. The aboveground biomass was harvested, and the belowground biomass was obtained by coring (30 cm deep). Then all the biomass samples were dried within (80 + 5) °C in oven and weighted. The conclusion is as follows: 1) From the northeast to the southwest in Inner Mongolia, along the gradient of grassland type, meadow steppe-*typical steppe-*desert steppe-*steppification desert-* desert, the cover degree of vegetation community reduces. 2) By unitary regression analysis, biomass is negatively correlated with mean annual temperature, s^CTC accumulated temperature, ^10°C accumulated temperature and annual hours of sunshine, among which mean annual temperature is crucial, and positively with mean annual precipitation and mean annual relative humidity, and the correlation coefficient between biomass and mean annual relative humidity is higher. Altitude doesn't act on it evidently. Result of multiple regression analysis indicates that as the primary restrictive factor, precipitation affects biomass through complicated way on large scale, and its impaction is certainly important. Along the gradient of grassland type, total biomass reduces. The proportion of aboveground biomass to total biomass reduces and then increases after desert steppe. The trend of below ground biomass's proportion to total biomass is adverse to that of aboveground biomass. 3) Precipitation is not always the only driving factor along the transect for below-/aboveground biomass ratio of different vegetation type composed by different species, and distribution of temperature and precipitation is more important, which is much different among climatic regions, so that the trend of below-/aboveground biomass ratio along the grassland transect may change much through the circumscription of semiarid region and arid region. 4) Among reproductive allocation of aboveground biomass, only the proportion of stem in total biomass notably correlates to the given parameters. Stem/leaf biomass ratio decreases when longitude and latitude increase, caloric variables decrease, and variables about water increase from desert to meadow steppe. The change trends are good modeled by logarithm or binomial equations. 5) 0'-10 cm belowground biomass highly correlates to environmental parameters, whose proportion to total biomass changes most distinctly and increases along the gradient from the west to the east. The deeper belowground biomass responses to the environmental change on the adverse trend but not so sensitively as the surface layer. Because the change value of 0~10 cm belowground biomass is always more than that of below 10 cm along the gradient, the deference between them is balanced by aboveground biomass's change by the resource allocation equilibrium hypothesis.

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

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

Clinical significance and prognostic value of urinary nucleosides in breast cancer patients

Objective: Thirteen urinary nucleosides, primarily degradation products of tRNA, were evaluated as potential tumor markers for breast cancer patients.