Origin of the Qinghai-Tibetan Plateau endemic Milula (Liliaceae): further insights from karyological comparisons with Allium

Milula, a monotypic genus endemic to the Qinghai-Tibetan Plateau, was found to be nested deeply within Allium by the molecular phylogeny despite the aberrant morphology. It remains unknown what had contributed to the rapid evolution of morphology and origin of this exceptional species. In contrast to a previous report of its karyotypes with 2n = 16 = 8M+8SM (2SAT), similar to most species of Allium, a rather different karyotype, 2n = 20 = 4M +10SM+6T (2SAT), was found in examined 31 individuals from 6 populations of M. spicata distributed in the central Tibet. Karyotypes of 7 Allium species occurring in the Qinghai-Tibetan Plateau were further reported. The basic number x = 8 was confirmed for all of them and their karyotypes consist mainly of metacentric and submetacentric chromosomes with rare subterminal and terminal chromosomes. The karyotype of M. spicata is distinctly different from that of most Allium species occurring in the plateau through a complete comparison of all available species in this region and adjacent areas. However, the same chromosome number and similar karyotypic structure were found in A. fasciculatum of Sect. Bromatorrhiza, indicating a possible close relationship between them. But this similarity is contradictory to the preliminary molecular phylogenetic analysis that Milula was closely related to A. cyathophorum of Sect. Bromatorrhiza with x=8, but the other species with x=10 and 11 in this section were clearly placed in the other clade. We therefore suggested that the paralleling evolution from x=8 to x=9, 10 and 11 with increasing asymmetry of karyotype possibly due to the chromosomal Robertsonian translocation might occur separately in the two recognized phylogenetic lineages of Allium. In addition to aneuploidy and following change of the chromosomal structures, the habitat isolation due to the recent uplift of the Qinghai-Tibetan Plateau and the Quaternary climatic oscillation, plays a greater role in origin of Milula and other endemic species (genera) with aberrant morphology from their progenitors.

Uniformity of karyotypes in Ligularia (Asteraceae : Senecioneae), a highly diversified genus of the eastern Qinghai-Tibet Plateau highlands and adjacent areas

Ligularia, a highly diversified genus in the eastern Qinghai-Tibet Plateau and adjacent areas, was chosen as a suitable subject in which to study speciation patterns in this 'hot spot' area at the chromosomal level. Chromosome numbers and karyotypes were studied in 23 populations of 14 species, most of which are endemic to this area. The basic number x = 29 was confirmed for all species. Ligularia virgaurea was found to have diploid and triploid cytotypes, 2n = 58 and 87. Other species are only diploid, with 2n = 58. The karyotypes of all populations within any species, and all species spanning most sections and covering most of the morphological range in Ligularia, are very similar to each other, belonging to type 2A according to Stebbin's classification. This karyotype was also found in its close allies, e.g. Cremanthodium, Ligulariopsis, Parasenecio, and Sinacalia. Aneuploid reduction of chromosome number from 2n = 60 to 58 and karyotypic variation was found in Ligularia and its allies. Such a chromosomal pattern with few polyploids infers that variation of karyotype structure at the diploid level seems to be the predominant feature of chromosomal evolution in this group and sympatric speciation via hybridization and polyploidization has played a minor role in its species diversity. (C) 2004 The Linnean Society of London

储层评价目标驱动的地震叠前叠后综合反演方法研究

The exploration and development of natural gas in the north of Ordos basin have been one important part in China’s energy stratagem. Reservoir in upper Palaeozoic group is of lithological trap and its prediction is a crux in a series of works. Based on foregoing seismic reservoir prediction, seismic data are re-processed with some optical methods and pre-stack information is used in corresponding inversions. Through the application of diverse methods, a series of techniques for reservoir prediction come into being. Several results are achieved as flowing: 1. A set of log processing and interpretation methods is developed. Porosity, permeability and gas saturation models are rebuilt. 2. Based on the petro-physics analysis of reservoirs in upper Palaeozoic group, the equations about lithology, property, hydrocarbon and elastic parameters are established. 3. Forward modeling based on elastic wave theory is first applied in the study area and increases the resolution of modeling results. 4. A series of techniques such as pre-stack time migration and others are combined to improve the data quality. 5. Pre-stack seismic inversion is first employed in the north of Ordos Basin and brings the results of EI, P-impedance, S-impedance and other elastic parameters. 6. In post-stack inversion, logs indicating reservoir parameters are rebuilt and boost the resolution of lithology inversion. 7. Amplitude, coherence, frequency-discomposed amplitude, waveform and other sensitive attributes are extracted to describe sands distribution. Seismic modes standing for sands of P1x3, P1x2 are established. 8. Among 9 proposed wells, 8 wells encountered sands and became production wells. The output of DK13 amounts to 510,000 m3 per day. Keywords：the north of Ordos Basin, reservoir prediction, pre-stack inversion, post-stack inversion, seismic attributes．

地震数据波动方程反演方法研究

Geophysical inversion is a theory that transforms the observation data into corresponding geophysical models. The goal of seismic inversion is not only wave velocity models, but also the fine structures and dynamic process of interior of the earth, expanding to more parameters such as density, aeolotropism, viscosity and so on. As is known to all, Inversion theory is divided to linear and non-linear inversion theories. In rencent 40 years linear inversion theory has formed into a complete and systematic theory and found extensive applications in practice. While there are still many urgent problems to be solved in non-linear inversion theory and practice. Based on wave equation, this dissertation has been mainly involved in the theoretical research of several non-linear inversion methods: waveform inversion, traveltime inversion and the joint inversion about two methods. The objective of gradient waveform inversion is to find a geologic model, thus synthetic seismograms generated by this geologic model are best fitted to observed seismograms. Contrasting with other inverse methods, waveform inversion uses all characteristics of waveform and has high resolution capacity. But waveform inversion is an interface by interface method. An artificial parameter limit should be provided in each inversion iteration. In addition, waveform information will tend to get stuck in local minima if the starting model is too far from the actual model. Based on velocity scanning in traditional seismic data processing, a layer-by-layer waveform inversion method is developed in this dissertation to deal with weaknesses of waveform inversion. Wave equation is used to calculate the traveltime and derivative (perturbation of traveltime with respect to velocity) in wave-equation traveltime inversion (WT). Unlike traditional ray-based travetime inversion, WT has many advantages. No ray tracing or traveltime picking and no high frequency assumption is necessary and good result can be got while starting model is far from real model. But, comparing with waveform inversion, WT has low resolution. Waveform inversion and WT have complementary advantages and similar algorithm, which proves that the joint inversion is a better inversion method. And another key point which this dissertation emphasizes is how to give fullest play to their complementary advantages on the premise of no increase of storage spaces and amount of calculation. Numerical tests are implemented to prove the feasibility of inversion methods mentioned above in this dissertation. Especially for gradient waveform inversion, field data are inversed. This field data are acquired by our group in Wali park and Shunyi district. Real data processing shows there are many problems for waveform inversion to deal with real data. The matching of synthetic seismograms with observed seismograms and noise cancellation are two primary problems. In conclusion, on the foundation of the former experiences, this dissertation has implemented waveform inversions on the basis of acoustic wave equation and elastic wave equation, traveltime inversion on the basis of acoustic wave equation and traditional combined waveform traveltime inversion. Besides the traditional analysis of inversion theory, there are two innovations: layer by layer inversion of seimic reflection data inversion and rapid method for acoustic wave-equation joint inversion.

高分辨率非线性岩性物性联合反演方法和应用

In the prediction of complex reservoir with high heterogeneities in lithologic and petrophysical properties, because of inexact data (e.g., information-overlapping, information-incomplete, and noise-contaminated) and ambiguous physical relationship, inversion results suffer from non-uniqueness, instability and uncertainty. Thus, the reservoir prediction technologies based on the linear assumptions are unsuited for these complex areas. Based on the limitations of conventional technologies, the thesis conducts a series of researches on various kernel problems such as inversions from band-limited seismic data, inversion resolution, inversion stability, and ambiguous physical relationship. The thesis combines deterministic, statistical and nonlinear theories of geophysics, and integrates geological information, rock physics, well data and seismic data to predict lithologic and petrophysical parameters. The joint inversion technology is suited for the areas with complex depositional environment and complex rock-physical relationship. Combining nonlinear multistage Robinson seismic convolution model with unconventional Caianiello neural network, the thesis implements the unification of the deterministic and statistical inversion. Through Robinson seismic convolution model and nonlinear self-affine transform, the deterministic inversion is implemented by establishing a deterministic relationship between seismic impedance and seismic responses. So, this can ensure inversion reliability. Furthermore, through multistage seismic wavelet (MSW)/seismic inverse wavelet (MSIW) and Caianiello neural network, the statistical inversion is implemented by establishing a statistical relationship between seismic impedance and seismic responses. Thus, this can ensure the anti-noise ability. In this thesis, direct and indirect inversion modes are alternately used to estimate and revise the impedance value. Direct inversion result is used as the initial value of indirect inversion and finally high-resolution impedance profile is achieved by indirect inversion. This largely enhances inversion precision. In the thesis, a nonlinear rock physics convolution model is adopted to establish a relationship between impedance and porosity/clay-content. Through multistage decomposition and bidirectional edge wavelet detection, it can depict more complex rock physical relationship. Moreover, it uses the Caianiello neural network to implement the combination of deterministic inversion, statistical inversion and nonlinear theory. Last, by combined applications of direct inversion based on vertical edge detection wavelet and indirect inversion based on lateral edge detection wavelet, it implements the integrative application of geological information, well data and seismic impedance for estimation of high-resolution petrophysical parameters (porosity/clay-content). These inversion results can be used to reservoir prediction and characterization. Multi-well constrains and separate-frequency inversion modes are adopted in the thesis. The analyses of these sections of lithologic and petrophysical properties show that the low-frequency sections reflect the macro structure of the strata, while the middle/high-frequency sections reflect the detailed structure of the strata. Therefore, the high-resolution sections can be used to recognize the boundary of sand body and to predict the hydrocarbon zones.

煤系地层三维三分量地震数据解释方法研究

Conventional 3D seismic exploration cannot meet the demand of high yield and high efficiency safe production in coal mine any more. Now it is urgent to improve the discovery degree of coal mine geological structures for coal production in China. Based on 3D3C seismic exploration data, multi-component seismic information is fully excavated. First systematic research on 3D3C seismic data interpretation of coal measure strata is carried out. Firstly, by analyzing the coal measure strata, the seismic-geologic model of coal measure strata is built. Shear wave logging is built by using regression analysis. Horizon calibration methods of PP-wave and PS-wave are studied and the multi-wave data are used together to interpret small faults. Using main amplitude analysis technology, small faults which cannot be found from PP-wave sections can be interpreted from the low frequency PS-wave sections. Thus, the purpose to applying PS-wave data to fine structure assistant interpretation is achieved. Secondly, PP- and PS-wave post-stack well constrained inversion methods of coal measure strata are studied. Joint PP- and PS-wave post-stack inversion flow is established. More attribute parameters, which are applied in fine lithology interpretation of coal measure strata, are obtained from combinations of the inversion results. Exploring the relation between rock with negative Poisson’s ratio and anisotropy, fracture development in coal seam are predicted. Petrophysical features of coal measure strata are studied, and the relations between elastic parameters and lithology, fluid and physical properties are established. Inversions of the physical parameters such as porosity, permeability and water saturation, which reflect lithology and fluid property, are obtained. Finally, the approaches of shear wave splitting and Thomsen parameters inversion, which provide new ideas for seismic anisotropy interpretation of coal measure strata, are studied to predict fracture development. The results of practical application indicate that the methods in this paper have good feasibility and applicability. They have positive significance for high yield and high efficiency safe production in coal mine.

地震资料叠前储层预测综合研究--以辽中凹陷东二下段浊积扇气藏为例

Reservoir prediction techniques from prestack seismic are among the most important ones for exploration of lithologic hydrocarbon reservoir. In this paper, we set the turbidite fan sandstone reservoir in Liao-Zhong depress as our researching target, and aims to solve the apllication difficulties on pre-stack inversion in the area, where the drilling data is scarce and the reservoir is lateral varied. Meanwhile, AVO analysis and pre-stack inversion for gas-bearing detection is systematically researched. The seismic reflection characters of gas-bearing sandstone in turbidite fan with different fluid content are defined, after analyzing results from AVO seismic simulation and porous fluid replacement of real log data, and under the guides of the seismic characters from classical gas-bearing sandstone reservoir and numerical simulation for complicate gas-bearing sandstone. It is confirmed that detecting gas-bearing sandstone in turbidite fan via AVO technologies is feasible. In terms of AVO analysis, two AVO characters, fluid detection factor and product of intercept and gradient, can effectively identify top and bottom boundaries and lateral range of tuibidite gas sand by comparing real drilling data. Cross-plotting of near and far angle stack data could avoid the correlation existing in P-G analysis. After comparing the acoustic impedance inversions with routine stacked data and AVO intercept, impedance derived from AVO intercept attribute could reduce the acoustic impedance estimating error which is caused by AVO. On the aspect of elastic impedance inversion, the AVO information in the pre-stack gathers is properly reserved by creating partial angle stack data. By the far angle elastic impedance alone, the gas sand, with abnormally low range of values, can be identified from the background rocks. The boundary of gas sand can also be clearly determined by cross-plotting of near and far angle elastic impedances. The accuracy of far angle elastic impedance is very sensitive to the parameter K, and by taking the statistical average of Vp/Vs on the targeted section in key wells, the accuracy of low frequency trends is gurranteed; the intensive absorsion within the area of the gas sand, which tends to push the spectral of seismic data to the lower end, will cause errors on the inversion result of elastic impedance. The solution is to confine the inversion on the interested area by improving the wavelet. On the aspect of prestack AVA simultaneous inversion, the constraint of local rock-physical trends between velocities of P-wave、S-wave and density successfully removes the instability of inversion, thus improves the precision of the resulting elastic parameters. Plenty of data on rock properties are derived via AVO analysis and prestack seismic data inversion. Based on them, the fluid anomaly is analysized and lithological interpretation are conducted. The distribution of gas sand can be consistently determined via various of ways, such as cross-plotting of P and G attributes, near and far partial angle stack data, near and far angle elastic impedances, λρ and Vp/Vs, etc. The shear modulo and density are also reliable enough to be used for lithological interpretation. We successfully applied the AVO analysis and pre-stack inversion techniques to gas detecting for turbidite fan sand reservoir in Liao-Zhong depression.

低感应数条件下大地电导率测量方法研究

A frequency domain electromagnetic (conductivity) method for near surface soundings at low frequencies is discussed in this thesis. Its elementary principle is to detect the conductivity of the earth by the secondary magnetic fields induced by a current dipole on the earth. According to the EM induction theory, a coil with alternating current on the earth will generate a magnetic field in whole space which is referred to as the primary field Hp. The primary field would induce secondary currents in the earth which go down to depth like a batch of smoking rings. These currents further produce secondary magnetic field Hs .The primary and secondary magnetic fields are collected together by a receiver coil. Generally speaking，the secondary magnetic field is a complicated function of coil spacing, transmitting frequency and earth conductivity. But at low induction numbers, the secondary field is deduced to as a simple function of frequency, spacing and conductivity. Especially the ratio of secondary to primary field shares a linear proportion to the apparent conductivity. The earth conductivity can be interpreted by proper inversions with the apparent conductivity. The method is discussed at three steps: (1)Derivation of primary and secondary magnetic fields arising from vertical and horizontal magnetic dipoles on the earth based on the basic EM induction theory. (2)Field techniques and equipment developed for the method. (3)An interpretation technique was introduced using a cumulative and relative response function. Finally a test example is presented for examining the effectiveness of the method.

粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

接收函数方法与中国及邻区上地幔间断面的研究

The receiver function method applied in researching the discontinuities in upper mantle was systematically studied in this paper. Using the theoretical receiver functions, the characteristics of P410S and P660S phases were analyzed, and the influencing factors for detection of these phases were discussed. The stability of receiver function was studied, and a new computational method of receiver function, RFSSMS (Receiver Function of Stack and Smooth of Multi seismic-records at a Single station), was put forward. We built initial reference velocity model for the media beneath each of 18 seismic stations respectively; then estimated the buried depths of 410-km and 660-km discontinuities(simply marked as '410' and '660') under the stations by using the arrive time differences of P410S and P660S with P. We developed a new receiver function inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion), to obtain the whole crust and upper mantle velocity structure and the depths of discontinuities beneath a station. The major works and results could be summarized as follows: (1) By analysis of the theoretical receiver functions with different velocity models and different ray parameters, we obtain the knowledge: The amplitudes of P410S and P660S phases are decreasing with the increasing of epicentral distance A , and the arrival time differences of these phases with P are shorter as A is longer. The multiple refracted and/or reflected waves yielded on Moho and the discontinuities in the crust interfere the identification of P410S. If existing LVZ under the lithosphere, some multiple waves caused by LVZ will interfere the identification of P410S. The multiple waves produced by discontinuity lied near 120km depth will mix with P410s phase in some range of epicentral distance; and the multiple waves concerned with the discontinuity lied near 210km depth will interfere the identification of P660S. The epicentral distance for P4i0s identification is limited, the upper limit is 80° . The identification of P660S is not restricted by the epicenter distance obviously. The identification of P410S and P6gos in the theoretical receiver functions is interfered weakly from the seismic wave attenuation caused by the media absorption if the Q value in a reasonable range. (2) The stability of receiver function was studied by using synthetic seismograms with different kind of noise. The results show that on the condition of high signal-noise-ratio of seismic records, the high frequency background noise and the low frequency microseism noise do not influence the calculating result of receiver function. But the media "scattering noise" influence the stability of receiver function. When the scattering effect reach some level, the identification of P4iOs and P66os is difficult in single receiver function which is yielded from only one seismic record. We provided a new method to calculate receiver function, that is, with a group of earthquake records, stacking the R and Z components respectively in the frequency domain, and weighted smooth the stacked Z component, then compute the complex spectrum ratio of R to Z. This method can improve the stability of receiver function and protrude the P4i0s and P66os in the receiver function curves. (3) 263 receiver functions were provided from 1364 three component broadband seismograms recorded at 18 stations in China and adjacent areas for the tele-earthquakes. The observed arrival time differences of P410S and P660S with P were obtained in these receiver functions. The initial velocity model for every station was built according to the prior research results. The buried depths of '410' and '660' under a station were acquired by the way of adjusting the depths of these two discontinuities in the initial velocity model until the theoretical arrival time differences of P410S and P660S with P well conformed to the observed. The results show an obvious lateral heterogeneity of buried depths of ' 410' and (660' . The depth of '410' is shallower beneath BJI, XAN, LZH and ENH, but deeper under QIZ and CHTO, and the average is 403km . The average depth of '660' is 663km, deeper under MDJ and MAJO, but shallower under QIZ and HYB. (4) For inversing the whole crust and upper mantle velocity structure, a new inversion method -PGARFI (Peeling-Genetic Algorithm of Receiver Function Inversion) has- been developed here. The media beneath a station is divided into segments, then the velocity structure is inversed from receiver function from surface to deep successively. Using PGARFI, the multi reflection / refraction phases of shallower discontinuities are isolated from the first order refraction transform phase of deep discontinuity. The genetic algorithm with floating-point coding was used hi the inversion of every segment, and arithmetical crossover and non-uniform mutation technologies were employed in the genetic optimization. 10 independent inversions are completed for every segment, and 50 most excellent velocity models are selected according to the priority of fitness from all models produced in the inversion process. The final velocity structure of every segment is obtained from the weighted average of these 50 models. Before inversion, a wide range of velocity variation with depth and depth range of the main discontinuities are given according to priori knowledge. PGARFI was verified with numerical test and applied in the inversion of the velocity structure beneath HIA station down to 700km depth.

中强地震震源的破裂过程

With the large developments of the seismic sources theory, computing technologies and survey instruments, we can model and rebuild the rupture process of earthquakes more realistically. On which earthquake sources' properties and tectonic activities law are realized more clearly. The researches in this domain have been done in this paper as follows. Based on the generalized ray method, expressions for displacement on the surface of a half-space due to an arbitrary oriented shear and tensile dislocation are also obtained. Kinematically, fault-normal motion is equivalent to tensile faulting. There is some evidence that such motion occurs in many earthquakes. The expressions for static displacements on the surface of a layered half-space due to static point moment tensor source are given in terms of the generalized reflection and transmission coefficient matrix method. The validity and precision of the new method is illustrated by comparing the consistency of our results with the analytical solution given by Okada's code employing same point source and homogenous half-space model. The computed vertical ground displacement using the moment tensor solution of the Lanchang_Gengma earthquake displays considerable difference with that of a double couple component .The effect of a soft layer at the top of the homogenous half-space on a shallow normal-faulting earthquake is also analyzed. Our results show that more seismic information would be obtained utilizing seismic moment tensor source and layered half-space model. The rupture process of 1999 Chi-Chi, Taiwan, earthquake investigated by using co-seismic surface displacement GPS observations and far field P-wave records. In according to the tectonic analysis and distributions of aftershock, we introduce a three-segment bending fault planes into our model. Both elastic half-space models and layered-earth models to invert the distribution of co-seismic slip along the Chi-Chi earthquake rupture. The results indicate that the shear slip model can not fit horizontal and vertical co-seismic displacements together, unless we add the fault-normal motion (tensile component) in inversions. And then, the Chi Chi earthquake rupture process was obtained by inversion using the seismograms and GPS observations. Fault normal motions determined by inversion, concentrate on the shallow northern bending fault from Fengyuan to Shuangji where the surface earthquake ruptures reveal more complexity and the developed flexural slip folding structures than the other portions of the rupture zone For understanding the perturbation of surface displacements caused by near-surface complex structures, We have taken a numeric test to synthesize and inverse the surface displacements for a pop-up structure that is composed of a main thrust and a back thrust. Our result indicates that the pop-up structure, the typical shallow complex rupture that occurred in the northern bending fault zone form Fengyuan to Shuangji, can be modeled better by a thrust fault added negative tensile component than by a simple thrust fault. We interpret the negative tensile distributions, that concentrate on the shallow northern bending fault from Fengyuan to Shuangji, as a the synthetic effect including the complexities of property and geometry of rupture. The earthquake rupture process also reveal the more spatial and temporal complexities form Fenyuan to SHuangji. According to the three-components teleseismic records, the S-wave velocity structure beneath the 59 teleseismic stations of Taiwan obtained by using the transform function method and the SA techniques. The integrated results, the 3D crustal structure of Taiwan reveal that the thickest part of crustal local in the western Central Range. This conclusion is consistent with the result form the Bouguer gravity anomaly. The orogenic evolution of Taiwan is young period, and the developing foot of Central Range dose not in static balancing. The crustal of Taiwan stays in the course of dynamic equilibrium. The rupture process of 2003)2,24,Jiashi, Xinjiang earthquake was estimated by the finite fault model using far field broadband P wave records of CDSN and IRIS. The results indicate that the earthquake focal is north dip trust fault including some left-lateral strike slip. The focal mechanism of this earthquake is different form that of earthquakes occurred in 1997 and 1998, but similar to that of 1996, Artux, Xinjiang earthquake. We interpreted that the earthquake caused trust fault due to the Tarim basin pushing northward and orogeny of Tianshan mountain. In the end, give a brief of future research subject: Building the Real Time Distribute System for rupture process of Large Earthquakes Based on Internet.