Donut-Shaped Fingerprint In Homologous Polypeptide Relationships-A Topological Feature Related To Pathogenic Structural Changes In Conformational Disease

Features of homologous relationship of proteins can provide us a general picture of protein universe, assist protein design and analysis, and further our comprehension of the evolution of organisms. Here we carried Out a Study of the evolution Of protein molecules by investigating homologous relationships among residue segments. The motive was to identify detailed topological features of homologous relationships for short residue segments in the whole protein universe. Based on the data of a large number of non-redundant Proteins, the universe of non-membrane polypeptide was analyzed by considering both residue mutations and structural conservation. By connecting homologous segments with edges, we obtained a homologous relationship network of the whole universe of short residue segments, which we named the graph of polypeptide relationships (GPR). Since the network is extremely complicated for topological transitions, to obtain an in-depth understanding, only subgraphs composed of vital nodes of the GPR were analyzed. Such analysis of vital subgraphs of the GPR revealed a donut-shaped fingerprint. Utilization of this topological feature revealed the switch sites (where the beginning of exposure Of previously hidden "hot spots" of fibril-forming happens, in consequence a further opportunity for protein aggregation is Provided; 188-202) of the conformational conversion of the normal alpha-helix-rich prion protein PrPC to the beta-sheet-rich PrPSc that is thought to be responsible for a group of fatal neurodegenerative diseases, transmissible spongiform encephalopathies. Efforts in analyzing other proteins related to various conformational diseases are also introduced. (C) 2009 Elsevier Ltd. All rights reserved.

Hematological and plasma biochemical responses of crucian carp (Carassius auratus) to intraperitoneal injection of extracted microcystins with the possible mechanisms of anemia

Alterations in hematological indices such as decreases in blood cell counts (RBC), hematocrit (Ht) and hemoglobin (Hb) concentrations are key symptoms of anemia. However, few experiments were conducted to examine changes in hematological indices of fish exposed to microcystins that are believed to be fatal to circulatory systems of vertebrates. An acute toxicological experiment was designed to study hematological changes of crucian carp injected intraperitoneally (i.p.) with extracted microcystins at two doses, 50 and 200 mu g MC-LReqkg(-1) body weight. After being i.p. injected with microcystins, the fish exhibited behavioral abnormity. There were significant decreases in RBC in the high-dose group, and in Ht and Hb concentrations in both dose groups, while erythrocte sedimentation rate (ESR) significantly increased, indicating the appearance of normocytic anemia. There were no prominent changes in the three red cell indices, mean corpuscular volume (MCV), mean corpuscular hemoglobin (MCH,), and mean corpuscular hemoglobin concentration (MCHC). Increases in blood urea nitrogen (BUN) and creatinine (CR) in both dose groups suggest the occurrence of kidney impairment. Alteration in blood indices was reversible at the low dose group. Conclusively, anemia induced by kidney impairment was a key factor to cause abnormity of swimming behaviors and high mortality of crucian carp. (c) 2007 Elsevier Ltd. All rights reserved.

Introduction of the human lactoferrin gene into grass carp (Ctenopharyngodon idellus) to increase resistance against GCH virus

Haemorrhage can be an epidemic and fatal condition in grass carp. It is known now that the Grass Carp Haemorrhage Virus (GCHV) triggers haemorrhage. Human lactoferrin (hLF) plays an important role in the non-specific immune system, making some organisms more resistant to some viruses. Sperm of grass carp was mixed with linearized pCAhLFc, which is a DNA construct containing an hLF cDNA and the promoter of common carp beta-actin gene, and then electroporated. Then, mature eggs were fertilized in vitro with the treated sperm cells. The fry were sampled and analyzed by polymerase chain reaction (PCR). Results indicated that the foreign gene had been transferred successfully into the cells of some fry. Under optimal electroporation conditions, the efficiency of gene transfer was as high as 46.8%. About 35.7% of treated 5-month-old grass carp contained foreign genes. Most transgenic fry demonstrated significant delays in onset of symptoms of haemerrhage after injection of GCHV, suggesting a significant positive relationship between hLF cDNA and levels of disease resistance (P < 0.01). Results suggest that transgenic grass carp could be bred for increased resistance to haemorrhage. (C) 2002 Elsevier Science B.V. All rights reserved.

活性炭纤维(ACF)吸附法处理河流突发污染事故模拟研究

环境突发污染事故给人民生活、经济发展和生态环境造成重大影响，研究污染物泄漏造成河流突发污染事故的应急处理方法十分必要。本论文选取苯酚、苯胺和亚甲基蓝等典型污染物为实验对象，采用吸附容量大、密度与水接近的活性炭纤维(ACF)为吸附剂。在自制的河流模型中，研究了ACF以苯酚、苯胺和亚甲基蓝为典型污染物的吸附过程，考察了吸附剂投加量、污染物浓度、吸附剂比表面积、吸附剂投加方式、水流速度与水质等对吸附速率与吸附效果的影响。实验结果表明，ACF能以较快的速率吸附苯酚、苯胺和亚甲基蓝，吸附率都在95％以上； ACF投加量是影响吸附速率最重要的因素，当一次性投加ACF质量之比为 1:2:4时，吸附速率常数之比近似为1:2:4；污染物浓度对吸附速率的影响显著，浓度较低时吸附速率较高。苯酚初始浓度为7mg·L-1时，经过86分钟的吸附，处理后的浓度可以达到地表水Ⅴ类水中挥发酚的限值要求(0.1mg·L-1)；在吸附11分钟左右追加适量的ACF，能够明显提高吸附速率；河水流速和河流中的天然有机物、浊度、河水硬度对ACF吸附都不产生显著影响，这说明ACF作为河流突发污染事故应急处理的吸附剂，有广泛的适应性。在实际河水中，ACF对苯酚的吸附过程与在模拟河水中相似，吸附效果显著。实验结果还表明，ACF对苯酚的吸附是放热反应，符合Freundlich模型和Langmuir模型。事故应急处理后，应该及时将吸附了污染物的ACF打捞上来，有利于进行后续处理。 Emergency environmental pollution accidents pose significant impacts on our living, economic development and ecological environment. The study on the approach of emergency control for the contingency caused by leakage of pollutants in rivers is very necessary. In the experiment, phenol, aniline and methylene blue were selected as representative pollutant and activated carbon fiber (ACF) was selected as adsorbent, which has strong adsorption capacity and similar density to water. In the self-made river model, the effects of ACF dosage, pollutant concentration, ACF surface area, ACF adding ways, water flow rate and water quality on adsorption courses were investigated. The experimental results showed that ACF could adsorb pollutant quickly and effectively. The ACF dosage was the most important factor that affected adsorption rate .When the ACF dosage rate was 1:2:4, the constants of adsorption rate was approximately 1:2:4. The effect of pollutant concentrations on the adsorption rate was notable. Faster adsorption rates were achieved at low pollutant concentrations. Phenol concentration reached the limits of volatile phenol in Category Ⅴ surface water (0.1 mg·L-1) after 86 minutes of adsorption with initial phenol concentration of 7 mg·L-1. After 11 minutes of adsorption, certain amount of ACF was added and the adsorption rate was improved significantly. River flow rate and water quality have little effect on the adsorption rate. The adsorption results obtained in actual river water were comparable with that in simulating river water. The results also showed that, ACF on the absorption of phenol is exothermic reaction, witch matched with the Freundlich model and the Langmuir model. After emergency treatment, the ACF absorbed pollutants should be promptly salvaged for follow-up treatment.

功能菌强化应急处理苯胺突发污染事故研究

近年来各种环境污染事故频发，据统计仅2001~2003年间，发生的各类环境污染事故就高达5606次，其中水污染事故3235次，占全部的57.7%。这些事故不仅给人民生命财产造成巨大损失，也给生态环境造成严重的破坏。因此开发安全高效的应急处理技术迫在眉睫。本研究以筛选高效苯胺降解菌为基础，通过对高效菌降解性能的研究指导将高效菌作为功能郡主投加到已有生物处理系统强化应急处理苯胺突发污染事故废液，取得了良好的效果。 苯胺高效降解菌AN-P1为红球菌（Rhodococcus sp.），其通过间位途径降解苯胺，AN-P1利用苯胺生长和降解的最佳pH为6，最适浓度为2000 mg/L，最适温度为30 ℃，最佳接种量为0.3‰。AN-P1降解含500 mg/L、1000 mg/L、2000 mg/L苯胺的培养物分别经过28 h、24 h、32 h降解，出水苯胺含量能达到《污水综合排放标准》（GB8978-1996）一级标准。但由于苯胺降解过程中释放了大量氨氮，出水氨氮仍较高未能达标排放。而常规SBR系统应急处理效果较差，苯胺和COD去除率均低于10%，出水未能达标排放。活性碳吸附后的回收和后续处理也会带来操作不变和二次污染问题，且处理后出水往往难于达标排放，尚需进行进一步处理。 生物处理系统应急处理后恢复运行处理效果监测和PCR-DGGE图谱分析显示，用AN-P1菌强化应急处理系统后不仅能快速高效的去除苯胺，而且可以有效保障处理系统对污染物的净化性能，有效的保护系统中的功能微生物免受苯胺毒害。 研究结果表明，从实际处理效果、对原有生物系统性能保护及实际应用操作等多方面考虑，用AN-P1菌强化应急处理苯胺突发污染事故在技术上都是可行的。本研究为应急处理苯胺突然污染事故废液提供了新的方法。 Recent years, environment pollution accidents happened frequently, the data showed that there are 5606 accidents between 2001 and 2003, including 3235 water environment accidents, which is 57.7% of all. These accedents not only caused money lost and life lost but also caused serious damage to the ecologicl environment. So exploring highly-effective and secure methods to solve these accidents is an urgent mission. We screened a highly-effective aniline-degrading bacterium and did some researches on its ability to degrade aniline, in order to guide the emergency treatment of aniline containing wastewater that caused by sudden accident pollution with bioaugmentation. A highly-effective aniline-degrading bacterium AN-P1 was isolate and characterized as Rhodococcus sp. It degrades aniline through meta-cleavage pathway. The optimal pH and temperature for cell growth and aniline degradation were 6 and 30 ℃, respectively, and the opitimal concentration of aniline was 2000 mg/L, the optimal inoculation amount was 0.3‰.It took bacterium AN-P1 only 18 h, 24 h and 32 h, respectively, for the treatment of MSB containing 500 mg/L, 1000 mg/L, 2000 mg/L aniline to meet the first grade of national some of the NH4+-N which caused by aniline degradation. It took bacterium AN-P1 only 10 h, 20 h and 32 h, respectively, for the treatment of wastewater containing 500 mg/L, 1000 mg/L, 2000 mg/L aniline to meet the first grade of national integrated wastewater discharge standard. The bacterium AN-P1 can also remove some of the NH4+-N which caused by aniline degradation. It took bacterium AN-P1 only 10 h, 20 h and 32 h, respectively, for the treatment of wastewater containing 500 mg/L, 1000 mg/L, 2000 mg/L aniline to meet the first grade of national integrated wastewater discharge standard. By combing AN-P1 with regular SBR system, it took only 36 h for the emergency treatment of wastewater containing 2000 mg/L aniline under simulating engineering conditions to meet the discharge standard. While the NH4+-N of effluent can not meet the standard because of the high amount NH4+-N caused by aniline degradation. The regular SBR system was not good at aniline and COD removal. The removal efficiency of which are less than 10%. It cost 67.8 g activated carbon to absorbed 1000 mg aniline. It is inconvenient to transport and use it for the emergency treatment of aniline when the sudden pollution accident happened. Meanwhile, it was complex ad hard to recycle the activated carbon and treat the aniline wastewater get from activated carbon recycling too. Hard to meet the effluent standard was also a problem of activated carbon absorption method. According to the PCR-DGGE profile and removal efficiency of pollutants and COD when the systerm recover from emergency treatment, AN-P1 can efficiently protect the microbial community of regular activated sludge system against the aniline. It proved that combing AN-P1 with regular biological system is a feasible strategy for emergency treatment of aniline sudden pollution accident. The research offered a new way for emergency treatment of aniline sudden pollution accident.

Prematurely condensed chromosome fragments in human lymphocytes induced by high doses of high-linear-energy-transfer irradiation

This study provides a useful biodosimetry protocol for radiation accidents that involve high doses of heavy particle radiation. Human peripheral blood lymphocytes (PBLs) were irradiated in vitro with high doses (5–50 Gy) of charged heavy-ion particles (carbon ions, at an effective linear-energy-transfer (LET) of 34.6 keV/ m), and were then stimulated to obtain dividing cells. PBLs were treated with 100nMcalyculin A to force chromosomes to condense prematurely, and chromosome spreads were obtained and stained with Giemsa. The G2 prematurely condensed chromosome (G2-PCC) index and the number of G2-PCC including fragments (G2-PCC-Fs) per cell for each radiation dose point were scored. Dose-effect relationships were obtained by plotting the G2-PCC indices or G2-PCC-Fs numbers against radiation doses. The G2-PCC index was greater than 5% up to doses of 15 Gy; even after a 30Gy radiation dose, the index was 1 to 2%. At doses higher than 30 Gy, however, the G2-PCC indices were close to zero. The number of G2-PCC-Fs increased steeply for radiation doses up to 30 Gy at a rate of 1.07 Gy−1. At doses higher than 30 Gy, the numbers of G2-PCC-Fs could not be accurately indexed because of the limited numbers of cells for analysis. Therefore, the number of G2-PCC-Fs could be used to estimate radiation doses up to 30 Gy. In addition, a G2-PCC index close to zero could be used as an indicator for radiation doses greater than 40 Gy.

Biodosimetry estimate for high-LET irradiation

The purpose of this paper is to prepare for an easy and reliable biodosimeter protocol for radiation accidents involving high-linear energy transfer (LET) exposure. Human peripheral blood lymphocytes were irradiated using carbon ions (LET: 34.6 keV mu m(-1)), and the chromosome aberrations induced were analyzed using both a conventional colcemid block method and a calyculin A induced premature chromosome condensation (PCC) method. At a lower dose range (0-4 Gy), the measured dicentric (dics) and centric ring chromosomes (cRings) provided reasonable dose information. At higher doses (8 Gy), however, the frequency of dics and cRings was not suitable for dose estimation. Instead, we found that the number of Giemsa-stained drug-induced G2 prematurely condensed chromosomes (G2-PCC) can be used for dose estimation, since the total chromosome number (including fragments) was linearly correlated with radiation dose (r = 0.99). The ratio of the longest and the shortest chromosome length of the drug-induced G2-PCCs increased with radiation dose in a linear-quadratic manner (r = 0.96), which indicates that this ratio can also be used to estimate radiation doses. Obviously, it is easier to establish the dose response curve using the PCC technique than using the conventional metaphase chromosome method. It is assumed that combining the ratio of the longest and the shortest chromosome length with analysis of the total chromosome number might be a valuable tool for rapid and precise dose estimation for victims of radiation accidents.

三峡库区巴东组泥灰岩溶蚀作用与岩体力学性质弱化研究

Abstract The karsrt erosion engineering geology became a highlight problem in recent years, in particularly, the karst erosion of marlite of Badong formation made the rock mechanics weaken in Three Gorges Reservoir area, which reduces the safety of slope. During the immigrant construction, many high slopes have been formed, whose instabilities problems pose serious threats to the safety of the people and properties. The accidents of the slope failure take place now and then. By testing, it has been found that the karst erosion pattern and dissolution rate of marlite are not weaker than that of the pure limestone. Furthermore, owning to the weathering and unloading, the karst erosion of the marlite will reach certain depth of the slope, which is named infiltrated karst erosion. The karst erosion made the rock mass quality of slope or foundation worse in a large scale. The karst erosion geological disasters, taken place or not, has become the main restrictive factors to the social stability and economic development. Thus the karst erosion process and mechanism of marlite of Badong formation are studied as the main content of this dissertation. The weakening characteristic of rock mass mechanics parameters are studied along with the rock mass structure deformation and failure processes in the course of the karst erosion. At first, the conditions and influencing factors of the karst erosion are analyzed in the investigative region, on the basis of different karst erosion phenomenon of the marlite and different failure modes of slope. Then via indoor the karst erosion tests, it is analyzed that the karst erosion will change the rock mass composition and its structure. Through test, the different karst erosion phenomena between micro and macro have been observed, and the karst erosion mechanism of the marlite has been summarized. Damage theory is introduced to explain the feature of dissolution pore and the law of crack propagation in the marlite. By microscope and the references data, it can be concluded that the karst erosion process can be divided into rock minerals damage and rock structural damage. And the percent of karst erosion volume is named damage factor, which can be used to describe the quantify karst erosion degree of marlite. Through test, the rock mechanical properties in the different period of karst erosion are studied. Based on the damage mechanics theory and the test result, the relation between the karst erosion degree of marlite and weakening degree of mechanical properties is summarized. By numerical simulations, the karst erosive rock mass mechanics is verified. The conclusion is drawn as below: to the rock mass of marlite, the karst erosion damage made mechanics parameters variation, the deformation modulus, cohesion, and inter friction angle reduce as the negative exponent with the increasing of the karst erosion volume, however, the Poisson ratio increases as the positive exponent with the karst erosion volume increasing. It should be noticed that the deduced formulations are limited to the test data and certain conditions. It is suitable to the rock mass parametric weakening process after the karst erosion of marlite in Three Gorges Reservoir area. Based on the failure types of marlite slope in the field, the karst erosion and weathering process of rock mass are analyzed. And the evolution law of deformation and failure of the marlite mass is studied. The main failure feature of the marlite slope is the karst erosive structure subsidence mode in Three Gorges Reservoir area. The karst erosive structure subsidence mode is explained as follows: the rock mass undergoes the synthetic influence, such as weathering, unloading, corrosion, and so on, many pores and cavities have been formed in the rock mass interior, the rock mass quality is worsen and the rock mass structure is changed, and then the inherent structure of rock mass is collapsed under its gravity, therefore, the failure mode of compaction and subsidence take place. Finally, two examples are used to verify the rock mass parameters in Three Gorges Reservoir area, and the relationship between the marlite slope stability and the time of karst erosion is proposed.

基于精细结构描述及数值试验的断续节理岩体力学参数智能选取

Evaluating the mechanical properties of rock masses is the base of rock engineering design and construction. It has great influence on the safety and cost of rock project. The recognition is inevitable consequence of new engineering activities in rock, including high-rise building, super bridge, complex underground installations, hydraulic project and etc. During the constructions, lots of engineering accidents happened, which bring great damage to people. According to the investigation, many failures are due to choosing improper mechanical properties. ‘Can’t give the proper properties’ becomes one of big problems for theoretic analysis and numerical simulation. Selecting the properties reasonably and effectively is very significant for the planning, design and construction of rock engineering works. A multiple method based on site investigation, theoretic analysis, model test, numerical test and back analysis by artificial neural network is conducted to determine and optimize the mechanical properties for engineering design. The following outcomes are obtained: (1) Mapping of the rock mass structure Detailed geological investigation is the soul of the fine structure description. Based on statistical window，geological sketch and digital photography，a new method for rock mass fine structure in-situ mapping is developed. It has already been taken into practice and received good comments in Baihetan Hydropower Station. (2) Theoretic analysis of rock mass containing intermittent joints The shear strength mechanisms of joint and rock bridge are analyzed respectively. And the multiple modes of failure on different stress condition are summarized and supplied. Then, through introducing deformation compatibility equation in normal direction, the direct shear strength formulation and compression shear strength formulation for coplanar intermittent joints, as well as compression shear strength formulation for ladderlike intermittent joints are deducted respectively. In order to apply the deducted formulation conveniently in the real projects, a relationship between these formulations and Mohr-Coulomb hypothesis is built up. (3) Model test of rock mass containing intermittent joints Model tests are adopted to study the mechanical mechanism of joints to rock masses. The failure modes of rock mass containing intermittent joints are summarized from the model test. Six typical failure modes are found in the test, and brittle failures are the main failure mode. The evolvement processes of shear stress, shear displacement, normal stress and normal displacement are monitored by using rigid servo test machine. And the deformation and failure character during the loading process is analyzed. According to the model test, the failure modes quite depend on the joint distribution, connectivity and stress states. According to the contrastive analysis of complete stress strain curve, different failure developing stages are found in the intact rock, across jointed rock mass and intermittent jointed rock mass. There are four typical stages in the stress strain curve of intact rock, namely shear contraction stage, linear elastic stage, failure stage and residual strength stage. There are three typical stages in the across jointed rock mass, namely linear elastic stage, transition zone and sliding failure stage. Correspondingly, five typical stages are found in the intermittent jointed rock mass, namely linear elastic stage, sliding of joint, steady growth of post-crack, joint coalescence failure, and residual strength. According to strength analysis, the failure envelopes of intact rock and across jointed rock mass are the upper bound and lower bound separately. The strength of intermittent jointed rock mass can be evaluated by reducing the bandwidth of the failure envelope with geo-mechanics analysis. (4) Numerical test of rock mass Two sets of methods, i.e. the distinct element method (DEC) based on in-situ geology mapping and the realistic failure process analysis (RFPA) based on high-definition digital imaging, are developed and introduced. The operation process and analysis results are demonstrated detailedly from the research on parameters of rock mass based on numerical test in the Jinping First Stage Hydropower Station and Baihetan Hydropower Station. By comparison，the advantages and disadvantages are discussed. Then the applicable fields are figured out respectively. (5) Intelligent evaluation based on artificial neural network (ANN) The characters of both ANN and parameter evaluation of rock mass are discussed and summarized. According to the investigations, ANN has a bright application future in the field of parameter evaluation of rock mass. Intelligent evaluation of mechanical parameters in the Jinping First Stage Hydropower Station is taken as an example to demonstrate the analysis process. The problems in five aspects, i. e. sample selection, network design, initial value selection, learning rate and expected error, are discussed detailedly.

大型反倾岩体变形破坏机制及稳定性分析———以金沙江宽谷坝址龙蟠边坡为例

A large number of catastrophic accidents were aroused by the instability and destruction of anti-dip rock masses in the worldwide engineering projects, such as hydropower station, mine, railways and so on. Problems in relation to deformation and failure about anti-dip rock slopes are significant for engineering geology research. This dissertation takes the Longpan slope in the Jinsha River as a case to study the deformation mechanism of large-scale anti-dip rock masses and the slope stability analysis method. The primary conclusions are as follows. The Dale Reach of Jinsha River, from Longpan to the debouchment of Chongjiang tributary, is located in the southeastern margin of the Qinghai-Tibet Plateau. Longpan slope is the right embankment of Dale dam, it is only 26 km to the Shigu and 18 km to Tiger Leaping Gorge. The areal geology tectonic structures here area are complicated and blurry. Base on the information of geophysical exploration (CSAMT and seismology) and engineering geological investigation, the perdue tectonic pattern of Dale Reach is put forward for the first time in this paper. Due to the reverse slip of Longpan fault and normal left-rotation of Baihanchang fault, the old faulted valley came into being. The thick riverbed sediments have layered characters of different components and corresponding causes, which attribute to the sedimentary environments according with the new tectonic movements such as periodic mountain uplifting in middle Pleistocene. Longpan slope consists of anti-dip alternate sandstone and slate stratums, and the deformable volume is 6.5×107m3 approximately. It was taken for an ancient landslide or toppling failure in the past so that Dale dam became a vexed question. Through the latest field surveying, displacement monitoring and rock masses deforming characters analyses, the geological mechanism is actually a deep-seated gravitational bending deformation. And then the discrete element method is used to simulate the deforming evolution process, the conclusion accords very well with the geo-mechanical patterns analyses. In addition strength reduction method based on DEM is introduced to evaluate the factor of safety of anti-dip rock slope, and in accordance with the expansion way of the shear yielding zones, the progressive shear failure mechanism of large-scale anti-dip rock masses is proposed for the first time. As an embankment or a close reservoir bank to the lower dam, the stability of Longpan slope especially whether or not resulting in sliding with high velocity and activating water waves is a key question for engineering design. In fact it is difficult to decide the unified slip surface of anti-dip rock slope for traditional methods. The author takes the shear yielding zones acquired form the discrete element strength reduction calculation as the potential sliding surface and then evaluates the change of excess pore pressure and factor of stability of the slope generated by rapid drawdown of ponded water. At the same time the dynamic response of the slope under seismic loading is simulated through DEM numerical modeling, the following results are obtained. Firstly the effective effect of seismic inertia force is resulting in accumulation of shear stresses. Secondly the discontinuous structures are crucial to wave transmission. Thirdly the ultimate dynamic response of slope system takes place at the initial period of seismic loading. Lastly but essentially the effect of earthquake load to bringing on deformation and failure of rock slope is the coupling effect of shear stresses and excess pore water pressure accumulation. In view of limitations in searching the critical slip surface of rock slope of the existing domestic and international software for limit equilibrium slope stability analyses, this article proposes a new method named GA-Sarma Algorithm for rock slope stability analyses. Just as its name implies, GA-Sarma Algorithm bases on Genetic Algorithm and Sarma method. GA-Sarma Algorithm assumes the morphology of slip surface to be a broken line with traceability to extend along the discontinuous surface structures, and the slice boundaries is consistent with rock mass discontinuities such as rock layers, faults, cracks, and so on. GA-Sarma Algorithm is revolutionary method that is suitable for global optimization of the critical slip surface for rock slopes. The topics and contents including in this dissertation are closely related to the difficulties in practice, the main conclusions have been authorized by the engineering design institute. The research work is very meaningful and useful for the engineering construction of Longpan hydropower station.

地震波场模拟及叠前深度偏移的李群方法

The processes of seismic wave propagation in phase space and one way wave extrapolation in frequency-space domain, if without dissipation, are essentially transformation under the action of one parameter Lie groups. Consequently, the numerical calculation methods of the propagation ought to be Lie group transformation too, which is known as Lie group method. After a fruitful study on the fast methods in matrix inversion, some of the Lie group methods in seismic numerical modeling and depth migration are presented here. Firstly the Lie group description and method of seismic wave propagation in phase space is proposed, which is, in other words, symplectic group description and method for seismic wave propagation, since symplectic group is a Lie subgroup and symplectic method is a special Lie group method. Under the frame of Hamiltonian, the propagation of seismic wave is a symplectic group transformation with one parameter and consequently, the numerical calculation methods of the propagation ought to be symplectic method. After discrete the wave field in time and phase space, many explicit, implicit and leap-frog symplectic schemes are deduced for numerical modeling. Compared to symplectic schemes, Finite difference (FD) method is an approximate of symplectic method. Consequently, explicit, implicit and leap-frog symplectic schemes and FD method are applied in the same conditions to get a wave field in constant velocity model, a synthetic model and Marmousi model. The result illustrates the potential power of the symplectic methods. As an application, symplectic method is employed to give synthetic seismic record of Qinghai foothills model. Another application is the development of Ray+symplectic reverse-time migration method. To make a reasonable balance between the computational efficiency and accuracy, we combine the multi-valued wave field & Green function algorithm with symplectic reverse time migration and thus develop a new ray+wave equation prestack depth migration method. Marmousi model data and Qinghai foothills model data are processed here. The result shows that our method is a better alternative to ray migration for complex structure imaging. Similarly, the extrapolation of one way wave in frequency-space domain is a Lie group transformation with one parameter Z and consequently, the numerical calculation methods of the extrapolation ought to be Lie group methods. After discrete the wave field in depth and space, the Lie group transformation has the form of matrix exponential and each approximation of it gives a Lie group algorithm. Though Pade symmetrical series approximation of matrix exponential gives a extrapolation method which is traditionally regarded as implicit FD migration, it benefits the theoretic and applying study of seismic imaging for it represent the depth extrapolation and migration method in a entirely different way. While, the technique of coordinates of second kind for the approximation of the matrix exponential begins a new way to develop migration operator. The inversion of matrix plays a vital role in the numerical migration method given by Pade symmetrical series approximation. The matrix has a Toepelitz structure with a helical boundary condition and is easy to inverse with LU decomposition. A efficient LU decomposition method is spectral factorization. That is, after the minimum phase correlative function of each array of matrix had be given by a spectral factorization method, all of the functions are arranged in a position according to its former location to get a lower triangular matrix. The major merit of LU decomposition with spectral factorization (SF Decomposition) is its efficiency in dealing with a large number of matrixes. After the setup of a table of the spectral factorization results of each array of matrix, the SF decomposition can give the lower triangular matrix by reading the table. However, the relationship among arrays is ignored in this method, which brings errors in decomposition method. Especially for numerical calculation in complex model, the errors is fatal. Direct elimination method can give the exact LU decomposition But even it is simplified in our case, the large number of decomposition cost unendurable computer time. A hybrid method is proposed here, which combines spectral factorization with direct elimination. Its decomposition errors is 10 times little than that of spectral factorization, and its decomposition speed is quite faster than that of direct elimination, especially in dealing with a large number of matrix. With the hybrid method, the 3D implicit migration can be expected to apply on real seismic data. Finally, the impulse response of 3D implicit migration operator is presented.