栉孔扇贝大规模死亡病原学研究

用平板画线法从患病栉孔扇贝（Chlamys farreri）体内分离到了一种原核生物（简称QDP）。QDP可以在改进的液体培养基MEM（含2.2%NaCl，5%小牛血清）和脑心浸液（含2.2% NaCl）中生长；菌落在显微镜下（150×）为无色、透明的小点状；革兰氏染色阴性；菌体为圆形或近似圆形。QDP在发育过程中有两种状态，一种为未成熟阶段，直径小于100nm；另一种为成熟阶段，直径变化很大，最小约60nm，最大可达4µm以上。较小的个体有拟核、核糖体和新月状的空泡，未见细胞壁；较大的个体有细胞壁，胞内大部分被空泡充满，未见拟核和核糖体。栉孔扇贝组织超簿切片电镜观查证实QDP的存在。QDP的密度随着生长发育时间的不同而有所变化，繁殖高峰期密度较大。 建立了密度梯度离心结合滤膜过滤分离技术，优化人工培养条件。最适生长温度为23℃，最适生长pH值为7.4，最适生长盐度相当于细胞培养液所需的盐浓度（0.85%NaCl）。 提取的QDP核酸能被RNase A 降解，且没有检测到DNA。以PCR、RT-PCR扩增其16SrRNA基因序列片段，PCR反应没有扩增出扩增子，而RT-PCR则扩增出了16S rRNA基因序列片段，经测定其序列全长度为1430bp，经与GENEBANK中的16S rRNA片段比较分析，与6种不同科的微生物的同源率最高的为95%-95.47%。 采用温度梯度和病原浓度梯度回归感染实验方法，较为系统地研究了QDP的致病性。研究结果表明：QDP对栉孔扇贝有强烈的致病作用，高温（23℃以上）是其致病的必要条件，证实DQP是栉孔扇贝大规模死亡的病原体之一。

Classification of common oysters from North China

Oysters are commonly found on rocky shores along China's northern coast, although there is considerable confusion as to what species they are. To determine the taxonomic status of these oysters, we collected specimens from nine locations north of the Yangtze River and conducted genetic identification using DNA sequences. Fragments from three genes, mitochondrial 165 rRNA, mitochondria! cytochrome oxidase I (COI), and nuclear 285 rRNA, were sequenced in six oysters from each of the nine sites. Phylogenetic analysis of all three gene fragments clearly demonstrated that the small oysters commonly found on intertidal rocks in north China are Crassostrea gigas (Thunberg, 1793), not C. plicatula (the zhe oyster) as widely assumed. Their small size and irregular shell characteristics are reflections of the stressful intertidal environment they live in and not reliable characters for classification. Our study confirms that the oysters from Weifang, referred to as Jinjiang oysters or C. rivularis (Gould, 1861), are C. ariakensis (Wakiya, 1929). We found no evidence for the existence of C. talienwhanensis (Crosse, 1862) and other Crassostrea species in north China. Our study highlights the need for reclassifying oysters of China with molecular data.

轮-履复合可变形机器人的移动机构参数分析

针对非结构环境地形特点,结合轮式、履带式移动机构在复杂地形环境中的运动特点,设计了一种适应于非结构环境的轮-履互换、履带几何形状可变化的复合型移动机器人。该轮-履复合型移动机器人的机构系统由一个控制箱体和两个相同的轮-履复合移动单元模块组成。每个轮-履复合移动单元模块由运动轮和几何形状可变的履带轮构成,在一个驱动力的作用下根据不同的约束环境而表现为轮式和履带式两种运动特征,同时此单元模块具有轮-履互换、履带几何形状可变等特点。运用理论推导、参数优化等方法对轮-履复合移动单元模块的机构参数进行了分析综合,并经仿真分析验证了机构参数选择的合理性,为轮-履复合、履带可变形移动机器人的机构参数化设计提供了理论依据。

模块化变形机器人非同构构形表达与计数

提出一种新型链式可重构模块机器人平台,该机器人平台具有手动可重构和自动变形的特点,介绍一种三模块变形机器人样机。组成机器人的单个模块可以简化为由模块本体、连接臂和偏置关节组成。模块的数量可以根据实际工作的需要进行选择,模块间的连接具有规则连接和非规则连接两种方式;同时,由连接模块的偏置关节的运动,模块间的相对位置可以改变。由于模块连接方式的不同和模块间相对位置的变化,变形机器人具有多种非同构构形;为此,根据模块的物理结构和邻接关系提出了用构形矩阵来表达机器人结构的拓扑信息,并在仿真环境下进行等效描述;提出基于组合计数原理的递归算法,用于多模块变形机器人的非同构构形的计数,并根据构形矩阵的拓扑信息对构形进行评价。最后根据仿真结果给出了一种三模块变形机器人样机对称构形的设计示例,验证了算法的可行性。

AUV自主导航航位推算算法的分析研究

地球形状的不规则性,各种导航传感器本身的误差,以及仪器的安装偏差等,使得AUV(自治水下机器人)在进行远距离自主航行时,自主导航的精度大大下降。针对以上问题及实际工程需要,论文对AUV自主导航的航位推算算法做了进一步研究并加以改进,以提高其自主导航精度。最后,利用2004年中国科学院沈阳自动化所水下机器人研究中心进行AUV湖试所获得的数据,对文中提出的算法进行了验证。结果表明,AUV的自主导航精度得到大大提高,可以用于修正原来的自主导航算法。

碳酸盐岩缝洞型储层地球物理响应特征及预测方法研究

Proven by the petroleum exploration activities, the karsts-fissure reservoir in carbonate rocks is significant to find out the large scale oil & gas field. They are made up of the four reservoir types: karsts-cave, karsts-crack, crack-cave and fracture-pore-cave. Each reservoir space and each reservoir bed has different features of reservoir heterogeneity and small scale of pore-crack-cave. The fracture-cave reservoir in carbonate rocks is characteristic by multi-types and long oiliness well. The reservoir shape is controlled by the irregular pore-crack-cave. The development level of fracture and karst-cave is the key element of hydrocarbon enriching, high productivity and stable production. However, most of Carbonate formation are buried deeply and the signal-ration-noise of seismic reflection are very low. It is reason why the fracture-cave reservoir are difficult to be predicted effectively. In terms of surveyed and studied lots of the former research outcome, The author applied the methods of synthetical reservoir geophysical prediction from two ways including macrosopic and microcomic technics in terms of the reservoir-cap condition, geophysics and geology feature and difficulty of prediction in carbonate rocks. It is guiden by the new ideas of stratigraphy, sedimentology, sedimentography, reservoir geology and karst geology. The geophysics technology is key technics. In aspects of macroscopic studies, starting off the three efficiencies of controlling the reservoir distribution including sedimental facies, karst and fracture, by means of comprehensive utilization of geology, geophysics, boring well and well log, the study of reservoir features and karst inside story are developed in terms of data of individual well and multiple well. Through establishing the carbonate deposition model, karstic model and fracture model, the macro-distribution laws of carbonatite are carried out by the study of coherence analysis, seismic reflection feature analysis and palaeotectonics analysis. In aspects of microcosmic studies, starting off analysis in reservoir geophysical response feature of fracture and karst-cave model according to guidance of the macroscopic geological model in carbonate reservoir, the methods of the carbonate reservoir prediction are developed by comprehensively utilization of seismic multi-attribution intersection analysis, seismic inversion restricted by log, seismic discontinuity analysis, seimic spectrum attenuation gradient, moniliform reflection feature analysis and multiparameter karst reservoir appraisement.Through application of carbonate reservoir synthetical geophysics prediction, the author r successfully develops the beneficial reservoir distribution province in Ordovician of Katake block 1in middle Tarim basin. The fracture-cave reservoir distributions are delineated. The prospect direction and favorable aims are demonstrated. There are a set of carbonate reservoir prediction methods in middle Tarim basin. It is the favorable basic technique in predicting reservoir of the Ordovician carbonate in middle Tarim. Proven by exploration drilling, the favorable region of moniliform reflection fracture and pore-cave and cave-fracture in lower-middle Ordovician are coincidence with the region of hydrocarbon show. It’s indicated that the reservoir prediction methods described in the study of Ordovician carbonate formation are feasible practicably.

重庆老龙洞二叠系-三叠系界线地层微生物岩研究

The largest mass extinction in the Phanerozoic happened at the end of the Permian. The microbialites formed in the extreme environments after the mass extinction has become a hotspot for geologists and paleontologists throughout the world. The dendroid microbialites that were described for the first time in 1999 from the Permian-Triassic boundary section at Laolongdong, Chongqing, have been studied by many geologists from China and overseas. Two important viewpoints about their origin have been proposed. Some researchers believed that they resemble Quaternary travertine shrubs in form, and may belong to microbialites. Some other researchers proposed that the dendroid structure is composed of clots formed by coccoidal cynaobacteria, and is microbialite. Our detailed survey on the section reveals that: (1) there is an interval of speckled “microbialite” in the section, and it underlies the dendroid “microbialite”, (2) the dendroid “microbialite” does not always have dendroid appearance; they are dendroid only in very local places; they are not dendroid in most places; for this reason, they are not comparable to recent tufa; (3) the volume of the dendroid structure greatly increases toward the top of the dendroid microbialite interval: accounting to 70% of the whole rock in the top part. This distribution pattern implies that the formation of this structure may be related to downward migration of the diagenetic fluid. Examination of thin sections reveals that the dendroid structure or point-like structure in the “microbialite” look as lighter areas in the thin sections and are composed of large blocky clear calcites containing scattered yellow dirty small calcite rhombi and irregular “points” of relict lime mudstone or wackestone or packstone. Their formation is by any one of the following two processes: (1) dissolution → filling of large blocky calcite; (2) dolomitization → dedolomitization → dissolution by meteoric fresh water → filling by large blocky calcites. It has been found that there are at least two sea-level falls during the P-T transition. As the sea level fall, the carbonate deposits came into supratidal environment, and suffered dolomitization caused by evaporative fluid or mixing water of sea water and meteoric water. Since the fluid migrated downward from the top of the deposits and in random pathway, the dolomitization formed dendroid or speckled dolomitic areas. As the deposits came into subaerial environments, the meteoric fresh water migrated along the dendroid or speckled dolomitic area with higher porosity, and dissolution happened, which caused the rock became spongy or alveolate. In later time, after the strata came into phreatic zone, large clear blocky calcites grew in and filled the pores in the spongy areas. The dendroid and speckled structure were formed in this way, rather than composed of clots formed by coccoid cyanobecteria. The microbial fossils in Laolongdong section include two types. The first is the tube-like cyanobecteria in middle Bed 3, which are generally less than 1 mm in length, taper toward one end, and are internally filled by microspars. They are straight or sinuous, with micritic wall 0.005~0.01 mm thick. Since this kind of microbial fossils are abundant in middle Bed 3, this rock belongs to microbialite. The second type occurs in Bed 5 and lower and middle Bed 6. They are irregular globular in shape, generally 0.2 ~ 0.5 mm in size, with several outward progresses, and internally filled by one layer of needle-like calcite cements on the wall and the large blocky calcite in the inner space. According to their shape and preservation way, it is inferred that this kind of fossils were formed from some kind of bacterial colony. The bacterial colony may be cuticle in composition, since it has some hardness as it is indicated by its resistance to deposit loading. These organisms discomposed during diagenetic time, and formed good porosity. In later diagenetic time, these pores were firstly cemented by needle-like calcites and later filled by large blocky calcites. So, the bacterial colony promoted the formation of dendroid and speckled structures. However, they did not always form such structures. On the other hand, even though no bacterial colony or other microbes or any kind of fossils were present, dendroid or speckled structures can form. Bed 4 of Laolongdong section contains abundant gastropods but no microbial fossils, and is not microbialite, even though it is speckled. The top of Bed 6 is dendroid, but contain no microbial fossils, and is not micrbialite.

金川镍矿岩体移动与开采稳定性研究

The two major issues in mining industry are work safety and protection of ground environment when carrying on underground mining activities. Cut-and-fill mining method is increasingly applied in China owing to its advantages of controlling ground pressure and protecting the ground environment effectively. However, some cut-and-fill mines such as Jinchuan nickel mine which has big ore body, broken rock mass and high geostress have unique characteristics on the law of ground pressure and rock mass movement that distinguish from other mining methods. There are still many problems unknown and it is necessary for the further analysis. In this dissertation, vast field survey, geology trenching and relative data analysis are carried out. The distribution of ground fissures and the correlation of the fissures with the location of underground ore body is presented. Using of monitoring data by three-dimension fissure meter and GPS in Jinchuan Deposit Ⅱ, the rule of the surface deformation and the reason of ground fissures generation are analyzed. It is shown that the stress redistribution in surrounding rocks resulting from mining, the existence of the void space underground and the influence of on-going mining activities are three main reasons for the occurrence of ground fissures. Based on actual section planes of No.1 ore body, a large-scale 3D model is established. By this model, the complete process of excavation and filling is simulated and the law of rock mass movement and stability caused by Cut-and-fill Mining is studied. According to simulation results, it is concluded that the deformation of ground surface is still going on developing; the region of subsidence on the ground surface is similar with a circle; the area on the hanging wall side is larger than one on the lower wall side; the contour plots show the centre of subsidence lay on the hanging wall side and the position is near the ore body boundary of 1150m and 1250m where ore body is the thickest. Along strike-line of Jinchuan Deposit Ⅱ, the deformation at the middle of filling body is larger than that in the two sides. Because of the irregular ore body, stress concentrates at the boundary of ore body. With the process of excavation and filling, the high stress release and the stress focus disappear on the hanging wall side. The cut-and-fill mechanism is studied based on monitoring data and numerical simulation. The functions of filling body are discussed. In this dissertation, it is concluded that the stress of filling body is just 2MPa, but the stress of surrounding rock mass is 20MPa. We study the surface movement influenced by the elastic modulus of backfill. The minimal value of the elastic modulus of backfill which can guarantee the safety production of cut-and-fill mine is obtained. Finally, based on the real survey results of the horizontal ore layer and numerical simulation, it is indicated that the horizontal ore layer has destroyed. Key words: cut-and-filling mining, 3D numerical simulation, field monitoring, rock mass movement, cut-and-filling mechanism, the elastic modulus of backfill, the horizontal ore layer

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

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.

积分偏移的计算几何与并行实现

The Second Round of Oil & Gas Exploration needs more precision imaging method, velocity vs. depth model and geometry description on Complicated Geological Mass. Prestack time migration on inhomogeneous media was the technical basic of velocity analysis, prestack time migration on Rugged surface, angle gather and multi-domain noise suppression. In order to realize this technique, several critical technical problems need to be solved, such as parallel computation, velocity algorithm on ununiform grid and visualization. The key problem is organic combination theories of migration and computational geometry. Based on technical problems of 3-D prestack time migration existing in inhomogeneous media and requirements from nonuniform grid, parallel process and visualization, the thesis was studied systematically on three aspects: Infrastructure of velocity varies laterally Green function traveltime computation on ununiform grid, parallel computational of kirchhoff integral migration and 3D visualization, by combining integral migration theory and Computational Geometry. The results will provide powerful technical support to the implement of prestack time migration and convenient compute infrastructure of wave number domain simulation in inhomogeneous media. The main results were obtained as follows: 1. Symbol of one way wave Lie algebra integral, phase and green function traveltime expressions were analyzed, and simple 2-D expression of Lie algebra integral symbol phase and green function traveltime in time domain were given in inhomogeneous media by using pseudo-differential operators’ exponential map and Lie group algorithm preserving geometry structure. Infrastructure calculation of five parts, including derivative, commutating operator, Lie algebra root tree, exponential map root tree and traveltime coefficients , was brought forward when calculating asymmetry traveltime equation containing lateral differential in 3-D by this method. 2. By studying the infrastructure calculation of asymmetry traveltime in 3-D based on lateral velocity differential and combining computational geometry, a method to build velocity library and interpolate on velocity library using triangulate was obtained, which fit traveltime calculate requirements of parallel time migration and velocity estimate. 3. Combining velocity library triangulate and computational geometry, a structure which was convenient to calculate differential in horizontal, commutating operator and integral in vertical was built. Furthermore, recursive algorithm, for calculating architecture on lie algebra integral and exponential map root tree (Magnus in Math), was build and asymmetry traveltime based on lateral differential algorithm was also realized. 4. Based on graph theory and computational geometry, a minimum cycle method to decompose area into polygon blocks, which can be used as topological representation of migration result was proposed, which provided a practical method to block representation and research to migration interpretation results. 5. Based on MPI library, a process of bringing parallel migration algorithm at arbitrary sequence traces into practical was realized by using asymmetry traveltime based on lateral differential calculation and Kirchhoff integral method. 6. Visualization of geological data and seismic data were studied by the tools of OpenGL and Open Inventor, based on computational geometry theory, and a 3D visualize system on seismic imaging data was designed.

消除采集脚印的处理方法研究

Seismic technique is in the leading position for discovering oil and gas trap and searching for reserves throughout the course of oil and gas exploration. It needs high quality of seismic processed data, not only required exact spatial position, but also the true information of amplitude and AVO attribute and velocity. Acquisition footprint has an impact on highly precision and best quality of imaging and analysis of AVO attribute and velocity. Acquisition footprint is a new conception of describing seismic noise in 3-D exploration. It is not easy to understand the acquisition footprint. This paper begins with forward modeling seismic data from the simple sound wave model, then processes it and discusses the cause for producing the acquisition footprint. It agreed that the recording geometry is the main cause which leads to the distribution asymmetry of coverage and offset and azimuth in different grid cells. It summarizes the characters and description methods and analysis acquisition footprint’s influence on data geology interpretation and the analysis of seismic attribute and velocity. The data reconstruct based on Fourier transform is the main method at present for non uniform data interpolation and extrapolate, but this method always is an inverse problem with bad condition. Tikhonov regularization strategy which includes a priori information on class of solution in search can reduce the computation difficulty duo to discrete kernel condition disadvantage and scarcity of the number of observations. The method is quiet statistical, which does not require the selection of regularization parameter; and hence it has appropriate inversion coefficient. The result of programming and tentat-ive calculation verifies the acquisition footprint can be removed through prestack data reconstruct. This paper applies migration to the processing method of removing the acquisition footprint. The fundamental principle and algorithms are surveyed, seismic traces are weighted according to the area which occupied by seismic trace in different source-receiver distances. Adopting grid method in stead of accounting the area of Voroni map can reduce difficulty of calculation the weight. The result of processing the model data and actual seismic demonstrate, incorporating a weighting scheme based on the relative area that is associated with each input trace with respect to its neighbors acts to minimize the artifacts caused by irregular acquisition geometry.

地震波模拟：非规则网格离散的数值解法、吸收边界及其实际应用

This dissertation presents a series of irregular-grid based numerical technique for modeling seismic wave propagation in heterogeneous media. The study involves the generation of the irregular numerical mesh corresponding to the irregular grid scheme, the discretized version of motion equations under the unstructured mesh, and irregular-grid absorbing boundary conditions. The resulting numerical technique has been used in generating the synthetic data sets on the realistic complex geologic models that can examine the migration schemes. The motion equation discretization and modeling are based on Grid Method. The key idea is to use the integral equilibrium principle to replace the operator at each grid in Finite Difference scheme and variational formulation in Finite Element Method. The irregular grids of complex geologic model is generated by the Paving Method, which allow varying grid spacing according to meshing constraints. The grids have great quality at domain boundaries and contain equal quantities of nodes at interfaces, which avoids the interpolation of parameters and variables. The irregular grid absorbing boundary conditions is developed by extending the Perfectly Matched Layer method to the rotated local coordinates. The splitted PML equations of the first-order system is derived by using integral equilibrium principle. The proposed scheme can build PML boundary of arbitrary geometry in the computational domain, avoiding the special treatment at corners in a standard PML method and saving considerable memory and computation cost. The numerical implementation demonstrates the desired qualities of irregular grid based modeling technique. In particular, (1) smaller memory requirements and computational time are needed by changing the grid spacing according to local velocity; (2) Arbitrary surfaces and interface topographies are described accurately, thus removing the artificial reflection resulting from the stair approximation of the curved or dipping interfaces; (3) computational domain is significantly reduced by flexibly building the curved artificial boundaries using the irregular-grid absorbing boundary conditions. The proposed irregular grid approach is apply to reverse time migration as the extrapolation algorithm. It can discretize the smoothed velocity model by irregular grid of variable scale, which contributes to reduce the computation cost. The topography. It can also handle data set of arbitrary topography and no field correction is needed.

土石混合体细观结构力学及其边坡稳定性研究

Soil-rock mixture (S-RM) refers to one extremely uneven loose rock and soil materials system with certain stone content. Its formation has started since Quaternary and it is composed of block stone, fine grained soil and pore with certain project scale and high strength. S-RM has extensive distribution in nature, especially in southwest China where the geotectonic background is complicated, the fracture activity is developed and the geomorphological characteristics of high mountain and steep gorge area are protuberant. This kind of complicated geologic body has developed wider in these areas. S-RM has obvious difference with the general soil or rock (rock mass) in physical and mechanical properties because its two components-“soil” and “rock-block” has extreme differences in physical and mechanical properties. The proposition of S-RM and its deep research are needed in the modern engineering construction. It is also the necessity in the modern development of rock and soil mechanics. The dissertation starts from the meso-structural characteristics of soil-rock and takes a systematic research on its meso-structural mechanics, deformation and failure mechanism and the stability of S-RM slope. In summary, it achieves the following innovative results and conclusions. There are various views on the conception of S-RM and its classification system. Based on the large number of field tests, the dissertation makes the conception and classification of S-RM more systematic. It systematically proposed the conception of meso-structural mechanics of S-RM. Thus the dissertation has laid a foundation for its deep study. With the fast development of the computer technology and digital image processing theory, digital image processing technology has been successfully applied in many fields and provided reliable technology support for the quantitative description of the structural characteristics of S-RM. Based on the digital image processing technology, the dissertation systematically proposes and developed the quantitative analysis method and quantitative index for the meso-structure of S-RM. The results indicate that the meso-structure such as its internal soil-rock granularity composition, the soil-rock shape and the orientability has obvious self-organization in the macro statistical level. The dissertation makes a systematic research on the physical mechanical properties, deformation and failure mechanism of S-RM based on large field test. It proposes the field test for the underwater S-RM and deduces the 3D data analysis method of in-situ horizontal push-shear test. The result indicates that S-RM has significant phenomenon of shear dilatancy in the shearing process, and its dilatancy will be more obvious with the increased proportion of rock or the decreased confining pressure. The proportion of rock has great effect on the strength of S-RM and rock-block, especially the spatial position of particles with comparatively big size has great effect on the shape and spatial position of the sample shear zone. The dissertation makes some improvements in the single ring infiltration test equipment and its application on the permeability of S-RM. The results indicate that the increasing of rock-block would make it more difficult for the soil to fill in the vacuity between the rock-block and the proportion would increase which would result in the increased permeability coefficient. The dissertation builds the real meso-structural model of S-RM based on the digital image processing technology. By using geometric reconstruction technology, it transfers the structural mode represented by Binary image into CAD format, which makes it possible to introduce the present finite element analysis software to take research on numerical experimental investigation. It systematically realizes leaping research from the image，geometric mode, to meso-structural mechanics numerical experiment. By using this method, the dissertation takes large scale numerical direct-shear test on the section of S-RM. From the mesoscopic perspective, it reveals three extended modes about the shear failure plane of S-RM. Based on the real meso-structural model and by using the numerical simulation test, the character and mechanics of seepage failure of S-RM are studied. At the same time, it builds the real structural mode of the slope based on the analysis about the slope crosssection of S-RM. By using the strength reduction method, it takes the research on the stability of S-RM and gets great achievements. The three dimensional geometric reconstruction technology of rock block is proposed, which provides technical support for the reconstruction of the 3D meso-structural model of S-RM. For the first time, the dissertation builds the stochastic structure model of two-dimensional and three-dimensional polygons or polyhedron based on the stochastic simulation technique of monte carlo method. It breaks the traditional research which restricted to the random generation method of regular polygon and develops the relevant software system (R-SRM2D/3D) which has great effect on meso-structural mechanics of S-RM. Based on the R-SRM software system which randomly generates the meso-structural mode of S-RM according to the different meso-structural characteristics, the dissertation takes a series of research on numerical test of dual axis and real three-axis, systematically analyses the meso destroy system, the effects of meso-structural characteristics such as on the stone content, size composition and block directionality on the macro mechanical behavior and macro-permeability. Then it proposes the expression of the upper and lower limit for the macro-permeability coefficient of the inhomogeneous geomaterials, such as S-RM. By using the strength reduction FEM, the dissertation takes the research on the stability of the slope structural mode of the randomly formed S-RM. The results indicate that generally, the stability coefficient of S-RM slope increases with the increasing of stone content; on the condition of the same stone content, the stability coefficient of slope will be different with different size composition and the space position of large block at the internal slop has great effect on the stability. It suggests that meso-structural characteristics, especially the space position of large block should be considered when analyzing the stability of this kind of slope and strengthening design. Taking Xiazanri S-RM slope as an example, the dissertation proposes the fine modeling of complicated geologic body based on reverse engineering and the generation method of FLAC3D mode. It resolves the bottleneck problem about building the fine structural mode of three-dimensional geological body. By using FLAC3D, the dissertation takes research on the seepage field and the displacement field of Xiazanri S-RM slope in the process of reservoir water level rising and decreasing. By using strength reduction method, it analyses the three-dimension stability in the process of reservoir water level rising and decreasing. The results indicate that the slope stability firstly show downward trend in the process of reservoir water level rising and then rebound to increase; the sudden drawdown of reservoir water level has great effect on the slope stability and this effect will increase with the sudden drawdown amplitude rising. Based on the result of the rock block size analysis of S-RM, and using R-SRM2D the stochastic structure model of Xiazanri S-RM slope is built. By using strength reduction method, the stability of the stochastic structure model is analysis, the results shows that the stability factor increases significantly after considering the block.

基于广义Lippmann-Schwinger方程的地震波场模拟算法与应用

In the last several decades, due to the fast development of computer, numerical simulation has been an indispensable tool in scientific research. Numerical simulation methods which based on partial difference operators such as Finite Difference Method (FDM) and Finite Element Method (FEM) have been widely used. However, in the realm of seismology and seismic prospecting, one usually meets with geological models which have piece-wise heterogeneous structures as well as volume heterogeneities between layers, the continuity of displacement and stress across the irregular layers and seismic wave scattering induced by the perturbation of the volume usually bring in error when using conventional methods based on difference operators. The method discussed in this paper is based on elastic theory and integral theory. Seismic wave equation in the frequency domain is transformed into a generalized Lippmann-Schwinger equation, in which the seismic wavefield contributed by the background is expressed by the boundary integral equation and the scattering by the volume heterogeneities is considered. Boundary element-volume integral method based on this equation has advantages of Boundary Element Method (BEM), such as reducing one dimension of the model, explicit use the displacement and stress continuity across irregular interfaces, high precision, satisfying the boundary at infinite, etc. Also, this method could accurately simulate the seismic scattering by the volume heterogeneities. In this paper, the concrete Lippmann-Schwinger equation is specifically given according to the real geological models. Also, the complete coefficients of the non-smooth point for the integral equation are introduced. Because Boundary Element-Volume integral equation method uses fundamental solutions which are singular when the source point and the field are very close，both in the two dimensional and the three dimensional case, the treatment of the singular kernel affects the precision of this method. The method based on integral transform and integration by parts could treat the points on the boundary and inside the domain. It could transform the singular integral into an analytical one both in two dimensional and in three dimensional cases and thus it could eliminate the singularity. In order to analyze the elastic seismic wave scattering due to regional irregular topographies, the analytical solution for problems of this type is discussed and the analytical solution of P waves by multiple canyons is given. For the boundary reflection, the method used here is infinite boundary element absorbing boundary developed by a pervious researcher. The comparison between the analytical solutions and concrete numerical examples validate the efficiency of this method. We thoroughly discussed the sampling frequency in elastic wave simulation and find that, for a general case, three elements per wavelength is sufficient, however, when the problem is too complex, more elements per wavelength are necessary. Also, the seismic response in the frequency domain of the canyons with different types of random heterogeneities is illustrated. We analyzed the model of the random media, the horizontal and vertical correlation length, the standard deviation, and the dimensionless frequency how to affect the seismic wave amplification on the ground, and thus provide a basis for the choice of the parameter of random media during numerical simulation.

湘西地区上震旦统上部－下寒武统底部黑色岩系沉积学及地球化学

Widespread black chert-shales occur in the Ediacaran-Cambrian(E-C) boundary successions along the flank of Yangtze Platform, South China, remarkable changes in sedimentology, geochemistry and biology were recorded. Although extensive studies were carried out upon this boundary succession, the origin of black chert-shales still remain controversial. This paper focuses on the E-C black chert-shales in western Hunan, South China, upon which detailed depositional and geochemical changes are documented, accordingly a depositional model for black chert-shales is proposed. Stratigraphic anatomy across the depositional strike demonstrates that the shallow-water Dengying dolostone along the platform margin sharply pass basinward into the Liuchapo chert successions, which indicate syndepositional extensional faulting at depth could have occurred along the platform margin. The deep-water Niutitang phosphorite-rich black shales are either underlain by the Dengying dolostones on the platform margin toward platform interior or directly by the Liuchaopo chert successions farther basinwards. By detailed investigation, silica chimneys are firsly identified approximately in the chert along platform margin; two types of silica chimneys, including mounded and splayed/funnelized chert(generally brecciated) bodies are further sorted out. The mounded chert are exitbited by domed or hummocky surfaces on the top and irregular spongy to digitiform internal fabrics; within the silica mounds, abundant original vesicles/voids and/or channels were mostly plugged by initial chalcedony, quartze crystals with minor dolomite and bladed barite crystals. Splayed/funnelized brecciated chert “intrusion” cross-cut the uppermost dolostones capping to the horizon underneath, and are directly overlain by the Niutitang phosphorite-rich black shales. Their similarities to the silica chimneys reported from the oceanic spreading centres suggest a similar origin responsible for these unique silica bodies which is also supported by the microthermonmetric data and element geochemistry. High P, Ba, Fe contents and positive correlation between Fe and TOC concentrations in the Niutitang black shales indicate a high palaeo-productivity in the Early Cambrian ocean. The low Th/U and the high V/Cr, V/Sc, V/(V+Ni) ratios in the black shales suggest an anoxic water condition during this interval. Furthermore, Positive Eu anomalies and high Ba contents in the sediments also imply a hydrothermal influence on the formation of Niutitang black shales. To better constrain the placement of deep-water successions straddling the E-C boundary and the timing of hydrothermal silica chimneys, sensitive high-resoluton ion microprobe(SHRIMP) U-Pb dating of zircon grains from tuffs within the chert succession of Liuchapo Formation at Ganziping was conducted and yields a weighted-mean 206Pb/238Pb age at 536.6±5.5Ma, younger than E-C boundary age(542.0±0.3Ma). This age combined with carbon isotopic data is then proposed to correspond to the U-Pb age of zircons(538.2±1.5Ma) from the Zhongyicun member of Meishucun Formation at Meishucun in eastern Yunna, thus, the E-C boundary in Gazngziping was placed between the Dengying formations and Liuchapo formatioms. therefore, the silica chimneys took place at the beginning of the Cambrian period. The temporal coincidence of silica chimneys and negative excursions of δ13C and δ34Spy pairs suggest hydrothermal activities were likely responsible for the isotopic changes. Under such a circumstance, vast amounts of greenhouse gases(CO2, CH4, H2S), with highly 13C-depleted carbon and 34S-depleted sulfur would be released into the ocean and atmosphere. A positive shift in δ34Scas and Δ34S values from the late Ediacaran to the Early Cambrian could be a reflection of enhanced bacterial sulfate reduction(BSR), strengthened by the intensified oceanic anoxia stimulated by hydrothermal activities. Based on the analyses of sedimentology and geochemistry, a model- “oceanic anoxia induced by hydrothermal–volcanic activies” was proposed to responsible for the formation of black chert-shales during this E-C transition. Under this case, hydrothermal-volcanic activies could release large large amount of greenhouse into atmosphere and metal micronutrients into the ocean, which may lead to global warming, stratified ocean, thereby a high palaeoproductivity; on the other hand, the massive releasing of reduced hydrothermal fluids with abundant H2S, could have in turn enhanced the ocean anoxia. All of these were favourable the for preservation of organic matter, and subsequent extensive deposition of black silica-shales.