Diel vertical migration of zooplankton following optimal food intake under predation

A model is developed to investigate the trade-offs between benefits and costs involved in zooplanktonic diel vertical migration (DVM) strategies. The 'venturous revenue' (VR) is used as the criterion for optimal trade-offs. It is a function of environmental factors and the age of zooplankter. During vertical migration, animals are assumed to check instantaneously the variations of environmental parameters and thereby select the optimal behavioral strategy to maximize the value of VR, i.e. taking up as much food as possible with a certain risk of mortality. The model is run on a diel time scale (24 h) in four possible scenarios during the animal's life history. The results show that zooplankton can perform normal DVM balancing optimal food intake against predation risk, with the profile of DVM largely modified by the age of zooplankter.

Control mechanisms of diel vertical migration: Theoretical assumptions

We explore control mechanisms underlying the vertical migration of zooplankton in the water column under the predator-avoidance hypothesis. Two groups of assumptions in which the organisms are assumed to migrate vertically in order to minimize realized or effective predation pressure (type-I) and to minimize changes in realized or effective predation pressure (type-II), respectively, are investigated. Realized predation pressure is defined as the product of light intensity and relative predation abundance and the part of realized predation pressure that really affects organisms is termed as effective predation pressure. Although both types of assumptions can lead to the migration of zooplankton to avoid the mortality from predators, only the mechanisms based on type-II assumptions permit zooplankton to undergo a normal diel vertical migration (morning descent and evening ascent). The assumption of minimizing changes in realized predation pressure is based on consideration of DVM induction only by light intensity and predators. The assumption of minimizing changes in effective predation pressure takes into account, apart from light and predators also the effects of food and temperature. The latter assumption results in the same expression of migration velocity as the former one when both food and temperature are constant over water depth. A significant characteristic of the two type-II assumptions is that the relative change in light intensity plays a primary role in determining the migration velocity. The photoresponse is modified by other environmental variables: predation pressure, food and temperature. Both light and predation pressure are necessary for organisms to undertake DVM. We analyse the effect of each single variable. The modification of the phototaxis of migratory organisms depends on the vertical distribution of these variables. (C) 2001 Academic Press.

Feeding activities of zooplankton in the Bohai Sea

The Bohai Sea was the site of the Chinese national GLOBEC programme. During the June 1997 cruises of R/V Science No.1, observations and experiments on zooplankton feeding were conducted. At five 48 h time-series stations the following observations and measurements on zooplankton were carried out: (1) diurnal vertical migration, by collecting samples at different layers every 3 h with a closing net; (2) diurnal feeding rhythms, by gut pigment analysis; and (3) ingestion rate, by both gut pigment analysis and the dilution method. A classification by body size was used to deal with the diversity of species and developmental stages of zooplankton assemblages. Samples were separated into three size groups: small (200-500 mu m), medium (500-1000 mu m) and large (> 1000 mu m). The results showed that the copepods (Calanus sinicus, Paracalanus parvus, Acartia bifilosa and Centropages mcmurrichi) performed clear diurnal vertical migrations. However, their behaviour was different at different stations. The variation in gut pigment content over the 24 h cycle showed strong diurnal feeding rhythms, particularly for the large size group. Gut pigment contents reached their daily maximum during the time from dusk to midnight (18:00-24:00). The peak value was about 10 times the minimum observed in the daytime. The in situ daily grazing rate, based on gut pigment contents and evacuation experiments, was 4.00-12.65 ng chla ind(-1) day(-1) for the small size group, 5.99-66.58 ng chla ind(-1) day(-1) for the medium size group and 31.31-237.13 ng chla ind(-1) day(-1) for the large size group. The copepods consumed only a small part (2.90-13.52%) of the phytoplankton biomass hut about 77% of the daily production. The grazing mortality rate of phytoplankton by microzooplankton (<200 mu m) measured by the dilution method ranged from 0.43 to 0.69 day(-1) The calculated daily consumption of phytoplankton biomass was 35-50%, and 85-319% of the potential production.

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.

高阶广义屏波动方程叠前深度偏移

Attaining sufficient accuracy and efficiency of generalized screen propagator and improving the quality of input gathers are often problems　of wave equation presack depth migration,　in this paper,a high order formula of generalized screen propagator　for one-way wave equation is proposed by using the asymptotic expansion of single-square-root operator. Based on the formula，a new generalized screen propagator is developed ,which is composed of split-step Fourier propagator and high order correction terms，the new generalized screen propagator not only improving calculation precision without sharply increasing the quantity of computation，facilitates the suitability of generalized screen propagator to the media with strong lateral velocity variation.　As wave-equation prestack depth migration is sensitive to the quality of input gathers, which greatly affect the output，and the available seismic data processing system has inability to obtain traveltimes corresponding to the multiple arrivals, to estimate of great residual statics， to merge seismic datum from different projects and to design inverse Q filter, we establish　difference equations with an embodiment of Huygens’s principle for obtaining traveltimes corresponding to the multiple arrivals,bring forward a time variable matching filter for seismic datum merging by using the fast algorithm called Mallat tree for wavelet transformations， put forward a method for estimation of residual statics by applying the optimum model parameters estimated by iterative inversion with three organized algorithm,i.e,the CMP intertrace cross-correlation algorithm,the Laplacian image edge extraction algorithm,and the DFP algorithm, and present phase-shift inverse Q filter based on Futterman’s amplitude and phase-velocity dispersion formula and wave field extrapolation theory. All of their numerical and real data calculating results shows that our theory and method are practical and efficient. Key words: prestack depth migration, generalized screen propagator, residual statics，inverse Q filter ,traveltime，3D seismic datum mergence

高邮凹陷陈堡及陈堡东地区油气成藏规律研究

With the development of petroleum exploration in Gaoyou Depression, both old and new areas have been the active exploration targets, so the study of petroleum accumulation is significant to the petroleum exploration in the study area and the integrated oil and gas accumulation theory. Based on hydrocarbon accumulation theory and systematical research methods and combined with the structural characteristics of Gaoyou Depression, Chenbao and East of Chenbao were selected as the study areas in this dissertation, oil and gas migration pathways, accumulation periods, as well as accumulation models were studied, and favorable exploration targets were proposed. There develop three sets source rocks, which are Tai-2 Member, Fu-2 Member and Fu-4 Member respectively. Tai-2 Member is the predominant source rock in the eastern part. Fu-2 Member mainly occurs in the northern slope, while Fu-4 Member develops in the deep depression. In the study area, oil mainly comes from Fu-2 Member of Liuwushe subsag. The lower limit of TOC is 0.4%, and active source rock mostly distributed in the south fault-step zone. The source rock in Liuwushe subsag began to generate hydrocarbon in the late of Dainan depositional stage and the threshold was 2300m. The macro and micro characteristics of reservoirs and the reservoir heterogeneity characteristics of the Fu-1 Member were studied systematicly. The results show that Fu-1 Member, which has better reservoir properties, are medium porosity-medium permeability reservoir. The reservoir permeability has good correlation with porosity connectivity. The reservoirs have strong dissolution, pores are mainly thin to medium throat, and throat radii are distributed concentratedly, the sorting is good and pore structures are homogeneous. Sandstone reservoirs whether in the plan view, interlayer or in layers have a certain degree of heterogeneity, in particular, the heterogeneity in layers directly affect and control the oil and gas migration and accumulation. By analyzing the lithology correlation of the fault walls, shale smear, cross section stress, the configuration of fracture active periods and hydrocarbon generation and expulsion periods and fuzzy comprehensive evaluation, the main faults sealing were evaluated. The results show that the faults in Chenbao and East of Chenbao had poor sealing properties in Sanduo period and could be used as the migration pathways at that time. After Sanduo period, the tectonic stress fields in the area changed largely, and, consequently, the fault properties converted from tensional shear to compressive shear, the faults changed progressively from close to open, so the faults sealing became better and were conducive to the preservation of oil and gas reservoirs. According to the seismic event suspension modes and profile configurations above and under the unconformities, combined with tectonic evolutions of the study areas, the unconformity types can be classified into truncation unconformity, overlapped unconformity and parallel unconformity and the distribution characteristics of unconformities in the plan view was also studied. The unconformity structure was divided into basal conglomerate, weathered clay and semi-weathered layer vertically in the study area and this kind of structure make unconformities to be effective oil and gas migration pathways and is significant to hydrocarbon accumulation in a parts of areas. With the analyses of typical oil and gas reservoirs in the study area, combined with the research results of pathway systems, hydrocarbon accumulation models were established and the oil and gas accumulation laws in Chenbao and East of Chenbao analyzed. The oil and gas came from Liuwushe subsag and Liuliushe subsag. The oil and gas from Liuwushe subsag mainly migrated from the structural high parts into the fault-step zone along strata in northeast direction, a part of them migrated upward into the fault-step zone and the Wubao low uplift along Wu-1 Fault in northeast direction. The oil and gas from Liuliushe subsag mainly migrated into the upper reservoirs through Wu-2 fault, and lesser oil and gas migrated into the fault-step zone because of the controls of cross-section orientation, depression center and the hydrocarbon formation tendency. The favorable exploration targets in Chenbao and East of Chenbao have been concluded: the southern fault-step zone is a favorable oil and gas accumulation zone of Liuwushe subsag, and they are fault block reservoirs where fault acted as the barriers, the main target intervals are Fu-1 Member and Fu-3 Member in palaeocene; Oil and gas in the middle and northern fault-step zone mainly laterally migrated from the south areas, and the main target interval is Fu-3 Member in palaeocene; Fu-1 Member and the reserviors above the Wubao subsag are still the focuses in future explorations. The results of this study have important guiding significance for the future oil and gas exploration.

叠前地震数据规则化、重建及噪音压制

I address of reconstruction of spatial irregular sampling seismic data to regular grids. Spatial irregular sampling data impairs results of prestack migration, multiple attenuations, spectra estimation. Prestack 5-D volumes are often divided into sub-sections for further processing. Shot gathers are easy to obtain from irregular sampling volumes. My strategy for reconstruction is as follows: I resort irregular sampling gathers into a form of easy to bin and perform bin regularization, then utilize F-K inversion to reconstruct seismic data. In consideration of poor ability of F-K regularization to fill in large gaps, I sort regular sampling gathers to CMP and proposed high-resolution parabolic Radon transform to interpolate data and extrapolate offsets. To strong interfering noise--multiples, I use hybrid-domain high-resolution parabolic Radon transform to attenuate it. F-K regularization demand ultimately for lower computing costs. I proposed several methods to further improve efficiency of F-K inversion: first I introduce 1D and 2D NFFT algorithm for a rapid calculation of DFT operators; then develop fast 1D and 2D CG method to solve least-square equations, and utilize preconditioner to accelerate convergence of CG iterations; what’s more, I use Delaunay triangulation for weight calculation and use bandlimit frequency and varying bandwidth technique for competitive computation. Numerical 2D and 3D examples are offered to verify reasonable results and more efficiency. F-K regularization has poor ability to fill in large gaps, so I rearrange data as CMP gathers and develop hybrid-domain high-resolution parabolic Radon transforms which be used ether to interpolate null traces and extrapolate near and far offsets or suppress a strong interfere noise: multiples. I use it to attenuate multiples to verify performances of our algorithm and proposed routines for industrial application. Numerical examples and field data examples show a nice performance of our method.

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 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.

利用接收函数方位变化研究地壳各向异性

Because of its sensitivity to the velocity discontinuity of the earth, receiver function technique has become a routine procedure used to probe interior structure of the earth. Receiver functions contain anisotropic information of the earth’s interior, however, traditional receiver function techniques such as migration imaging and waveform inversion method, which are based on isotropic media assumption, can not effectively extract the anisotropy information contained in the azimuth variation pattern. Only by using the anisotropic media, e.g. a model with symmetric axis of arbitrary orientation, computing the response, can we obtain the detailed anisotropy information hidden in the radial and transversal receiver function. Focusing on the receiver function variation pattern changing wtih different back azimuths, we introduced different kinds of symmetric systems of seismic anisotropy used often, and summarized some possible causes of anisotropy formation. We show details about how to calculate the response of a stratified anisotropy model with symmetric axis of arbitrary orientation. We also simulated receiver functions among different models and analyzed how the changing of anisotropic parameters influence the azimuth variation pattern of receiver functions. The anisotropy study by receiver function analysis was applied to Taihang Mountain Range (TMR) in North China in this thesis. The maximum entropy spectrum deconvolution technique was used to extract radial and transversal receiver functions from the waveforms of 20 portable seismic stations deployed in TMR. Considering the signal-to-noise ratio and the azimuth coverage, we got the variation pattern of receiver functions for 11 stations. After carefully analyzing the pattern of the receiver functions that we got, we obtained the reliable evidence on the existence of anisotropy in the shallow crust in TMR. Our results show that, although the thickness of the upper crustal layer is only about 1 km, the layer shows a strong anisotropy with magnitude of 8~15%; in the deeper of crust, the magnitudes of anisotropy is about 3%~5%, showing a pattern with fast-symmetric-axis. The crust anisotropy beneath TMR in North China obtained in this study also shows a significant difference in both the lateral and vertical scale, which might imply a regional anisotropy characteristic in the studied region.

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

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.

乍得西南部中新生代火山作用与盆地演化

Located in the Central and West African, Chad, which is not well geological explored, is characterized by Mesozoic- Cenozoic intra-continental rift basins. The boreholes exposed that, during Mesozoic-Cenozoic times, volcanic activities were intense in these basins, but study on volcanic rocks is very weak, especially on those embedded in rift basins, and so far systematic and detailed work has still no carried out. Based on the project of China National Oil and Gas Exploration and Development Corporation, “The analysis of reservoir condition and the evaluation of exploration targets of seven basins in block H in Chad”, and the cooperative project between Institute of Geology and Geophysics, CAS and CNPC International (Chad) Co. Ltd., “Chronology and geochemistry studies on Mesozoic-Cenozoic volcanic rocks from southwestern Chad Basins”, systematic geochronology, geochemistry and Sr-Nd-Pb isotopic geochemistry studies on volcanic rocks from southwestern Chad basins have been done in the thesis for the first time. Detailed geochronological study using whole-rock K-Ar and Ar-Ar methods shows the mainly eruption ages of these volcanic rocks are Late Cretaceous- Paleogene. Volcanic rocks in the well Nere-1 and Figuier-1 from Doba basin are products of the Late Cretaceous which majority of the K-Ar (Ar-Ar) ages fall in the interval 95-75 Ma, whereas volcanic rocks in the well Ronier-1 from Bongor Basin and the Well Acacia-1 from Lake Chad Basin formed in the Paleogene which the ages concentrated in 66-52Ma. Two main periods of volcanic activity can be recognized in the study area, namely, the Late Cretaceous period and the Paleogene period. Volcanic activities have a general trend of south to north migration, but this may be only a local expression, and farther future studies should be carried on. Petrology study exhibits these volcanic rocks from southwestern Chad basins are mainly tholeiitic basalt. Major- and trace elements as well as Sr-Nd-Pb isotopic geochemistry studies show that the late Cretaceous and the Paleogene basalts have a definitely genetic relationship, and magmas which the basalts in southwestern Chad basins derived from were produced by fractional crystallization of olivine and clinopyroxene and had not do suffered from crustal contamination. These basalts are prominently enriched light rare earth elements (LREE), large-ion lithophile elements (LILE) and high field strength elements (HFSE) and depleted compatible elements. They have positive Ba, Pb, Sr, Nb, Ta, Zr, Hf anomalies and negative Th, U, P,Y anomalies. It is possible that the basalts from southwestern Chad basins mainly formed by mixing of depleted mantle (DM) and enriched mantle (EMⅡ) sources. The late Cretaceous basalts have higher (87Sr/86Sr)i ratios than the Paleogene basalts’, whereas have lower (143Nd/144Nd)i ratios than the latter, showing a significant temporal evolution. The mantle sources of the Late Cretaceous basalts may have more enriched mantle(EMⅡ) compositions, whereas those of the Paleogene basalts are relatively more asthenospheric mantle (DM) components. The mantle components with temporal change observed in basalts from Chad basins were probably correlated with the asthenospheric mantle upwelling and lithospheric thinning in Central and Western Africa since Mesozoic. Mesozoic- Cenozoic Volcanism in Chad basins probably is a product of intra- plate extensional stress regime, corresponded to the tectonic setting of the whole West and Central African during Cretaceous. Volcanism is closely correlated with rifting. As time passed from early period to late, the basaltic magma of Chad basins, characterized with shallower genetic depth, higher density and smaller viscosity, probably indicates the gradual strengthening evolution of the rifting. In the initial rife stage, volcanic activities are absent in the study area. Volcanic activities are basiccally corresponded with the strong extensional period of Chad basins, and the eruption of basalts was slightly lagged behind the extensional period. In the post-rift stage (30-0Ma), these basins shifted to the thermal sag phase, volcanic activities in the study area significantly decreased and then terminated.

成矿流体构造动力学研究—以冀西北后沟金矿和山东七宝山金矿为例

Tectonic dynamics of metallogenetic fluids is a new crossed subjects among fluid geology, mineral deposit geology and structural geology, and is one of the major current projects of geosciences. It is mainly focused on structures and tectonic dynamic induced by fluid motion, variation of physical condition of fluids (such as temperature and pressure), and interaction between chemical component of fluids and wall rocks in the crust. It takes features of deformation and metamorphysim, which formed during interaction between fluids and rocks and have been perserved in rocks, as basic research objects. After studying types, orders, distributions and fabrics of these features, and analyzing and testing physical and chemical information from these features by some techniques, it is intended to reconstruct moving process of fluids, dynamics of interaction between fluids and rocks, and dynamics of mineralizations. Three problems of tectonic dynamics of metallogenetic fluids, which have not been paid much attentions before, have been studied and discussed in this report. Three relative topics are including: 1)Double-fracturing induced by thermal stress and pressure of fluids and mineralization of Gold-copper in Breccia Pipe at the Qibaoshan in Shandong Province; 2)Parting structures induced by K-metasomatism in the Hougou area, northwestern Heibei province; 3)Migration mechanism of dissolved mass in Fe&S-rich fluids in Hougou gold deposit in Heibei province. After a synthetical study of two years, the author has made some new processes and progresses. The main new advances can be summaried as the following: 1)Thermal stress of fluids formed by temperature difference between fluids and country rock, during upword migration process of fluids with high temperature and pressure, can make rock to break, and some new fractures, which surfaces were uasally dry, formed. The breccia pipe at the Qibaoshan area in Shandong province has some distinct texture of fluidogenous tectonics, the breccia pipe is caused by double-fracturing induced by thermal stress and pressure, distribution of gold-corpper ore bodies are controlled powerfully by fluidogenous tectonics in the breccia pipe. 2)The author discovered a new kind of parting structures in K-alterated rocks in the northwestern part of Hebei province. The parting structures have some distinct geometry and fabrics, it is originated from the acting and reacting fores caused by K-metasomatism. Namely, the crystallizations of metasomatic K-feldspars are a volume expansion process, it would compress the relict fluid bodies, and the pressures in the relict fluid bodies gathered and increased, when the increased pressure of the fluid relict bodies is bigger than the strength of K-feldspars, the K-feldspars were broken with the strong compression, and the parting structures formed. 3)Space position replacing is a important transport pattern of dissolved mass in Fe&S-rich fluid. In addition, basing on views of tectonic dynamics of metallogenic fluids, and time-space texture of fluid-tectonic-lithogenetic-mineralization of the known gold-corpper mineral deposit and the subvolcanic complex at Qibaoshan area in Shandong province, this report does a detail prodict of position-shape-size of two concealed ore-bearing breccia pipe.