Circadian rhythms in the growth and reproduction of the brown alga Undaria pinnatifida and gametogenesis under different photoperiods

Circadian growth rhythm of the juvenile sporophyte of the brown alga Undaria pirznatifida was measured with the computer-aided image analysis system in constant florescent white light under constant temperature (10 degrees C). The growth rhythm persisted for 4 d in constant light with a free-running period of 25. 6 h. Egg release from filamentous gametophytes pre-cultured in the light - dark regime was evaluated for six consecutive days at fixed time intervals in constant white light and 12 h light per day. Egg release rhythm persisted for 3 d in both regimes, indicating the endogenous nature. Temporal scale of egg release and gametogenesis in 18, 16, 12 and 8 h light per day were evaluated respectively using vegetatively propagated filamentous gametophytes. Egg release occurred 2 h after the onset of dark phase and peaked at midnight. Evaluation of the rates of oogonium formation, egg release or fertilization revealed no significant differences in four light-dark regimes, indicating; the great plasticity of sexual reproduction. No photoperiodic effect in gametogenesis in terms of oogonium formation and egg release was found, but fertilization in short days was significantly higher than in long days. Results of this investigation further confirmed the general occurrence of circadian rhythms in inter-tidal seaweed species.

Settling abalone veliger larvae in a free-swimming microalgal culture

In the current abalone hatchery in China, insufficient diatoms on vertically placed corrugated pvc plates at later stage often could not support the growth of postlarvae up to the stage that they can feed on live macroalgae. As a result, stripping the spats (35 mm) off by anaesthetization and switching the diet from live diatoms to artificial powdered diet in combination has to be performed in most of the abalone farms. This manipulation normally leads to more than 50% mortality. Here we report the direct use of the unicellular green alga Platymonas helgolandica Kylin var. tsingtaoensis as a potential alga to be used to settle the veliger larvae of the Pacific abalone Haliotis discus hannai and to feed the postlarvae. Settlement rate of 2-day-old veliger larvae in mono culture of P helgolandica could be as high as 92% ( +/- 4.2%) on day 10 in small scale trials, higher than that in the selected benthic diatom strain (53.6% +/- 12.7%) when settled in the water in which bacteria propagation was controlled by treatment of 2 ppm of benzylpenicillinum calcium and streptomycin sulfate. Postlarvae fed solely on P. helgolandica or the selected benthic diatom Navicula-2005-A grew at rates of 40.1 ( +/- 1.9) and 45.8 (+/- 13.4) mu m day(-1), respectively, when raised at 22 degrees C until day 50 postfertilization. P. helgolandica was shown to have distinct diurnal settling rhythm characterized with a peak of settled cells in the middle of the night for cell division and a peak of free-swimming cells in the middle of the day. High density of attached P. helgolandica cells on the inner surface of the culture facility in the night fits the nocturnal feeding behavior of the abalone spats. Judged by the promising larvae settling rate, growth and survival rates of the postlarvae fed with this alga, the free-swimming micro-green alga P. helgolandica constitutes a potential species for settling the veliger larvae and for supporting the growth of postlarvae as well. (c) 2006 Published by Elsevier B.V.

Observations on the nocturnal activity and feeding behavior of Anguilla japonica glass eels under laboratory conditions

Glass eels of the temperate anguillid species, Anguilla japonica, clearly showed a nocturnal activity rhythm under laboratory conditions. Light-dark cycle was a determinant factor affecting their photonegative behavior, nocturnal locomotor activity, and feeding behavior. Under natural light conditions, glass eels remained in shelters with little daytime feeding, but came out to forage during darkness. They moved and foraged actively in the following dark, and then their activity gradually declined possibly because of food satiation. They finally buried in the sand or stayed in tubes immediately after the lights came on. Under constant light, glass eels often came out of the shelters to forage in the lights but spent little time moving outside the shelters (e.g. swimming or crawling on the sand). Glass eels took shelter to avoid light and preferred tubes to sand for shelter possibly because tubes were much easier for them to take refuge in than sand. Feeding and locomotor activities of the glass eels were nocturnal and well synchronized. They appeared to depend on olfaction rather than vision to detect and capture prey in darkness. Feeding was the driving force for glass eels to come out of sand under constant light. However, in the dark, some glass eels swam or crept actively on sand even when they were fully fed. The lunar cycles of activity rhythms of glass eels that have been observed in some estuarine areas were not detected under these laboratory conditions.

Summer diet of two sympatric species of raptors Upland Buzzard (Buteo hemilasius) and Eurasian Eagle Owl (Bubo bubo) in alpine meadow: problem of coexistence

Summer diets of two sympatric raptors Upland Buzzards (Buteo hemilasius Temminck et Schlegel) and Eurasian Eagle Owls (Bubo bubo L. subsp. Hemachalana Hume) were studied in an alpine meadow (3250 m a.s.l.) on Qinghai-Tibet Plateau, China. Root voles Microtus oeconomus Pallas, plateau pikas Ochotona curzoniae Hodgson, Gansu pikas O. cansus Lyon and plateau zokors Myospalax baileyi Thomas were the main diet components of Upland Buzzards as identified through the pellets analysis with the frequency of 57, 20, 19 and 4%, respectively. The four rodent species also were the main diet components of Eurasian Eagle Owls basing on the pellets and prey leftovers analysis with the frequency of 53, 26, 13 and 5%, respectively. The food niche breadth indexes of Upland Buzzards and Eurasian Eagle Owls were 1.60 and 1.77 respectively (higher value of the index means the food niche of the raptor is broader), and the diet overlap index of the two raptors was larger (C-ue = 0.90) (the index range from 0 - no overlap - to I - complete overlap). It means that the diets of Upland Buzzards and Eurasian Eagle Owls were similar (Two Related Samples Test, Z = -0.752, P = 0.452). The classical resource partitioning theory can not explain the coexistence of Upland Buzzards and Eurasian Eagle Owls in alpine meadows of Qinghai-Tibet Plateau. However, differences in body size, predation mode and activity rhythm between Upland Buzzards and Eurasian Eagle Owls may explain the coexistence of these two sympatric raptors.

黄骅坳陷港西开发区河流相储层构型与剩余油分布模式研究

In the production tail of oilfield, water-cut is very high in thick channel sand oil reservoir, but recovery efficiency is relative low, and recoverable remaining oil reserves is more abundant, so these reserves is potential target of additional development. The remaining oil generally distributed with accumulation in certain areas, controlled by the reservoir architecture that mainly is the lateral accretion shale beddings in the point bar, so the study of reservoir architecture and the remaining oil distribution patterns controlled by architecture are very significant. In this paper, taking the Minghuazhen formation of Gangxi oilfield as a case, using the method of hierarchy analysis, pattern fitting and multidimensional interaction, the architecture of the meandering river reservoir is precisely anatomized, and the remaining oil distribution patterns controlled by the different hierarchy architecture are summarized, which will help to guide the additional development of oil fields. Not only is the study significant to the remaining oil forecasting, but also it is important for the theory development of reservoir geology. With the knowledge of sequence correlation and fluvial correlation model, taking many factors into account, such as combination of well and seismic data, hierarchical controlling, sedimentary facies restraint, performance verification and 3-D closure, an accurate sequence frame of the study area was established. On the basis of high-resolution stratigraphic correlation, single layer and oil sand body are correlated within this frame, and four architecture hierarchies, composite channel, single channels, point bars and lateral accretion sandbody are identified, The result indicates that Minghuazhen Formation of Gangxi oilfield are dominated by meandering river deposition, including two types of channel sandbodies, narrow band and wide band channel sandbody, and each of them has different characteristics of facies variation laterally. Based on the identification of composite channel, according to the spatial combination patterns and identified signs of single channel, combined with channel sandbody distribution and tracer material data, single channel sandbodies are identified. According to empirical formula, point-bar scales of the study area are predicted, and three identification signs are summarized, that is, positive rhythm in depositional sequence, the maximum thick sand and near close to the abandoned channel, and point bars are identified. On the basis of point bar recognition, quantitative architecture models inner point bar are ascertained, taking the lateral accretion sand body and lateral accretion shale beddings in single well as foundation, and quantitative architecture models inner point bar as guidance, and result of tracer material data as controlling, the the lateral accretion sand body and lateral accretion shale beddings are forecasted interwell, so inner architecture of point bar is anatomied. 3-D structural model, 3-D facies model and 3-D petrophysical properties models are set up, spatial distribution characteristics of sedimentary facies and petrophysical properties is reappeared. On the basis of reservoir architecture analysis and performance production data, remaining oil distribution patterns controlled by different hierarchy architecture units, stacked channel, single channel and inner architecture of point bar, are summarized, which will help to guide the additional development of oil fields.

黄骅坳陷马东地区沙河街组重力流水道储层评价及剩余油研究

In this paper, taking Madong district of Huanghua depression as a case, based on the theory of sequence stratigraphy, sedimentology, reservoir geology and geophysics, according to core analysis, seismic attribute analysis, logging constrained inversion, multi-data correlation of strata, reservoir modeling, etc. the lower and middle first member of Shahejie formation of the study area was forecasted and evaluated. As a result, a number of reservoir prediction and remaining oil distribution methods suitable to oil exploitation of gravity flow channel reservoir are presented. Scientific foundation is provided to the next adjustment of development program and exploitation of the remaining oil. According to high resolution sequence stratigraphy theory, precise stratigraphic framework was founded, the facies types and facies distribution were studied under the control of stratigraphic framework, the technologies of seismic attribute abstraction and logging constrained inversion. Result shows that gravity flow channel, as the main facies, developed in the rising period of base-level cycle, and it was formed during the phase of contemporaneous fault growth. As the channel extends, channel width was gradually widened but thickness thined. The single channels were in possession of a great variety of integrated modes, such as isolated, branching off, merging and paralleling, forming a kind of sand-mud interblending complex sedimentary units. Reservoir quality differs greatly in vertical and horizontal direction, and sedimentary microfacies is main controlling factor of the reservoir quality. In major channel, deposition thickness is great, and petrophysical property is well. While in marginal channel, reservoir is thinner, and petrophysical property is unfavorable. Structure and reservoir quality are main factors which control the oil and gas distribution in the study area. On the basis of the research about the reservoir quality, internal, planar and 3-D reservoir heterogeneities are characterized, and the reservoir quality was sorted rationally. At last, based on the research of reservoir numerical simulation of key well group, combined with reservoir performance analysis and geological analysis above, remaining oil distribution patterns controlled by internal rhythm of gravity flow channel were set up. Through this research, a facies-restrained reservoir prediction method integrating multi-information was presented, and potential orientation of remaining oil distribution in gravity flow channel reservoir is clarified.

复杂条件地层滤波预测的理论与关键技术研究

The theory researches of prediction about stratigraphic filtering in complex condition are carried out, and three key techniques are put forward in this dissertation. Theoretical aspects: The prediction equations for both slant incidence in horizontally layered medium and that in laterally variant velocity medium are expressed appropriately. Solving the equations, the linear prediction operator of overlaid layers, then corresponding reflection/transmission operators, can be obtained. The properties of linear prediction operator are elucidated followed by putting forward the event model for generalized Goupillaud layers. Key technique 1: Spectral factorization is introduced to solve the prediction equations in complex condition and numerical results are illustrated. Key technique 2: So-called large-step wavefield extrapolation of one-way wave under laterally variant velocity circumstance is studied. Based on Lie algebraic integral and structure preserving algorithm, large-step wavefield depth extrapolation scheme is set forth. In this method, the complex phase of wavefield extrapolation operator’s symbol is expressed as a linear combination of wavenumbers with the coefficients of this linear combination in the form of the integral of interval velocity and its derivatives over depth. The exponential transform of the complex phase is implemented through phase shifting, BCH splitting and orthogonal polynomial expansion. The results of numerical test show that large-step scheme takes on a great number of advantages as low accumulating error, cheapness, well adaptability to laterally variant velocity, small dispersive, etc. Key technique 3: Utilizing large-step wavefield extrapolation scheme and based on the idea of local harmonic decomposition, the technique generating angle gathers for 2D case is generalized to 3D case so as to solve the problems generating and storing 3D prestack angle gathers. Shot domain parallel scheme is adopted by which main duty for servant-nodes is to compute trigonometric expansion coefficients, while that for host-node is to reclaim them with which object-oriented angle gathers yield. In theoretical research, many efforts have been made in probing into the traits of uncertainties within macro-dynamic procedures.

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.

砂泥岩薄互层储层预测方法与应用研究

With the deeply development of exploration and development in petroleum in China, new increasing reserves are found in old oil fields and the verge of the old ones through re-study of geological property. It is more and more important to discovery and develop thin layer or thin inter-bedded layers reservoirs. All of the targets are thin sand-shale inter-bedded reservoirs and the core technology is reservoir predictions between wells in thin sand-shale inter-bedded layers. The continuity of the thin sand-shale inter-bedded layers in space or separating and heterogeneity is the key of reservoir geology research. The seismic reflection, high resolution analysis method and inversion method to thin sand-shale inter-bedded layers are thorough discussed and deeply studied in this paper to try to find the methods and resolutions of reservoir geology research. The below is followed. 1. Based on the pre-research of other people, five models are created: the sand sphenoid body, interlay sandstone and interlay shale of the equal thickness, interlay sandstone of the equal thickness and interlay shale of the unequal thickness, interlay sandstone of the unequal thickness and interlay shale of the unequal thickness, interlay sandstone of the changing thickness in sequence and interlay shale of the changing thickness in sequence. Then the study of the forward modeling are conducted on the thin layer and thin inter-bedded layers geological characters and seismic reflections including amplitude, frequency, phase, wave shape and time-frequency responding in the domains of time and frequency. The affect of petro-physics difference of layers, single thin layer thickness, thickness of inter-bedded, layer number of inter-bedded, incident wavelet domain frequency and types, sample interval to seismic reflection characters, frequency spectrum and time-frequency respond of reflectivity is theoretically discussed. 2. Qualitatively analyzing the sedimentary rhythm of the thin inter-bedded layers in vertical orientation and computing the single layer thickness or the average thickness with the method of generalized S transform. Identifying the reflecting interface or lithology interface using the amplitude value of amplitude spectrum domain frequency. 3. Based on the seismic respond of thin sand-shale inter-bedded layers, bring out the high resolution analysis method of seismic data in thin sand-shale inter-bedded layers using wavelet analysis and the idea of affecting low and high frequency with middle frequency. Then analyzing the effect to the method and testing some wavelets in the method. This method is applied to the theoretical models and the field data. 4. Bring forward one improved very fast simulated annealing method (IVFSA) to resolve the problem nonlinearity and multi-parameters of the inversion in thin inter-bedded layers. And IVFSA is more productive and higher precision than general ways. 5. New target constrained function is used in the inversion based on the property of the inversion in thin inter-bedded layers. 6. Making the full use of geological and logging information, IVFSA and the new function are applied in the non-linear inversion to improve reservoir prediction and evaluation in thin inter-bedded formations combined with the idea of logging and seismic inversion. This method was applied to the field data and got good results.

煤系地层三维三分量地震数据解释方法研究

Conventional 3D seismic exploration cannot meet the demand of high yield and high efficiency safe production in coal mine any more. Now it is urgent to improve the discovery degree of coal mine geological structures for coal production in China. Based on 3D3C seismic exploration data, multi-component seismic information is fully excavated. First systematic research on 3D3C seismic data interpretation of coal measure strata is carried out. Firstly, by analyzing the coal measure strata, the seismic-geologic model of coal measure strata is built. Shear wave logging is built by using regression analysis. Horizon calibration methods of PP-wave and PS-wave are studied and the multi-wave data are used together to interpret small faults. Using main amplitude analysis technology, small faults which cannot be found from PP-wave sections can be interpreted from the low frequency PS-wave sections. Thus, the purpose to applying PS-wave data to fine structure assistant interpretation is achieved. Secondly, PP- and PS-wave post-stack well constrained inversion methods of coal measure strata are studied. Joint PP- and PS-wave post-stack inversion flow is established. More attribute parameters, which are applied in fine lithology interpretation of coal measure strata, are obtained from combinations of the inversion results. Exploring the relation between rock with negative Poisson’s ratio and anisotropy, fracture development in coal seam are predicted. Petrophysical features of coal measure strata are studied, and the relations between elastic parameters and lithology, fluid and physical properties are established. Inversions of the physical parameters such as porosity, permeability and water saturation, which reflect lithology and fluid property, are obtained. Finally, the approaches of shear wave splitting and Thomsen parameters inversion, which provide new ideas for seismic anisotropy interpretation of coal measure strata, are studied to predict fracture development. The results of practical application indicate that the methods in this paper have good feasibility and applicability. They have positive significance for high yield and high efficiency safe production in coal mine.

三维非均匀各向异性介质中的纵波走时反演研究

Many observations show that seismic anisotropy is very common in the crust and upper mantle of the Earth. Seismic anisotropy can provide some clue about the changing and transporting process inside the earth. in recent years, abundant earthquake travel time data are accumulated, computers become more powerful, and these make the inversion of earthquake travel time data practical. In this thesis we studied the theory of elastic wave in anisotropic media, some formule for travel time inversion were derived. We present an iterative procedure to determine 21 elastic parameters from qP wave travel times. No a priori assumptions about heterogeneity and anisotropy of the model are made. The procedure is suitable for the case when we know nothing about the symmetry of anisotropy of the media, as well as for the case of earthquake travel time inversion which may contain various symmetry of anisotropy. The procedure is tested with a synthetic multiple-source offset VSP experiment. The results proved that the formulae are correct, and the procedure is practical. The results and the related theory indicate that the anisotropic inversion needs more rays than isotropic case. For a 2-D weak anisotropic (WA) medium, we need at least 5 rays in different directions to retrieve the elastic parameters on one grid point, and for a 3-D WA medium we need at least 15 rays in different directions to retrieve the elastic parameters on one grid point. The results also indicate that the starting background velocity has no influence on the final results, at least for the model we specified. Our results also show that insufficient illumination coverage will slow down the convergence rate, and make the results more sensitive to noise. We apply the procedure to a set of field travel time data. The data is from an artificial seismic observation. This observation is for locating micro-seismic events around a tunnel, its purpose is to find out if the digging process and the stress condition around the tunnel can generate micro-cracks. The size of this area is around 100m. The anisotropy derived from qP travel times is the same as the anisotropy showed by apparent velocities, and is also consistent with the anisotropy derived from S-wave splitting phenomena.

河流成因储层沉积微相与剩余油分布及形成机理研究

Based on the study of fluvial sandstone reservoir in upper of Guantao group in Gudao and Gudong oilfields, this paper first introduces A.D.Miall's(1996a) architectural-element analysis method that was summarized from ground outcrop scale into the reservoir formation research of the study area, more subtly divides sedimentary microfacies and establishes sedimentary model of research area.on this base, this paper summarizes the laws of residual oil distribution of fluvial formation and the control effect of sedimentary microfacies to residual oil distribution, and reveals residual oil formation mechanism. These results have been applied to residual oil production, and the economic effect is good. This paper will be useful for residual oil research and production and enhancement of oil recovery in similar reservoir. The major conclusions of this paper are as follows. 1. Using the architectural-element analysis method to the core data, a interfacial division scheme of the first to the dixth scale is established for the studied fluvial formation. 2.Seven architectural-elements are divided in upper of Guantao group of study area. The sandstone group 5～1+2 of Neogene upper Gutao group belongs to high sinuous fine grain meandering river, and the sandstone group 6 is sandy braided river. 3. Inter layer, the residual oil saturation of "non-main layer" is higher than "main layer", but the residual recoverable reserve of former is larger. Therefore, "main layer" is the main body of residual oil distribution. The upper and middle part of inner layer has lower permeability and strong seeping resistance. Addition to gravity effect in process of driving, its driving efficiency is low; residual oil saturation is high. Because of controlling of inside non-permeable interlayer or sedimentary construction, the residual oil saturation of non-driving or lower driving efficiency position also is high. On plane, the position of high residual oil saturation mostly is at element LV, CS, CH (FF), FF etc, Which has lower porosity and permeability, as well as lens sand-body and sand-body edge that is not controlled by well-net, non-perfect area of injection and production, lower press difference resort area of inter-well diffiuent-line and shelter from fault, local high position of small structure. 4.Microscopic residual oil mainly includes the non-moved oil in the structure of fine pore network, oil in fine pore and path, oil segment in pore and path vertical to flow direction, oil spot or oil film in big pore, residual oil in non-connective pore. 5.The most essential and internal controlling factor of fluvial formation residual oil distribution is sedimentary microfacies. Status of injection and production is the exterior controlling factor of residual oil distribution. 6. The controlling effect of formation sedimentary microfacies to residual oil distribution indicates inter-layer vertical sedimentary facies change in scale of injection and production layer-series, planar sedimentary face change and inner-layer vertical sedimentary rhythm and interbed in single layer to residual oil distribution. 7. It is difficult to clear up the inter-layer difference in scale of injection and production layer-series. The using status of minor layer is not good and its residual oil saturation is high relatively. It is obvious that inter-layer vertical sedimentary facies changes control inter-layer residual oil distribution at the same or similar conditions of injection and production. For fluvial formation, this vertical sedimentary facies change mainly is positive gyration. Namely, from down to top, channel sediment (element CHL, LA) changes into over-bank sediment (element LV, CR, CS). 8. In water-injection developing process of transverse connecting fluvial sandstone oil formation, injection water always comes into channel nearby, and breaks through along channel and orientation of high pressure gradient, does not expand into side of channel until pressure gradient of channel orientation changes into low. It brings about that water-driving status of over-bank sedimentary element formation (LV, CR, CS) is not good, residual oil saturation is high. In non-connective abandoned channel element (CH) formation with channel, because this reverse is difficult to control by injection and production well-series, its using status is not good, even terribly not good, residual oil is enrichment. 9. The rhythm and sedimentary structure, sedimentary facies change in single sand body brings about vertical changes of formation character, growth character of inner layer interbed. These are important factor of controlling and affecting vertical water spread volume and inner layer residual oil forming and distribution in single sand body. Positive rhythm, is the principal part of fluvial sandstone inner layer sedimentary rhythm. Namely, from down to upside, rock grain granularity changes from coarse to fine, seeping ability changes from strong to feebleness. It brings about that water-driving status of inner layer upside is not good, residual oil saturation is high. Inner layer interbed has different degree affecting and controlling effect to seeping of oil and water. Its affecting degree lies on interbed thickness, extending scale, position, and jeted segment of production or injection well. The effect of interbed at upside of oil formation to oil and water seeping is less; the effect of interbed at middle of oil formation to oil and water seeping is more. 10. Indoor experiment and research indicate that wettability, permeability step, vertical permeability, position of Kmax and ratio of oil viscousity and water viscousity all have great effect on the water-driving recovery ratio. 11. Microscopic residual oil distribution is affected and controlled by formation pore network structure, pressure field distribution, and oil characteristic. 12.The residual oil forming mechanism: the over-bank sedimentary element and upper part of a positive rhythm sandstone have fine pore and throat network, permeability is low, displacement pressure of pore and throat is high. The water-driving power usually falls short of displacement pressure that brings about injection water does not spread into these pore and throat network, thereby immovable oil area, namely residual oil, is formed. At underside of channel sedimentary element and positive rhythm sandstone, porosity and permeability is relatively high, connecting degree of pore and throat is high, displacement pressure of pore and throat is low. Thereby injection water is easy to enter into pore and throat, driving oil in them. Because the pore space is irregular, the surface of pore wall is coarse and non-flat. That the oil locate on concave hole of pore wall and the dead angle of pore, and the oil attaches on surface of pore wall by surface tension, are difficult to be peeled off, becoming water-driving residual oil (remaining oil). On the other hand, Because flowing section lessens, flowing resistance increase, action of capillary fore, or seeping speed decreases at process of transfer at pass narrow throat path in the course carried by driving water. The "oil drop", "oil bead", or "oil segment" peeled off by driving water is difficult to carry and to drive out by water at less pressure difference. Thereby they are enclosed in pore to form discontinuous residual oil. 13.This results described above have been applied in nine develop blocks of Gudao and Gudong oilfield. Its applying effect is marked through local injection production adjustment, deploying replacement well, repair hole, replacement envelop, block off water and profile control etc. Relative method and technology can be applied to other oil production area of Shengli oilfield, and obtain better economic and societal effect.

攀西地区红格、新街岩体的岩石地球化学特征

The mafic-ultramafic layered intrusions in the Panxi, China contain large V-Ti-magnetite deposits. These layered intrusions are related with the Emeishan continental flood basalts in space and time. Two layered intrusions, Hongge and Xinjie have clear PGE mineralization at the base of the intrusions. Thus the detailed investigations of these two intrusions not only have a geological but also have an economic significance. This thesis aims to characterize the elemental and Sr-Nd isotopic features of diverse rock zones within the intrusion on the basis of systematic studies of the major, trace element and isotope ratios, therefore to constrain the petrogenesis, mantle source and evolution of the Hongge and Xinjie intrusions. Generally, both Hongge and Xinjie intrusions show the same Fe-Ti-rich and Si-M-poor characteristics. They are also enriched in rare-earth elements (REE) and large-ion lithophile elements (LILE) as well as in Sr-Nd isotope ratios (Hongge: initial Sr = 0.7056-0.7076, ε_(Nd)(t) and (Nd/Sm)_N-ε_(Nd)(t) plots, the Hongge intrusion has a similar elemental and isotopic features to the Emeishan low-Ti (LT) basalts, whereas the Xinjie intrusion was close to the Emeishan high-Ti (HT) basalt. Therefore, the Hongge intrusion may be co-genetic with the LT basalt, formed by the partial melting of the spinel-garnet transition mantle that had a slight enriched isotope character. In contrast, the Xinjie intrusion and the HT basalts are probably derived from the garnet-phases mantle with a primitive isotope character. The involvement of the components of mantle wedge into the source is considered to be the major reason of the REE and LILE enrichment and Nd isotope depletion in the Xinjie intrusion. In contrast with the systematic variations in TiO_2 content, Mg#, transition elements (Ni, Cu, Co), REE concentrations, and La/Yb, La/Sm ratios from the lower zone to upper zone, the different rock zones of the Hongge intrusion have no clear Sr-Nd isotope variations. This suggests that the Hongge intrusions were formed by the crystal fractionation from the same magma source. The rhythm may be formed by slow injection of the co-genetic magma during the crystal fractionation. The increase in K_2O and Al_2O_3 contents, REE abundance, and the degree of the REE fractionation in the base of the intrusion, together with the relatively low ε_(Nd)(t) value, may imply that the base of the Hongge intrusion was contaminated with the local crust rocks. Xinjie intrusion shows the clearly elemental and isotopic differences in diverse cumulus cycles. The observation of the systematic variations in TiO_2 content, Mg# value, transition elements (Ni, Cu, Co), REE concentrations, and La/Yb, La/Sm ratios in first cycle was not occurred in second cumulus cycle. In addition, the ε_(Nd)(t) value in second cumulus cycle is apparently higher than that of the first one. Thus the abruptly elemental and isotopic changes at the base of second cycle demonstrate that there is considerable new and depleted magma addition to the residue magma after the crystallization of the first cycle. These features are very similar to those of the well-known PGE-rich Bushveld and Stillwater layered intrusions. The PGE mineralization in Xinjie intrusion is much better than in Hongge intrusion. Therefore, the layered intrusion similar to the Xinjie in Panxi area posses the better prospects for the PGE deposits.

各向异性介质中的地震面波和S波分裂

Surface wave propagation in the anisotropic media and S-wave splitting in China mainland are focused in this M.S. dissertation. We firstly introduced Anderson parameters in the research of surface wave propagation in the anisotropic media were deduced, respectively. By applying the given initial model to the forward calculation of Love wave, we compared dispersion curves of Love wave in the anisotropic media with the one in the isotropic media. the results show that, although the two kind of results are similar with each other, the effect of anisotropy can not be neglected. Furthermore, the variation of anisotropy factors will result in the variation of dispersion curves, especially for high-mode one. The method of grid dispersion inversion was then described for further tectonic inversion. We also deduced inversion equation on the condition that the layered media is anisotropic, and calculated the phase-velocity partial derivatives with respect to the model parameters, P- and S-wave velocities, density, anisotropic parameters for Rayleigh wave and Love wave. Having analyzed the results of phase-velocity partial derivatives, we concluded that the derivatives within each period decreased with the depth increasing, the phase-velocity of surface wave is sensitive to the S-wave velocities and anisotropic factors and is not sensitive to the densities of layers. Dispersion data of Love wave from the events occurred during the period from 1991 to 1998 around the Qinghai and Tibet Plateau, which magnitudes are more than 5.5, have been used in the grid dispersion inversion. Those data have been preprocessed and analyzed in the F-T domain. Then the results of 1°*1° grid dispersion inversion, the pure path dispersion data, in the area of Qianghai and Tibet Plateau were obtained. As an example, dispersion data have been input for the tectonic inversion in the anisotropic media, and the results of anisotropic factors under the region of Qianghai and Tibet Plateau were initially discussed. As for the other part of this dissertation. We first introduced the phenomena of S-wave splitting and the methods for calculation the splitting parameters. Then, We applied Butterworth band-pass filter to S-wave data recorded at 8 stations in China mainland, and analyzed S-wave splitting at different frequency bands. The results show the delay time and the fast polarization directions of S-wave splitting depend upon the frequency bands. There is an absence of S-wave splitting at the station of Wulumuqi (WMQ) for the band of 0.1-0.2Hz. With the frequency band broaden, the delay time of S-wave splitting decreases at the stations of Beijing (BJI), Enshi (ENH), Kunming (KMI) and Mudanjiang (MDJ); the fast polarization direction at Enshi (ENH) changes from westward to eastward, and eastward to westward at Hailaer (HIA). The variations of delay time with bands at Lanzhou (LZH) and qiongzhong (QIZ) are similar, and there is a coherent trend of fast polarization directions at BJI, KMI and MDJ respectively. Initial interpretations to the results of frequency band-dependence of S-wave splitting were also presented.

边坡动力响应分析及应用研究

As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum VHlhro) point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fN