Molecular cloning, expression of a big defensin gene from bay scallop Argopecten irradians and the antimicrobial activity of its recombinant protein

Antimicrobial peptides are important components of the host innate immune responses by exerting broad-spectrum microbicidal activity against pathogenic microbes. The first mollusk big defensin (designated AiBD) cDNA was cloned from bay scallop Argopecten irradians by expressed sequence tag (EST) and rapid amplification of cDNA ends (RACE) techniques. The scallop AiBD consisted of 531 nucleotides with a canonical polyadenylation signal sequence AATAAA and a poly(A) tail, encoding a polypeptide of 122 amino acids. The high similarity of AiBD deduced amino acid sequence with big defensin from Tachypleus tridentatus and Branchiostoma belcheri tsingtaunese indicated that AiBD should be a member of big defensin family. The expression of AiBD in various tissues was measured by using Northern blotting analysis. mRNA transcripts of AiBD could be detected in haemocytes of unchallenged scallops. The temporal expression of AiBD in haemolymph after Vibrio anguilarum challenge was recorded by quantitative real time PCR. The relative expression level of AiBD in haemolymph was up-regulated evenly in the first 8 h, followed by a drastic increase, and increased 131.1-fold at 32 h post-injection. These results indicated that AiBD could be induced by bacterial challenge, and it should participate in the immune responses of A. irradians. Biological activity assay revealed that recombinant AiBD could inhibit the growth of both Gram-positive and Gram-negative bacteria, and also showed strong fungicidal activity towards the expression host. Recombinant expression of AiBD made it possible to further characterize its functions involved in immune responses, and also provided a potential therapeutic agent for disease control in aquaculture. (c) 2006 Elsevier Ltd. All rights reserved.

cDNA cloning and mRNA expression of a selenium-dependent glutathione peroxidase from Zhikong scallop Chlamys farreri

The glutathione peroxidases are essential enzymes of the cellular antioxidant defence system. In the present study, the full-length cDNA sequence encoding an extracellular glutathione peroxidase (designated CfGPx3) was isolated from Zhikong scallop Chlamys farreri. The complete cDNA was of 1194 bp, containing a 5' untranslated region (UTR) of 50 bp, a 3' UTR of 490 bp and an open reading frame (ORF) of 654 bp encoding a polypeptide of 217 amino acids. CfGPx3 possessed all the conserved features critical for the fundamental structure and function of glutathione peroxidase, such as the selenocysteine encoded by stop codon UGA, the GPx signature motif ((96)LGVPCNQFI(103)) and the active site motif ((WNFEKF184)-W-179). The high similarity of CfGPx3 with GPx from other organisms indicated that CfGPx3 should be a new member of the glutathione peroxidase family. By fluorescent quantitative real-time PCR, the CfGPx3 mRNA was universally detected in the tissues of haemocytes, gill, gonad, muscle and hepatopancreas with the highest expression in hepatopancreas. After scallops were challenged by Listonella anguillarum, the expression level of CfGPx3 transcript in haemocytes was significantly up-regulated (P<0.05) at 8 h post challenge. These results suggested that CfGPx3 was potentially involved in the immune response of scallops and perhaps contributed to the protective effects against oxidative stress. (C) 2010 Elsevier Inc. All rights reserved.

Peptidoglycan recognition protein of Chlamys farreri (CfPGRP-S1) mediates immune defenses against bacterial infection

Peptidoglycan recognition protein (PGRP) is an essential molecule in innate immunity for both invertebrates and vertebrates, owing to its prominent ability in detecting and eliminating the invading bacteria. Several PGRPs have been identified from mollusk, but their functions and the underlined mechanism are still unclear. In the present study, the mRNA expression profiles, location, and possible functions of PGRP-S1 from Zhikong scallop Chlamys farreri (CfPG RP-St) were analyzed. The CfPGRP-S1 protein located in the mantle, gill, kidney and gonad of the scallops. Its mRNA expression in hemocytes was up-regulated extremely after PGN stimulation (P < 0.01), while moderately after the stimulations of LPS (P < 0.01) and beta-glucan (P < 0.05). The recombinant protein of CfPGRP-S1 (designated as rCfPGRP-S1) exhibited high affinity to PGN and moderate affinity to LPS, but it did not bind beta-glucan. Meanwhile, rCfPGRP-S1 also exhibited strong agglutination activity to Gram-positive bacteria Micrococcus luteus and Bacillus subtilis and weak activity to Gram-negative bacteria Escherichia coli. More importantly, rCfPGRP-S1 functioned as a bactericidal amidase to degrade PGN and strongly inhibit the growth of E. coli and Staphyloccocus aureus in the presence of Zn2+. These results indicated that CfPGRP-S1 could not only serve as a pattern recognition receptor recognizing bacterial PGN and LPS, but also function as a scavenger involved in eliminating response against the invaders. (C) 2010 Elsevier Ltd. All rights reserved.

The cDNA cloning and mRNA expression of heat shock protein 70 gene in the haemocytes of bay scallop (Argopecten irradians, Lamarck 1819) responding to bacteria challenge and naphthalin stress

Heat shock protein 70 (HSP70) is an important member of the heat shock protein superfamily, and it plays a key role in the process of protecting cells, facilitating the folding of nascent peptides and responding to stress. The cDNA of bay scallop Argopecten irradians HSP70 (designated AIHSP70) was cloned by the techniques of homological cloning and rapid amplification of cDNA end (RACE). The full length of AIHSP70 cDNA was 2651 bp in length, having a 5' untranslated region (UTR) of 96 bp, a 3' UTR of 575 bp, and an open reading frame (ORF) of 1980 bp encoding a polypeptide of 659 amino acids with an estimated molecular mass of 71.80 kDa and an estimated isoelectric point of 5.26. BLAST analysis revealed that the AIHSP70 gene shared high identity with other known HSP70 genes. Three classical HSP signature motifs were detected in AIHSP70 by InterPro, analysis. 3-D structural prediction of AIHSP70 showed that its N terminal ATPase activity domain and,C terminal substrate-binding domain shared high similarity with that in human heat shock protein 70. The results indicated that the AIHSP70 was a member of the heat shock protein 70 family. A semi-quantitive RT-PCR method was used to analyse the expression of AIHSP70 gene after the treatment of naphthalin which is one kind of polycyclic aromatic hydrocarbon (PAH) and the challenge of bacteria. mRNA expression of AIHSP70 in scallop was up-regulated significantly after the stimulation of naphthalin and increased with increasing naphthalin concentration. A clearly time-dependent expression pattern of AIHSP70 was observed after the scallops were infected by Vibrio anguillarum, and the mRNA expression reached a maximum level at 8 h and lasted to 16 h, and then dropped progressively. The results indicated that AIHSP70 could play an important role in mediating the environmental stress and immune response in scallop. (c) 2006 Elsevier Ltd. All rights reserved.

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.

Jiangella gansuensis gen. nov., sp nov., a novel actinomycete from a desert soil in north-west China

A novel actinomycete strain, designated YIM 002(T), was isolated from a desert soil sample in Gansu Province, north-west China. This actinomycete isolate formed well-differentiated aerial and substrate mycelia. In the early stages of growth, the substrate mycelia fragmented into short or elongated rods. Chemotaxonomically, it contained LL-2,6-diaminopimelic acid in the cell wall. The cell-wall sugars contained ribose and glucose. Phospholipids present were phosphatidylinositol mannosides, phosphatidylinositol and diphosphatidylglycerol. MK-9(H-4) was the predominant menaquinone. The major fatty acids were anteiso C-15:0 (35.92%), anteiso C-17:0 (15.84%), iso C-15:0 (10.40%), iso C-16:0 (7.07%) and C(17:10)w8c (9.37%). The G+C content of the DNA was 70 mol%. Phylogenetic analysis and signature nucleotide data based on 16S rRNA gene sequences showed that strain YIM 002(T) is distinct from all recognized genera of the family Nocardioidaceae in the suborder Propionibacterineae. On the basis of the phenotypic and genotypic characteristics, it is proposed that isolate YIM 002(T) be classified as a novel species in a new genus, Jiangella gansuensis gen. nov., sp. nov. The type strain is YIM 002(T) (= DSM 44835(T) = CCTCC AA 204001(T) = KCTC 19044(T)).

Pg震相走时的各向异性影响及各向异性结构成像

As the first arrival of seismic phase in deep seismic sounding, Pg is the important data for studying the attributes of the sedimentary layers and the shape of crystalline basement because of its high intensity and reliable detection. Conventionally, the sedimentary cover is expressed as isotropic, linear increasing model in the interpretation of Pg event. Actually, the sedimentary medium should be anisotropic as preferred cracks or fractures and thin layers are common features in the upper crust, so the interpretation of Pg event needs to be taken account of seismic velocity anisotropy. Traveltime calculation is the base of data processing and interpretation. Here, we only study the type of elliptical anisotropy for the poor quality and insufficiency of DSS data. In this thesis, we first investigate the meaning of elliptical anisotropy in the study of crustal structure and attribute, then derive Pg event’s traveltime-offset relationship by assuming a linear increasing velocity model with elliptical anisotropy and present the invert scheme from Pg traveltime-offset dataset to seismic velocity and its anisotropy of shallow crustal structure. We compare the Pg traveltime calculated by our analytic formula with numerical calculating method to test the accuracy. To get the lateral variation of elliptical anisotropy along the profiling, a tomography inversion method with the derived formula is presented, where the profile is divided into rectangles. Anisotropic imaging of crustal structure and attribute is efficient method for crust study. The imaging result can help us interprete the seismic data and discover the attribute of the rock to analyze the interaction between layers. Traveltime calculation is the base of image. Base on the ray tracing equations, the paper present a realization of three dimension of layer model with arbitrary anisotropic type and an example of Pg traveltime calculation in arbitrary anisotropic type is presented. The traveltime calculation method is complex and it only adapts to nonlinear inversion. Perturbation method of travel-time calculation in anisotropy is the linearization approach. It establishes the direct relation between seismic parameters and travetime and it is fit for inversion in anisotropic structural imaging. The thesis presents a P-wave imaging method of layer media for TTI. Southeastern China is an important part of the tectonic framework concerning the continental margin of eastern China and is commonly assumed to comprise the Yangtze block and the Cathaysia block, the two major tectonic units in the region. It’s a typical geological and geophysical zone. In this part, we fit the traveltime of Pg phase by the raytracing numerical method. But the method is not suitable here because the inefficiency of numerical method and the method itself. By the analytic method, we fit the Pg and Sg and get the lateral variation of elliptical anisotropy and then discuss its implication. The northeastern margin of Qinghai-Tibetan plateau is typical because it is the joint area of Eurasian plate and Indian plate and many strong earthquakes have occurred there in recent years．We use the Pg data to get elliptical anisotropic variation and discuss the possible meaning.

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.

空间电流体系地磁反演及理论研究

The space currents definitely take effects on electromagnetic environment and also are scientific highlight in the space research. Space currents as a momentum and energy provider to Geospace Storm, disturb the varied part of geomagnetic field, distort magnetospheric configuration and furthermore take control of the coupling between magnetosphere and ionosphere. Due to both academic and commercial objectives above, we carry on geomagnetic inverse and theoretical studies about the space currents by using geomagnetic data from INTERMAGNET. At first, we apply a method of Natural Orthogonal Components (NOC) to decomposition the solar daily variation, especially for (solar quiet variation). NOC is just one of eign mode analysis, the most advantage of this method is that the basic functions (BFs) were not previously designated, but naturally came from the original data so that there are several BFs usually corresponding to the process really happened and have more physical meaning than the traditional spectrum analysis with the fixed BFs like Fourier trigonometric functions. The first two eign modes are corresponding to the and daily variation and their amplitudes both have the seasonal and day-to-day trend, that will be useful for evaluating geomagnetic activity indices. Because of the too strict constraints of orthogonality, we try to extend orthogonal contraints to the non-orthogonal ones in order to give more suitable and appropriate decomposition of the real processes when the most components did not satisfy orthogonality. We introduce a mapping matrix which can transform the real physical space to a new mathematical space, after that process, the modified components which associated with the physical processes have satisfied the orthogonality in the new mathematical space, furthermore, we can continue to use the NOC decomposition in the new mathematical space, and then all the components inversely transform back to original physical space, so that we would have finished the non-orthogonal decomposition which more generally in the real world. Secondly, geomagnetic inverse of the ring current’s topology is conducted. Configurational changes of the ring current in the magnetosphere lead to different patterns of disturbed ground field, so that the global configuration of ring current can be inferred from its geomagnetic perturbations. We took advantages of worldwide geomagnetic observatories network to investigate the disturbed geomagnetic field which produced by ring current. It was found that the ring current was not always centered at geomagnetic equator, and significantly deviated off the equator during several intense magnetic storms. The deviation owing to the tilting and latitudinal shifting of the ring current with respect to the earth’s dipole can be estimated from global geomagnetic survey. Furthermore those two configurational factors which gave a quantitative description of the ring current configuration, will be helpful to improve the Dst calibration and understand the dependence of ring current’s configuration on the plasma sheet location relative to the equator when magnetotail field warped. Thirdly, the energization and physical acceleration process of ring current during magnetic storm has been proposed. When IMF Bz component increase, the enhanced convection electric field drive the plasma injection into the inner magnetosphere. During the transport process, a dynamic heating is happened which make the particles more ‘hot’ when the injection is more deeply inward. The energy gradient along the injection path is equivalent to a kind of force, which resist the plasma more earthward injection, as a diamagnetic effect of the magnetosphere anti and repellent action to the exotically injected plasma. The acceleration efficiency has a power law form. We use analytical way to quantitatively describe the dynamical process by introducing a physical parameter: energization index, which will be useful to understand how the particle is heated. At the end, we give a scheme of how to get the from storm time geomagnetic data. During intense magnetic storms, the lognormal trend of geomagnetic Dst decreases depend on the heating dynamic of magnetosphere controlling ring current. The descending pattern of main phase is governed by the magnetospheric configuration, which can be describled by the energization index. The amplitude of Dst correlated with convection electric field or south component of the solar wind. Finally, the Dst index is predicted by upstream solar wind parameter. As we known space weather have posed many chanllenges and impacts on techinal system, the geomagnetic index for evaluating the activity space weather. We review the most popular Dst prediction method and repeat the Dst forecasting model works. A concise and convnient Key Points model of the polar region is also introduced to space weather. In summary, this paper contains some new quantitative and physical description of the space currents with special focus on the ring current. Whatever we do is just to gain a better understanding of the natural world, particularly the space environment around Earth through analytical deduction, algorithm designing and physical analysis, to quantitative interpretation. Applications of theoretical physics in conjunction with data analysis help us to understand the basic physical process govering the universe.

数据驱动型多次波衰减方法的研究

In the practical seismic profile multiple reflections tend to impede the task of even the experienced interpreter in deducing information from the reflection data. Surface multiples are usually much stronger, more broadband, and more of a problem than internal multiples because the reflection coefficient at the water surface is much larger than the reflection coefficients found in the subsurface. For this reason most attempts to remove multiples from marine data focus on surface multiples, as will I. A surface-related multiple attenuation method can be formulated as an iterative procedure. In this essay a fully data-driven approach which is called MPI —multiple prediction through inversion (Wang, 2003) is applied to a real marine seismic data example. This is a pretty promising scheme for predicting a relative accurate multiple model by updating the multiple model iteratively, as we usually do in a linearized inverse problem. The prominent characteristic of MPI method lie in that it eliminate the need for an explicit surface operator which means it can model the multiple wavefield without any knowledge of surface and subsurface structures even a source signature. Another key feature of this scheme is that it can predict multiples not only in time but also in phase and in amplitude domain. According to the real data experiments it is shown that this scheme for multiple prediction can be made very efficient if a good initial estimate of the multiple-free data set can be provided in the first iteration. In the other core step which is multiple subtraction we use an expanded multi-channel matching filter to fulfil this aim. Compared to a normal multichannel matching filter where an original seismic trace is matched by a group of multiple-model traces, in EMCM filter a seismic trace is matched by not only a group of the ordinary multiple-model traces but also their adjoints generated mathematically. The adjoints of a multiple-model trace include its first derivative, its Hilbert transform and the derivative of the Hilbert transform. The third chapter of the thesis is the application for the real data using the previous methods we put forward from which we can obviously find the effectivity and prospect of the value in use. For this specific case I have done three group experiments to test the effectiveness of MPI method, compare different subtraction results with fixed filter length but different window length, invest the influence of the initial subtraction result for MPI method. In terms of the real data application, we do fine that the initial demultiple estimate take on a great deal of influence for the MPI method. Then two approaches are introduced to refine the intial demultiple estimate which are first arrival and masking filter respectively. In the last part some conclusions are drawn in terms of the previous results I have got.

离散裂缝诱导各向异性粘弹性地震波模拟与特征分析

The real media always attenuate and distort seismic waves as they propagate in the earth. This behavior can be modeled with a viscoelastic and anisotropic wave equation. The real media can be described as fractured media. In this thesis, we present a high-order staggered grid finite-difference scheme for 2-D viscoelastic wave propagation in a medium containing a large number of small finite length fractures. We use the effective medium approach to compute the anisotropic parameters in each grid cell. By comparing our synthetic seismogram by staggered-grid finite-difference with that by complex-ray parameter ray tracing method, we conclude that the high-order staggered-grid finite-difference technique can effectively used to simulate seismic propagation in viscoelastic-anisotropic media. Synthetic seismograms demonstrate that strong attenuation and significant frequency dispersion due to viscosity are important factors of reducing amplitude and delaying arrival time varying with incidence angle or offset. On the other hand, the amount of scattered energy not only provides an indicator of orientation of fracture sets, but can also provide information about the fracture spacing. Analysis of synthetic seismograms from dry- and fluid-filled fractures indicates that dry-filled fractures show more significant scattering on seismic wavefields than fluid-filled ones, and offset-variations in P-wave amplitude are observable. We also analyze seismic response of an anticlinal trap model that includes a gas-filled fractured reservoir with high attenuation, which attenuates and distorts the so-called bright spot.

鄂尔多斯盆地志丹地区延长组储层综合评价

Sedimentary provenance direction，sedimentary facies，reservoir geological characteristic，pore structure； physical property characteristic，reservoir classification and evaluation ，forthermore，favorable area area are forecasted of Yanchang formation in ZhiDan region， by mainly study on the data of field section observation ，core observation, well logging explaination and routine microscope slice identification，scanning Electron Microscope and reservoir analysis of lithology and physical property , Under the guidance of such advanced theories and methods as sedimentology，reservoir sedimentology，lithological oil pool and so on，in the thesis. The following fruits are mainly achieved in this paper： Yanchang formation stratum is divided and correlated in this entire region， and the characteristic of oil layer unit is detailed discussed , respectively. According to main marker bed and supplementary ones.and research result shows that the source of provenance direction during Yanchang Formation mianly is north-east. Delta and lake are mainly developed in study area ，sub-facies and micro-facies are divided，distribution of sedimentary micro-facies in plane and palaeogeographic evolution are described，and gentle slope type- shallow water delta depositional model is established. Fine-grain arkose sandstone is the main reservoir，and which have experienced such different degree diagenesis as compaction, cementation, replacement and dissolution, and in which compaction and cementation are mainly factors to reduce sandstone physical property and dissolution effectively improved physical property during burial diagenesis procedure. All reservoirs of Yanchang Formation have entered A period of late diagenetic stage according to scheme of diagenesis period division . Intergranular porosity，dissolution porosity，fissure porosity are main pore types. And porosity structure are analyse by mercury penetration capillary pressure parameter，fine-shortness type and fine- length throat type are mainly developed. as a whole，the reservoir， with the characteristic of porosity and permeability altering apparently，strong inhomogeneity ， is a medium- porosity and medium permeability one. In plane，higher- porosity and higher-permeability are corresponded well with distributary channel area， physical property and inhomogeneity are affected by both deposition and diagenesis，and distributary channel and underwater distributary channel are favorable facies . According to such characteristic as lithology，physical property，pore structure ，diagenesis and sandstone distribution, the sandy reservoir can be classified 4 types, and the main sandy in every oil layer unit are evaluated according to the standard. The analysis result of petroleum concentration rule shows that Yanchang Formation are with not only favourable oil source rock，reservoir，covering combination ，but also good entrapment condition in study area. Lithology and structure-lithology oil pool are mainly developed ，based on condition of favorable reservoir developments，accounting for deliverability and sandstone superface elevation，zone of profitabilitis are forecasted.

二氧化碳在孔隙咸水层中多相流动数值模拟研究

There has been a growing concern about the use of fossil fuels and its adverse effects on the atmospheric greenhouse and ecological environment. A reduction in the release rate of CO2 into the atmosphere poses a major challenge to the land ecology of China. The most promising way of achieving CO2 reduction is to dispose of CO2 in deep saline aquifers. Deep aquifers have a large potential for CO2 sequestration in geological medium in terms of volume and duration. Through the numerical simulation of multiphase flow in a porous media, the transformation and motion of CO2 in saline aquifers has been implemented under various temperature and hydrostatic pressure conditions, which plays an important role to the assessment of the reliability and safety of CO2 geological storage. As expected, the calculated results can provide meaningful and scientific information for management purposes. The key problem to the numerical simulation of multiphase flow in a porous media is to accurately capture the mass interface and to deal with the geological heterogeneity. In this study, the updated CE/SE (Space and time conservation element and solution element) method has been proposed, and the Hybrid Particle Level Set method (HPLS) has extended for multiphase flows in porous medium, which can accurately trace the transformation of the mass interface. The benchmark problems have been applied to evaluate and validate the proposed method. In this study, the reliability of CO2 storage in saline aquifers in Daqingzi oil field in Sunlong basin has been discussed. The simulation code developed in this study takes into account the state for CO2 covering the triple point temperature and pressure to the supercritical region. The geological heterogeneity has been implemented, using the well known geostatistical model (GSLIB) on the base of the hard data. The 2D and 3D model have been set up to simulate the CO2 multiphase flow in the porous saline aquifer, applying the CE/SE method and the HPLS method .The main contents and results are summarized as followings. (1) The 2D CE/SE method with first and second –order accuracy has been extended to simulate the multiphase flow in porous medium, which takes into account the contribution of source and sink in the momentum equation. The 3D CE/SE method with the first accuracy has been deduced. The accuracy and efficiency of the proposed CE/SE method have been investigated, using the benchmark problems. (2) The hybrid particle level set method has been made appropriate and extended for capturing the mass interface of multiphase flows in porous media, and the numerical method for level set function calculated has been formulated. (3) The closed equations for multiphase flow in porous medium has been developed, adept to both the Darcy flow and non-Darcy flow, getting over the limitation of Reynolds number to the calculation. It is found that Darcy number has a decisive influence on pressure as well as velocity given the Darcy number. (4) The new Euler scheme for numerical simulations of multiphase flows in porous medium has been proposed, which is efficient and can accurately capture the mass interface. The artificial compressibility method has been used to couple the velocities and pressure. It is found that the Darcy number has determinant effects on the numerical convergence and stability. In terms of the different Darcy numbers, the coefficient of artificial compressibility and the time step have been obtained. (5) The time scale of the critical instability for critical CO2 in the saline aquifer has been found, which is comparable with that of completely CO2 dissolved saline aquifer. (6) The concept model for CO2 multiphase flows in the saline aquifer has been configured, based on the temperature, pressure, porosity as well as permeability of the field site .Numerical simulation of CO2 hydrodynamic trapping in saline aquifers has been performed, applying the proposed CE/SE method. The state for CO2 has been employed to take into account realistic reservoir conditions for CO2 geological sequestration. The geological heterogeneity has been sufficiently treated , using the geostatistical model. (7) It is found that the Rayleigh-Taylor instability phenomenon, which is associated with the penetration of saline fluid into CO2 fluid in the direction of gravity, has been observed in CO2 multiphase flows in the saline aquifer. Development of a mushroom-type spike is a strong indication of the formation of Kelvin-Helmholtz instability due to the developed short wavelength perturbations present along the interface and parallel to the bulk flow. Additional key findings: the geological heterogeneity can distort the flow convection. The ascending of CO2 can induce the persistent flow cycling effects. The results show that boundary conditions of the field site have determinant effects on the transformation and motion of CO2 in saline aquifers. It is confirmed that the proposed method and numerical model has the reliability to simulate the process of the hydrodynamic trapping, which is the controlling mechanism for the initial period of CO2 storage at time scale of 100 years.

东喜马拉雅构造结新生代碰撞变形与岩石圈演化

The Eastern Himalayan Syntaxis (EHS) is one of the strongest deformation area along the Himalayan belt resulted from the collision between Indian plate and the Eurasian Plate since the 50～60Ma, and has sensitivity tracked and preserved the whole collisional processes. It should depend on the detail geological investigations to establish the deformational accommodate mode, and the uplift history, to elucidate the deep structure and the crust-mantle interaction of the Tibet Plateau of the EHS. The deep-seated (Main Mantle Thrusts) structures were exhumed in the EHS. The MMT juxtapose the Gangdese metamorphic basement and some relic of Gangdese mantle on the high Himalayan crystalline series. The Namjagbawa group which is 1200～1500Ma dated by U/Pb age of zircon and the Namla group which is 550Ma dated by U/Pb age of zircon is belong to High Himalayan crystalline series and Gangdese basement respectively. There is some ophiolitic relic along the MMT, such as metamorphic ocean mantle peridotite and metamorphic tholeiite of the upper part of ocean-crust. The metamorphic ocean mantle peridotites (spinel-orthopyroxene peridotite) show U type REE patterns. The ~(87)Sr/~(86)Sr ratios were, 0.709314～0.720788, and the ~(143)Nd/~(144)Nd ratios were 0.512073～0.512395, plotting in the forth quadrant on the ~(87)Sr/~(86)Sr-~(143)Nd/~(144)Nd isotope diagram. Some metamorphic basalt (garnet amphibolite) enclosures have been found in the HP garnet-kynite granulite. The garnet amphibolites can be divided two groups, the first group is deplete of LREE, and the second group is flat or rich LREE, and their ~(87)Sr/~(86)Sr, ~(143)Nd/~(144)Nd ratios were 0.70563～0.705381 and 0.512468～0.51263 respectively. Trace element and isotopic characteristics of the garnet amphibolites display that they formed in the E-MORB environment. Some phlogolite amphibole harzburgites, which exhibit extensive replacement by Phl, Amp, Tc and Dol etc, were exhumed along the MMT. The Phl-Amp harzburgites are rich in LREE and LILE, such as Rb, K etc, and depletes Eu (Eu~* = 0.36 ～ 0.68) and HFSE, such as Nb, Ta, Zr, Hf, P, Ti etc. The trace element indicate that the Phl-Amp harzburgites have island arc signature. Their ~(87)Sr/~(86)Sr are varied from 0.708912 to 0.879839, ~(143)Nd/~(144)Nd from 0.511993 to 0.512164, ε Nd from- 9.2 to - 12.6. Rb/Sr isochrone age of the phlogolite amphibole harzburgite shows the metasomatism took place at 41Ma, and the Amp ~(40)Ar/~(39)Ar cooling age indcate the Phl-Amp harzburgite raising at 16Ma. There is an intense crust shortening resulted from the thrust faults and folds in the Cayu block which is shortened more 120km than that of the Lasha block in 35～90Ma. With the NE corner of the India plate squash into the Gangdese arc, the sinistral Pai shear fault and the dextral Aniqiao shear fault on the both sides of the Great bent of Yalun Zangbu river come into active in 21～26Ma. On the other hand, the right-lateral Gongrigabu strike-slip faults come into activity at the same period, a lower age bound for the Gongrigabu strike-slip fault is estimated to be 23～24Ma from zircon of ion-probe U/Pb thermochronology. The Gongrigabu strike-slip faults connect with the Lhari strike-slip fault in the northwestern direction and with the Saganing strike-slip at the southeastern direction. Another important structure in the EHS is the Gangdese detachment fault system (GDS) which occurs between the sedimental cover and the metamorphic basement. The lower age of the GDS is to be 16Ma from the preliminary 40Ar/39Ar thermochronology of white mica. The GDS is thought to be related to the reverse of the subducted Indian crust and the fast uplift of the EHS. Structural and thermochronology investigation of the EHS suggest that the eastern Tibet and the western Yunnan rotated clockwise around the EHS in the period of 35～60Ma. Later, the large-scale strike-slip faults (RRD, Gaoligong and Saganing fault) prolongate into the EHS, and connect with the Guyu fault and Gongrigabu fault, which suggest that the Indianchia block escape along these faults. Two kind of magmatic rocks in the EHS have been investigated, one is the mantle-derived amphibole gabbro, dioposide diorite and amphibole diorite, another is crust origin biotit-garnet adamellite, biotit-garnet granodiorite and garnet-amphibole-biotite granite. The amphibole gabbro dioposite diorite and amphibole diorite are rich in LREE, and LILE, such as Ba, Rb, Th, K, Sr etc, depleted in HFSE, such as Nb, Ta, Zr, Hf, Ti etc. The ratio of ~(87)Sr/~(86)Sr are from 0.7044 to 0.7048, ~(143)Nd/~(144)Nd are from 0.5126 to 0.5127. The age of the mantle origin magamatic rocks, which result from the partial melt of the raising and decompression anthenosphere, is 8Ma by ~(40)Ar/~(39)Ar dating of amphibole from the diorite. The later crust origin biotite-garnet adamellite, biotite-garnet granodiorite and garnet-amphibole-biotite granite are characterized by aboudance in LREE, and strong depletion of Eu. The ratios of ~(87)Sr-~(86)Sr are from 0.795035 to 0.812028, ~(143)Nd/~(144)Nd from 0.51187 to 0.511901. The ~(40)Ar/~(39)Ar plateau age of the amphibole from the garnet-amphibole-biotite granite is 17.5±0.3Ma, and the isochrone age is 16.8±0.6Ma. Their geochemical characteristics show that the crust-derived magmatic rocks formed from partial melting of the lower curst in the post-collisional environment. A group of high-pressure kaynite-garnet granulites and enclave of high-pressure garnet-clinopyroxene grnulites and calc-silicate grnulites are outcroped along the MMT. The peak metamorphic condition of the high-pressure granulites yields T=800～960 ℃, P=1.4～1.8Gpa, corresponding the condition of 60km depth. The retrograde assemblages of the high-pressure grnulites occur at the condition of T=772.3～803.3 ℃, P=0.63～0.64Gpa. The age of the peak metamorphic assemblages are 45 ～ 69Ma indicated by the zircon U/Pb ion-plobe thermochronology, and the retrograde assemblage ages are 13～26Ma by U/Pb, ~(40)Ar/~(39)Ar thermochronology. The ITD paths of the high-pressure granulites show that they were generated during the tectonic thickening and more rapid tectonic exhumation caused by the subducting of the Indian plate and subsequent break-off of the subducted slab. A great deal of apatite, zircon and sphene fission-track ages, isotopic thermochronology of the rocks in the EHS show that its rapid raising processes of the EHS can be divided into three main periods. There are 35～60Ma, 13～25Ma, 0～3Ma. 3Ma is a turn in the course of raising in the EHS which is characterized by abruptly acceleration of uplifting. The uplift ratios are lower than 1mm .a~(-1) before 3Ma, and higher than 1mm .a~(-1) with a maximum ratio of 30mm .a~(-1) since 3Ma. The bottom (knick point) of the partial anneal belt is 3.8km above sea level in the EHS, and correspond to age of 3Ma determined by fission-track age of apatite. The average uplift ratio is about 1.4 mm .a~(-1) below the knick point. The EHS has raised 4.3km from the surface of 2.36km above sea level since 3Ma estimated by the fossil partial anneal belt of the EHS. We propose a two-stage subduction model (B+A model) basing on Structural, thermochronological, magmatical, metamorphic and geophysical investigations of the EHS. The first stage is the subduction of the Indian continental margin following after the subduction of the Tethys Ocean crust and subsequent collision with the Gangdese arc, and the second stage is the Indian crust injecting into the lower crust and upper mantle of the Tibet plateau. Slab break-off seems to be occurred between these two stages.