黄土高原旱区长期施肥条件下土壤钾素形态空间分布特征及有效性研究

在黄土高原南部旱地长期肥料定位试验的基础上研究了土壤钾素空间分布特征及其有效性。结果表明:长期施肥后土壤中特殊吸附性钾(SAK)和非特殊吸附性钾(NSAK)储量增加,但水溶性钾(WSK)和非交换性钾(NEK)则有明显的下降,单施N水溶性钾下降了48.24%,单施P下降32.32%,NP配施和NPK配施分别下降10.61%和17.93%,非交换性钾降幅为8.56%~24.91%。增施钾肥可以缓解因长期施肥作物生长所携出的钾素,增加耕层土壤中的水溶性钾、非特殊吸附性钾及特殊吸附性钾。相关分析结果表明,土壤不同形态钾素对速效钾的重要性依次为WSK>NSAK>SAK>NEK,土壤速效钾与水溶性钾、非特殊吸附性钾呈显著相关,与特殊吸附性钾和非交换性钾无显著相关性。

九寨沟湖泊湿地植物群落与植物地理研究

九寨沟湖泊湿地在维持九寨沟的生态平衡中起着重要的作用，在旅游产业的发展下，湿地生态系统及生物多样性面临着较大的威胁。尽管九寨沟湿地具有重要的生态价值，但目前对其研究尚比较薄弱。湿地植物群落和植物地理研究可以为湿地资源的可持续利用和监测保护提供科学依据。作者从2004年8月到2007年11月对九寨沟湿地的植物物种组成、地理分布、优势植物群落的结构、生长动态、湿地土壤种子库进行了调查研究。主要结果如下： 1. 九寨沟湿地物种组成、地理分布特点及湿地植物群落特点 九寨沟湿地共有苔藓植物8科13属16种，维管植物为48科107属199种。九寨沟湿地植物的地理成份较为丰富，维管植物在科级水平上有7种地理分布型（变型），在属级水平上有13种地理分布型（变型）， 在种级水平上共有29种地理分布型（变型）。九寨沟湿地植物以温带成份和我国特有成份为主，同时兼有热带、亚热带成份和环极—高山成份。九寨沟湿地植物的分布表现出明显的垂直地带性和水平地带性。湿地植物群落可划分21个群落类型，不同植物群落类型的物种多样性及物种组成存在较大的差异。九寨沟湿地植物的物种多样性和群落多样性以及较高的生产力特征，是维持其湿地生态景观多样性和稳定性的基础。 2. 土壤、水环境、海拔等对湿地植物的分布及生物多样性的影响 九寨沟湿地土壤、水等环境因子存在较大的差异。帕米尔苔草和宽叶香蒲等群落的凋落物较多，土壤有机碳、土壤总磷较高，可能是九寨沟湿地的重要土壤碳库。 九寨沟湿地植物沿水环境梯度的分布规律表现为：沉水植物（轮藻—篦齿眼子菜，水苦荬，杉叶藻）——挺水植物（水木贼，芦苇，宽叶香蒲）——湿生草本（苔草、节节草、披散木贼）——湿生灌木（柳灌丛，小檗灌丛）等。海拔也影响湿地植物的物种组成。 水深对物种多样性有影响，水深与物种丰富度负相关。随着水深的增加，水木贼、芦苇、杉叶藻、宽叶香蒲等群落的物种多样性下降；在长期淹水和季节性淹水的地方，水木贼群落物种多样性存在显著差异。土壤总氮与水木贼群落物种丰富度正相关。 3. 土壤营养元素、水环境对植物生长的影响 水深影响湿地植物生物量的分配。芦苇无性系分株在47 cm水深的环境中单株平均生物量最大；在干滩地中（地面水深0 cm），叶生物量百分比最大，而茎生物量百分比最小，茎的生物量百分比和生长速率随水深的增加而增加；在较干的滩地生境中，开花率、花序的生物量百分比明显大于水较深的生境。 水深与水木贼地上生物量负相关，但水木贼地上生物量在长期淹水和季节性淹水的地方没有显著的差异。在水浅的地方，杉叶藻、水木贼、芦苇等植物群落中，其他伴生物种的生物量占样方总生物量的百分比较大。 土壤有机碳、土壤总氮、土壤总磷等对湿地植物生物量的影响比较大：宽叶香蒲地上生物量与土壤总磷正相关；水木贼地上生物量与土壤总氮正相关；杉叶藻地上生物量与土壤有机碳正相关。 水深、土壤营养成分对湿地植物高度、密度等有影响。水木贼的平均高度在季节性淹水的地方比长期淹水的地方低，平均密度在长期淹水的地方比季节性淹水的地方低；除了5月份，其他观察月份水木贼的密度都与水深负相关，同时与土壤有机碳正相关。另外，芦苇密度与土壤有机碳含量正相关，宽叶香蒲密度与水深负相关，帕米尔苔草高度与土壤有机碳负相关。 4. 优势植物群落的动态变化 在优势植物群落中，优势种的高度、密度、盖度、生物量等在群落中占绝对优势。除五花海，水木贼群落的物种组成、高度、生物量在两年间没有显著的变化。芦苇群落的物种丰富度在近两年有所增加。 湿地植物生长表现为明显的季节动态，生长的峰值大多在7月－8月。优势植物群落的物候与水文周期有关。湿地植物群落的物种组成和密度，可以作为对湿地监测和保护的生物指示。 5. 九寨沟湿地土壤种子库特征及其在湿地生物多样性恢复中的作用 水深和现存植被物种丰富度可以解释湿地土壤种子库的变化。水深可以解释表层物种丰富度45%的变化。现存植被物种丰富度可以分别解释10 cm土层、2－5 cm土层及5－10 cm土层土壤种子库45%、48%和25%的变化。 湿地土壤种子库的密度为0－15945粒m-2, 种子库中共发现23个物种。现存植被优势物种和种子库优势物种不同。各层土壤种子库密度和物种丰富度并不存在显著的差异，但第二层土壤种子库密度最大。海拔、现存植被优势种盖度、土壤总磷、土壤总氮、土壤有机碳对湿地土壤种子库的密度和垂直结构没有影响。土壤种子库物种丰富度小于地上植被物种丰富度。湿地土壤种子库与地上植被的相关性不大。在浅水区域，湿地土壤种子库在湿地植被恢复中有一定作用。但在深水区域，保护现存植被更重要。 The lakeshore wetlands are valuable ecological units of the Jiuzhaigou lakes. Pressure for travel industry development pose a continuing and severe threat to the biodiversity-support function of the wetland system. Despite the ecological importance of wetlands in Jiuzhaigou, they are so far poorly studied. Both general plant communties and biogeographical studies are needed in order to attain basis for sustainable use the wetland resources and adequate protection of these areas. The present study was undertaken to examine aquatic plants distribution and the species compositon, structure and growth dynamics of their communities with variations of environmental factors along altitudes, water depth and soil properities gradients in Jiuzhaigou. Analysis of field survey data collected during August 2004 and November 2007 in lakeshore wetlands in Jiuzhaigou National Nature Reserve, Sichuan, China. The results were as following: (i) Species composition and biogeography in wetland vegetation 8 families, 13 genus, 16 species of moss and 48 families, 107 genus and 199 species of vascular plants in Jiuzhaigou wetlands were found. The floristic compositions were abundunt. Ten geographical distribution types at family level, 13 geographical distributions types at generic level and 29 geographical distribution types at specific level in vascular plants were found. Most species in Jiuzhaigou wetlands are temperate elements and Chinese endemic elements, with a few of tropical and subtropical and some circumarctic elements. And the plant distributions show clear vertical and horizontal patterns. There were 21 major wetland plant community types. Species composition and species richness in different plant communities are different. The species diversity and plant community diversity and their high biomass are the basis for the diversity and stability of wetland landscapes in Jiuzhaigou. (ii) Water depth, soil nutrients and altitudes influence on the species diversity and plant distribution. Total phosphorous and organic cabon in soil were higher in C. pamiernensis and T. latifolia communities, where are important cabon reservoirs in Jiuzhaigou wetlands. Along gradients of water depth, among populations of the dominant plant species present: submerged macrophytes (Chara vulgaris, Potagemonton pectinatus, Veronica anagalis-aquatica，Hippuris vulgaris), emergent macrophytes (Equisetum fluviatile, Phragamites australis, Typha latifolia), helophytes (Carex pamirensis )and shrubs (Salix sp., Berberis sp. ). Altitudes influence on the assemblage of plant communities. Water depth negatively correlated with species richness. Specie richness showed differences between permanently flooded sites and seasonally flooded sites in E. fluvatile communities. And total nitrogen in soil was negatively correlated with species richness in E. fluviatile communities. Altitudes show no significant influence on species richness, but in fact, through our analyses, they do have influence on the assemblage of wetland plants. (iii) Water depth, soil nutrients influence on the plant growth Water depth influences the biomass allocation in Phragmities australis. The average aboveground biomass of a single ramet (4.2 g) was the largest in the habitat with water level 47 cm above the soil surface. At the habitat with water level under soil surface 15 cm (－15 cm), the leaf biomass percentage (of the total ramet biomass) was the largest (46.1%), and the height and percentage of ramose ramets ( with branches on stem )(of the total ramets in a plot) were found obviously different. The deeper in water, the larger the biomass percentage and growth rate of stems were. The flowering rate and biomass of panicles were greater in shallow water than those in deep water. Water depth negatively correlated with aboveground biomass of E. fluviatile. However, above-ground biomass of E. fluviatile showed no significant difference between permanently flooded sites and seasonally flooded sites. But in shallow water, more biomasses of accompanying species were found in dominant plant communities such as H. vulgaris communities, E. fluviatile communities and P. australis communities. Water depth, soil nutrients influence on shoot density and shoot length of wetland plants. The shoot density of E. fluviatile was correlated to water depth in all growth months. Annual average density was significantly lower at permanently flooded sites than at seasonally flooded sites. But the annual average shoot length was significantly lower at seasonally flooded sites than at permanently flooded sites. (iv) Growth dynamics of dominant communities in Jiuzhaigou wetland The shoot length and shoot density, coverage and biomass of domiant species were dominated in plant communities. The species composition increased in P. australis communities in recent two years. The species richness in E. fluviatile communities showed no difference between 2005 and 2007. The above-ground biomass and shoot density in Five－flower Lake from July 2005 to July 2007 were significantly different, while in other sites, the differences were not significant. Shoot height, shoot density and above-ground biomass showed significant seasonal changes in all sites. Growth dynamics correlated with the cycle of water levels in lakes. Most plants growth parameters peaked at July or August. The biomass of T. latifolia peaked in August. But the shoot length of T. latifolia in deeper water peaked in July. The shoot length of E. fluviatile increased significantly from May to August except in seasonally flooded sites in Arrow-bamboo Lake. The species composition of communities and shoot density can be used as bioindicators in Jiuzhaigou wetland. (v) Soil seed bank in Jiuzhaigou wetland and its role in vegetation restoration Seed density in all soil layer samples was negatively correlated to water depth. Water depth can explain 45% variance of species richness in surface layer in sediment. Species richness in extant vegetation can explain 45%, 48%, 25% variance of species richness in total 10 cm and in 2－5 cm and 5－10 cm layer sediment respectively. Mean seed densities in wetlands ranged from 0 to 15945 m–2. A total of 23 species germinated in seed bank. The dominant species in seed bank and extant vegetation showed great difference. The total number of species and seedlings that germinated in different layers was not significantly different. But the second layer had the greatest seed density. In shallow water, seed bank can contribute to vegetation restoration, while in deeper water, protection of extant vegetation may be a better strategy.

川西北高寒草甸土壤呼吸研究：季节动态和放牧效应

除植被冠层的光合作用之外，土壤的呼吸作用是陆地生态系统碳收支中最大的通量。土壤呼吸即使发生较小的变化也能显著地减缓或加剧大气中CO2浓度的增加，从而明显影响到全球气候变化。土壤呼吸速率变化与否以及变化的方向可以反映生态系统对环境变化的敏感程度和响应模式。尽管如此，土壤呼吸仍是一个为人们了解不多的生态系统过程。 草地生态系统是陆地生态系统的一个重要组成部分。针对草地土壤呼吸进行野外实验研究和相应方法论的探讨将对区域乃至全球碳源汇性质的准确估算具有重要的科学意义。然而，近几年来关于草地土壤呼吸的主要研究工作都集中在温带草原和部分热带草原，而针对高寒草甸生态系统土壤呼吸的研究报道还很少。 2008年4月至2009年4月期间，我分别在2008年6、8、10、12月和2009年2月和4月分6次对川西北的典型高寒草甸群落的土壤呼吸进行观测，分析了不同类型高寒草甸群落土壤呼吸的季节变化特征以及环境因子和放牧模式对其影响。主要研究结果如下： 1）该地区高寒草甸生态系统在生长季（6月～8月）土壤呼吸速率较大（6.07～9.30μmolCO2¡m-2¡s-1 ） ， 在非生长季（ 12 月～ 2 月） 较小（ 0.16 ～0.49μmolCO2¡m-2¡s-1 ） 。土壤CO2 年累积最大释放量为3963 ～ 5730gCO2¡m-2¡yr-1，其中，生长季土壤CO2的释放量占年总释放量的85%～90%。非生长季占10%～15%。非生长季所占比例略小于冬季积雪覆盖地区的冬季土壤呼吸占年土壤呼吸量的比例（14%～30%）。温度，尤其地温，是影响该地区高寒草甸生态系统土壤呼吸速率的最主要环境因子。土壤呼吸速率与地上生物量和土壤水分之间没有显著相关性，但是土壤含水量过大会导致土壤呼吸速率下降。 2）在观测期内，草丘区的土壤呼吸显著高于对照区的土壤呼吸，其最大土壤呼吸速率为16.77μmolCO2¡m-2¡s-1，土壤CO2 年累积最大释放量为8145gCO2¡m-2¡yr-1，是对照区的近2 倍。由于草丘在高寒草甸中占有较大的面积比例（近30%），因此，它将对高寒草甸生态系统的碳循环起着重要的作用。 3）放牧模式不仅可以影响高寒草甸群落的土壤CO2 排放，而且还可以改变土壤呼吸的温度敏感性（Q10）。本研究表明，在生长季有长期放牧活动干扰时将会增加土壤向大气中释放二氧化碳的速度，促使土壤碳库中碳的流失。禁牧样地的土壤呼吸速率在刚禁牧时先迅速增大，随着禁牧时间的延长土壤呼吸速率将会下降。此外，与其它放牧模式相比，冬季放牧将高寒草甸群落土壤呼吸速率在生长季达到最大值的时间明显向后推迟。不同放牧模式下高寒草甸群落土壤呼吸的Q10 值大小顺序为：禁牧一年群落>冬季放牧群落>禁牧三年群落>夏季放牧群落>自由放牧群落。 4）基于呼吸室技术的观测方法中，测量前的剪草处理可以明显改变该地区高寒草甸群落的土壤温度和土壤呼吸速率。在生长季，剪草处理将使土壤呼吸速率的瞬时响应增加90%左右。由于剪草处理明显增加了剪草样方白天的土壤温度，而土壤温度与土壤呼吸之间存在着极显著的指数相关关系，因而剪草处理导致土壤呼吸速率迅速增加。因此，在高寒地区基于呼吸室技术观测的土壤呼吸应当进行校正。 综上所述，川西北高寒草甸生态系统土壤呼吸速率在生长季较高，而在非生长季较低。土壤温度是影响该地区土壤呼吸的最主要环境因子。在实验观测期，草丘区土壤呼吸速率显著高于对照区的，是对照区土壤呼吸速率的近2倍。由于测量前的剪草处理可以明显改变待测点的土壤呼吸速率，因此，应对在高寒地区基于呼吸室技术观测的土壤呼吸进行校正。 Soil respiration is the second largest component (less than plant phtotosynthesis) of carbon dioxide flux between terrestrial ecosystems and the atmosphere. A minor change in soil respiration rate can significantly slow down or accelerate the increase of atmospheric CO2 concentration that is closely related to global climatic change. In turn, the change in the flux direction and rate of soil respiration may indicate the elasticity and stability of ecosystems to global changes and human disturbance. However, soil respiration is still an ecosystem process that has been poorly understood. Grassland ecosystem is an important component of the terrestrial ecosystem. Accurately estimating the CO2 flux from soil to atmosphere in situ is the key to evaluating the carbon resource and sink regionally or globally. Despite of extensive studies on the temperate and tropic grasslands, the soil respiration of alpine meadows has not substantially been measured. In the current study, soil respiration was measured for an annual cycle from April, 2008 to April, 2009 for the alpine meadow in northwestern Sichuan Province of China to determine the seasonal variation of soil respiration for the typical plant communities. The results are shown as follows: 1) Large seasonal variation of soil respiration was observed in the alpine meadow. The rate of soil respiration was the greatest (6.07～9.30μmolCO2¡m-2¡s-1) in June and the smallest (0.16 ～ 0.49μmolCO2¡m-2¡s-1) occurred from December to February in the non-growing season. The total emission of soil CO2 was estimated as 3963～5730 gCO2¡m-2¡yr-1, 85%~90% of which was released during the growing season, and 10%~15% during the non-growing season which was slightly less than the ratio of winter and annual CO2 flux from soil. Temperature, particularly the soil temperature, was the major environmental factor regulating the soil respiration. Significant and positive relationships were not found between soil respiration and soil moisture and between soil respiration and plant above-ground biomass, but excessive soil water content would decrease in the rate of soil respiration. 2) The rate of soil respiration in grass hummock communities was up to 16.77μmolCO2¡m-2¡s-1, which was about twice as great as in the controls (communities located in low and even sites). Considering the large proportion (about 30% on average) of hummock area in the meadow, it can be concluded that the hummocks played an important role in the carbon cycling of the study ecosystem. 3) Grazing patterns affected the flux of CO2 emission and the temperature sensitivity of soil respiration (Q10) in the alpine meadow. Grazing during growing season increased the rate of soil respiration. The rate of soil respiration increased significantly immediately after the alpine meadow being fenced, but thereafter decreased. In addition, grazing in winter delayed the peak respiration rate relative to the non-grazing mode. The Q10 value was the largest in the non-grazed area for one year, and next came the area with grazing in winter, followed by the non-grazed area for three years, the area with grazing in summer, and the non-limited grazed area. 4) In the chamber-based techniques, clipping manipulation before each measurement increased the transient rate of soil respiration by about 90% in the summer of the alpine meadow. As increase in soil temperature at daytime in the clipped plots by clipping and the exponential relationship between soil respiration and temperature, clipping manipulation led to increase in the rate of soil respiration. This suggested that a correction should be done for the techniques if employed in alpine and cold regions. In summary, the rate of soil respiration in the alpine meadow was the greatest in June and the smallest occurred from ecember to February in the non-growing season. Soil temperature was the major environmental factor regulating the soil respiration. The rate of soil respiration in grass hummock communities was up to 16.77μmolCO2¡m-2¡s-1, which was about twice as great as in the controls. A correction should be done for the techniques if employed in alpine and cold regions, because of the effect of clipping manipulation on soil temperature and respiration.

Methods for physical separation and characterization of soil organic matter fractions

Turnover of soil organic matter (SOM) is coupled to the cycling of nutrients in soil through the activity of soil microorganisms. Biological availability of organic substrate in soil is related to the chemical quality of the organic material and to its degree of physical protection. SOM fractions can provide information on the turnover of organic matter (OM), provided the fractions can be related to functional or structural components in soil. Ultrasonication is commonly used to disrupt the soil structure prior to physical fractionation according to particle size, but may cause redistribution of OM among size fractions. The presence of mineral particles in size fractions can complicate estimations of OM turnover time within the fractions. Densiometric separation allows one to physically separate OM found within a specific size class from the heavier-density mineral particles. Nutrient contents and mineralization potential were determined for discrete size/density OM fractions isolated from within the macroaggregate structure of cultivated grassland soils. Eighteen percent of the total soil C and 25% of the total soil N in no-till soil was associated with fine-silt size particles having a density of 2.07-2.21 g/cm3 isolated from inside macroaggregates (enriched labile fraction or ELF). The amount of C and N sequestered in the ELF fraction decreased as the intensity of tillage increased. The specific rate of mineralization (mug net mineral N/mug total N in the fraction) for macroaggregate-derived ELF was not different for the three tillage treatments but was greater than for intact macroaggregates. The methods described here have improved our ability to quantitatively estimate SOM fractions, which in turn has increased our understanding of SOM dynamics in cultivated grassland systems.

黄土丘陵沟壑区退耕地人工柠条林土壤养分特征及其空间变异

应用时空互代方法,以柠条林为例,对黄土丘陵沟壑区不同利用年限人工林土壤养分特征、空间变异及其演变进行了系统的研究。结果表明,该区人工林土壤肥力处于较低水平;人工林表层土壤养分中速效磷和速效钾的空间变异性较大;各环境因子对土壤养分有一定的影响,海拔、坡度和坡向等环境因子与人工林地土壤养分间呈负相关关系,坡位和地形与人工林地土壤养分间呈正相关关系。随着利用年限的增加,人工林土壤养分各指标含量均增加,与利用年限有显著的相关性。有机质、全氮、有效氮和速效磷的增加量不明显,全磷含量保持相对稳定的水平,速效钾经过多年积累有明显的增加。从土壤养分指数模型可以算出该区的人工林土壤有机质、全氮、有效氮约需35a能达到中上等养分水平,速效钾则约需27a才能达到中上等养分水平。

Probabilistic modelling of soil moisture dynamics of irrigated cropland in the North China Plain

A probabilistic soil moisture dynamic model is used to estimate the soil moisture probability distribution and plant water stress of irrigated cropland in the North China Plain. Soil moisture and meteorological data during the period of 1998 to 2003 were obtained from an irrigated cropland ecosystem with winter wheat and maize in the North China Plain to test the probabilistic soil moisture dynamic model. Results showed that the model was able to capture the soil moisture dynamics and estimate long-term water balance reasonably well when little soil water deficit existed. The prediction of mean plant water stress during winter wheat and maize growing season quantified the suitability of the wheat-maize rotation to the soil and climate environmental conditions in North China Plain under the impact of irrigation. Under the impact of precipitation fluctuations, there is no significant bimodality of the average soil moisture probability density function.

Particle-attached and free-living bacterial communities in a contrasting marine environment: Victoria Harbor, Hong Kong

Diversity of particle-attached and free-living marine bacteria in Victoria Harbor, Hong Kong, and its adjacent coastal and estuarial environments was investigated using DNA fingerprinting and clone library analysis. Denaturing gradient gel electrophoresis (DGGE) analysis of 16S rRNA genes showed that bacterial communities in three stations of Victoria Harbor were similar, but differed from those in adjacent coastal and estuarine stations. Particle-attached and free-living bacterial community composition differed in the Victoria Harbor area. DNA sequencing of 28 bands from DGGE gel showed Alphaproteobacteria was the most abundant group, followed by the Bacteroidetes, and other Proteobacteria. Bacterial species richness (number of DGGE bands) differed among stations and populations (particle-attached and free-living; bottom and surface). BIOENV analysis indicated that the concentrations of suspended solids were the major contributing parameter for the spatial variation of total bacterial community structure. Samples from representative stations were selected for clone library (548 clones) construction and their phylogenetic distributions were similar to those of sequences from DGGE. Approximately 80% of clones were affiliated to Proteobacteria, Bacteroidetes and Cyanobacteria. The possible influences of dynamic pollution and hydrological conditions in the Victoria Harbor area on the particle-attached and free-living bacterial community structures were discussed.

Particle-attached and free-living bacterial communities in a contrasting marine environment: Victoria Harbor, Hong Kong

Diversity of particle-attached and free-living marine bacteria in Victoria Harbor, Hong Kong, and its adjacent coastal and estuarial environments was investigated using DNA fingerprinting and clone library analysis. Denaturing gradient gel electrophoresis (DGGE) analysis of 16S rRNA genes showed that bacterial communities in three stations of Victoria Harbor were similar, but differed from those in adjacent coastal and estuarine stations. Particle-attached and free-living bacterial community composition differed in the Victoria Harbor area. DNA sequencing of 28 bands from DGGE gel showed Alphaproteobacteria was the most abundant group, followed by the Bacteroidetes, and other Proteobacteria. Bacterial species richness (number of DGGE bands) differed among stations and populations (particle-attached and free-living; bottom and surface). BIOENV analysis indicated that the concentrations of suspended solids were the major contributing parameter for the spatial variation of total bacterial community structure. Samples from representative stations were selected for clone library (548 clones) construction and their phylogenetic distributions were similar to those of sequences from DGGE. Approximately 80% of clones were affiliated to Proteobacteria, Bacteroidetes and Cyanobacteria. The possible influences of dynamic pollution and hydrological conditions in the Victoria Harbor area on the particle-attached and free-living bacterial community structures were discussed.

Geophysical evaluation methods for buried hill reservoirs in the Jiyang superdepression of the Bohai Bay basin, eastern China

The Jiyang superdepression is one of the richest hydrocarbon accumulations in the Bohai Bay basin, eastern China. Comprehensive seismic methods have been used in buried hill exploration in Jiyang to describe these fractured reservoirs better. Accurate seismic stratigraphic demarcation and variable-velocity mapping were applied to reveal the inner structure of the buried hills and determine the nature of the structural traps more precisely. Based on the analysis of rock properties and the characteristics of well-developed buried hill reservoirs, we have successfully linked the geology and seismic response by applying seismic forward technology. Log-constrained inversion, absorption coefficient analysis and tectonic forward-inversion with FMI loggings were applied to analyse and evaluate the buried hill reservoirs and gave satisfying results. The reservoir prediction was successful, which confirmed that the comprehensive utilization of these methods can be helpful in the exploration of buried hill reservoirs.

基于正演方法的井间地震响应特征研究

地震勘探的基本任务是测量由人工震源激发的各种地震波，通过对所得到的地震波场的数据记录进行必要的处理和解释而提取所需要的地质信息。不同的地震波场特征反应了不同的地质构造等信息，因此，对地震波场特征进行研究就显得非常重要。井间地震是在一口井中激发、另一口井中接收的一种特殊的地震勘探方法，其实际记录要比常规的地震记录复杂的多，其波场中的各种波也更丰富复杂，如何识别不同类型的波并研究它们的响应特征是一个非常必要的任务。正演模拟可以为地震数据采集、处理、解释提供理论依据，是井间地震波场特征分析的一种有效方法。 本文主要研究了井间地震的数值正演模拟问题，系统地介绍了射线追踪正演模拟和波动方程正演模拟方法，利用Visual Basic 6.0语言编写了相应的射线追踪正演模拟程序，而且分两种方法详细分析了井间地震数值模型的地震响应特征，为直观的井间地震波场识别和勘探地震数据反演提供了必要的依据。本文取得的主要成果如下： 1)在前人的VSP正演模拟的基础上，改进并总结出了一套井间地震射线追踪正演模型的制作方法。 2)编制了井间地震射线追踪正演模拟程序V1.0。 3)利用编写的井间地震正演模拟程序，作了两种射线追踪数值模拟试算；利用Tesseral声波方程正演模拟试算了含有楔形体和丘状体的理论模型；取得了较好的模拟数据。

Capacity of soil to protect organic carbon and biochemical characteristics of density fractions in Ziwulin Haplic Greyxems soil

Physical protection is one of the important ways to stabilize organic carbon in soils. In order to understand the role of soils as a carbon sink or source in global climatic change and carbon cycles and properly manage soils as a carbon sink, we ought to know how many organic carbon (OC) in a given soil could be protected. By a density fractionation approach and ultrasonic technique, each soil sample was divided into three fractions: free light fraction (free-LF), occluded fraction (occluded-LF) and heavy fraction (HF). The obtained fractions were analyzed for total OC content, carbohydrate content and recalcitrant OC content. The results showed: (i) In the whole soil profile, dominance of OC consistently decreased in the following order: HF, free-LF, occluded-LF. This suggested that OC in soils were mostly protected. From 0-10 to 60-80 cm horizons, the OC in free-LF decreased from 25.27% to 3.72%, while OC in HF they were increased from 72.57% to 95.39%. The OC in occluded-LF was between 2.16% and 0.89%. (ii) Organic carbon recalcitrance in free-LF was similar to that in HF, and was even higher than that in HF below the surface horizon. This suggested that free-LF was not always the most fresh and non-decomposed fraction. OM quality of HF was higher than that of free-LF in the surface 10 cm below, namely the protected OM had higher quality than free OM in these horizons.

裂缝性油气藏储层预测方法研究

In this paper, based on the E & P situation in the oilfield and the theory of geophysical exploration, a series researches are conducted on fracture reservoir prediction technology in general，and it especially focus on some difficult points. The technological series which integrated amplitude preserved data processing、interpretation and its comprehensive application research as a whole were developed and this new method can be applied to the other similar oilfield exploration and development. The contents and results in this paper are listed as follows: 1. An overview was given on the status and development of fracture reservoir estimation technique, compare and analyze those geophysical prediction methods. This will be very helpful to the similar reservoir researches. 2. Analyze and conclude the characters of geologies and well logging response of burial hills fracture reservoir, those conclusions are used to steer the geophysical research and get satisfying results. 3. Forward modeling anisotropy seismic response of fracture reservoir. Quantitatively describe the azimuthal amplitude variation. Amplitude ellipse at each incidence angle is used to identify the fracture orientation. 4. Numerical simulation of structure stress based on finite difference method is carried out. Quantitatively describe and analyze the direction and intensity of fracture. 5. Conventional attributes extraction of amplitude preserved seismic data、attributes with different azimuthal angle and different offset are used to determine the relationship between the results and fracture distribution. 6. With spectrum decomposition method based on wavelet transform， the author disclose the reservoir distribution in space. It is a powerful tool to display its anisotropy. 7. Integrated seismic wave impendence、elastic impendence、spectrum decomposition、attribute extraction、fracture analysis result as a whole to identify and evaluate the fracture reservoir. An optimum workflow is constructed. It is used to practical oil&gas production and good results are obtained. This can indicate the wide foreground of this technique series.

辽东湾隐蔽油气藏地球物理识别技术研究

In China and world, more than half the recent basin discovered reserves involve lithologic hydrocarbon reservoir reserves. The major target for further hydrocarbon basin exploration is the subtle reservoir. The Liaodong Bay prospect is much important in Bohai Sea, which includes Liaoxi low uplift, Liaodong uplift, Liaoxi sag and Liaozhong sag. After dozens years’ exploration in Liaodong Bay, few unexplored big-and-middle-sized favorable structural traps are remained and most of the stock structure targets are bad for fragmentary. Thus seeking for new prospect area and making a breakthrough, have become the unique way to relieve the severe exploration condition in Liaodong Bay. Technique Route Based on the petrophysical property of target area, the seismic forward inference of typical subtle trap model is expanded with analysis of logging, seismic and geologic data. According to petrophysical characteristics and forward inference and research on seismic response of actual seismic data in target area, the optimization of geophysical technique is used in subtle trap identification and the geophysical identification technique system of subtle reservoir is formed. The Key Research ① Petrophysical Model The petrophysical parameter is the basic parameter for seismic wave simulation. The seismic response difference of rocks bearing different fluids is required. With the crossplot of log data, the influence of petrophysical parameters on rock elastic properties of target area is analyzed, such as porosity, shale index, fluid property and saturation. Based on the current research on Biot-Gassmann and Kuster-Toksoz model, the petrophysical parameter calculator program which can be used for fluid substitution is established. ② S-wave evaluation based on conventional log data The shear velocity is needed during forward inference of AVO or other elastic wave field. But most of the recent conventional log data is lack of shear wave. Thus according to the research on petrophysical model, the rock S-wave parameter can be evaluated from conventional log data with probability inverse method. ③ AVO forward modeling based on well data For 6 wells in JZ31-6 block and 9 wells in LD22-1 block, the AVO forward modeling recording is made by log curve. The classification of AVO characteristics in objective interval is made by the lithologic information. ④ The 2D parameter model building and forward modeling of subtle hydrocarbon trap in target area. According to the formation interpretation of ESS03D seismic area, the 2D parameter model building and seismic wave field forward modeling are carried on the given and predicted subtle hydrocarbon trap with log curve. ⑤ The lithology and fluid identification of subtle trap in target area After study the seismic response characteristics of lithology and fluid in given target area, the optimization of geophysical technique is used for lithology identification and fluid forecast. ⑥The geophysical identification technique system of subtle reservoir The Innovative Points of this Paper ① Based on laboratory measurement and petrophysical model theory, the rock S-wave parameter can be evaluated from conventional log data with probability inverse method. Then the fluid substitution method based on B-G and K-T theory is provided. ② The method and workflow for simulating seismic wave field property of subtle hydrocarbon trap are established based on the petrophysical model building and forward modeling of wave equation. ③ The description of subtle trap structural feature is launched. According to the different reflection of frequency wave field structural attribute, the fluid property of subtle trap can be identified by wave field attenuation attribute and absorption analysis. ④ It’s the first time to identify subtle trap by geophysical technique and provide exploration drilling well location. ⑤ The technique system of subtle reservoir geophysical identification is formed to provide available workflow and research ideas for other region of interest.

基于广义正交多项式迭积微分算子的地震波场数值模拟研究

Seismic Numerical Modeling is one of bases of the Exploratory Seismology and Academic Seismology, also is a research field in great demand. Essence of seismic numerical modeling is to assume that structure and parameters of the underground media model are known, simulate the wave-field and calculate the numerical seismic record that should be observed. Seismic numerical modeling is not only a means to know the seismic wave-field in complex inhomogeneous media, but also a test to the application effect by all kinds of methods. There are many seismic numerical modeling methods, each method has its own merits and drawbacks. During the forward modeling, the computation precision and the efficiency are two pivotal questions to evaluate the validity and superiority of the method. The target of my dissertation is to find a new method to possibly improve the computation precision and efficiency, and apply the new forward method to modeling the wave-field in the complex inhomogeneous media. Convolutional Forsyte polynomial differentiator (CFPD) approach developed in this dissertation is robust and efficient, it shares some of the advantages of the high precision of generalized orthogonal polynomial and the high speed of the short operator finite-difference. By adjusting the operator length and optimizing the operator coefficient, the method can involve whole and local information of the wave-field. One of main tasks of the dissertation is to develop a creative, generalized and high precision method. The author introduce convolutional Forsyte polynomial differentiator to calculate the spatial derivative of seismic wave equation, and apply the time staggered grid finite-difference which can better meet the high precision of the convolutional differentiator to substitute the conventional finite-difference to calculate the time derivative of seismic wave equation, then creating a new forward method to modeling the wave-field in complex inhomogeneous media. Comparing with Fourier pseudo-spectral method, Chebyshev pseudo-spectral method, staggered- grid finite difference method and finite element method, convolutional Forsyte polynomial differentiator (CFPD) method has many advantages: 1. Comparing with Fourier pseudo-spectral method. Fourier pseudo-spectral method (FPS) is a local operator, its results have Gibbs effects when the media parameters change, then arose great errors. Therefore, Fourier pseudo-spectral method can not deal with special complex and random heterogeneous media. But convolutional Forsyte polynomial differentiator method can cover global and local information. So for complex inhomogeneous media, CFPD is more efficient. 2. Comparing with staggered-grid high-order finite-difference method, CFPD takes less dots than FD at single wave length, and the number does not increase with the widening of the studying area. 3. Comparing with Chebyshev pseudo-spectral method (CPS). The calculation region of Chebyshev pseudo-spectral method is fixed in , under the condition of unchangeable precision, the augmentation of calculation is unacceptable. Thus Chebyshev pseudo-spectral method is inapplicable to large area. CFPD method is more applicable to large area. 4. Comparing with finite element method (FE), CFPD can use lager grids. The other task of this dissertation is to study 2.5 dimension (2.5D) seismic wave-field. The author reviews the development and present situation of 2.5D problem, expatiates the essentiality of studying the 2.5D problem, apply CFPD method to simulate the seismic wave-field in 2.5D inhomogeneous media. The results indicate that 2.5D numerical modeling is efficient to simulate one of the sections of 3D media, 2.5D calculation is much less time-consuming than 3D calculation, and the wave dispersion of 2.5D modeling is obviously less than that of 3D modeling. Question on applying time staggered-grid convolutional differentiator based on CFPD to modeling 2.5D complex inhomogeneous media was not studied by any geophysicists before, it is a fire-new creation absolutely. The theory and practices prove that the new method can efficiently model the seismic wave-field in complex media. Proposing and developing this new method can provide more choices to study the seismic wave-field modeling, seismic wave migration, seismic inversion, and seismic wave imaging.

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

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.