辐状肋柱花的苷类成分

利用硅胶柱层析、Sephadex LH-20及反相硅胶RP-18分离及纯化技术，从辐状肋柱花全草乙醇提取物的正丁醇萃取部分得到7个水溶性成分，经^1H NMR、^13C NMR等波谱技术鉴定为异荭草苷、芒果苷、Swertipunicoside、当药醇苷、异牡荆苷、当药黄素、和7-O-[α-L-吡喃鼠李糖-（1-2）-β-D-吡喃木糖]-1，8-二羟基-3-甲氧基[口山]酮。除异荭草苷外，其余化合物均首次从该植物中得到。

15种獐牙菜属植物中主要药用成分的高效液相色谱测定

对青藏高原和云贵高原的15种獐牙菜属植物进行了3种苦味苷,即獐牙菜苦苷(swertiamarin)、龙胆苦苷(gentiopicroside)、苦龙苷(amarogentin)、一种黄酮苷-当药黄素(swertisin)、及5种(口山) 山酮苷-芒果苷(mangiferin)、当药醇苷(swertianolin)、7-O-[α-L-吡喃鼠李糖-(1→2)-β-D-吡喃木糖]-1,8-二羟基-3-甲氧基(口山) 山酮(7-O-[α-L-rhamnopyranosyl-(1→2)-β-D-xylopyranosyl]-1,8-dihydroxy-3-methoxyxanthone)、7-O-β-D-吡喃木糖-1,8-二羟基-3-甲氧基(口山) 山酮(7-O-β-D-xylopyranosyl-l,8-dihydroxy-3-methoxyxanthone)、3-O-β-D-吡喃葡萄糖-1,8-二羟基-5-甲氧基(口山) 山酮(3-O-β-D-glucopyranosyl-1,8- dihydroxy-5-methoxyxanth-one)等9种主要药效成分同时进行了高效液相色谱的含量测定(Kromasil C18柱,甲醇-水梯度洗脱,二级管阵列检测);并对其主要药效成分的分布进行了比较。

青海南部七种翠雀属植物的核型

报道了采自青海南部７ 种翠雀属（Delphiniun L.） 植物的染色体数目和核型。大通翠雀花D.pylzowii 的核型公式为２ｎ＝ １６ ＝ ２ｍ ＋ ４ｓｍ ＋ １０st； 毛翠雀花D.trichophorum的核型公式为２ｎ＝ １６ ＝ ２ｍ ＋ ４ｓｍ ＋ １０st； 蓝白翠雀花D.albocoerelum 的核型公式为２ｎ ＝ １６ ＝ ２ｍ ＋ ４ｓｍ（２ＳＡＴ） ＋ １０st； 囊谦翠雀花D.nangchienense的核型公式为２ｎ ＝ １６ ＝ ２ｍ ＋ ６ｓｍ ＋ ８ｓｔ； 唐古拉翠雀花D.tangkulaense 的核型公式为2ｎ＝１６＝２ｍ（２ＳＡＴ）＋６ｓｍ＋８ｓｔ；单花翠雀花D.candelabrum var.monanthum 的核型2n=16=2m+6sm+8st;展毛翠雀花D.kamaoense var.glabrescens的核型公式为2n=16=2m+6sm+8st；前5种植物的染色体数目和核型为首次报道。

青藏高原高寒灌丛生态系统二氧化碳通量研究

本论文以青藏高原东北部海北地区高寒灌丛（Alpine Shrub）生态系统为研究对象，利用微气象观测系统及涡度相关（Eddy Covariance）技术，自2003年1月1日至2005年12月31日对该类广布于青藏高原的典型高寒草地类型进行长期连续观测。在对生态系统CO2净交换（NEE）以及群落叶面积指数（LAI）、生物量等生物学指标和光合有效辐射（PAR）、温度、土壤水分、脉冲性降水事件等主要环境因子进行连续监测的基础上，重点分析和探讨了海北地区高寒灌丛生态系统净生态系统CO2交换(NEE)在时、日、月及年际尺度上的变化模式，生长季与非生长季高寒灌丛生态系统CO2净交换特征，高寒灌丛生态系统大气CO2源/汇年际差异，土壤温度、昼夜温差、光合有效辐射、脉冲性降水事件等主要环境因子影响。从而，揭示了不同时间尺度下的高寒灌丛生态系统NEE变化规律，阐明主要环境因子对生态系统NEE的影响，明确了该生态系统大气CO2源/汇状况及其季节分布模式；同时，也为青藏高原区域尺度的高寒草地生态系统CO2通量研究和碳收支的估算提供科学依据和基础数据，对进一步揭示我国乃至亚洲陆地生态系统的碳收支状况有着重要意义。主要研究结果概括为以下几个方面： 1、海北地区高寒灌丛生态系统净生态系统CO2交换时动态特征存在很大的季节性差异，暖季小时NEE变化振幅大，CO2净吸收的极值一般出现在午间，最大吸收量为1.7 g CO2 m-2 h-1左右。夜间为CO2净释放，净生态系统交换值较为稳定（0.5~ 0.9 g CO2 m-2 h-1）；冷季日变化振幅极小，除14:00~18:00时一定量CO2释放外，其余时段通量均很小。 2、从日平均净生态系统CO2交换来看，6~9月日平均NEE一般为负值（CO2净吸收），2003~2005年6~9 月间日平均NEE分别为-5.65 g CO2 m-2 d-1、-6.08 g CO2 m-2 d-1和-4.81 g CO2 m-2 d-1；而10~12月及翌年1~5月期间日平均NEE通常为正值（CO2净释放），该时段3年高寒灌丛日平均净生态系统CO2交换分别为1.91 g CO2 m-2 d-1、1.90 g CO2 m-2 d-1和2.19 g CO2 m-2 d-1。2003~2004年高寒灌丛生态系统CO2净释放维持天数分别为249 d、 254 d和264 d，2003年净释放维持天数最少，而净吸收维持天数2005年最少（101d）。2003、2004和2005年全年日平均CO2净吸收分别为0.611 g CO2 m-2 d-1、0.759 g CO2 m-2 d-1和0.167 g CO2 m-2 d-1。 3、就季节差异而言，2003、2004和2005年整个生长季节高寒灌丛平均CO2日净生态系统交换分别为-3.99 g CO2 m-2 d-1、-4.59 g CO2 m-2 d-1、-3.27 g CO2 m-2 d-1。7、8月生长季节CO2净吸收的最高，2003、2004、2005年7月和8月份高寒灌丛生态系统CO2净吸收分别为222 g CO2 m-2 和224 g CO2 m-2、355 g CO2 m-2和216 g CO2 m-2、263 g CO2 m-2和186 g CO2 m-2。在相对短暂的生长季节海北地区高寒灌丛生态系统表现出显著的大气CO2净吸收能力，2003、2004和2005年生长季节高寒灌丛生态系统CO2净吸收量分别为610 g CO2 m-2、701 g CO2 m-2和500 g CO2 m-2。相对于温度等环境因子，高寒灌丛生态系统生长季白昼NEE小时变化规律更受光合有效辐射变化的影响。 4、2003～2005年非生长季节日平均NEE分别为1.83 g CO2 m-2、2.01 g CO2 m-2和2.07 g CO2 m-2。4月和10月是非生长季节CO2净释放的最高月份，2003、2004和2005年全月净释放量为105 g CO2 m-2和77 g CO2 m-2、105 g CO2 m-2和117 g CO2 m-2及105 g CO2 m-2和138 g CO2 m-2，2003～2005年整个非生长季CO2净释放分别为CO2为388 g CO2 m-2、425 g CO2 m-2和439 g CO2 m-2。非生长季节海北地区高寒灌丛生态系统NEE小时变化与5 cm土壤温度存在极显著的正相关关联，表明在非生长季节土壤温度是影响青藏高原高寒灌丛生态系统NEE的重要环境因子。 5、从生态系统CO2源/汇特征来看，海北地区高寒灌丛生态系统2003、2004和2005年全年净CO2固定总量分别为223 g CO2 m-2 a-1、277 g CO2 m-2 a-1和61 g CO2 m-2 a-1，3年平均CO2值为187 g CO2 m-2 a-1。在为期3年的研究时段海北地区高寒灌丛生态系统表现为弱的大气二氧化碳的汇。 6、高寒灌丛群落表观光合量子产额（a）和表观最大光合速率（Pmax）受叶面积指数的影响。在6～9月份期间，由于LAI的不同，a和Pmax值差异明显，7、8月份较高而6月和9月明显较低。海北地区高寒灌丛生态系统a和Pmax值高于西藏当雄地区高寒草甸生态系统，但低于平原地区相关生态系统。 维持天数2005年最少（101d）。2003、2004和2005年全年日平均CO2净吸收分别为0.611 g CO2 m-2 d-1、0.759 g CO2 m-2 d-1和0.167 g CO2 m-2 d-1。 3、就季节差异而言，2003、2004和2005年整个生长季节高寒灌丛平均CO2日净生态系统交换分别为-3.99 g CO2 m-2 d-1、-4.59 g CO2 m-2 d-1、-3.27 g CO2 m-2 d-1。7、8月生长季节CO2净吸收的最高，2003、2004、2005年7月和8月份高寒灌丛生态系统CO2净吸收分别为222 g CO2 m-2 和224 g CO2 m-2、355 g CO2 m-2和216 g CO2 m-2、263 g CO2 m-2和186 g CO2 m-2。在相对短暂的生长季节海北地区高寒灌丛生态系统表现出显著的大气CO2净吸收能力，2003、2004和2005年生长季节高寒灌丛生态系统CO2净吸收量分别为610 g CO2 m-2、701 g CO2 m-2和500 g CO2 m-2。相对于温度等环境因子，高寒灌丛生态系统生长季白昼NEE小时变化规律更受光合有效辐射变化的影响。 4、2003～2005年非生长季节日平均NEE分别为1.83 g CO2 m-2、2.01 g CO2 m-2和2.07 g CO2 m-2。4月和10月是非生长季节CO2净释放的最高月份，2003、2004和2005年全月净释放量为105 g CO2 m-2和77 g CO2 m-2、105 g CO2 m-2和117 g CO2 m-2及105 g CO2 m-2和138 g CO2 m-2，2003～2005年整个非生长季CO2净释放分别为CO2为388 g CO2 m-2、425 g CO2 m-2和439 g CO2 m-2。非生长季节海北地区高寒灌丛生态系统NEE小时变化与5 cm土壤温度存在极显著的正相关关联，表明在非生长季节土壤温度是影响青藏高原高寒灌丛生态系统NEE的重要环境因子。 5、从生态系统CO2源/汇特征来看，海北地区高寒灌丛生态系统2003、2004和2005年全年净CO2固定总量分别为223 g CO2 m-2 a-1、277 g CO2 m-2 a-1和61 g CO2 m-2 a-1，3年平均CO2值为187 g CO2 m-2 a-1。在为期3年的研究时段海北地区高寒灌丛生态系统表现为弱的大气二氧化碳的汇。 6、高寒灌丛群落表观光合量子产额（a）和表观最大光合速率（Pmax）受叶面积指数的影响。在6～9月份期间，由于LAI的不同，a和Pmax值差异明显，7、8月份较高而6月和9月明显较低。海北地区高寒灌丛生态系统a和Pmax值高于西藏当雄地区高寒草甸生态系统，但低于平原地区相关生态系统。

移动机器人基于激光测距和单目视觉的室内同时定位和地图构建

该文研究了部分结构化室内环境中自主移动机器人同时定位和地图构建问题.基于激光和视觉传感器模型的不同,加权最小二乘拟合方法和非局部最大抑制算法被分别用于提取二维水平环境特征和垂直物体边缘.为完成移动机器人在缺少先验地图支持的室内环境中的自主导航任务,该文提出了同时进行扩展卡尔曼滤波定位和构建具有不确定性描述的二维几何地图的具体方法.通过对应用于SmartROB-2移动机器人平台所获得的实验结果和数据的分析讨论,论证了所提出方法的有效性和实用性.

一个简单的机器人视觉引导装配实验系统

本文对视觉控制下的一个简单实验室装配系统作了介绍,讨论了系统组成、机器人控制、二维图象特征的提取、对物体自动识别、定位定向、系统标定、实现垒积木装配工作.本实验系统用的是我所研制的国内第一台示教再现机器人.

改进叠前深度偏移的计算效率＆叠前高分辨率成像

An high-resolution prestack imaging technique of seismic data is developed in this thesis. By using this technique, the reflected coefficients of sheet sands can be gained in order to understand and identify thin oil reservoirs. One-way wave equation based migration methods can more accurately model seismic wave propagation effect such as multi-arrivals and obtain almost correct reflected energy in the presence of complex inhomogeneous media, and therefore, achieve more superiorities in imaging complex structure. So it is a good choice to apply the proposed high-resolution imaging to the presatck depth migration gathers. But one of the main shorting of one-way wave equation based migration methods is the low computational efficiency, thus the improvement on computational efficiency is first carried out. The method to improve the computational efficiency of prestack depth migration is first presented in this thesis, that is frequency-dependent varying-step depth exploration scheme plus a table-driven, one-point wavefield interpolation technology for wave equation based migration methods; The frequency-dependent varying-step depth exploration scheme reduces the computational cost of wavefield depth extrapolation, and the a table-driven, one-point wavefield interpolation technology reconstructs the extrapolated wavefield with an equal, desired vertical step with high computational efficiency. The proposed varying-step depth extrapolation plus one-point interpolation scheme results in 2/3 reduction in computational cost when compared to the equal-step depth extrapolation of wavefield, but gives the almost same imaging. The frequency-dependent varying-step depth exploration scheme is presented in theory by using the optimum split-step Fourier. But the proposed scheme can also be used by other wave equation based migration methods of the frequency domain. The proposed method is demonstrated by using impulse response, 2-D Marmousi dataset, 3-D salt dataset and the 3-D field dataset. A method of high-resolution prestack imaging is presented in the 2nd part of this thesis. The seismic interference method to solve the relative reflected coefficients is presented. The high-resolution imaging is obtained by introducing a sparseness- constrained least-square inversion into the reflected coefficient imaging. Gaussian regularization is first imposed and a smoothed solution is obtained by solving equation derived from the least-square inversion. Then the Cauchy regularization is introducing to the least-square inversion , the sparse solution of relative reflected coefficients can be obtained, that is high-resolution solution. The proposed scheme can be used together with other prestack imaging if the higher resolution is needed in a target zone. The seismic interference method in theory and the solution to sparseness-constrained least-square inversion are presented. The proposed method is demonstrated by synthetic examples and filed data.

基于随机介质理论的多尺度储层非均质性定量描述

Application of long-term exploration for oil and gas shows that the reservoir technology of prediction is one of the most valuable methods. Quantitative analysis of reservoir complexity is also a key technology of reservoir prediction. The current reservoir technologies of prediction are based on the linear assumption of various physical relationships. Therefore, these technologies cannot handle complex reservoirs with thin sands, high heterogeneities in lithological composition and strong varieties in petrophysical properties. Based on the above-mentioned complex reservoir, this paper conducts a series of researches. Both the comprehending and the quantitative analysis of reservoir heterogeneities have been implemented using statistical and non-linear theories of geophysics. At the beginning, the research of random media theories about reservoir heterogeneities was researched in this thesis. One-dimensional (1-D) and two-dimensional (2-D) random medium models were constructed. The autocorrelation lengths of random medium described the mean scale of heterogeneous anomaly in horizontal and deep directions, respectively. The characteristic of random medium models were analyzed. We also studied the corresponding relationship between the reservoir heterogeneities and autocorrelation lengths. Because heterogeneity of reservoir has fractal nature, we described heterogeneity of reservoir by fractal theory based on analyzing of the one-dimensional (1-D) and two-dimensional (2-D) random medium models. We simulated two-dimensional (2-D) random fluctuation medium in different parameters. From the simulated results, we can know that the main features of the two-dimensional (2-D) random medium mode. With autocorrelation lengths becoming larger, scales of heterogeneous geologic bodies in models became bigger. In addition, with the autocorrelation lengths becoming very larger, the layer characteristic of the models is very obvious. It would be difficult to identify sandstone such as gritstone, clay, dense sandstone and gas sandstone and so on in the reservoir with traditional impedance inversion. According to the obvious difference between different lithologic and petrophysical impedance, we studied multi-scale reservoir heterogeneities and developed new technologies. The distribution features of reservoir lithological and petrophysical heterogeneities along vertical and transverse directions were described quantitatively using multi-scale power spectrum and heterogeneity spectrum methods in this paper. Power spectrum (P spectrum) describes the manner of the vertical distribution of reservoir lithologic and petrophysical parameters and the large-scale and small-scale heterogeneities along vertical direction. Heterogeneity spectrum (H spectrum) describes the structure of the reservoir lithologic and petrophysical parameters mainly, that is to say, proportional composition of each lithological and petrophysical heterogeneities are calculated in this formation. The method is more reasonable to describe the degree of transverse multi-scale heterogeneities in reservoir lithological and petrophysical parameters. Using information of sonic logs in Sulige oil field, two spectral methods have been applied to the oil field, and good analytic results have been obtained. In order to contrast the former researches, the last part is the multi-scale character analysis of reservoir based on the transmission character of wave using the wavelet transform. We discussed the method applied to demarcate sequence stratigraphy and also analyzed the reservoir interlayer heterogeneity.

黄海海域前新生代残留盆地综合地质地球物理研究

China’s annual oil import volume has been increasing in recent years, but the oil price in the international market fluctuates and poses a severe threat to China’s economic development and national security. Therefore, it is of great importance to study the gas and oil exploration of Pre-Cenozoic Residual Basins in Yellow Sea. Yellow Sea has widespread and thick Mesozoic and Paleozoic strata that contain multilayer source rock. Hence, Yellow Sea Mesozoic and Paleozoic strata have good conditions of forming Pre-Cenozoic hydrocarbon reservoirs. Pre-Cenozoic Residual Basins are usually buried deep and then transformed many times in its long evolutional history. These characteristics make it difficult to apply a single method in exploring Pre-Cenozoic Residual Basins. On the other hand, it is highly effective to solve key problems of gas and oil exploration of Pre-Cenozoic Residual Basins in Yellow Sea by using integrated geological and geophysical methods which make full use of the advantages of various exploring techniques. Based on the principle of “the region controls the local; the deep restricts the shallow,” this study focuses on Pre-Cenozoic Residual Basins in Yellow Sea to describe the structure frame of its distribution, with gravity, magnetic, seismic, drill-hole and geological data and previous research findings. In addition, the distribution characteristics of Pre-Cenozoic Residual Basins in Yellow Sea are also analyzed. This paper explores the characteristics of error between gravity forward with constant density and gravity forward with variable density through the study on 2-D and 3-D gravity forward in frequency domain. The result shows that there is a linear relationship between error and depth of 2-D geological model but there is a nonlinear relationship between error and depth of 3-D geological model. The error can be removed according to its linear characteristics or statistical nature of nonlinear characteristics. There is also error between gravity inversion with constant density and gravity inversion with variable density due to variable density and edge-effect. Since there are not noticeable rules between the error and the two causes as variable density and edge-effect, this study adopts gravity inversion with variable density and methods to eliminate the edge-effect in basement inversion to improve inversion accuracy. Based on the study on the rock physical properties and strata distribution of Yellow Sea and adjacent regions, this study finds that there is a big density contrast between Cretaceous-Jurassic strata and their substratum. The magnetic basement of south Yellow Sea is regarded as top of Archeozoic-Proterozoic early strata, and there are double magnetic basements in north Yellow Sea. Gravity and magnetic data are used to inverse the gravity basement and magnetic basement of Yellow Sea, with seismic and drill-hole data as constrains. According to data of gravity and magnetic basement distribution, the depth of Cenozoic strata and previous research findings, this paper calculates the thickness of the Mesozoic and Pre-Mesozoic Residual Basins, draws the distribution outline of Pre-Cenozoic Residual Basins in Yellow Sea, and analyzes its macro-distribution characteristics. Gravity inversion is applied on a typical geological profile in Yellow Sea to analyze the characteristics of its fractures and magnetic basements. The characteristics of Pre-Cenozoic Residual Basins distribution outline in Yellow Sea and the fractures and magnetic basements of its typical profile shown by profile inversion provides new geophysical evidence for these structure views such as “the South Yellow Sea and the North Yellow Sea belong to different structural units” and “Sino-Korea and Yangtze blocks combine along Yellow Sea East Fractured Zone in Yellow Sea”.

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

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

青藏高原东部长周期大地电磁测深研究

As a complement to conventional MT, Long-period Magnetotellurics (LMT) has been developed at low frequency for soundings of deep electric structures. Eastern Himalayan Syntaxis (EHS) and surrounding area is a key place for the study of dynamics of the uplift of Tibetan plateau. Experiments in the pioneer studies for EHS3D project showed that the study area shares an unusual low resistive crust and upper mantle. Conventional MT could not provide sufficient information about the deep structures of the study area that requested long period MT measurement to be complemented. This thesis presents the LMT studies in eastern Tibet along the EHS3D-3 Profile from Xiachayu to Yushu including data acquisition, processing, inversion and interpretation. The effective period of the measured LMT signals extend from 10s up to 30000s for the duration more than one week measurement. The resulting model shows that the LMT sounding coincides with the MT data in overlapped periods. Especially the induction arrows and tippers derived from LMT data provide more information about the base of the conductors beneath the plateau with higher resolution. Anomalous induction coefficients and 2-D model suggest extensive conductive bodies beneath Lhasa block and Qiangtang terrain which would be a possible evidence for partial melt and fluids at depth.

青藏高原东部电性结构特征及其构造意义

Eastern Himalayan Syntaxis (EHS) and its surroundings (eastern margin of Tibet) is one of the most complicated tectonic areas in the world. As the exhaust opening of the balanced materials of the Tibetan Plateau during the collision of Indan and Eurasian plates, the deep structure beneath EHS surrounding region is referred to as the key to the study of the dynamics of the plateau. EHS3D project, sponsored by NSFC, has been proposed to explore the deep electric features of the area. During the first stage of EHS3D(2006-2008), MT+LMT measurements have been conducted along two lines from Chayu to Qingshuihe (EHS3D-3) and Chayu to Ruoergai (EHS3D-2). This paper will discuss the MT models of EHS3D-3 line. By the data procrssing, including distortion analysis, Robust estimation and strike decomposition, rotated apparent resitivities and phases have been obtained for each station. Then conventional 2-D inversion algorithms (NLCG and RRI) were employed to produce 2-D models. The final preferred 2-D model suggests that the upper crust consists of resistive blocks while in mid-lower crust there are two extensive conductive bodies beneath Lhasa block and Qiangtang terrain respectively. Jinshajiang suture is a gradient belt and Bangong-Nujiang suture appear a conductive belt dipping to the north. . We concluded that the formation of the two conductive bodies attributed to the partial melt and fluids in the lower crust. The regional electric strike derived from decomposition analysis indicates that the crust and upper mantle move in different manners. The upper crust moves like slips of rigid blocks along major slip faults while the lower crust creeps as a flow in the conductive channels.

中高层大气对低层大气声重波响应的分析

Acoustic Gravity waves (AGW) play an important role in balancing the atmospheric energy and momentum budget. Propagation of gravity wave in the atmosphere is one of the important factors of changing middle and upper atmosphere and ionosphere. The purpose of this dissertation is to study the propagation of gravity wave in a compression atmosphere whit means of numerical simulation and to analyze the response of middle and upper atmosphere to pulse disturbance from lower atmosphere. This work begins with the establishment of 2-D fully nonlinear compressible atmospheric dynamic model in polar coordinate, which is used ton numerically study gravity wave propagation. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. We also simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model and analyze the data we obtained by using Fourier Transform (FT), Short-time Fourier Transform (STFT) and Empirical Mode Decomposition (EMD) method which is an important part of Hilbert-Huang Transform (HHT). The research content is summarized in the following: 1. By using a two-dimensional full-implicit-continuous-Eulerian (FICE) scheme and taking the atmospheric basic motion equations as the governing equations, a numerical model for nonlinear propagation of acoustic gravity wave disturbance in two-dimensional polar coordinates is solved. 2. Then the propagation characteristics of acoustic gravity wave packets are investigated and discussed. Results of numerical simulation show that the acoustic gravity wave packets propagate steadily upward and keep its shape well after several periods. 3. We simulate the response of middle and upper atmosphere to pulse disturbance of lower atmosphere in background winds or without background winds by using this model, and obtain the distribution of a certain physical quantity in time and space from earth’s surface to 300km above. The results reveal that the response of ionosphere occurs at a large horizontal distance from the source and the disturbance becomes greater with increasing of height. The situation when the direction of the background wind is opposite to or the same as the direction of disturbed velocity of gravity-wave is studied. The results show that gravity wave propagating against winds is easier than those propagating along winds and the background wind can accelerate gravity wave propagation. Just upon the source, an acoustic wave component with period of 6 min can be found. These images of simulation are similar to observations of the total electron content (TEC) disturbances caused by the great Sumatra-Andaman earthquake on December 26 in 2004. 4. Using the EMD method the disturbed velocity data of a certain physical quantity in time and space can be decomposed into a series of intrinsic mode function (IMF) and a trend mode respectively. The results of EMD reveal impact of the gravity wave frequency under the background winds.

离散裂缝诱导各向异性粘弹性地震波模拟与特征分析

The real media always attenuate and distort seismic waves as they propagate in the earth. This behavior can be modeled with a viscoelastic and anisotropic wave equation. The real media can be described as fractured media. In this thesis, we present a high-order staggered grid finite-difference scheme for 2-D viscoelastic wave propagation in a medium containing a large number of small finite length fractures. We use the effective medium approach to compute the anisotropic parameters in each grid cell. By comparing our synthetic seismogram by staggered-grid finite-difference with that by complex-ray parameter ray tracing method, we conclude that the high-order staggered-grid finite-difference technique can effectively used to simulate seismic propagation in viscoelastic-anisotropic media. Synthetic seismograms demonstrate that strong attenuation and significant frequency dispersion due to viscosity are important factors of reducing amplitude and delaying arrival time varying with incidence angle or offset. On the other hand, the amount of scattered energy not only provides an indicator of orientation of fracture sets, but can also provide information about the fracture spacing. Analysis of synthetic seismograms from dry- and fluid-filled fractures indicates that dry-filled fractures show more significant scattering on seismic wavefields than fluid-filled ones, and offset-variations in P-wave amplitude are observable. We also analyze seismic response of an anticlinal trap model that includes a gas-filled fractured reservoir with high attenuation, which attenuates and distorts the so-called bright spot.

波前重建法射线追踪及走时表压缩

As a fast and effective method for approximate calculation of seismic numerical simulation, ray tracing method, which has important theory and practical application value, in terms of seismic theory and seismic simulation, inversion, migration, imaging, simplified from seismic theory according to geometric seismic, means that the main energy of seismic wave field propagates along ray paths in condition of high-frequency asymptotic approximation. Calculation of ray paths and traveltimes is one of key steps in seismic simulation, inversion, migration, and imaging. Integrated triangular grids layout on wavefront with wavefront reconstruction ray tracing method, the thesis puts forward wavefront reconstruction ray tracing method based on triangular grids layout on wavefront, achieves accurate and fast calculation of ray paths and traveltimes. This method has stable and reasonable ray distribution, and overcomes problems caused by shadows in conventional ray tracing methods. The application of triangular grids layout on wavefront, keeps all the triangular grids stable, and makes the division of grids and interpolation of a new ray convenient. This technology reduces grids and memory, and then improves calculation efficiency. It enhances calculation accuracy by accurate and effective description and division on wavefront. Ray tracing traveltime table, which shares the character of 2-D or 3-D scatter data, has great amount of data points in process of seismic simulation, inversion, migration, and imaging. Therefore the traveltime table file will be frequently read, and the calculation efficiency is very low. Due to these reasons, reasonable traveltime table compression will be very necessary. This thesis proposes surface fitting and scattered data compression with B-spline function method, applies to 2-D and 3-D traveltime table compression. In order to compress 2-D (3-D) traveltime table, first we need construct a smallest rectangular (cuboidal) region with regular grids to cover all the traveltime data points, through the coordinate range of them in 2-D surface (3-D space). Then the value of finite regular grids, which are stored in memory, can be calculated using least square method. The traveltime table can be decompressed when necessary, according to liner interpolation method of 2-D (3-D) B-spline function. In the above calculation, the coefficient matrix is stored using sparse method and the liner system equations are solved using LU decomposition based on the multi-frontal method according to the sparse character of the least square method matrix. This method is practiced successfully in several models, and the cubic B-spline function can be the best basal function for surface fitting. It make the construction surface smooth, has stable and effective compression with high approximate accuracy using regular grids. In this way, through constructing reasonable regular grids to insure the calculation efficiency and accuracy of compression and surface fitting, we achieved the aim of traveltime table compression. This greatly improves calculation efficiency in process of seismic simulation, inversion, migration, and imaging.