Thermally induced evolution of phase transformations in gas hydrate sediment

Thermally induced evolution of phase transformations is a basic physical-chemical process in the dissociation of gas hydrate in sediment (GHS). Heat transfer leads to the weakening of the bed soil and the simultaneous establishment of a time varying stress field accompanied by seepage of fluids and deformation of the soil. As a consequence, ground failure could occur causing engineering damage or/and environmental disaster. This paper presents a simplified analysis of the thermal process by assuming that thermal conduction can be decoupled from the flow and deformation process. It is further assumed that phase transformations take place instantaneously. Analytical and numerical results are given for several examples of simplified geometry. Experiments using Tetra-hydro-furan hydrate sediments were carried out in our laboratory to check the theory. By comparison, the theoretical, numerical and experimental results on the evolution of dissociation fronts and temperature in the sediment are found to be in good agreement.

功能菌剂的混合发酵研究

随着化工行业的发展，大量有毒有害难降解有机物随工业废水的排放进入环境，这些物质能够在环境中长期存在、积累和扩散，通过食物链对动植物的生存及人类的健康造成不良影响。本文以苯酚、对氯硝基苯、氯苯和十六烷为模拟污染物，以前期研制的功能菌剂为对象，经过紫外线线诱变筛选出优于出发菌株的功能菌，对诱变后功能菌的理化性能进行了研究，对菌种进行了鉴定，在此基础上，就其相互之间的微生态关系进行研究，为混合发酵提供理论基础，并就其最佳发酵条件及发酵参数进行了研究，最后对发酵产品的性能进行了检测。目前，国内外有关功能菌剂的研究还存在多方面的不足，主要包括：①由于多菌种混合发酵过程较为复杂，各菌之间存在复杂的相互作用，影响因素较多，关于菌种之间的相互关系研究得很少，环境功能菌剂的发酵方法大多采用单独发酵后混合的方式。单独发酵对原材料、设备和能源的利用率较低，对于多菌种制剂发酵，在设备、能源和原材料的方面造成的浪费更大，将会大幅增加菌剂的生产成本，影响多菌种功能菌剂的发展；②功能菌剂生产过程的质量控制方面研究得较少；③功能菌剂产品的稳定性、抗冲击性能研究得较少，对环境微生物制剂的研究主要集中在菌种选育和培养条件优化方面。 通过本论文研究，得到以下主要结论。 （1）在紫外线诱变处理中，用紫外线对发生一定程度退化的出发菌株进行诱变处理后，六株具有高效降解性能的菌株被筛选出来，诱变筛选出的菌株形态和ERIC-PCR指纹图谱与出发菌株相比发生了明显改变；而且诱变后的菌株对目标难降解底物的降解能力均得到改善，其中，FPN、FCB、F14、FEm对目标底物的降解率提高了20%以上；诱变后菌株经过7次连续传代接种后，对目标难降解底物的降解率无显著变化，具有一定的遗传稳定性。并对诱变后的功能菌进行了初步的鉴定，这6株菌都分别是芽孢杆菌。 （2）对诱变后的功能菌相互之间的微生态关系进行了研究，通过抑菌实验、生长量以及基质消耗量的比较，确定它们之间的生长关系是无害共栖关系，可以进行混合发酵。 （3）对该功能菌剂进行发酵培养条件研究，结果表明发酵培养基的最佳成分（g/L）：葡萄糖 31.0g/L、玉米粉10.0g/L、磷酸氢二钾1.0g/L、硫酸铵1.1g/L、硫酸镁0.55g/L。通过研究不同的培养条件对菌体生长和降解性能的影响，确定了最佳培养条件：培养基初始pH7.5；最适温度32℃；培养基装液量125mL（250 mL三角瓶），以及培养时间对降解性能的影响，培养20 h的产物对降解最为有利。通过研究添加不同目标污染物对菌体生长和降解性能的影响，确定了添加目标污染物的最佳量以及最佳时间：苯酚投加量：1.125 g/L，对氯硝基苯投加量：0.1 g/L；最佳投加时间为发酵培养开始后4 h。 （4）以摇瓶分批发酵最优条件为基础，对FPN、F10、FCB、FNa、F14 和 FEm进行了摇瓶分批发酵试验。以摇瓶分批发酵试验数据为依据，对功能菌剂分批发酵动力学进行了研究，建立了菌体生长和基质消耗的动力学模型，拟合模型能较好的反映功能菌剂分批发酵过程。 （5）功能菌剂和活性污泥协同作用，可以提高系统的生物降解能力，功能菌剂投加量为2%，新鲜活性污泥3500 mg/L，降解24 h条件下，功能菌剂和活性污泥的协同作用对COD的去除率和对照组相比，最多的提高了36.8%。功能菌剂和活性污泥协同作用以及活性污泥的单独作用，其生物降解过程均符合一级反应动力学过程，功能菌剂和活性污泥协同作用的生物降解动力学方程为：，相关系数97%。采用SBR运行方式，引入功能菌剂的SBR系统明显能够改善和提高生物降解的效率。与仅有活性污泥的系统相比，系统对COD的平均去除率可以提高27.1%，同时，系统的耐负荷冲击以及耐毒害冲击的性能比仅有活性污泥的SBR系统强，特别是负荷冲击对引入功能菌剂的SBR系统影响很小。仅有活性污泥的SBR系统经过负荷冲击和毒害冲击之后，不能恢复到冲击之前的水平，而且系统有效作用时间的周期比引入功能菌剂的SBR系统相比大大缩短，而引入功能菌剂的SBR系统处理效果较为稳定，恢复能力很强。 Along with the development of industries, many recalcitrant organic chemicals have been discharged into natural environments together with wastewaters and can exist in waters, soil and sediments for a long time without degradation. These haz-ardous substances, their byporducts and metabolizabilities can be highly toxic, mu-tagenic and carcinogenic, thereby threatening animals, plants and human health through food chain. Consequently the removal of these compounds is of significant interest in the area of wastewater treatment. In this dissertation, the phenol, hydro-quinone, chlorobenzene and hexadecane treated as the model pollutants, the func-tional microorganism agent was used as the starting strains, they treated with ultra-violet light, and then the mutant strains with high degradation ability were screened out and identified primarily, the relationship between these stains were studied, the medium composition and fermentation conditions were optimized, the degradation ability of the fermented production was tested. The literature survey indicates that the study of the microorganism agent is far from complete and more information is re-quired on following problems. 1, Because of the complexity of relationship in mixed fermentation and the complicated factors, the study is hardly to process.2, There is a lack of information on the quality control of the producing process .3, And there is a lack of information on the stability about the microorganism agent. In this dissertation, the main results of the present study could be summarized as follows: (1)The degenerate starting strains were treated with the ultraviolet light, and six mutant strains with high biodegradation ability were screened out by using the me-dium with selective pressure of model pollutants. The mutant strains had great changes in colonialmorphology and ERIC-PCR fingerprinting. And the mutant strains got obvious advantages over the starting strains in degradation ability and over 20% improvement of removal rates was achieved for FPN、FCB、F14 and FEm. The de-gradation ability of the mutant strains was stable after seven generations. After that, the mutant strains were primarily identified as bacillus respectively. (2) The relationship between these mutant strains was studied. By the compari-son of antibiosis effect, biomass and consumption of substrate, the relationships were neutralism and they could be mixed fermented. (3) The optimized cultivation conditions were as follows: glucose 31.0 g/L, corn power 10 g/L, K2HPO4 1.0 g/L, (NH4)2SO4 1.1 g/L, MgSO4 0.55 g/L, initial pH7.5, temperature 32℃, working volume 125 mL/250 mL, and cultivation time 20h (con-sidering the time effect on degradation ability), adding pollutants phenol (1.125 g/L) and hydroquinone (0.1 g/L) into the broth at 4 h after cultivation. (4) Based on the above optimum condition, the batch fermentation was per-formed with strains FPN, F10, FCB, FNa, F14 and FEm in shake flask. The batch fermentation kinetics was studied based on the experimental data. Two kinetic models were constructed which could reflect the regularity of growth and substrate consump-tion in the process of batch fermentation. (5) The co-operation of functional microorganism agent and activated sludge could raise biodegradation of system by adding some microorganism agent and 3500 mg/L fresh activated sludge. Bioaugumentation by the addition of high effective deg-radation culture enhanced the treatment effect of SBR system and the COD removal rate was increased by 20%-36.8%. Its biodegradation matched first-order dynamical reaction equation, and the reaction equation was ln0.2327.391ct=−+. The micro-organism agent had the effect of optimization to activated sludge micro-ecosystem. The SBR system adding 2% microorganism agent, the average COD removal rate of that was increased by 27.1% and stronger anti-shock ability to load and toxicant were achieved (compared with SBR system just adding activated sludge). Especially the load-shock has barely effect to the SBR system adding microorganism agent. After the load and toxicant shock, the SBR system just adding activated sludge couldn’t come back to original level and the activated sludge micro-ecosystem was frustrated. The applying of microorganism agent increased biological activity and system’s re-sistance ability to load shock and toxicant shock.

Constrain intergalactic medium from the SZ effect map

In this paper, we try to detect the SZ effect in the 2MASS DWT clusters and less bound objects in order to constrain the warm-hot intergalactic medium distribution on large scales by cross-correlation analysis. The results of both observed WMAP and mock SZ effect map indicate that the hot gas distributes from inside as well as outside of the high density regions of galaxy clusters, which is consistent with the results of both observation and hydro simulation. Therefore, the DWT measurement of the cross-correlation would be a powerful tool to probe the missing of baryons in the universe.

黄土丘陵沟壑区林地水文生态效应

根据安塞水土保持试验站1993～2002年林地径流小区的降雨产流产沙的定位观测资料及2002年土壤含水量资料,分析了不同树种对坡面尺度降雨产流产沙及土壤水分的影响。结果表明:场降雨径流小区的产流量、产沙量与降雨量具有较好的相关性;多元回归分析表明,场降雨产流量和产沙量与降雨量和最大30min雨强的乘积呈正相关,与植被覆盖度呈负相关,场降雨产沙量回归方程复相关系数为0.253,各处理场降雨产流量回归方程复相关系数的变化范围为0.465～0.723,均达到了极显著的水平(P<0.01)。同时,各树种均具有良好的减流减沙功能,与农地相比,年均产流量和产沙量分别减少4.8%～52.9%和26.8%～86.0%;沙棘纯林及其混交林的减流减沙效益优于油松纯林。同时,沙棘纯林及其各混交林在造林初期就表现出良好的减流减沙效益,随着树龄的增长,其作用更加明显;而油松纯林在造林初期作用不明显,甚至出现产流量和产沙量大于农地的现象,但随着树龄的增长,减流减沙作用逐渐呈现并增大。沙棘纯林及其混交林30cm以下土壤含水量在整个生长季中均呈递减趋势,生长季初(4月份)土壤含水量最高,而生长季末(10月份)降到最低值。2002年沙棘纯林的...

小流域水分行为、生态效应及其优化调控研究方面的若干问题

把小流域水文生态研究 ,分解为小流域水分行为、生态效应及其优化调控研究 3个方面。在分析国内外研究现状与问题的基础上 ,指出了其需要系统加以探讨的若干内容 ,认为由此可构建流域水文生态学的基本框架 ,进而发展流域水文生态学。

黄土高原森林水文生态效应和林草适宜覆盖指标

概述了黄土高原森林各层次的水文生态功能和水土保持效益。在此基础上 ,根据自然条件、地形地貌和国民经济发展要求 ,提出了为保护生态环境和协调经济发展的适宜林草覆盖指标

ANUDEM和TIN两种建立DEM方法的对比研究

针对流域水文和土壤侵蚀定量模拟分析需要,对基于TIN和Hutchinson方法建立的DEM表现地形形态和起伏的能力进行了对比研究。结果表明,基于TIN建立的DEM始终存在一些平顶现象,一些较小的侵蚀沟被忽略,其上提取的河流不完全连续,多处出现多重线条河流,因而不能如实地反映地形起伏的细部特征。而基于ANUDEM建立的DEM,其派生等高线的形状与输入等高线吻合较好,较好地表现了地形的形态和起伏,对地形和坡度的反映更加连续和光滑,其上提取的河流信息基本与地形图上的河流一致。

Supramolecular isomerism in zinc hydroxide coordination polymers with pyridine-2,4-dicarboxylic acid: Two polymorphs with centrosymmetric two-dimensional and acentric three-dimensional coordination networks

Reactions of Zn(BF4)(2) and pyridine-2,4-dicarboxylic acid (2,4-pydcH(2)) in the presence of 1,2-bis( 4-pyridyl) ethylene or 1,3-bis(4-pyridyl) propane under hydro(solvo) thermal conditions yielded two polymorphic metal-organic coordination polymers formulated as Zn-2(OH)(2)(2,4-pydc) (1 and 2). Polymorph 1 features a two-dimensional (2-D) layer-like structure that is constructed by 2,4-pydc ligands bridging between the Zn-OH-Zn double-chain units. Each single Zn-OH-Zn chain is composed of mu(2)-OH groups connecting trigonal bipyramidal and tetrahedral Zn centers. Polymorph 2 is a 3-D coordination polymer containing 2-D Zn-OH-Zn sheets that consist of mu(2)- and mu(3)-OH groups and trigonal bipyramidal Zn centers. The sheets are pillared by 2,4-pydc ligands to form an acentric structural architecture. 1 and 2 are rare examples that the two polymorphs exhibit a centrosymmetric 2-D coordination network and an acentric 3-D coordination network, respectively. The different structures lead to differences in photoluminescent properties and thermal stabilities for 1 and 2.

Fabrication and characterization of hybrid Langmuir-Blodgett films of polyquinoline/stearic acid/rare earth-substituted heteropolymolybdate

Three kinds of hybrid Langmuir-Blodgett films are obtained by the organization of poly(1-hydro-2,2,4-trimethyl)quinoline (PQ), stearic acid(SA) and rare earth-substituted heteropolymolybdates (RE(PMo11)(2), RE = Ce-III, Eu-III, La-III) using the Langmuir-Blodgett technique. They are characterized by pi-A isotherms, absorption spectra, fluorescence spectra, IR and atomic force microscope. The absorption spectra indicate that the molecules of PQ and heteropolymolybdates are incorporated into the LB films. The atomic force microscope reveals that heteropolymolybdates aggregate at the surface of the LB film.

降雨条件下碎石土滑坡稳定性研究

Debris Landslide is one of the types of landslides with the widest distribution, largest quantity, and the closest relationship with engineering construction. It is also one of the most important types of landslides that can cause disaster. This kind of landslide often occurs in the loose slopes which are made up of loose congeries formed by earth filling, residual soil, slope wash, dilapidation, landslide or full weathered material of hard rock. Rainfall is always the chief inducing factor of debris Landslide. Therefore, to research stability of debris Landslide during rainfall not only has important theoretical significance for understanding developing law and deformation and failure mechanism of debris landslide, but also has important practical significance for investigating, appraising, forecasting, preventing and controlling debris landslides. This thesis systematically summarized the relationships between rainfall and landslide, the method to survey water table in the landslides, the deformation and failure mechanism of debris landslide, and the progress in the stability analysis of landslides based on the analyses of data collected widely at home and abroad. The problems in the study of the stability of debris landslide during rainfall was reviewed and discussed. Due to the complicated geological conditions and the random rainfall conditions, the research on the landslides' stability must be based on engineering geological qualitative analysis. Through the collection of the data about the Panxi region and the Three Gorges Reservoir region, the author systematically summarized the engineering geological conditions, hydro-geological condition, distribution characteristics of stress field in the slope, physical and mechanical properties and hydro-mechanical properties of debris. In the viewpoint of dynamics of soil water and hydromechanics, physical process of rainfall to supply groundwater of debris landslides can be divided into two phases, i.e. non-saturated steady infiltrating phase and saturated unsteady supplying phase. The former can be described by mathematical model of surface water infiltration while the latter can be described by equivalent continuous medium model of groundwater seepage. With regard to specific hydrological geology system, we can obtain the dynamic variation law of water content, water table, landslide stability of rock and soil mass, along with quantity and duration of rainfall after the boundary condition on hydrological geology has been ascertained. This is a new way to study the response law of groundwater in the landslides during rainfall. After wet face of rock and soil mass connects with ground water table, the raising of water table will occur due to the supply of rainfall. Then interaction between ground water and rock and soil mass will occur, such as the action of physics, water, chemistry and mechanics, which caused the decrease of shearing strength of sliding zone. According to the action of groundwater on rock and soil mass, a concise mechanical model of debris landslide’s deformation was established during rainfall. The static equilibrium condition of landslide mass system was achieved according to the concise mechanical model, and then the typical deformation and failure process and failure mode of debris landslide during rainfall were discussed. In this thesis, the former limiting equilibrium slice method was modified and improved based on shearing strength theory of , a stability analysis program of debris landslide was established and developed taking account of the saturated-unsaturated seepage, by introducing the shearing strength theory of unsaturated soil mass made by (1978). The program has reasonable data storage and simple interface and is easy to operate, and can be perfectly used to carry out sensitivity analysis of influencing factors of landslides' stability, integrated with the program of Office Excel. The design of drainage engineering are always bases on empirical methods and is short of effective quantitative analysis and appraise, therefore, the conception of critical water table of debris landslide was put forward. For debris landslides with different kinds of slide face in the engineering practice, a program to search the critical water table of debris landslide was developed based on native groundwater table. And groundwater table in the slope should be declined below the critical water table in the drainage works, so the program can be directly used to guide drainage works in the debris landslide. Taking the slope deformation body in the back of former factory building of Muli Shawan hydroelectric power station as an example, a systematic and detailed research on debris landslides' stability during rainfall was researched systematically, the relationship among quantity of rainfall, water table and stability of slope was established, the debris landslides' stability in process of rainfall from dynamic viewpoint was analyzed and researched.

典型人类活动对边坡变形及稳定性的影响研究

China is a mountainous country in which geological hazards occurred frequently, especially in the east of China. Except the geology, topography and extreme climate, the large scale human activities have become a major factor to landslides. Typical human activities which induced landslides are fill, cut and underground mining. On the topic of the deformation mechanism and slope stability, taking three different man-made slopes as examples, deformation mechanism and slope stability were studied by several methods, such as field work, numerical modeling and monitor. The details are as following: (1) The numerical modeling approach advantages over other conventional methods such as limit methods, so the numerical modeling is the major tool in this thesis. So far, there is no uniform failure criterion for numerical simulation. The failure criterion were summarized and analyzed firstly, subsequently the appropriate criterion was determinated. (2) Taking 220kV Yanjin transformation substation fill slope as example, the deformable characteristic, unstable mode and laboratory tests were studied systematically. The results show: the slope deformation was probably caused by a combination effect of unfavorable topographic, geological and hydro geological conditions, and external loading due to filling. It was concluded that the creep deformation of the slope was triggered by external loading applied at the back of the slope. In order to define the calculating parameters, a set of consolidated drained (CD) tests, consolidated undrained (CU) tests, repeated direct shear tests and UCS tests were carried out. The stability of the slope before and after reinforcement was assessed using 3D numerical modeling and shear strength reduction technique. The numerical modeling results showed: the factor of safety (FOS) of the slope was 1.10 in the natural state, and reduced to 1.03 after fill, which was close to the critical state and it caused creeping slip or deformation under rainfall. The failure surface in the slope is in active shear failure, whereas tensile failure occurs at the slope crest. After the site was reinforced with piles, the FOS was 1.27. Therefore, the slope is stable after reinforcement measures were taken. (3) The cut slope stability is a complex problem. Taking the left cut slope of Xiangjiaba as example in this thesis, the deformation and slope stability were studied systematically by numerical modeling and monitor methods. The numerical results show: the displacement is gradually increasing along with the cutting, and the largest displacement is 27.5mm which located at the bench between the elevation 340 and 380. Some failure state units distribute near the undermining part and there is no linked failure state occurred from crest to bottom during cutting. After cutting, some failure units appeared at the ground surface between elevation 340 and 360. The increasing tense stress made the disturbed rock failed. The slope is stable after cutting by the monitor method, such as surface monitor, multipoint displacement meter, inclinometer and anchor cable tensometer. (4) The interaction between underground mining and slope stability is a common situation in mountainous. The slope deformation mechanism induced by underground mining may contributed significantly to slope destabilization. The Mabukan slope in xiangjiaba was analyzed to illustrate this. Failure mechanism and the slope stability were presented by numerical modeling and residual deformation monitor. The results show: the roof deformed to the free face and the floor uplift lightly to the free face. The subsidence basin is formed, but the subsidence and the horizontal movement is small, and there is no failure zone occurred. When the underground mining is going on, the roof deformation, subsidence and the horizontal movements begin increasing. The rock deformation near the free face is larger than the ground surface, and the interaction between these coal seams appeared. There are some tensile failures and shear failures occurred on the roof and floor, and a majority of failure is tensile failure. The roof deformation, subsidence and the horizontal movements increased obviously along with the underground mining. The failure characteristic is shear failure which means the tensile stress transformed to the compressive stress. So the underground mining will induced tensile stress first which lead to structure crack, subsequently the compressive stress appeared which result in slippage. The crest was subjected to horizontal tension which made the rock crack along with the joint. The long term residual deformation monitor demonstrates that the slope is stable after the underground mining stopped.

海底金属矿床充水条件及矿坑突水的断层效应研究

The Xinli mine area of Sanshandao mine is adjacent to the Bohai Sea and its main exploitable ore deposit occurs in the undersea rock mass. The mine is the biggest undersea gold mine of China after production. The mine area faces a latent danger of water bursting, even sudden seawater inrush. There is no mature experience in undersea mining in China so far. The vein ore deposit is located in the lower wall of a fault; its possible groundwater sources mainly include bittern, Quaternary pore water and modern seawater. To ensure the safety of undersea mining, to survey the flooding conditions of the ore deposit using proper measures and study the potential seawater inrush pattern are the key technical problems. With the Xinli mine area as a case study, the engineering geological conditions of the Xinli mine area are surveyed in situ, the regional structural pattern and rock mass framework characteristics are found out, the distribution of the structural planes are modeled by a Monte Carlo method and the connectivity coefficients of rock mass structural planes are calculated. The regional hydro-geological conditions are analyzed and the in-situ hydro-geological investigation and sampling are performed in detail, the hydrochemistry and isotopes testing and groundwater dynamic monitoring are conducted, the recharge, runoff, discharge conditions are specified and the sources of flooding are distinguished. Some indices are selected from the testing results to calculate the proportion of each source in some water discharge points and in the whole water discharge of the Xinli mine area. The temporal and spatial variations of each water source of the whole ore deposit flooding are analyzed. According to the special project conditions in the Xinli mine area, the permeability coefficient tensors of the rock mass in Xinli mine area are calculated based on a fracture geometry measurement method, in terms of the connectivity and a few hydraulic testing results, a modified synthetic permeability coefficient are calculated. The hydro-geological conceptual and mathematical model are established，the water yield of mine is predicted using Visual Modflow code. The spreading law of surrounding rock mass deformation and secondary stress are studied by numerical analysis; the intrinsic mechanism of the faults slip caused by the excavation of ore deposit is analyzed. The results show that the development of surrounding rock mass deformation and secondary stress of vein ore deposit in the lower wall of a fault, is different from that in a thick-big ore deposit. The secondary stress caused by the excavation of vein ore deposit in the lower wall of a fault, is mainly distributed in the upper wall of the fault, one surface subsidence center will occur. The influences of fault on the rock mass movement, secondary stress and hydro-geological structures are analyzed; the secondary stress is blocked by the fault and the tensile stress concentration occurs in the rock mass near the fault, the original water blocking structure is destructed and the permeable structure is reconstructed, the primary structural planes begin to expand and newborn fissures occur, so the permeability of the original permeable structure is greatly enhanced, so the water bursting will probably occur. Based on this knowledge, the possible water inrush pattern and position of the Xinli mine area are predicted. Some computer programs are developed using object-oriented design method under the development platform Visual Studio.Net. These programs include a Monte Carlo simulation procedure, a joint diagrammatizing procedure, a structural planes connectivity coefficient calculating procedure, a permeability tensor calculating procedure, a water chemical formula edit and water source fixture conditions calculating procedure. A new computer mapping algorithm of joint iso-density diagram is raised. Based on the powerful spatial data management and icon functions of Geographic Information System, the pit water discharge dynamic monitoring data management information systems are established with ArcView.

地下水流动数值模拟的并行计算方法研究

Numerical modeling of groundwater is very important for understanding groundwater flow and solving hydrogeological problem. Today, groundwater studies require massive model cells and high calculation accuracy, which are beyond single-CPU computer’s capabilities. With the development of high performance parallel computing technologies, application of parallel computing method on numerical modeling of groundwater flow becomes necessary and important. Using parallel computing can improve the ability to resolve various hydro-geological and environmental problems. In this study, parallel computing method on two main types of modern parallel computer architecture, shared memory parallel systems and distributed shared memory parallel systems, are discussed. OpenMP and MPI (PETSc) are both used to parallelize the most widely used groundwater simulator, MODFLOW. Two parallel solvers, P-PCG and P-MODFLOW, were developed for MODFLOW. The parallelized MODFLOW was used to simulate regional groundwater flow in Beishan, Gansu Province, which is a potential high-level radioactive waste geological disposal area in China. 1. The OpenMP programming paradigm was used to parallelize the PCG (preconditioned conjugate-gradient method) solver, which is one of the main solver for MODFLOW. The parallel PCG solver, P-PCG, is verified using an 8-processor computer. Both the impact of compilers and different model domain sizes were considered in the numerical experiments. The largest test model has 1000 columns, 1000 rows and 1000 layers. Based on the timing results, execution times using the P-PCG solver are typically about 1.40 to 5.31 times faster than those using the serial one. In addition, the simulation results are the exact same as the original PCG solver, because the majority of serial codes were not changed. It is worth noting that this parallelizing approach reduces cost in terms of software maintenance because only a single source PCG solver code needs to be maintained in the MODFLOW source tree. 2. P-MODFLOW, a domain decomposition–based model implemented in a parallel computing environment is developed, which allows efficient simulation of a regional-scale groundwater flow. The basic approach partitions a large model domain into any number of sub-domains. Parallel processors are used to solve the model equations within each sub-domain. The use of domain decomposition method to achieve the MODFLOW program distributed shared memory parallel computing system will process the application of MODFLOW be extended to the fleet of the most popular systems, so that a large-scale simulation could take full advantage of hundreds or even thousands parallel processors. P-MODFLOW has a good parallel performance, with the maximum speedup of 18.32 (14 processors). Super linear speedups have been achieved in the parallel tests, indicating the efficiency and scalability of the code. Parallel program design, load balancing and full use of the PETSc were considered to achieve a highly efficient parallel program. 3. The characterization of regional ground water flow system is very important for high-level radioactive waste geological disposal. The Beishan area, located in northwestern Gansu Province, China, is selected as a potential site for disposal repository. The area includes about 80000 km2 and has complicated hydrogeological conditions, which greatly increase the computational effort of regional ground water flow models. In order to reduce computing time, parallel computing scheme was applied to regional ground water flow modeling. Models with over 10 million cells were used to simulate how the faults and different recharge conditions impact regional ground water flow pattern. The results of this study provide regional ground water flow information for the site characterization of the potential high-level radioactive waste disposal.

东海陆架盆地西南部烟囱构造与油气运聚关系研究

In recent years, chimney structure has been proved one of important indicators and a useful guide to major petroleum fields exploration through their exploration history both at home and abroad. Chimney structure, which has been called "gas chimney" or "seismic chimney", is the special fluid-filled fracture swarm, which results from the boiling of active thermal fluid caused by abruptly decreasing of high pressure and high temperature in sedimentary layers of upper lithosphere. Chimney structure is well developed in continental shelf basin of East China Sea, which indicates the great perspectives of petroleum resources there. However, the chimney structure also complicated the petroleum accumulation. So the study of chimney structure on its formation, its effect on occurrence and distribution of petroleum fields is very important not only on theoretical, but also on its applied research. It is for the first time to make a clear definition of chimney structure in this paper, and the existence and practical meaning of chimney structure are illustrated. Firstly, on the viewpoint of exploration, this will amplify exploration area or field, not only in marine, but also on continent. Secondly, this is very important to step-by-step exploration and development of petroleum fields with overpressure. Thirdly, this will provide reference for the study on complex petroleum system with multi-sources, commingled sources and accumulation, multi-stage accumulations, and multi-suits petroleum system in the overlay basin. Fourthly, when the thermal fluid enters the oceanic shallow layer, it can help form gas hydrate under favorable low-temperature and high-pressure conditions. Meanwhile, the thermal fluid with its particular component and thermal content will affect the physical, chemical and ecological environments, which will help solving the problem of global resources and environment. Beginning from the regional tectonic evolution characteristics, this paper discussed the tectonic evolution history of the Taibei depression, then made an dynamical analysis of the tectonic-sedimentary evolution during the Mesozoic and Cenozoic for the East China Sea basin. A numerical model of the tectonic-thermal evolution of the basin via the Basin-Mod technique was carried out and the subsidence-buried history and thermal history of the Taibei depression were inverse calculated: it had undergone a early rapid rift and sag, then three times of uplift and erosion, and finally depressed and been buried. The Taibei depression contains a huge thick clastic sedimentary rock of marine facies, transitional facies and continental facies on the complex basement of ante-Jurassic. It is a part of the back-arc rifting basins occurred during the Mesozoic and Cenozoic. The author analyzed the diagenesis and thermal fluid evolution of this area via the observation of cathodoluminescence, scanning electron microscope and thin section, taking advantage of the evidences of magma activities, paleo-geothermics and structural movement, the author concluded that there were at least three tectonic-thermal events and three epochs of thermal-fluid activities; and the three epochs of thermal-fluid activities were directly relative to the first two tectonic-thermal events and were controlled by the generation and expulsion of hydrocarbon in the source rock simultaneously. Based on these, this paper established the corresponding model between the tectonic-thermal events and the thermal-fluid evolution of the Taibei Depression, which becomes the base for the study on the chimney structures. According to the analyses of the gas-isotope, LAM spectrum component of fluid inclusion, geneses of CO_2 components and geneses of hydrocarbon gases, the author preliminarily verified four sources of the thermal fluid in the Taibei Depression: ① dehydration of mud shale compaction, ② expulsion of hydrocarbon in the source rock; ③ CO_2 gas hydro-thermal decomposition of carbonatite; ④magma-derived thermal fluid including the mantle magma water and volatile components (such as H_2O, CO_2, H_2S, SO_2, N_2 and He etc.). On the basis of the vitrinite reflectance (Ro), homogenization temperature of fluid inclusion, interval transit time of major well-logging, mud density of the wells, measured pressure data and the results of previous studies, this paper analyzed the characteristics of the geothermal fields and geo-pressure fields for the various parts in this area, and discussed the transversal distribution of fluid pressure. The Taibei depression on the whole underwent a temperature-loss process from hot basin to cold basin; and locally high thermal anomalies occurred on the regional background of moderate thermal structure. The seal was primarily formed during the middle and late Paleocene. The overpressured system was formed during the middle and late Eocene. The formation of overpressured system in Lishui Sag underwent such an evolutionary process as "form-weaken-strengthen-weaken". Namely, it was formed during the middle and late Eocene, then was weakened in the Oligocene, even partly broken, then strengthened after the Miocene, and finally weakened. The existence of the thermal fluid rich in volatile gas is a physical foundation for the boiling of the fluid, and sharply pressure depletion was the major cause for the boiling of the fluid, which suggests that there exists the condition for thermal fluid to boil. According to the results of the photoelastic simulation and similarity physical experiments, the geological condition and the formation mechanism of chimnestructures are summarized: well compartment is the prerequisite for chimney formation; the boiling of active thermal fluid is the original physical condition for chimney formation; The local place with low stress by tension fault is easy for chimney formation; The way that thermal fluid migrates is one of the important factors which control the types of chimney structures. Based on where the thermal fluid come from and geometrical characteristics of the chimney structures, this paper classified the genetic types of chimney structures, and concluded that there existed three types and six subtypes chimney structures: organic chimney structures generated by the hydrocarbon-bearing thermal fluid in middle-shallow layers, inorganic and commingling-genetic chimney structures generated by thermal fluid in middle-deep layers. According to the seismic profiles interpretations, well logging response analysis and mineralogical and petrological characteristics in the study area, the author summarized the comprehensive identification marks for chimney structures. Especially the horizon velocity analysis method that is established in this paper and takes advantage of interval velocity anomaly is a semi-quantitative and reliable method of chimney structure s identification. It was pointed out in this paper that the occurrence of the chimney structures in the Taibei depression made the mechanism of accumulation complicated. The author provided proof of episodic accumulation of hydrocarbon in this area: The organic component in the boiling inclusion is the trail of petroleum migration, showing the causality between the boiling of thermal fluid and the chimney structures, meanwhile showing the paroxysmal accumulation is an important petroleum accumulation model. Based on the evolutionary characteristics of various types of chimney structures, this paper discussed their relationships with the migration-accumulation of petroleum respectively. At the same time, the author summarized the accumulating-dynamical models associated with chimney structures. The author analyzed such accumulation mechanisms as the facies state, direction, power of petroleum migration, the conditions of trap, the accumulation, leakage and reservation of petroleum, and the distribution rule of petroleum. The author also provides explanation for such practical problems the existence of a lot of mantle-derived CO_2, and its heterogeneous distribution on plane. By study on and recognition for chimney structure, the existence and distribution of much mantle-derived CO_2 found in this area are explained. Caused by tectonic thermal activities, the deep magma with much CO_2-bearing thermal fluid migrate upward along deep fault and chimney structures, which makes two wells within relatively short distance different gas composition, such as in well LF-1 and well LS36-1-1. Meanwhile, the author predicted the distribution of petroleum accumulation belt in middle-shallow layer for this area, pointed out the three favorable exploration areas in future, and provided the scientific and deciding references for future study on the commingling-genetic accumulation of petroleum in middle-deep layer and the new energy-gas hydrate.

面向决策支持的综合地质信息系统功能开发和GIS应用

Synthetic Geology Information System(SGIS) is a part of the theory of Engineering Geomechanics-mate-Synthetics(EGMS), is also a development of its technical methodology. SGIS includes ways of geology engineering investigation, design, and construction. Although SGIS has an integrate theory frame, and some parts of it have gained great progress, the completion of SGIS is a continuous and accumulative process. This paper analyses the ways and principle of building knowledge database and model database, summarizes the experts' experience on exploration methods selection and the characters of exploration models, combining with the application of Decision Support System(DSS) in Decision support of Synthetic Exploration Methods for Railway engineering Geology. By the analysis of hierarchy structure of the model database, the effects of geology engineering factors on the selection of exploration methods are expressed. By the usage of fuzzy patterns recognize, hierarchy structure analysis, fuzzy collection closement analysis etc, the software of DSS for engineering design and construction are developed. At same time, by the development of Monitoring Data Analysis System and experiment data management system of Hydro-power project, this paper discussed the data management of science experiment of Hydro-power project by the usage of synthetic database and the usage of Geography Information System(GIS) and DSS technics. The technic of visual operation of data process and project monitoring system are presented. The intelligence algorithm of self-adoption is carried out to improve the data process and analysis of monitoring. Items of the project theoretical analysis and data process are designed in detail. All the theory and technical methods presented in this paper are one part of SGIS, in which the application of DSS and GIS, is an important step of the progress and completion of SGIS.