Preliminary study of the dynamic origin of the distribution pattern of bottom sediments on the continental shelves of the Bohai Sea, Yellow Sea and East China Sea

The bottom sediment types in the Bohai Sea, Yellow Sea and East China Sea (BYECS) are diversified, and their distribution pattern is very complicated. However, the bottom sediment types can be simplified to be sandy sediment, clayey sediment and mixed sediment, which comprise the complicated distribution pattern of bottom sediment in the BYECS. The continental shelves of the BYECS are broad, with shallow water depths and tidal currents which are permanent and dominate the marine dynamics in the BYECS. Based on numerical simulation of tidal elevations and currents in the BYECS, the rates of suspended load transport and bed load transport during a single tidal cycle for sediments of eight different grain size ranges are calculated. The results show that any sediment, whose threshold velocity is less than that of tidal current, has the same transport trend. Suspended load transport rare, bed load transport rate, and the ratio of the former to the latter decrease with grain size becoming coarser and coarser. The erosion/accretion patterns of sediments with different grain sizes are determined by the sediment transport rate divergences, and the results show that the patterns are the same for sediments with different grain sizes. Three main bottom sediment types, i.e. sandy sediment mainly composed of fine sand, clayey sediment mainly composed of silty clay, and mixed sediment mainly composed of fine sand, silt, and clay, are obtained by computation. The three bottom sediment types and their distribution pattern are consistent not only with sediment transport field and the sea bed erosion/accretion pattern obtained by simulation, but also with field data of bottom sediment types and divisions. In the BYECS, sand ridges form mainly in the areas with strong rectilinear tidal currents, sand sheets form mainly in the areas dominated by strong rotatory tidal currents, and clayey sediments, i.e. mud patches, form mainly in the areas with weak tidal currents. Hence, not only the sandy sediments but also the clayey sediments in the BYECS are formed under the control of the whole tidal current field of the BYECS. The three main bottom sediment types are not isolated respectively-in fact, they constitute a whole tidal depositional system. Under the condition with no cyclonic cold eddy, the clayey sediments in the BYECS can form in weak tidal current environments. Therefore, a cold eddy is not necessary for the deposition of clayey sediments in the BYECS. (C) 2000 Academic Press.

Characteristics of nitrogen forms in the surface sediments of southwestern Nansha Trough, South China Sea

The area of the southwestern Nansha Trough is one of the most productive areas of the southern South China Sea. It is a typical semi-deep sea area of transition from shoal to abyssal zone. To understand distributions and roles of nitrogen forms involved in biogeochemical cycling in this area, contents of nitrogen in four extractable forms: nitrogen in ion exchangeable form (IEF-N), nitrogen in weak acid extractable form (WAEF-N), nitrogen in strong alkali extractable form (SAEF-N) and nitrogen in strong oxidation extractable form (SOEF-N), as well as in total nitrogen content (TN) in surface sediments were determined from samples collected from the cruise in April-May 1999. The study area was divided into three regions (A, B and C) in terms of clay sediment (< 4 mu m) content at < 40%, 40%-60% and > 60%, respectively. Generally, region C was the richest in the nitrogen of all forms and region A the poorest, indicating that the finer the grain size is, the richer the contents of various nitrogen are. The burial efficiency of total nitrogen in surface sediments was 28.79%, indicating that more than 70% of nitrogen had been released and participated in biogeochemical recycling through sediment-water interface.

Intra-habitat heterogeneity of microbial food web structure under the regime of eutrophication and sediment resuspension in the large subtropical shallow Lake Taihu, China

Planktonic microbial community structure and classical food web were investigated in the large shallow eutrophic Lake Taihu (2338 km(2), mean depth 1.9 m) located in subtropical Southeast China. The water column of the lake was sampled biweekly at two sites located 22 km apart over a period of twelve month. Site 1 is under the regime of heavy eutrophication while Site 2 is governed by wind-driven sediment resuspension. Within-lake comparison indicates that phosphorus enrichment resulted in increased abundance of microbial components. However, the coupling between total phosphorus and abundance of microbial components was different between the two sites. Much stronger coupling was observed at Site 1 than at Site 2. The weak coupling at Site 2 was mainly caused by strong sediment resuspension, which limited growth of phytoplankton and, consequently, growth of bacterioplankton and other microbial components. High percentages of attached bacteria, which were strongly correlated with the biomass of phytoplankton, especially Microcystis spp., were found at Site 1 during summer and early autumn, but no such correlation was observed at Site 2. This potentially leads to differences in carbon flow through microbial food web at different locations. Overall, significant heterogeneity of microbial food web structure between the two sites was observed. Site-specific differences in nutrient enrichment (i.e. nitrogen and phosphorus) and sediment resuspension were identified as driving forces of the observed intra-habitat differences in food web structure.

The preparation and antioxidant activity of glucosamine sulfate

Glucosamine sulfate was prepared from glucosamine hydrochloride that was produced by acidic hydrolysis of chitin by ion-exchange method. Optical rotation and elemental analysis characterized the degree of its purity. In addition, the antioxidant potency of chitosan derivative-glucosamine sulfate was investigated in various established in vitro systems, such as superoxide (O (2) (-) )/hydroxyl (center dot OH) radicals scavenging, reducing power, iron ion chelating. The following results are obtained: first, glucosamine sulfate had pronounced scavenging effect on superoxide radical. For example the O (2) (-) scavenging activity of glucosamine sulfate was 92.11% at 0.8 mg/mL. Second, the center dot OH scavenging activity of glucosamine sulfate was also strong, and was about 50% at 3.2 mg/mL. Third, the reducing power of glucosamine sulfate was more pronounced. The reducing power of glucosamine sulfate was 0.643 at 0.75 mg/mL. However, its potency for ferrous ion chelating was weak. Furthermore, except for ferrous ion chelating potency, the scavenging rate of radical and reducing power of glucosamine sulfate were concentration-dependent and increased with their increasing concentrations, but its ferrous ion chelating potency decreased with the increasing concentration. The multiple antioxidant activities of glucosamine sulfate were evidents of reducing power and superoxide/hydroxyl radicals scavenging ability. These in vitro results suggest the possibility that glucosamine sulfate could be used effectively as an ingredient in health or functional food, to alleviate oxidative stress.

Studies on two coordination polymers [M(mu(4)-pz25dc)](n) (M = Cd or Zn, pz25dc=pyrazine-2,5-dicarboxylato) with three-dimensional pillared-layer three-nodal framework: Synthesis, structural characterization, strong optical non-linearities and optical limiting properties

Two isomorphous new candidates [M(mu(4)-pz25dc)](n) (M = Cd, 1; Zn, 2; pz25dc = pyrazine-2,5-dicarboxylato)for nonlinear optical (NLO) materials have been synthesized hydrothermally and characterized crystallographically as pillared-layer three-nodal frameworks with one four-connected metal nodes and two crystallographically different four-connected ligand nodes. Their optical non-linearities are measured by the Z-Scan technique with an 8 ns pulsed laser at 532 nm. These two coordination polymers both exhibit strong NLO absorptive abilities [alpha(2) = (63 +/- 6) x 10 (12) mW (1) 1, ( 46 +/- 6) x 10 (11) mW (1) 2] and effective self-focusing performance [n(2) = (67 +/- 5) x 10 (18) 1, (13 +/- 3) x 10 (18) m(2) W (1) 2] in 1.02 x 10 (4) 1 and 1.05 x 10 (4) mol dm (3) 2 DMF solution separately. The values of the limiting threshold are also measured from the optical limiting experimental data. The heavy atom effect plays important role in the enhancement of optical non-linearities and optical limiting properties. (C) 2009 Elsevier B. V. All rights reserved.

Strong temperature dependence and no moss photosynthesis in winter CO2 flux for a Kobresia meadow on the Qinghai-Tibetan plateau

We examined the CO2 exchange of a Kobresia meadow ecosystem on the Qinghai-Tibetan plateau using a chamber system. CO2 efflux from the ecosystem was strongly dependence on soil surface temperature. The COZ efflux-temperature relationship was identical under both light and dark conditions, indicating that no photosynthesis could be detected under light conditions during the measurement period. The temperature sensitivity (Q(10)) of the COZ efflux showed a marked transition around -1.0 degrees C; Q(10) was 2.14 at soil surface temperatures above and equal to -1.0 degrees C but was 15.3 at temperatures below -1.0 degrees C. Our findings suggest that soil surface temperature was the major factor controlling winter COZ flux for the alpine meadow ecosystem and that freeze-thaw cycles at the soil surface layer play an important role in the temperature dependence of winter CO2 flux. (c) 2005 Elsevier Ltd. All rights reserved.

Three phases of dark-recovery course from photoinhibition resolved by the chlorophyll fluorescence analysis in soybean leaves under field conditions

The chlorophyll fluorescence in soybean leaves was observed by a portable fluorometer CF-1000 under field conditions. On clear days, F-0 increased while F, and F-v/F-m decreased gradually in the morning. At midday F-O reached its maximum while F-v and F-v/F-m reached their minimum. The reverse changes occurred in the afternoon. At dusk these parameters could return to levels near those at dawn. Following exposure to a strong sunlight for more than 3 h, the dark-recovery process displayed three phases: (1) slow increases in F-0, F-v and F-v/F-m within the first hour; (2) a faster decrease in F-0 and faster increases in F-v and F-v/F-m within subsequent two hours; (3) a slow decrease in F-0 and slow increases in F-v and F-v/F-m within the fourth hour. In comparison with darkness, weak irradiance had no stimulating effect on the recovery from photoinhibition. Hence the photoinhibition in soybean leaves is mainly the reflection of reversible inactivation of some photosystem 2 reaction centres, but not the result of D1 protein loss.

三塘湖盆地牛圈湖区储层特征研究—以齐古组、头屯河组为例

In view of few researches into pore textures and anisotropy characteristics of Qiguzu-toutunhezu reservoir in Niuquanhu block in Santanghu Basin, in order to enhance the hydrocarbon recovery of the region of interest and improve the reservoir development effect, with the employment of the experiments like cast thin slice, scanning electron microscope, conventional high pressure mercury penetration, constant speed mercury penetration and nuclear magnetic resonance, the thesis makes a thorough research into pore texture characteristics and anisotropy characteristics of the sandstone reservoir. The research shows that the microscopic pore textures are complicated, the anisotropy is high and waterflood development water/oil displacement efficiency is low, which are mainly caused by the high microscopic anisotropy of the reservoir. Specially, the research shows that Qiguzu-toutunhezu reservoir belongs to braided delta front intrafacies, the intergranular pore is the main type of pores, which take up 65.50 percent of total pores, intergranular dissolved pores, feldspar dissolved pores and lithic dissolved pores are on the second place, and there are few carbonate dissolved pores. The reservoir belongs to mesopore-fine throat and mesopore-medium throat. The pore distribution of the sandstone reservoir is comparatively centralized. The reservoir inhomogeneity is mainly caused by the throat inhomogeneity. Diagenesis mainly includes compaction, cementation, metasomatism and dissolution. Compared with compaction， cementation exerts more direct impact on the decline of the porosity of the sandstone reservoir ( pore loss factor is 63.75 percent in the cementing process). Based on the classification of diagenetic stages, the reservoir diagenesis is on the stage of the late period of early diagenetic stage to early period of late diagenetic stage. The study area of the small sandstone reservoir layer of the anisotropy of the relatively strong; plane, the anisotropy of the relatively weak. microscopic anisotropy of its relatively strong.

Fourier有限差分法深度偏移与正演模拟

With the development of seismic exploration, the target becomes more and more complex, which leads to a higher demand for the accuracy and efficiency in 3D exploration. Fourier finite-difference (FFD) method is one of the most valuable methods in complex structure exploration, which keeps the ability of finite-differenc method in dealing with laterally varing media and inherits the predominance of the phase-screen method in stablility and efficiency. In this thesis, the accuracy of the FFD operator is highly improved by using simulated annealing algorithm. This method takes the extrapolation step and band width into account, which is more suitable to various band width and discrete scale than the commonely-used optimized method based on velocity contrast alone. In this thesis, the FFD method is extended to viscoacoustic modeling. Based on one-way wave equation, the presented method is implemented in frequency domain; thus, it is more efficient than two-way methods, and is more convenient than time domain methods in handling attenuation and dispersion effects. The proposed method can handle large velocity contrast and has a high efficiency, which is helpful to further research on earth absorption and seismic resolution. Starting from the frequency dispersion of the acoustic VTI wave equation, this thesis extends the FFD migration method to the acoustic VTI media. Compared with the convetional FFD method, the presented method has a similar computational efficiency, and keeps the abilities of dealing with large velocity contrasts and steep dips. The numerical experiments based on the SEG salt model show that the presented method is a practical migration method for complex acoustical VTI media, because it can handle both large velocity contrasts and large anisotropy variations, and its accuracy is relatively high even in strong anisotropic media. In 3D case, the two-way splitting technique of FFD operator causes artificial azimuthal anisotropy. These artifacts become apparent with increasing dip angles and velocity contrasts, which prevent the application of the FFD method in 3D complex media. The current methods proposed to reduce the azimuthal anisotropy significantly increase the computational cost. In this thesis, the alternating-direction-implicit plus interpolation scheme is incorporated into the 3D FFD method to reduce the azimuthal anisotropy. By subtly utilizing the Fourier based scheme of the FFD method, the improved fast algorithm takes approximately no extra computation time. The resulting operator keeps both the accuracy and the efficiency of the FFD method, which is helpful to the inhancements of both the accuracy and the efficiency for prestack depth migration. The general comparison is presented between the FFD operator and the generalized-screen operator, which is valuable to choose the suitable method in practice. The percentage relative error curves and migration impulse responses show that the generalized-screen operator is much sensiutive to the velocity contrasts than the FFD operator. The FFD operator can handle various velocity contrasts, while the generalized-screen operator can only handle some range of the velocity contrasts. Both in large and weak velocity contrasts, the higher order term of the generalized-screen operator has little effect on improving accuracy. The FFD operator is more suitable to large velocity contrasts, while the generalized-screen operator is more suitable to middle velocity contrasts. Both the one-way implicit finite-difference migration and the two-way explicit finite-differenc modeling have been implemented, and then they are compared with the corresponding FFD methods respectively. This work gives a reference to the choosen of proper method. The FFD migration is illustrated to be more attractive in accuracy, efficiency and frequency dispertion than the widely-used implicit finite-difference migration. The FFD modeling can handle relatively coarse grids than the commonly-used explicit finite-differenc modeling, thus it is much faster in 3D modeling, especially for large-scale complex media.

济阳坳陷深层地质构造与演化

Based on geophysical and geological data in Jiyang depression, the paper has identified main unconformity surfaces (main movement surfaces) and tectonic sequences and established tectonic and strata framework for correlation between different sags. Based on different sorts of structural styles and characteristics of typical structures, the paper summarized characteristics and distribution of deep structures, discussed evolution sequence of structure, analyzed the relation between tectonic evolution and generation of petroleum. The major developments are as following: Six tectonic sequences could be divided from bottom to top in the deep zone of Jiyang depression. These tectonic sequences are Cambrian to Ordovician, Carboniferous to Permian, lower to middle Jurassic, upper Jurassic to lower Cretaceous, upper Cretaceous and Kongdian formation to the fourth member of Shahejie formation. The center of sedimentation and subsidence of tectonic sequences distinguished from each other in seismic profiles is controlled by tectonic movements. Six tectonic evolution stages could be summarized in the deep zone in Jiyang depression. Among these stages, Paleozoic stage is croton sedimentation basin; Indosinian stage, open folds of EW direction are controlled by compression of nearly SN direction in early Indosinian (early to middle Triassic) while fold thrust fault of EW – NWW direction and arch protruding to NNE direction are controlled by strong compression in late Indosinian (latter Triassic); early Yanshanian stage (early to middle Jurassic), in relatively weak movement after Indosinian compressional orogeny, fluviolacustrine is deposited in intermontane basins in the beginning of early Yanshanian and then extensively denudated in the main orogenic phase; middle Yanshanian (late Jurassic to early Cretaceous), strike-slipping basins are wide distribution with extension (negative reversion) of NW – SE direction; latter Yanshanian (late Cretaceous), fold and thrust of NE – NNE direction and positive reversion structure of late Jurassic to early Cretaceous strike-slipping basin are formed by strong compression of NW–SE direction; sedimentation stage of Kongdian formation to the fourth member of Shahejie formation of Cenozoic, half graben basins are formed by extension of SN direction early while uplift is resulted from compression of nearly EW direction latterly. Compression system, extension system and strike-slip system are formed in deep zone of Jiyang depression. According to identifying flower structure of seismic profiles and analysis of leveling layer slice of 3D seismic data and tectonic map of deep tectonic interface, strike-slip structures of deep zone in Jiyang depression are distinguished. In the middle of the Jiyang depression, strike-slip structures extend as SN direction, NNW direction in Huimin sag, but NNE in Zhandong area. Based on map of relict strata thickness, main faults activity and regional tectonic setting, dynamic mechanisms of deep structure are preliminary determination. The main reason is the difference of direction and character of the plate’s movement. Development and rework of multi-stage tectonic effects are benefit for favorable reservoir and structural trap. Based on tectonic development, accumulation conditions of deep sub-sags and exploration achievements in recent years, potential zones of oil-gas reservoir are put forward, such as Dongying sag and Bonan sag.

大陆地壳岩石结构重建：方法与实例

The composition of the continental crust has long been a subject of interest to earth scientists as it can provide key information about the crustal growth and evolution of the continents. In this paper we make a comparative study on the lithological discrimination schemes featuring with the use of different seismic attributes, such as P-wave velocity, P- to S-wave velocity ratio, acoustic or elastic impedances, Lame impedances and high-sensitive identification factors. The results demonstrate that Lame impedances have more powerful constrains than other seismic attributes. In order to fully take the advantage of make the best of the different seismic response of crustal rock, we firstly use seismic attribute that have weak distinguish power to construct loose constrained lithological model, then use seismic attributes that have stronger distinguish power to tighten the constrains of lithological discrimination. We propose a joint scheme (chain constrain technique) by combing all available constrains to reduce the non-uniqueness in mapping rock distribution. We adopt chain constrain technique to construct lithological model beneath Tunxi-Wenzhou transect, Southeastern China, Manzhouli-Suifenhe transect, Northeastern China, and geophysical profile in Bohai Bay Basin, North China. The results can be suumarized as the follows: (1) We compare the sensitivity of different seismic factor constraints on rock types, and conclude that Lame impedances have tighter constrains than seismic velocity, Vp/Vs, density. (2) We propose chain constrains to construct lithological model from integrated geophysical data, and reduce the non-uniqueness in mapping rock distribution. (3) We reconstruct crustal lithological model beneath Tunxi-Wenzhou transect, Southeastern China. The results suggested that Jiangshan-Shaoxing fault is a crust-scale, and it is the boundary between Cathaysia and Yanthze blocks. (4) We construct crustal lithological model beneath Manzhouli-Suifenhe transect, Northeastern China. (5) We map the petrologic distribution along a geophysical profile in Bohai Bay Basin, North China, and construct a three-layered petrology model from the depth 2 km to about 10 km, consisted of basalt (the first layer), pelitic siltstone (the second layer), and silty mudstone and fine sandstone (the third layer).

埕北断阶带油气运聚动力学研究

The problem of oil and gas migration and accumulation have been investigated for many years in petroleum geology field. However, it is still the most weak link. It is a challenge task to research the question about the dynamics of hydrocarbon migration and accumulation. The research area of this article,Chengbei step-fault zone is the important exploration area of Dagang oil field.The oil distribution is complicated in this area because of abundant faults and rock-reservoir-cap assemblage.In recent years, oil shows is often discovered, but no large-scale pool is found. The most important problem influencing exolore decision is lake of kowning about accumulation process of oil and resources potential. According to the geology characteristic and exolore difficult, the analysis principles of dynamics is used in this paper. The course from source to reservoir is considered as main research line, and relation of valid source rcok, migration dynamic and heterogeneous distribution of carrier is discussed especially in key time. By use of numerial model the couling of migration and passage is realized and dynamic process of oil migration is analysed quantitatively. On the basis of other research about structure and sendiment, basin model is built and parameters are choiced. The author has reconstructed characteristic and distribution of fluid dynamical in main pool-forming time by numerical model. The systems of oil migration and acuumulaiton are divided according to distribution of fluid potential. Furthermore, the scope of valid sourece rock and scale of discharging hydrocarbon is studied in geology history by the method of producting hydrocarbon poential. In carrier research, it is outstanding to analyse the function that fault controls the oil-gas migration and accumulation. According to the mechanism of fault sealing, the paper author puts forward a new quantitative method evaluating fault opening and sealing properties-fault connective probability by using the oil and gas shows in footwall and hangwall reservoir as the index of identifying fault sealing or non-sealing. In this method, many influencing factors are considered synthetically. Then the faut sealing propery of different position in third deimention of faults controlling hydrocarbon acummulation are quantitative evaluated, and it laies a foundation for building compex carrier systems. Ten models of carrier and dynamical are establishe by analysis of matching relation of all kinds of carriers in main pool-forming period. The forming process and distribution of main pathway has been studied quantitatively by Buoyancy-Percolation mode, which can conbine valid source rock, migration dynamical and carrier. On the basis of oil-gas migration and accumulation model, the author computes the loss of hydrocarbon in secondary migration, ahead of cap formation, and the quantity of valueless accumulation according to the stage of migration and accumulation and the losing mechamism. At the same time, resource potential is evaluated in every migration and accumulation system. It shows that the quanlity of middle systems arrive to 5.67×108t, which has a huge explore potential prospect. Finally, according to the result of quantitve analysis above mentioned, the favorable explore aims are forcasted by the way of overlapping migration pathway and valid trap and considering factors of pool-forming. The drilling of actual wells proved that the study result is credible. It would offer strong support to optimize explore project in Chengbei step-fault zone.

有源电磁扩散场克希霍夫积分偏移

As active electromagnetic method, field data of CSAMT method follow the equation of diffusion. Propagting in solid earth media, diffusion EM signal has strong attenuation and dispersion, otherwise seismic wave shows weak attenuation and dispersion, therefore the resolution power of CSAMT method is not better than seismic reflection method. However, there is consistence and similarity between EM signal and seismic wave in wave equation, we can apply Kirchhoff integral migration technique, a proven one in seismic method in time domain, to carry out seduo-seismic processing for CSAMT signal in frequency domain so that the attenuation and dispersion could be made compensated in some extent, and the resolution power and interpretation precision of active EM wave could be improved. Satisfying passive homogeneous Helmholtz quation, we proceed with Green theorem and combine the active inhomogenous Helmholtz quation, the Kirchhoff integral formula could be derived. Given practical problems, if we only consider the surface integral value, and assume that the intergral value in other interface is zero, combined with Green theorem in uniform half space, the expression could be simplified, and we can obtain frequency-domain Kirchhoff integral formula in surface, which is also called downward continuation of EM field in frequency domain. With image conditions and energy compensation considered, in order to get image conditions in time domain Fourier inverse transformation in frequency domain can be performed, so we can formulate the active Kirchhoff integral migration expression. At first, we construct relative stratified model, with different frequency series taken into account, then we change the distances between transmitter and reciever, the EM response can be obtained. Analyzing the EM properties, we can clarify near and far zone that can instruct us to carry out transmitter layout in practical application. Combined with field data surveyed in far zone, We perform Kirchhoff integral migration and compare the results with model to interpret. Secondly, with far field EM data, we apply TM mode to get EM response of given 2D model, then apply Kirchhoff integral migration on modelling data and interpret the results.

旋转轴对称弱各向异性介质模型的射线追踪方法

Ray tracing is a rapid and effective method for wave field calculation. Not only in the field of seismic-wave theory, but also in the field of seismic inversion and migration imaging,the seismic ray tracing method has become one of the most important methods. In anisotropic media, group velocity and phase velocity have different propagation directions. The seismic wave propagates along the direction of group velocity , it does not depend on the direction of phase velocity. Ray angle is a complex function with respect to phase angle, it is difficult to measure and calculate. But most rocks are weak anisotropic, so the expression of phase velocity can be simplified greatly. Based on the approximate expression of phase velocity this thesis for rotating axisymmetric weak anisotropic media deduces an expression of the partial derivative of phase velocity and an expression of group velocity with the method of linear approximation. This paper uses the fourth order Runge-Kutta method together with the two-dimensional interpolation and linear interpolation to obtain the parameters of the physical locations. At last the paths of seismic wave in rotating axisymmetric weak anisotropic media are computed. According to the analysis of the computational results, it indicates that the method developed in this paper has strong adaptability, high computational efficiency and high accuracy for rotating axisymmetric weak anisotropic media.

印度与欧亚板块强烈相互作用区深部结构特征研究

The central-south Tibet is a part of the products of the continental plate collision between Eurasia and India. To study the deep structure of the study area is significant for understanding the dynamics of the continental-continental collision. A 3-D density model matched well with the observations in the central-south Tibet was proposed in this study. In addition, this study has also used numerical simulation method to prove that Quasi-Love (QL) wave is deduced by anisotropy variation but not by lateral heterogeneity. Meanwhile, anisotropy variation in the upper mantle of the Qiangtang terrane and Lhasa terrane is detected by the QL waves observed in recorded seismograms. Based on the gravity modeling, some results are summarized as follows: 1) Under the constrain of geometrical structure detected by seismic data, a 3-D density model and Moho interface are proposed by gravity inversion of the central-south Tibet. 2) The fact that the lower crustal densities are smaller than 3.2 g/cm3, suggests absence of eclogite or partial eclogitization due to delamination under the central-south Tibet. 3) Seismicity will be strong or weak in the most negative Bouguer gravity anomaly. So there is no a certain relationship between seismicity and Bouguer gravity anomaly. 4) Crustal composition are determined after temperature-pressure calibration of seismic P wave velocity. The composition of lower crust might be one or a mixture of: 1. amphibolite and greenschist facies basalt beneath the Qiangtang terrane; 2. gabbro-norite-troctolite and mafic granulite beneath the Lhasa terrane. Because the composition of the middle crust cannot be well constrained by the above data set, the data set published by Rudnick & Fountain (1995) is used for comparison. It indicated the composition of the middle crust is granulite facies and might be pelitic gneisses.Granulite facies used to be interpreted as residues of partial melting, which coincidences with the previous study on partial melting middle crust. Amphibolite facies are thought to be produced after delamination, when underplating works in the rebound of the lower crust and lithospheric mantle. From the seismology study, I have made several followed conclusions: 1) Through the numerical simulation experiment of surface wave propagating in heterogeneity media, we can find that amplitude and polarization of surface wave only change a little when considering heterogeneity. Furthermore, it is proved that QL waves, generated by surface wave scattering, are caused by lateral variation of anisotropy but not by heterogeneity. 2) QL waves are utilized to determine the variation of uppermost mantle anisotropy of the Tibetan plateau. QL waves are identified from the seismograms of the selected paths recorded by the CAD station. The location of azimuth anisotropy gradient is estimated from the group velocities of Rayleigh wave, Love wave and QL wave. It suggests that south-north lateral variation of azimuthal anisotropy locates in Tanggula mountain, and east-west lateral variation in the north of Gandese mountain with 85°E longitude and near the Jinsha river fault with 85°E longitude.