一种轮式移动机器人实时速度估计方法

环境和机器人自身的不确定性影响轮式移动机器人的轨迹跟踪控制性能,此时仅仅使用里程计往往不能正确表达机器人的状态信息。在无速度传感器的情况下,讨论了使用加速度传感器和位置传感器的输出实时估计轮式移动机器人的速度。首先使用滑模观测器进行里程计信号处理,然后对车体加速度信号进行带通滤波提取车体扰动信息,通过频域融合信号表达轮式移动机器人的速度,并针对正交轮式全方位移动机器人进行了轨迹跟踪控制研究。试验结果表明采用融合数据可以更准确提供机器人的状态信息并得到更好的控制性能。

Theoretical and experimental research on time-optimal trajectory planning and control of industrial robots

提出了一种用于工业机器人时间最优轨迹规划及轨迹控制的新方法，它可以确保在关节位移、速度、加速度以及二阶加速度边界值的约束下，机器人手部沿笛卡尔空间中规定路径运动的时间阳短。在这种方法中，所规划的关节轨迹都采用二次多项式加余弦函数的形式，不仅可以保证各关节运动的位移、速度 、加速度连续而且还可以保证各关节运动的二阶加速度连续。采用这种方法，既可以提高机器人的工作效率又可以延长机器人的工作寿命以PUMA560机器人为对象进行了计算机仿真和机器人实验，结果表明这种方法是正确的有效的。它为工业机器人在非线性运动学约束条件下的时间最优轨迹规划及控制问题提供了一种较好的解决方案。

鄂尔多斯盆地上古生界低渗透砂岩有效储层识别与预测

Sulige Gasfield, with a basically proven reserve as high as one trillion cubic meters, is one giant gas field discovered in China. The major gas -bearing layers are Upper Paleozoic strata with fluvial-lacustrine sedimentary facies. Generally, gas reservoirs in this field are characteristic by "five low" properties, namely low porosity, low permeability, low formation pressure, low productivity and low gas abundance. Reservoirs in this field also feature in a large distribution area, thin single sandbody thickness, poor reservoir physical properties, thin effective reservoir thickness, sharp horizontal and/or vertical changes in reservoir properties as well as poor connectivity between different reservoirs. Although outstanding achievements have been acquired in this field, there are still several problems in the evaluation and development of the reservoirs, such as: the relation between seismic attributes and reservoir property parameters is not exclusive, which yields more than one solution in using seismic attributes to predict reservoir parameters; the wave impedance distribution ranges of sandstone and mudstone are overlapped, means it is impossible to distinguish them through the application of post-stack impedance inversion; studies on seismic petrophysics, reservoir geophysical properties, wave reflection models and AVO features have a poor foundation, makes it difficult to recognize the specific differences between tight sandstone and gas-bearing sandstone and their distribution laws. These are the main reasons causing the low well drilling success rate and poor economic returns, which usually result in ineffective development and utilization of the field. Therefore, it is of great importance to perform studies on identification and prediction of effective reservoirs in low permeable sandstone strata. Taking the 2D and 3D multiwave-multicomponent seismic exploration block in Su6-Su5 area of Sulige field as a study area and He 8 member as target bed, analysis of the target bed sedimentary characteristics and logging data properties are performed, while criteria to identify effective reservoirs are determined. Then, techniques and technologies such as pre-stack seismic information (AVO, elastic impedance, wave-let absorption attenuation) and Gamma inversion, reservoir litological and geophysical properties prediction are used to increase the precision in identifying and predicting effective reservoirs; while P-wave and S-wave impedance, ratio of P/S wave velocities, rock elastic parameters and elastic impedance are used to perform sandstone gas-bearing property identification and gas reservoir thickness prediction. Innovative achievements are summarized as follows: 1. The study of this thesis is the first time that multiwave-multicomponent seismic data are used to identify and predict non-marine classic reservoirs in China. Through the application of multiwave-multicomponents seismic data and integration of both pre-stack and post-stack seismic data, a set of workflows and methods to perform high-precision prediction of effective reservoirs in low permeable sandstone is established systematically. 2. Four key techniques to perform effective reservoir prediction including AVO analysis, pre-stack elastic wave impedance inversion, elastic parameters inversion, and absorption attenuation analysis are developed, utilizing pre-stack seismic data to the utmost and increasing the correct rate for effective reservoir prediction to 83% from the former 67% with routine methods. 3. This thesis summarizes techniques and technologies used in the identification reservoir gas-bearing properties using multiwave-multicomponent seismic data. And for the first time, quantitative analysis on reservoir fluids such as oil, gas, and/or water are carried out, and characteristic lithology prediction techniques through the integration of pre-stack and post-stack seismic prediction techniques, common seismic inversion and rock elastic parameters inversion, as well as P-wave inversion and converted wave inversion is put forward, further increasing the correct rate of effective reservoir prediction in this area to 90%. 4. Ten seismic attribute parameters are selected in the 3D multi-wave area to perform a comprehensive evaluation on effective reservoirs using weighted-factor method. The results show that the first class effective reservoir covers an area of 10.08% of the study area, while the second and the third class reservoirs take 43.8% and 46% respectively, sharply increasing the success rate for appraisal and development wells.

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

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 anchorages are unparalleled structures only in a suspension bridge, and as main bearing facilities, play an important role in connecting the superstructures and the ground. The tunnel anchorage, as one alternative type of the anchorages, has more advantages over its counterpart, the gravity anchorage. With the tunnel anchorages adopted, not only can surface excavation be reduced to protect the environment, and natural condition of the rock be utilized and potential bearing capacity of surrounding rock be mobilized to save engineering cost, but also the technological predominance of auxiliary engineering measures, such as prestressed concrete, anchoring piles, rock anchors and collar beam between the two separated anchorages, can be easily cooperated to work together harmoniously under the circumstances of poor rock quality. There are plentiful high mountains and deep canyons in west part of China, and long-span bridge construction is inevitably encountered in order to realize leapfrogging development of the transportation infrastructure. Western mountainous areas usually possess the conditions for constructing tunnel anchorages, and therefore, the tunnel anchorages, which are conformed to the conception of resource conservative and sustainable society, extremely have application and popularization value in western underdeveloped region. The scientific and technological problem about the design, construction and operation of tunnel anchorages should be further investigated. Combining the engineering of western tunnel anchorages for the Balinghe Suspension Bridge, this paper probed into the survey method and in-situ test method for tunnel anchorages, scientific rock quality evaluation of surrounding rock to provide reasonable physical and mechanical parameters for design, construction and operation of tunnel anchorages, bearing capacity estimation for tunnel anchorage, deformation prediction of the anchorage-rockmass system, tunnel-anchorage slope stability analysis and the evaluation of excavation stability and degree of safety of the anchorage tunnel. The following outcomes were obtained: 1. Materials of tunnel anchorages of suspension bridge built (and in progress) at home and abroad were systematically sorted out, with the engineering geological condition and geomechanical property of surrounding rock around the anchorage tunnel, the design size of anchorages and the construction method of anchorage tunnel paid more emphasis on, to unveil the internal relationship between the engineering geological conditions of surrounding rock and the design size and axis angle of anchorages and provide references for future design, construction and study of tunnel anchorages. 2. Physical and mechanical parameters were recommended based on three domestic and foreign methods of rock quality evaluation. 3. In-situ tests, adopting the back-thrust method, of two kinds of reduced scale model, 1/30 and 1/20, for the tunnel anchorages were conducted in the declining exploration drift with rock mass at the test depth being the same as surrounding rock around real anchorages, and reliable field rockmass displacement data were acquired. Attenuation relation between the increment of distance from the anchorage and the decrement of rockmass displacement under maximum test load, and influential scope suffered by anchorage load were obtained. 4. Using similarity theory, the magnitude of real anchorage and rockmass displacement under design load and degree of safety of the anchorage system were deduced. Furthermore, inversion analysis to deformation modulus of slightly weathered dolomite rock, the surrounding rock of anchorage tunnel, was performed by the means of numerical simulation. 5. The influential law of the geometrical size to the limit bearing capacity of tunnel anchorage was studied. 6. Based on engineering geological survey data, accounting for the combination of strata layer and adverse discontinuities, the failure patterns of tunnel anchorage slope were divided into three modes: sliding of splay saddle pier slope, superficial-layer slippage, and deep-layer slippage. Using virtual work principle and taking anchorage load in account, the stability of the three kinds of failure patterns were analyzed in detail. 7. The step-by-step excavation of anchorage tunnel, the numerical overload and the staged decrement of rock strength parameters were numerically simulated to evaluate the excavation stability of surrounding rock around anchorage tunnel, the overload performance of tunnel anchorage, and the safety margin of strength parameters of the surrounding rock.

藏东南巨厚松散堆积体发育条件、成灾模式及防治对策——以川藏公路然乌至鲁朗段为例

The unique geologic, geomorphic and climatic conditions of southeast Tibet have made the region to develop the multi-style and frequently occurring geologic hazards, especially the collapses and landslides and debris flows along the section of Ranwu-Lulang in Sichuan-Tibet highway. However, most of those geologic hazards have close relationship with the loose accumulations. That is, the loose accumulations are the main carrier of most geologic hazards. Thereof, the huge-thick accumulations along the highway is regarded as the objective in the thesis to study the geologic background, hazarding model and mitigation methods comprehensively, based on the multi-disciplinary theories and former materials. First of all, in the paper, based on field engineering geologic investigations, the genetic type and the characteristics of spatiotemporal distribution of the huge-thick loose accumulations along the highway, have been analysized from the factors of regional geology and geomorphy and climate, as well as the coupling acting of those factors with inoculation and eruption of the loose accumulations geologic hazards. The huge-thick loose accumulations has complex genetic types and specific regulations of spatiotemporal distribution, closely controlled by the outer environment of the region. The accumulations are composed of earth and boulder, with disorder structure and poor sorting, specific forming environments and depositing conditions. And its physical and mechanic properties are greatly distinguished from rock and common earth inland. When Sichuan-Tibet highway was firstly constructed along the north bank of Purlung Tsangpo River, the huge-thick loose accumulations was cut into many high and steep slopes. Through the survey to the cut-slopes and systematic investigation to their failures, the combination of height and angle of the accumulations slope has been obtained. At the same time, the types of genetic structure of those cut-slopes are also analysized and concluded, as well as their failure models. It is studied in the paper that there are piaster, duality, multielement and complexity types in genetic structure, and rip-dump-repose, rip-shear-slip and weathering-flake types in failure models. Moreover, it is briefly introduced present engineering performance methods and techniques dealing with the deformation and failure of the accumulations cut-slope. It is also suggested that several new techniques of slope enforcement and the method of landslide and rockfall avoiding should be applied. The research of high and steep cut-slope along the highway has broadened the acknowledgement of the combination of cut-slope height and angle. Especially, the dissertation also has made the monographic studies about the geologic background and hazarding models and prevention methods of some classic but difficult accumulations geologic hazards. They are: (1) Research of the engineering geologic background of the 102 landslide group and key problems about the project of tunnel. The 102 landslide group is a famous accumulational one composed of glacial tills and glaciofuvial deposit. The tunnel project is a feasible and optional one which can solve the present plight of “sliding after just harnessing” in the 102 section. Based on the glacial geomorphy and its depositing character, distribution of seepage line, a few drillhole materials and some surveying data, the position of contact surface between gneiss and accumulations has been recognized, and the retreating velocities of three different time scales (short, medium and long term) have been approximately calculated, and the weathering thickness of gneiss has also been estimated in the paper. On the basis of above acknowledgement, new engineering geomechnic mode is established. Numerical analysis about the stability of the No.2 landslide is done by way of FLAC program, which supplies the conclusion that the landslide there develops periodically. Thereof, 4 projects of tunnel going through the landslide have been put forwards. Safety distance of the tunnel from clinohefron has been numerically analysized. (2) Research of the geologic setting and disaster model and hazard mitigation of sliding-sand-slope. From the geologic setting of talus cone, it is indicated that the sliding-sand-slope is the process of the re-transportation and re-deposit of sand under the gravity action and from the talus cone. It is the failure of the talus cone essentially. The layering structure of the sliding-sand-slope is discovered. The models of movement and failure of the sliding-sand-slope has been put forwards. The technique, “abamurus+grass-bush fence+degradable culture pan”, is suggested to enforcement and green the sliding-sand-slope. (3) Characteristics and hazarding model and disaster mitigation of debris flow. The sources of solid material of three oversize debris flows have been analysized. It is found that a large amount of moraine existing in the glacial valley and large landslide dam-break are the two important features for oversize debris flow to be taken place. The disaster models of oversize and common debris flows have been generalized respectively. The former model better interpret the event of the Yigong super-large landslide-dam breaking. The features of common debris flow along the highway section, scouring and silting and burying and impacting, are formulated carefully. It is suggested that check dam is a better engineering structure to prevent valley from steeply scouring by debris flow. Moreover, the function of check dam in enforcing the slope is numerically calculated by FLAC program. (4) Songzong ancient ice-dammed lake and its slope stability. The lacustrine profile in Songzong landslide, more than 88 meters thick, is carefully described and measured. The Optical Simulated Luminescence (OSL) ages in the bottom and top of the silty clay layer are 22.5±3.3 kaB.P., 16.1±1.7 kaB.P., respectively. It is indicated by the ages that the lacustrine deposits formed during the Last Glacial Maximum ranging from 25ka B.P. to 15ka B.P. The special characteristics of the lacustrine sediment and the ancient lake line in Songzong basin indicated that the lacustrine sediment is related to the blocking of the Purlung Tsangpo River by the glacier in Last Glacial Maximum from Dongqu valley. The characteristics of the lacustrine profile also indicate that the Songzong ice-dammed lake might run through the Last Glacial Maximum. Two dimensional numerical modeling and analysis are done to simulate the slope stability under the conditions of nature and earthquake by FLAC program. The factor of safety of the lacusrtine slope is 1.04, but it will take place horizontal flow under earthquake activity due to the liquefaction of the 18.33 m silt layer. The realign to prevent the road from landslide is suggested.

青藏高原及邻区面波频散反演--岩石圈-软流层体系结构研究

The Qinghai-Tibet Plateau lies in the place of the continent-continent collision between Indian and Eurasian plates. Because of their interaction the shallow and deep structures are very complicated. The force system forming the tectonic patterns and driving tectonic movements is effected together by the deep part of the lithosphere and the asthenosphere. It is important to study the 3-D velocity structures, the spheres and layers structures, material properties and states of the lithosphere and the asthenosphere for getting knowledge of their formation and evolution, dynamic process, layers coupling and exchange of material and energy. Based on the Rayleigh wave dispersion theory, we study the 3-D velocity structures, the depths of interfaces and thicknesses of different layers, including the crust, the lithosphere and the asthenosphere, the lithosphere-asthenosphere system in the Qinghai-Tibet Plateau and its adjacent areas. The following tasks include: (1)The digital seismic records of 221 seismic events have been collected, whose magnitudes are larger than 5.0 over the Qinghai-Tibet Plateau and its adjacent areas. These records come from 31 digital seismic stations of GSN , CDSN、NCDSN and part of Indian stations. After making instrument response calibration and filtering, group velocities of fundamental mode of Rayleigh waves are measured using the frequency-time analysis (FTAN) to get the observed dispersions. Furthermore, we strike cluster average for those similar ray paths. Finally, 819 dispersion curves (8-150s) are ready for dispersion inversion. (2)From these dispersion curves, pure dispersion data in 2°×2° cells of the areas (18°N-42°N, 70°E-106°E) are calculated by using function expansion method, proposed by Yanovskaya. The average initial model has been constructed by taking account of global AK135 model along with geodetic, geological, geophysical, receiving function and wide-angle reflection data. Then, initial S-wave velocity structures of the crust and upper mantle in the research areas have been obtained by using linear inversion (SVD) method. (3)Taking the results of the linear inversion as the initial model, we simultaneously invert the S wave velocities and thicknesses by using non-linear inversion (improved Simulated Annealing algorithm). Moreover, during the temperature dropping the variable-scale models are used. Comparing with the linear results, the spheres and layers by the non-linear inversion can be recognized better from the velocity value and offset. (4)The Moho discontinuity and top interface of the asthenosphere are recognized from the velocity value and offset of the layers. The thicknesses of the crust, lithosphere and asthenosphere are gained. These thicknesses are helpful to studying the structural differentia between the Qinghai-Tibet Plateau and its adjacent areas and among geologic units of the plateau. The results of the inversion will provide deep geophysical evidences for studying deep dynamical mechanism and exploring metal mineral resource and oil and gas resources. The following conclusions are reached by the distributions of the S wave velocities and thicknesses of the crust, lithosphere and asthenosphere, combining with previous researches. (1)The crust is very thick in the Qinghai-Tibet Plateau, varying from 60 km to 80 km. The lithospheric thickness in the Qinghai-Tibet Plateau is thinner (130-160 km) than its adjacent areas. Its asthenosphere is relatively thicker, varies from 150 km to 230 km, and the thickest area lies in the western Qiangtang. India located in south of Main Boundary thrust has a thinner crust (32-38 km), a thicker lithosphere of about 190 km and a rather thin asthenosphere of only 60 km. Sichuan and Tarim basins have the crust thickness less than 50km. Their lithospheres are thicker than the Qinghai-Tibet Plateau, and their asthenospheres are thinner. (2)The S-wave velocity variation pattern in the lithosphere-asthenosphere system has band-belted distribution along east-westward. These variations correlate with geology structures sketched by sutures and major faults. These sutures include Main Boundary thrust (MBT), Yarlung-Zangbo River suture (YZS), Bangong Lake-Nujiang suture (BNS), Jinshajiang suture (JSJS), Kunlun edge suture (KL). In the velocity maps of the upper and middle crust, these sutures can be sketched. In velocity maps of 250-300 km depth, MBT, BNS and JSJS can be sketched. In maps of the crustal thickness, the lithospheric thickness and the asthenospheric thickness, these sutures can be still sketched. In particular, MBT can be obviously resolved in these velocity maps and thickness maps. (3)Since the collision between India and Eurasian plate, the “loss” of surface material arising from crustal shortening is caused not only by crustal thickening but also by lateral extrusion material. The source of lateral extrusion lies in the Qiangtang block. These materials extrude along the JSJS and BNS with both rotation and dispersion in Daguaiwan. Finally, it extends toward southeast direction. (4)There is the crust-mantle transition zone of no distinct velocity jump in the lithosphere beneath the Qiangtang Terrane. It has thinner lithosphere and developed thicker asthenosphere. It implies that the crust-mantle transition zone of partial melting is connected with the developed asthenosphere. The underplating of asthenosphere may thin the lithosphere. This buoyancy might be the main mechanism and deep dynamics of the uplift of the Qinghai-Tibet hinterland. At the same time, the transport of hot material with low velocity intrudes into the upper mantle and the lower crust along cracks and faults forming the crust-mantle transition zone.

青藏高原横波速度与各向异性结构及其壳幔形变印迹

As the most spectacular and youngest case of continental collision on the Earth, to investigate the crust and mantle of Tibetan plateau, and then to reveal its characters of structure and deformation, are most important to understand its deformation mechanism and deep process. A great number of surface wave data were initially collected from events occurred between 1980 and 2002, which were recorded by 13 broadband digital stations in Eurasia and India. Up to 1,525 source-station Rayleigh waveforms and 1,464 Love wave trains were analysed to obtain group velocity dispersions, accompanying with the detail and quantitative assessment of the fitness of the classic Ray Theory, errors from focal and measurements. Assuming the model region covered by a mesh of 2ox2o-sized grid-cells, we have used the damped least-squares approach and the SVD to carry out tomographic inversion, SV- and SH-wave velocity images of the crust and upper mantle beneath the Tibetan Plateau and surroundings are obtained, and then the radial anisotropy is computed from the Love-Rayleigh discrepancy. The main results demonstrate that follows, a) The Moho beneath the Tibetan Plateau presents an undulating shape that lies between 65 and 74 km, and a clear correlation between the elevations of the plateau and the Moho topography suggests that at least a great part of the highly raised plateau is isostatically compensated. b) The lithospheric root presents a depth that can be substantiated at ~140 km (Qiangtang Block) and exceptionally at ~180 km (Lhasa Block), and exhibits laterally varying fast velocity between 4.6 and 4.7 km/s, even ~4.8 km/s under northern Lhasa Block and Qiangtang Block, which may be correlated with the presence of a shield-like upper mantle beneath the Tibetan Plateau and therefore looked as one of the geophysical tests confirming the underthrusting of India, whose leading edge might have exceeded the Bangong-Nujiang Suture, even the Jinsha Suture. c) The asthenosphere is depicted by a low velocity channel at depths between 140 and 220 km with negative velocity gradient and velocities as low as 4.2 km/s; d) Areas in which transverse radial anisotropy is in excess of ~4% and 6% on the average anisotropy are found in the crust and upper mantle underlying most of the Plateau, and up to 8% in some places. The strength, spatial configuration and sign of radial anisotropy seem to indicate the existence of a regime of horizontal compressive forces in the frame of the convergent orogen at the same time that laterally varying lithospheric rheology and a differential movement as regards the compressive driving forces. e) Slow-velocity anomalies of 12% or more in southern Tibet and the eastern edge of the Plateau support the idea of a mechanically weak middle-to-lower crust and the existence of crustal flow in Tibet.

多波多分量地震资料走时近似方法的研究与应用

Compared with the conventional P wave, multi-component seismic data can markedly provide more information, thus improve the quality of reservoir evaluation like formation evaluation etc. With PS wave, better imaging result can be obtained especially in areas involved with gas chimney and high velocity formation. However, the signal-to-noise of multi-component seismic data is normally lower than that of the conventional P wave seismic data, while the frequency range of converted wave is always close to that of the surface wave which adds to the difficulty of removing surface wave. To realize common reflection point data stacking from extracted common conversion point data is a hard nut to crack. The s wave static correction of common receiver point PS wave data is not easy neither. In a word, the processing of multi-component seismic data is more complicated than P wave data. This paper shows some work that has been done, addressing those problems mentioned above. (1) Based on the AVO feature of converted wave, this paper has realized the velocity spectrum of converted waves by using Sarkar’s generalized semblance method taking into account of AVO factor in velocity analysis. (2)We achieve a method of smoothly offset division normal method.Firstly we scan the stacking velocities in different offset divisions for a t0, secondly obtain some hyperbolas using these stacking velocities, then get the travel time for every trace using these hyperbolas; in the end we interpolate the normal move out between two t0 for every trace. (3) Here realize a method of stepwise offset division normal moveout.It is similar to the method of smoothly offset division normal moveout.The main difference is using quadratic curve, sixth order curve or fraction curve to fit these hyperbolas. (4)In this paper, 4 types of travel time versus distance functions in inhomogeneous media whose velocity or slowness varies with depth and vertical travel time have been discussed and used to approximate reflection travel time. The errors of ray path and travel time based on those functions in four layered models were analyzed, and it has shown that effective results of NMO in synthetic or real data can be obtained. (5) The velocity model of converted PS-wave can be considered as that of P -wave based on the ghost source theory, thus the converted wave travel time can be approximated by calculation from 4 equivalent velocity functions: velocity or slowness vary linearly with depth or vertical travel time. Then combining with P wave velocity analysis, the converted wave data can be corrected directly to the P-wave vertical travel time. The improvements were shown in Normal Move out of converted waves with numerical examples and real data. (6) This paper introduces the methods to compute conversion point location in vertical inhomogeneous media based on linear functions of velocity or slowness versus depth or vertical travel time, and introduce three ways to choose appropriate equivalent velocity methods, which are velocity fitting, travel time approximation and semblance coefficient methods.

中低纬度变化磁场的暴时特征研究

Magnetic storm is a kind of severe disturbances in the whole solar-earth electromagnetic space. It has significant effects on communication, electric power, oil transport pipe and human activities in space. Therefore, magnetic storms are worth for applications systems, not only being a favorable issue for scientists. In this paper, the spatial and temporal distributions of the magnetic fields produced by the magnetosphere-ionosphere current systems during storms are studied. Four parts are included in this paper decomposion of different disturbances with different origins, topological structure of the ring current, the asymmetric characteristics of the ring current, and the statistic peculiarities of the day-to-day variability (DTD) of Sq. 1 The decomposition of magnetic disturbances at mid-low latitudes and its evolutions during storms Transient variations in the geomagnetic field recorded at mid-low latitudes mainly include the storm-time variation (Dst), solar quiet daily variation (Sq) and disturbance daily variation (SD). With the data of the geomagnetic meridian chain observatories in China, 25 storms during the period of 1997 to 1999 have been analyzed. According to the features of different variations, a method of “three-steps decomposition” is developed by using the method of Natural Orthogonal Components (NOC), Correlation Analysis and Fourier Analysis to separate those three components in turn. The results show that, the first eigenmode by the MNOC clearly describing the special distribution and temporal evolution of storm-time variation, in addition, Correlation Analysis and Fourier Analysis offer a useful method to extract the Sq and SD variations. The latitudinal shift of the Sq current focus seems to be the principal reason of the day-to-day variaitons in the daily range of Sq. The magnitude of SD reaches a maximum during the main phase, and then gradually decreases. 2 The topology structure of the ring current during storms Both the mechanism of the ring current and the geomagnetic data suggest that the central plane of the ring current is declining to the geomagnetic equator plane with a tilt angle δ. Using the H and Z component data at two stations in a meridian chain, we deduce a new parameter describing the invariable peculiarity of different storms. Then the δ angle is calculated by using the data from a meridian chain and tested with the ERC model. Finally the deduced tilt angles are used to modify Dst index. 3 The asymmetric characteristics of the ring current during storms The variations of the geomagnetic field at mid-low latitudes show a significant dawn-dusk asymmetry, resulting from the superposition of the fields from the symmetric ring current and the partial ring current. On the basis of the data from the 20°E, 30°E meridian chains and 30°N latitudinal chain, the dawn-dusk asymmetry is investigated by using three methods, namely, statistic analysis, ring current model calculation and typical event analysis. This characteristic implies the asymmetry of the spatial distribution of the ring current. In addition, during the main phase after the sudden commencement (SC), H field increases and reaches maximum around noontime, implying the effect of the Chapman-Ferraro current. 4 The statistic characteristics of the day-to-day variability and its mechanism The day-to-day variability of the geomagnetic Sq field is studied by using the magnetic data from a meridian chain of magnetometers along 120° E longitude. The method of NOC is applied to separate the Sq variation from complicated disturbances. The first eigenmode with the largest eigenvalue represents fairly well the Sq variation with a conspicuous day-to-day variability in the daily range. For the stations on the same north- or south-side of the Sq current system focus, the day-to-day variations show a positive correlation. In contrast, for the stations on the different sides of the Sq focus, they show a negative correlation, suggesting an important role of latitudinal shift of the Sq current system focus to the day-to-day variability of the Sq daily range. The Sq daily range is correlated with the magnetic indices Ap and Dst in a peculiar way: on some severe disturbed days, noticeably enhancements of the Sq are observed, implying increases of the ionospheric conductivities and/or tidal wind velocities; on other severe disturbed days, however, dramatically reduced Sq variations occur, suggesting dominant effects of the ‘disturbance dynamo’ process.

叠前时间偏移: 基于波动方程的架构、发展和应用

This dissertation starts from the point that the prestack time migration can been considered as an approximation of the prestack depth migration, giving a wave equation based prestack time migration approach. The new approach includes: analytically getting the travel time and amplitude based on the one way wave equation and the stationary-phase theory, using ‘spread’ imaging method and imaging following the prestack depth migration, updating the velocity model with respect to the flats of the events in CRP gathers. Based on this approach, we present a scheme that can image land seismic data without field static correction. We may determine the correct near surface velocities and stack velocities by picking up the residual correction of the events in the CRP gathers. We may get the rational migration section based on the updated velocities and correct the migration section from a floating datum plane to a universal datum plane. We may adaptively determine the migration aperture according to the dips of the imaging structures. This not only speed up the processing, but may suppress the migration noise produce by the extra aperture. We adopt the deconvolution imaging condition of wave equation migration. It may partially compensate the geometric divergence. In this scheme, we use the table-driven technique which may enhance the computational efficiency. If the subsurface is much more complicated, it may be impossible to distinguish the DTS curve. To solve this problem, we proposed a technique to determine the appropriate range of the DTS curve. We synthesize DTS panel in this range using different velocities and depths, and stack the amplitude around the zero time. Determine the correct velocity and location of the considered grid point by comparing the values.

利用台阵面波资料研究华北克拉通北缘燕山褶皱带及其邻区的地壳与上地幔结构

The continent of eastern China, especially the North China Craton (NCC), has endured intensive tectonic renovation during Mesozoic and Cenozoic, with the presence of widespread magmatism, high heat flow and development of large sedimentary basins and mountain ranges. The cratonic lithosphere of the region has been destroyed remarkably, which is characterized by not only a significant reduction in thickness but also complex modifications in physical and chemical properties of the lithosphere. As for the tectonic regime controlling the evolution of the NCC, various models have been put forward, including the impingement of mantle plumes (“mushroom cloud” model), the collision of south China block and north China block, the subduction of the Pacific plate, etc. Lithosphere delamination and thermal erosion were proposed as the two end-member mechanisms of the lithospheric thinning. However, given the paucity of the data, deep structural evidence is currently still scarce for distinguishing and testifying these models. To better understand the deep structure of the NCC, from 2000 to the present, temporary seismic array observations have been conducted in the NCC by the Seismological Laboratory of the Institute of the Geology and Geophysics, Chinese Academy of Sciences under the North China Interior Structure Project (NCISP). Many arrays extend from the North China Craton and the off-craton regions, and traverse a lot of main tectonic boundaries. A total of more than 300 broadband seismic stations have been deployed along several profiles that traversed the major tectonic units within the craton’s interior, at the boundary areas and in the neighboring off-craton regions. These stations recorded abundant high-quality data, which provides an unprecedented opportunity for us to unravel the deep structural features of the NCC using seismological methods. Among all the seismological methods, the surface wave method appears to be an efficient and widely adopted technique in studying the crustal and upper mantle structures. In particular, it can provide the absolute values of S-wave velocity that are difficult to obtain with other methods. Benefiting from the deployment of dense seismic arrays, progresses have been made in improving the spatial resolution of surface wave imaging, which makes it possible to resolve the fine-scale velocity structures of the crust and upper mantle based on surface wave analysis. Meanwhile, the differences in the S-wave velocities derived from Rayleigh and Love wave data can provide information on the radial anisotropy beneath the seismic arrays. In this thesis, using the NCISP-III broadband data and based on phase velocity dispersion analysis and inversion of fundamental mode Rayleigh and Love waves, I investigated the lateral variations in the S-wave velocity structure of the crust and uppermost mantle beneath the Yanshan Belt and adjacent regions at the northeastern boundary of the NCC. Based on the constructed structural images, I discussed possible deep processes of the craton destruction in the study region.

地震层析成像与天山地区的壳幔深部结构研究

The Tien Shan is the most prominent intracontinental mountain belt on the earth. The active crustal deformation and earthquake activities provide an excellent place to study the continental geodynamics of intracontinental mountain belt. The studies of deep structures in crust and upper mantle are significantly meaningful for understanding the geological evolution and geodynamics of global intracontinental mountain belts. This dissertation focuses on the deep structures and geodynamics in the crust and upper mantle in the Tien Shan mountain belt. With the arrival time data from permanent and temporal seismic stations located in the western and central Tien Shan, using seismic travel time tomographic method, we inversed the P-wave velocity and Vp/Vs structures in the crust and uppermost mantle, the Pn and Sn velocities and Pn anisotropic structures in the uppermost mantle, and the P-wave velocity structures in the crust and mantle deep to 690km depth beneath the Tien Shan. The tomographic results suggest that the deep structures and geodynamics have significant impacts not only on the deformations and earthquake activities in the crust, but also on the mountain building, collision, and dynamics of the whole Tien Shan mountain belt. With the strongly collision and deformations in the crust, the 3-D P-wave velocity and Vp/Vs ratio structures are highly complex. The Pn and Sn velocities in the uppermost mantle beneath the Tien Shan, specially beneath the central Tien Shan, are significantly lower than the seismic wavespeed beneath geological stable regions. We infer that the hot upper mantle from the small-scale convection could elevate the temperature in the lower crust and uppermost mantle, and partially melt the materials in the lower crust. The observations of low P-wave and S-wave velocities, high Vp/Vs ratios near the Moho and the absences of earthquake activities in the lower crust are consistent with this inference. Based on teleseismic tomography images of the upper mantle beneath the Tien Shan, we infer that the lithosphere beneath the Tarim basin has subducted under the Tien Shan to depths as great as 500 km. The lithosphere beneath the Kazakh shield may have subducted to similar depths in the opposite direction, but the limited resolution of this data set makes this inference less certain. These images support the plate boundary model of converge for the Tien Shan, as the lithospheres to the north and south of the range both appear to behave as plates.

岩石中地震波速度和衰减的应力变化响应研究

Stress change is one of key factors in seismic nucleating and triggering; therefore for understanding and forecasting earthquakes, it is necessary to research on stress status and its changes in rocks. Propagating in underground structures, wave velocity and attenuation contain information on stress changes of the Earth’s interior. For a better understanding of relationship between seismic data and stress changes, modeling and ultrasonic test supply significant references. In this article, acoustoelastic theory is introduced to explain nonlinear elastic characteristics of rocks. Based on the acoustoelastic theory, a solid-fluid coupled model is given to calculate velocity under different stress for porous and liquid fulfilled rocks. Except for the stress-velocity relationship, effects of pore pressure induced stress changes on ultrasonic coda attenuation are also studied. Intrinsic attenuation quality factors are calculated for a comparison purpose. Finally, the relationship between elastic constants and stress changes is thoroughly investigated, a mixture model from two phases of Hooke media is introduced to explain the differences between dynamic and static moduli, a relation among wave length, wave velocities and elastic moduli considering dimension of microstructure, dimension and state of surface between phases is presented. The most important aspect of this work is exploring and establishing relationships between the seismic properties of rocks and changes of their stress conditions, which will have its application in earthquake forecast and seismic hazard.

基于三重震相波形模拟研究亚洲东北部上地幔过渡带速度结构

Both the global and regional P wave tomographic studies have revealed significant deep structural heterogeneities in subduction zone regions. In particular, low-velocity anomalies have been observed beneath the descending high-velocity slabs in a number of subduction zones. The limited resolution at large depths and possible trade-off between the high and low velocities, however, make it difficult to substantiate this feature and evaluate the vertical extent of the low-velocity structure. From broadband waveform modeling of triplicated phases near the 660-km discontinuity for three deep events, we constrained both the P and SH wave velocity structures around the base of the upper mantle in northeast Asia. For the two events beneath the southern Kurile, the rays traveled through the lowermost transition zone and uppermost lower mantle under the descending Pacific slab. Our preferred models consistently suggest normal-to-lower P and significantly low SH velocities above and below the 660-km discontinuity extending to about 760-km depth compared with the global IASP91 model, corroborating previous observations for a slow structure underneath the slab. In contrast, both high P and SH velocity anomalies are shown in our preferred models for the Japan subduction zone region, likely reflecting the structural feature of a slab stagnant above the 660-km discontinuity. The velocity jumps across the 660-km discontinuity were found to be on average 4.5% and 7% for P and S waves under the south Kurile, and 3% and 6% under the Japan subduction zone. The respective velocity contrasts in the two regions are consistent with mineralogical models for colder slab interior and hotter under-slab areas. Based on mineral physics data, the depth-averaged ~1.5% P and ~2.5% SH velocity differences in the depth range of 560-760 km between the two regions could be primarily explained by a 350~450K temperature variation, although the presence of about 0.5wt%~1wt% water might also contribute to the subtle velocity variations near the base of the transition zone in the southern Kurile. From our modeling results, we speculate that the slow structure in the southern Kurile may be correlated to the low velocity zone observed previously around the 410-km discontinuity under Northern Honshu. Both are probably associated with a thermal anomaly rooted in the lower mantle beneath the subduction zone in northeast Asia.