A QUASI-SIMPLEX METHOD FOR STRICTLY CONVEX QUADRATIC PROGRAMMING

本文提出一个不用 Kuhn- Tucker条件而直接搜索严格凸二次规划最优目标点的鲁棒方法 .在搜索过程中 ,目标点沿约束多面体边界上的一条折线移动 .这种移动目标点的思想可以被认为是线性规划单纯形法的自然推广 ,在单纯形法中 ,目标点从一个顶点移到另一个顶点。

UUV推进系统容错控制分配算法研究

针对水下机器人(UUV)推进系统容错控制分配问题,本文提出了基于SVD分解(奇异值分解)与定点分配的混合算法。与传统的方法相比,它回避了求伪逆矩阵的问题,降低了计算量;能够满足推进器饱和约束限制。利用水下实验平台推进系统模型进行了仿真实验,验证了算法的正确性和有效性。

一种新型四自由度并联机器人运动学正问题的求解

对一种新型四自由度并联机器人运动学正问题进行了研究，利用３个变量构造 出求解正问题的３个约束方程，然后运用符号计算和析配消元法推导出了只含有一 个变元的３２次多项式方程，并且应用计算机软件系统Ｍａｔｈｅｍａｔｉｃａ进行了求正问 题实解的数值验证。

双机械臂协调系统的一类最优轨迹规划方法

对Allen的时间理论的某些改进

Ａｌｌｅｎ的时间理论因直观、易懂而倍受推崇，但它存在不能处理连续变化事件等缺欠。本文提出更为一般的时间理论框架，以扩展Ａｌｌｅｎ的理论，本框架的特点为：（ｌ）将Ａｌｌｅｎ的理论纳入其中；（２）可由时间点构造时区，并可处理时区和时间点；（３）以时间元素集的２Ｄ图形表示为基础的约束传播算法，既高效又可视化。

三维物体识别及姿态测定的推理系统

本文介绍了三维物体识别及姿态测定的一种新技术，从物体空间域模型出发，通过约束推理及几何推理，在物体三维信息部分给定的条件下，推断预测图象模型，并通过实测的图象数据反馈，推断出隐含在图象中未给定的三维信息，最终实现三维物体识别及姿态测定。整个系统在ＶＩＣＯＭ机上用Ｃ语言完成。

基于动力学约束的机械臂最优关节轨迹规划

本文提出了基于机械臂关节驱动力矩约束方程规划其关节最优运动轨迹的一种有效方法.该方法运用矩阵范数理论简化机械臂的动力学约束方程;在机械臂的关节空间内采用归一化的无因次量运用非线性规划法优化其运动轨迹.将所规划的无因次量轨迹方程作为机械臂产生实际运动轨迹的发生器,通过给定机械臂各运动段的起始和终止关节坐标,由系统的动力学约束方程计算出整个运动段所允许的最短运行时间,即生成所期望的运动轨迹.本文的轨迹规划方法计算效率高,可用于在线轨迹规划,文中通过算例证实了该方法的实用性.

柴达木盆地地震勘探关键技术研究

The rugged surface topography determined the seismic data acquisition construction conditions and the seismic wave explosive and receiver quality in Qaidam Basin. This dissertation systematically researched the seismic acquisition, imaging process and the attribute analysis techniques of complicated oil and gas reservoir. The main research achievements and cognitions are as follows: 1. Through the stimulation effects research and analysis from the aspect of lithologic water-containing differences, it’s specific that stable hydrous sand layer can effectively enhance the stimulation effects combined with the corresponding field tests. The seismic data S/N ratio has been improved due to the combination explosive stimulation. Through the fold number and maximum offset analyses of target horizon, the complicated geometry has been optimized and the S/N ratio of seismic data has been improved, which made an important basis for improvement of 3D seismic data. 2. It has been proved that the first arrival refraction static correction method under the model constraint of fine surface survey is suitable to the Qaidam Basin of western areas by the real seismic data processing. Although the refraction horizon of near surface has some changes in a certain extent, it’s steady basically. The refraction horizon can be continuously traced in sections, so it’s qualified for the refraction static correction method on the whole. 3. The research is based on the curved-ray pre-stack time migration techniques of rough topography, and improved the imaging precision of complex areas. This techniques adopted the constant and variable velocity scanning mode and enhanced the velocity analysis precision. The 3D pre-stack time migration techniques reasonably solved the imaging and velocity multiple solutions problems of steep-dip faults and the intersections of horizontal layers. What’s more, fine velocity analysis and mute are very important to enhance the imaging precision of the seismic data in complicated Wunan areas. 4. The 3D seismic data edge-preserving processing methods have been realized due to the image process techniques. Because this method uses the large range filter, it can attenuate the noise maximally. The faults, break points, lithologic pinchout points and lithologic body of small scale such as river will not be influenced by blur because of the edge-preserving characterization of the method which is really an effective assistant technique of low S/N ratio seismic data attribute analysis. 5. The use of spectral decomposition technique can effectively identify the reservoirs. The special geology body which will not be identified (or without obvious characters) in the seismic profile may be found through the details changes of different frequencies in the amplitude profiles.

层析静校正技术研究和应用

Static correction is one of the indispensable steps in the conventional onshore seismic data processing, particularly in the western part of China; it is theoretically and practically significant to resolve the issue of static correction. Conventional refraction static correction is put forward under the assumption that layered medium is horizontal and evenly distributed. The complicated nature of the near surface from western part of China is far from the assumption. Therefore, the essential way to resolve the static correction problem from the complex area is to develop a new theory. In this paper, a high-precision non-linear first arrival tomography is applied to solve the problem, it moved beyond the conventional refraction algorithm based on the layered medium and can be used to modeling the complex near surface. Some of the new and creative work done is as follows: One. In the process of first arrival tomographic image modeling, a fast high-order step algorithm is used to calculate the travel time for first arrival and ray path and various factors concerning the fast step ray tracing algorithm is analyzed. Then the second-order and third-order differential format is applied to the step algorithm which greatly increased the calculation precision of the ray tracing and there is no constraint to the velocity distribution from the complex areas. This method has very strong adaptability and it can meet the needs of great velocity variation from the complicated areas. Based on the numerical calculation, a fast high-order step is a fast, non-conditional and stable high-precision tomographic modeling algorithm. Two, in the tomographic inversion, due to the uneven fold coverage and insufficient information, the inversion result is unstable and less reliable. In the paper, wavelet transform is applied to the tomographic inversion which has achieved a good result. Based on the result of the inversion from the real data, wavelet tomographic inversion has increased the reliability and stability of the inversion. Three. Apply the constrained high-precision wavelet tomographic image to the static correction processing from the complex area. During tomographic imaging, by using uphole survey, refraction shooting or other weathering layer method, weathering layer can be identified before the image. Because the group interval for the shot first arrival is relatively big, there is a lack of precision for the near surface inversion. In this paper, an inversion method of the layer constraint and well constraint is put forward, which can be used to compensate the shallow velocity of the inversion for the shot first arrival and increase the precision of the tomographic inversion. Key words: Tomography ，Fast marching method，Wavelet transform, Static corrections, First break

南海大陆边缘地壳结构与盆地形成机制--以珠江口－琼东南盆地为例

South China Sea is located in the convergence of Eurasian plate, the Pacific Ocean plate and Indian Ocean-Australia plate. The total area is about 3,500,000 km2, the geologic structure is complicated, and the structure line cut off reciprocal is the marginal sea taking form by that the seafloor spreads during the middle Oligocene. South China Sea continental margin have developed more than 10 large oil-gas bearing basins and a number of medium-small sized basins. These basins contain abundant mineral resources such as oil & gas. The marginal deepwater area in the north part of South China Sea has become our country’s strategic energy prospecting frontier. The deepwater area of Zhujiangkou and Qiongdongnan basins is the research target in this thesis. The thesis studied deep structure and the earth dynamics of the north part of South China Sea margin, and these researches provide scientific basis for oil-gas resources strategic investigation and valuation in deepwater sea area of north part slope of South China Sea. In order to develop the research of rebuilding velocities and density architecture of earth shell in region of interest, in marginal deepwater area in the north part of South China, we adopted 14 long-cable seismic reflection profile data of 3556.41 kilometers in total, the gravity measurement data along profiles (3851.44 kilometers in total), the magnetic observation along profiles (3838.4 kilometers in total) and depth measurement along profile, the logging data of 11 wells in project, the interpreted fault parameter and preexisting geologic and geophysical research achievement. This thesis has carried out concretely studying research as follows: 1. Overlay-velocity data sampling and analysis, interval velocity calculation, time-depth conversion, model building of earth shell velocity and layering character of earth shell are studied on 14 deep sections. Velocity structure in region of interest has revealed: Changchang is the sag with thinnest crust in Qiongdongnan basin; the sedimentary thickness lowers gradually from north to south, and the thickness change from west to east is milder. The sags’ sedimentary velocities in Qiongdongnan basin have obvious demarcation. The velocity of the 8000 meters sedimentary rocks is 4700 m/s in Shunde sag and Baiyun sag, and is the lowest; at that depth, the velocity very different in Liwan sag and Baiyun sag, which is about 800m/s. 2. Extracting gravity data and building of initial crust density model along the section; With Bouguer gravity anomaly data as constraint, revising density distributes of initial model, and building the crust density model. 3. With crust velocity and density as constraint, correcting the effect of thermobaric field and constructing constitution structure of rock in region of interest. By this research, we known that rocks in Zhujiangkou upper crustal layer are chiefly granite-gneiss, quartzite, granodiorite and basalt, however, rocks in Qiongdongnan basin upper earth shell are chiefly composed of granite-gneiss, quartzite, granodiorite, diorite and basalt. 4. Synthetically crust velocity and density structure, gaining expanding factor on crust and entire crust along section. The result is indicated: the expanding factor in every sag rises from northwest to southeast, which have reflected thinning characteristic of crust from continent to ocean. Intra-crustal deformation degree in Changchang and Ledong-Lingshui sag is bigger than that in Songnan-Baodao sag. Entire crust extension factor in Changchang and Songnan-Baodao sag is greater than that in Ledong-Lingshui sag, which can make an explanation of frequently event and longer heating process in middle-east of Qiongdongnan basin. 5. Synthesize multidisciplinary information to discuss the earth dynamics significance of discordogenic seismic profile in deepwater area of Zhujiangkou and Qiongdongnan basins.

青岛晚白垩世岩石圈结构特征：劈石口晚白垩世基性脉岩及其捕虏体提供的证据

The occurrence of Late Cretaceous mafic dykes and their entrained peridotite and granulite xenoliths as well as clinopyroxene xenocrysts in the Qingdao region provide us a precious opportunity to unveil the nature and characteristics of the Late Mesozoic lithospheric mantle and lower crust beneath the Jiaodong region, and the change of the magma sources. These studies are of important and significant for understanding the lithospheric evolution in the eastern North China Craton. There were two periods of magma activities in Late Mesozoic in Qingdao Laoshan region, one was around 107Ma in the Early Cretaceous and the other around 86Ma in the Late Cretaceous according to the whole rock K-Ar age determination. The Early Cretaceous mafic dykes and the Late Cretaceous mafic dyke (i.e. Pishikou mafic dike) have completely different geochemical characteristics. The Early Cretaceous mafic dykes are enriched in LILE, strongly depleted in HFSE (Nb, Ta, Zr, Hf) and characterized by the highly radiogenic Sr and Nd isotopic compositions. These geochemical features indicate that the Early Cretaceous mafic dykes were derived from an enriched lithospheric mantle. In contrast, the Late Cretaceous mafic dyke is enriched in LILE, without HFSE depletion (Nb, Ta, Zr, Hf) and has less radiogenic Nd and Sr isotopic compositions. These geochemical features indicate that the Late Cretaceous mafic dyke was derived from the asthenosphere modified by subducted pelagic sediment contamination. The intrusive age of the Late Cretaceous mafic dyke provides further information for the termination of the lithosphere thinning for the eastern North China Crtaon. Pishikou Late Cretaceous mafic dyke contains abundant peridotitic xenoliths, granulite xenoliths and clinopyroxene xenocrysts. The peridotitic xenoliths can be divided into two types: high Mg# peridotites and low Mg# peridotites, according to their textural and mineral features. The high-Mg# peridotites have high Fo (up to 92.2) olivines and high Cr# (up to 55) spinels. The clinopyroxenes in the high# peridotites are rich in Cr2O3 and poor in Al2O3. The low-Mg# peridotites are typified by low Mg# (Fo <90) in olivines and low Cr# (Cr# <0.14) in spinels. The clinopyroxenes in the low-Mg# peridotites are rich in Al2O3 and Na2O and poor in Cr2O3. These two type peridotites have similar equilibrated temperatures of 950C-1100C. The Clinopyroxenes in the high-Mg# peridotites generally have high and variable REE contents (REE = 5.6-84 ppm) and LREE-enriched chondrite-normalized patterns ((La/Yb)N>1). In contrast, the clinopyroxenes in the low-Mg# peridotites have low REE contents (REE = 12 ppm) and LREE-depleted patterns ((La/Yb)N<1). The textural, mineral and elemental features of the low-Mg# peridotites are similar to those of the low-Mg peridotites from the Junan, representing the newly-accreted lithospheric mantle. However, the mineralogical and petrological features of the high-Mg# peridotites are similar to those of the high-Mg# peridotites from the Junan region, representing samples from the old refractory lithospheric mantle that was strongly and multiply affected by melts of different origins Late Cretaceous mafic dike in the Qingdao region also contains two types of granulite xenoliths according to the mineral constituents: the pyroxene-rich granulites and the plagioclase-rich granulites. Equilibrated temperatures calculated from the cpx-opx geothermometers are in a range of 861C - 910C for the pyroxene-rich granulites and of 847C - 890C for the plagioclase-rich granulites. The equilibrated pressure for the plagioclase-rich granulites is in a range of 9.9-11.7 kbar. Combined with the results of the peridotitic xenoliths, a 40C temperature gap exists between the peridotite and the granulite. The petrological Moho was 33~36 km at depths, broadly consistent with the seismic Moho estimated from the geophysical data. This indicates that there was no obvious crust-mantle transition zone in the Qingdao region in the Late Mesozoic. Pishikou Late Cretaceous mafic dyke entrained lots of clinopyroxene xenocrysts which are characterized by the chemical zoning. According to the zoning features, two types of clinopyroxene xenoliths can be classified, the normal zoning and the revise zoning. The normally-zoned clinopyroxene xenocrysts have LREE-depleted REE patterns in the cores. In contrast, the revisely-zoned clinopyroxenes have LREE-enriched REE patterns in the cores. According to the rim and core compositions of xenocrysts, all the rims are balanced with the host magma. Meanwhile, the origins of the cores were complicated, in which the normally-zoned clinopyroxenes were derived form the lithospheric mantle and some of the reversely-zoned clinopyroxnes were originated from the lower crust. Other revisely-zoned clinopyroxenes had experienced complex geological evolution and need to be further investigated. According to the above results, a simplified lithospheric profile has been established beneath the Qingdao region and a constraint on the nature and characteristics of the lithospheric mantle and lower crust has been made.

反假频重建理论及其在储层预测中的应用

To extend the cross-hole seismic 2D data to outside 3D seismic data, reconstructing the low frequency data to high frequency data is necessary. Blind deconvolution method is a key technology. In this paper, an implementation of Blind deconvolution is introduced. And optimized precondition conjugate gradient method is used to improve the stability of the algorithm and reduce the computation. Then high-frequency retrieved Seismic data and the cross-hole seismic data is combined for constraint inversion. Real data processing proved the method is effective. To solve the problem that the seismic data resolution can’t meet the request of reservoir prediction in the river face thin-layers in Chinese eastern oil fields, a high frequency data reconstruction method is proposed. The extrema of the seismic data are used to get the modulation function which operated with the original seismic data to get the high frequency part of the reconstruction data to rebuild the wide band data. This method greatly saves the computation, and easy to adjust the parameters. In the output profile, the original features of the seismic events are kept, the common feint that breaking the events and adding new zeros to produce alias is avoided. And the interbeded details are enhanced compared to the original profiles. The effective band of seismic data is expended and the method is approved by the processing of the field data. Aim to the problem in the exploration and development of Chinese eastern oil field that the high frequency log data and the relative low frequency seismic data can’t be merged, a workflow of log data extrapolation constrained by time-phase model based on local wave decomposition is raised. The seismic instantaneous phase is resolved by local wave decomposition to build time-phase model, the layers beside the well is matched to build the relation of log and seismic data, multiple log info is extrapolated constrained by seismic equiphase map, high precision attributes inverse sections are produced. In the course of resolve the instantaneous phase, a new method of local wave decomposition --Hilbert transform mean mode decomposition(HMMD) is raised to improve the computation speed and noise immunity. The method is applied in the high resolution reservoir prediction in Mao2 survey of Daqing oil field, Multiple attributes profiles of wave impedance, gamma-ray, electrical resistivity, sand membership degree are produced, of which the resolution is high and the horizontal continuous is good. It’s proved to be a effective method for reservoir prediction and estimation.

牡丹江碰撞带的组成及演化

Jiamusi Massif is an important tectonic unit in Northeast China. It’s significant for understanding the evolution of Paleo-Asian Ocean and reconstruction of the tectonic framework of Northeast China. Mudanjiang area is located in the southern margin of Jiamusi Massif and is the key to understand the evolution of Jiamusi Massif. However, the detailed geological research for Mudanjiang area has long been deficient in many important problems, such as the tectonic components of the Mudanjiang collision zone (MCZ), the age of collisional complexes and the scenario of tectonic evolution. Based on the lithology, geochemistry and the SHRIMP zircon U-Pb geochronology in Mudanjiang area, our new data and results come to some constraints for the tectonic reconstruction of MCZ as follows: 1） It is identified that the former suggestion, which the so-called “Heilongjiang Group” in Mudanjiang area is the vestige of oceanic crust, is correct. The oceanic relics represent the Neo-Proterozoic-Early Paleozoic oceanic basins based on the SHRIMP zircon U-Pb geochronology. 2） One sheet of gabbroic complex with oceanic island-type geochemical signature was discovered by this work in Mudanjiang area. 3） It is verified that the Proterozoic concordant U-Pb ages of the migmatites developed along the southern margin of Jiamusi massif, which represent the events of magmatic intrusion, as the direct evidence for the existence of the Proterozoic crystalline basements of the Jiamusi Massif. Based on geochronology, we suggest that the migmatization and coeval S-type granite magmatism of the southern margin of Jiamusi Massif took place about 490Ma. 4） The island arc complex has been found in the Heilongjiang Group, and the oceanic relics was found distributing on both sides, as provided important constraint for the tectonic reconstruction of the MCZ. 5） ~440Ma metamorphic event and coeval post-collisional granite magmatism have been firmly identified in the MCZ and its southern neighboring area. Together with previous data obtained by other researchers, our conclusions on the reconstruction of the tectonic architecture and evolution of the MCZ as follows: 1） The orogenic assemblages developed in the Mudanjiang collisional zone are featured by a sequence of ancient active continental margins and ensuing orogenic processing. The Mashan Group is the reworking basement of Jiamusi Massif, whereas the Heilongjiang Group represents arc and oceanic complexes characterized by imbricate deep-seated sliced and slivering sheets due to multi-phases of thrusting and nappe stacking. 2） The northern sub-belt of MCZ is probably the arc-continent collisional boundary related to the closure of main oceanic basin. The collisional age can be constrained by the events of syn-orogenic migmatization of migmatite, coeval S-type granite magmatism and the related granulite-facies metamorphism. Therefore, we suggested the collisional age of northern sub-belt is probably Cambrian-Early Ordovician. The extensive granulite-facies metamorphism of the Mashan Group in Jiamusi Massif, as affirmed by former works, was probably related with the collisional event. 3） The southern sub-belt of the MCZ was possibly related with the closure of back-arc basin. We presumed that the collisional age of southern sub-belt is about Ordovician-Early Silurian according to the ~440Ma extensive metamorphism and the occurrence of coeval post-collisional granite magmatism. 4） The extant structural architecture of the MCZ is related to the multi-phases of intra-continental superimposition, which is characterized by the Mesozoic nappe structure.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

预应力结构抗滑桩设计计算方法及其在滑坡治理中的应用研究

Landslide is a kind of serious geological hazards and its damage is very great. In recent years, landslides become more and more frequent along with increase of scale of engineering constructions and cause greater loss. Consequently, how to protect landslides has become important research subject in the engineering field. This paper improves the method how to compute landslide thrust and solves the irrational problem in the design of piles because of the irrational landslide thrust according to the theory and technology of existed anti-slide piles and pre-stressed cable anti-slide piles. Modern pre-stressing technology has been introduced and load balancing method has been used to improve the stressing behavior of anti-slide piles. Anchor cables, anti-slide piles and modern pre-stressing technology have been used to prevention complicated landslide. It is an important base to select values for the landslide thrust. An improved method to calculate design thrust of anti-slide piles has been presented in this paper on the base of residual thrust method by comparing existing methods to select values of landslide thrust in the design of anti-slide piles. In the method, residual landslide thrust behind the anti-slide piles and residual skid resistance before the piles has been analyzed, equitable distribution of residual landslide thrust behind the piles has been realized, and the method to select value of design thrust becomes more reasonable. The pre-stressed cable anti-slide piles are developed from the common anti-slide piles and are common method to prevent landslide. Their principle is that internal force of anti-slide piles is adjusted and size of section is diminished by changing constraint conditions of anti-slide piles. For landslides with deep slip surface and large scale of slopes, limitation of the method appears. Such landslides are in need of long piles and anchor cables which are not only non-economic but also can generate larger deformation and leave potential danger after prevention. For solving the problem, a new kind of anti-slide piles, inner pre-stressing force anti-slide piles, is presented in this paper, and its principle is that an additional force, which is generated in the inner anti-slide piles by arranging pre-stressed reinforcement or tight wire in a certain form in interior of anti-slide piles and stretching the steel reinforcement or tight wire, may balance out the internal force induced by landslide thrust whole or partly (load balancing method). The method will change bending moment which anti-slide piles are not good at bearing into compressive stress which piles are good at bearing, improve stressing performance of anti-slide piles greatly, diminish size of section, and make anti-slide piles not fissured in the natural service or postpone appearance of the fissures, and improve viability of anti-slide piles. Pre-stressed cable anti-slide piles and inner pre-stressing force anti-slide piles go by the general name of pre-stressed structure anti-slide piles in the paper, and their design and calculation method is also analyzed. A new calculation method is provided in the paper for design of anti-slide piles. For pre-stressed structure anti-slide piles, a new computation mode is firstly presented in the paper on the foundation of cantilever piles. In the mode, constraint form of load-bearing section of the anti-slide piles should be confirmed according to reservoir conditions in order to figure out amount of pre-stress of the anchor cables, and internal force should be analyzed for the load-bearing section of pre-stressed structure anti-slide piles so as to confirm anchorage section of anti-slide piles. Pre-stressed cables of the pre-stressed cable anti-slide piles can be arranged as required. This paper analyzes the load-bearing section of single-row and double-row pre-stressed cable anti-slide piles and provides a calculation method for design of the pre-stressed cable anti-slide piles. Inner pre-stressing force anti-slide piles are a new kind of structural style. Their load-bearing section is divided into four computation modes according to whether pre-stressed cables are applied for exterior of the anti-slide piles, and whether single-row or double-row exterior pre-stressed cables are applied. The load balancing method is used to analyze the computation modes for providing a method to design the inner pre-stressing force anti-slide piles rationally. Pre-stressed cable anti-slide piles and inner pre-stressing force anti-slide piles are applied to research on Mahe landfall in Yalong Lenggu hydropower station by the improved method to select value of design thrust of anti-slide piles. A good effect is obtained in the analysis.