伽师强震群区、西秦岭阿尼玛卿缝合带上部地壳精细结构
|
摘要
地壳结构,特别是上部地壳结构的探测与研究,对于地球学科本身的进步、防震减灾事业的长远发展及不可再生资源的进一步开发都有着十分重要的意义。上部地壳结构的知识是连接已有的大量地质学成果与较深部地壳构造的桥梁,同时上部地壳也是我们深入理解中下部地壳,乃至整个岩石圈结构的第一道屏障。上部地壳,特别是复杂构造地区上部地壳的结构,像“厚厚的云层”,阻碍了我们用地震方法这只“天文望远镜”对造成巨大地震灾害的陆内地震孕震区构造背景及地壳的演化历史和深部动力学过程的深入理解,只有对上部地壳这层在地质历史长河中永远不散的“厚厚云层”的结构认识清楚了,我们才有可能透过这“厚厚的云层”去更进一步深入认识中下地壳的孕震机理与环境,去更加科学地对地震危险区未来的震害作出评估,去为已有的大地构造理论提供更加可靠的深部地震学证据。而要做得这一点,首先要弄清楚这些地区的地壳结构,特别是上部地壳的精细结构特征。
复杂地区精细地壳结构的探测与研究主要用地震方法。大量的地震测深结果表明,地壳内部的结构,特别是造山带地区上部地壳的结构,是十分复杂的。一般说来,上部地壳的结构较中下地壳要复杂得多,而且在地表结构复杂的地区更是如此。出于对研究成本和成果产出及探测效率的考虑,宽角反射/折射和高分辨折射己成为精细地壳结构探测的最主要方法之一。上部地壳精细结构的知识则主要来自于宽角反射/折射和高分辨折射人工地震探测Pg波资料的处理结果。一般来说,探测的越细,要求的分辨就越高,所面对的研究目标就越复杂,传统地震数据处理方法的不足和缺陷就逐渐暴露出来。在复杂地壳结构地区,传统的基于程函方程有限差分解的Pg波资料处理方法就存在着以下缺陷:①有限差分反演结果只能给出速度分布,不能给出上部地壳内的可能的界面结构特征,而地壳内的界面是客观存在的,尤其是上部地壳内的折射界面较为发育,其结晶基底具有全球性,近几年在城市活断层探测中大量的高分辨折射地震资料表明,海相沉积盆地的结晶基底一般为一强折射界面;②在程函方程的有限差分解中假设波阵面处处可微,但在一些结构模型特别是在较复杂速度结构介质中波阵面可能自身相交,意味着其梯度并不存在,将导致数值计算上的不稳定性;③计算速度较慢,
The knowledge of earth crust, especially the upper crust, is very important for the progress of geoscience itself and the long term strategy of disaster reduction as well as further exploitation of non-regeneration resources. The learning of the upper crustal structures builds a bridge between the geological results and the deeper crustal tectonics, it is also a first barrier for us to deeply understand the lower crust and even the whole lithosphere. Especially, the upper crust in complex tectonic areas is just like the very thick "cloud" which prevents us to use the "astronomical telescope" of seismic methods to deeply understand the seisgenic background of interior continent earthquakes with great disaster and the crustal evolution history along with the deep geodynamic process. Only after clearly comprehending the upper crustal structure like very thick "cloud" which never disappears in geological history can we possibly further recognize the seisgenic background and its mechanics, make more scientific estimation for the future earthquake disaster and provide more reliable deep seismic evidence for geotectonic theories. For this reason, the first step should be to gain a clear idea about the crustal structures, particularly the fine structure of upper crust.Seismic methods are mainly used to investigate the fine crust structures in complex region. A large number of results from DSS (Deep Seismic Sounding) show that the interior crust structure, particularly the upper crust structure in orogenic area is very complex. Generally speaking, the interior structure of upper crust is more complex than the middle and lower crust, and even more just so in the areas with complex surface structures. Considering the balance among the study cost and achievements as well as investigation efficiency, the wide angle reflection/refraction and high resolution refraction profiles have been main methods for the fine crust structure survey. The learning of the fine upper crust structure mainly come from the results of Pg data of wide angle reflection/refraction and high resolution refraction lines. By and
large, the more fine structure to be investigated, the more high resolution to be required, and the more complex the survey objective to be faced, the inadequacy and defects of traditional processing method for Pg data are gradually exposed. So far, there are following drawbacks in traditional processing method for Pg data, which is based on the finite difference solution of eikonal equation:D The finite difference results only provide velocity distribution and can not supply the possible interface structure features in the upper crust. In fact, the interior crust interfaces do exist, especially in upper crust the refraction interface well developed and the crystalline basement exists globally. A large number of high resolution seismic refraction dada from urban active fault survey in recent several years present that in generally, the crystalline basement of the deposit basin with marine facies is a strong refracting interface;D It is presumed that everywhere on wavefront is derivable in solving the eikonal equation with finite difference method, but the wavefront may be crossed by itself in some structural model, especially in the comparatively complex media model, which means that its gradient does not exist and causes the algorithm to be unstable;□ The computing speed is not fast, which can not well suit the needs of the survey for high resolution and fine complex structures and results in very consuming cost;D The model is parameterized with square meshes, which enable the eikonal equation to be solved conveniently with the finite difference approximation,but it is not reasonable. The complicated level of real crust is different from place to place and even along the same profile, it has comparative difference among the different segments, maybe some part of which is simple and another one is very complex. Furthermore, the data amount along the different part at the same line may be different while denser observation is carried out at the complex structure part. Therefore, whether the model mesh lines are dense or sparse should correspond to both the model complex level and the data amount.In the wake of the quickly progressing of seismic observation
technique and continuously improving of survey system as well as accumulating with a great number of seismic data, especially the needs of high resolution and fine complex structure study in urban active fault investigation and basic researches for some important geoscience problems, the high resolution refraction method will be more and more widely used. For this reason, it is very necessary to improve the drawback in artificial seismic data processing discussed above and develop seismic data processing technique that can be applied to arbitrarily complex crustal structure, particularly the fine structure survey of complex upper crust.In this paper some problems mentioned above in processing seismic data have been settled and the complete process of Pg data interpretation along with some typically applied examples, such as in Jiashi strong earthquake swarm area and Animqing suit belt as well as its adjacent region, are given. Following contents are included in this thesis:(1) In this paper, the histories of seismic refraction investigation and wavefront imaging methods in artificial seismic survey were reviewed. The characteristics of wavefront tracing algorithm in complex crustal structures and its defects in processing the Pg data of artificial seismic investigation were expounded rather completely. Some drawbacks have been improved and a new method for Pg data processing, which had been used to process and to analyze the Pg data obtained from the fine crust structure investigation profiles for urban active fault survey and some key basic geoscience study projects, was proposed. Good results were achieved.( 2 ) The defects during the forward computing for the first arrival Pg wave of the artificial seismic survey in complex structure areas had been removed and the finite different wavefront imaging method based on Huygens principle for velocity structure was proposed. In comparatively complex structures, the ray trace of Pg wave is not just a single refracted path which is presumed in traditional data processing. It presents direct wave near shot points, and at the certain distances away from the shot points it maybe appears as refracted wave, head waved or
diffracted wave which depends on different upper crust and crystalline basement structure features. The improving scheme as follows: based on Hole (1992) finite difference algorithm, using the 5 operators of Lecomte et al.(2000) to simulate the advancing, expanding, evolving and propagating of Pg wave field and to achieve the tracking for the Pg wavefront, and applying backprojectkm of SIRT(Simultaneous Iterative Reconstruction Technique) to reconstruct the P wave velocity structures.(3) The RHA(Ray Hit Analysis) method for determining the interface patterns of upper crust basements was proposed and had been used to process the artificial seismic data obtained from the high resolving refraction profiles carried out in north Sichuan and south Gansu region(NSFC Project) and in Xinjian Jiashi earthquake swarm region. Based on the P velocity model determined by wavefront imaging method, the ray hit numbers for each mesh are computed and the ray hit distribution maps are drawn up, which can be used to determine possible interface patterns. This is because if the crustal interface exists, it is sure that there are denser seismic rays ,with emergent angles near to 90,close to the interface, and the larger velocity contrast between the two sides of the interface, the denser the rays to be distributed. So except the ray dense regions caused by shot point positions, the distribution features of comparatively dense ray hits, with emergent angles near to 90°, should reflect the crustal interface positions.(4) The drawbacks of HRIM(Hagedoorn Refractor Imaging Method) had been removed, and this method was applied in seismic data processing in Fuzhou, Zhangzhou, North Hainan and Jiashi strong earthquake swarm area(Xu Zhaofan, Zhang Xiankang et al.,2002,2005; Jia Shixue,Xu Zhaofan et al., 2005; Urban active fault survey reports of Geophysical Prospecting Center, CEA, 2002,2003, 2004, 2005).Hagedoorn proved that the interface is located at the point where the travel time is just equal to the sum of forward and inverse propagation times. According to this principle, the refractor can be reconstructed. Four conditions must be satisfied using Hagedoorn principle to image the
refractors, i.e. U measured reciprocal travel time data from both forward and reverse layout; D accurate reciprocal time; D known the velocity structure above the refractor; D reconstruction method for seismic wave field. Condition □ and □ may be satisfied during data collection. The velocity structures can be determined by finite difference inversion or wavefront imaging in terms of Huygens principle. An applied example was given by Aldridge et al.(1992) and Lecomte et al.(2000), respectively. In their algorithm, the velocity structures above the refractor need to be provided by other method and, it is difficult in real applications. Based on wavefront expanding in terms of Huygens principle and applying the Lecomte's (2000) five wavefront expanding operator to simulate the seismic wave propagating so as to realize forward computing of seismic wave field, and using Hole's (1992) back-project method for inversion along with visualizing of the results by Matlab algorithms, a complete refractor imaging method in terms of Hagedoorn principle can be completed(Xu Zhao fan & Zhang Xiangkang et al., 2002) and has been used to process the seismic data from high resolving refraction survey for the exploration of fine upper crust structures in urban active fault investigations in some cities and seismic active areas, and good results are achieved.(5) The data of the high resolving refraction profile in Jiashi strong earthquake swarm area were processed in detail with wavefront imaging method in terms of Huygens' principle, synthetical ray hit analysis and Hagedoorn wavefront refractor imaging principle. The fine P wave velocity structures of upper crust and the clear patterns of basement in Jiashi strong earthquake swarm area were obtained, and the possible faults were inferred. Based on the comparatively complete observation system in Jiashi strong earthquake swarm area, the relations among the inverse parameters (such as initial model, iterative times, smoothing scale etc.), wavefront imaging results and ray distribution characteristics were analyzed and discussed in detail. In Jiashi strong earthquake swarm area, the basement with the depth of about 3.0km , where the Pg velocity is about 4.6km/s, appears continuously and completely, which possibly is
the bottom boundary between the complete developed marine platform basin faces and carbonic acid rock build and under which it is light metamorphic rock in middle and late Proterozoic Group ; The interface located at the depth of about 9.0km shows comparatively large variation and its depth gradually deepens from NE to S W, the depthes of which are 9.0km at southwest end and 8.5km at northeast end. This interface is characterized by the crystalline basement where the propagating velocity of Pg is about 6.25km/s. The crystalline basement between the post numbers 38km and 65km is continuous and complete and its depth range is from 8.5km to 9.0km. From the post number 37km to 38km, its depth goes down abruptly and it is about 11.5km between the post numbers 25km and 37km, the location of which is consistent with the fault beneath Jiashi inferred from wide angle reflection/refraction profile JA1 (Zhang Xiankang et al. 2002).From the wavefront imaging results of the velocity structures we can find that the crustal velocity structure above the depth of 1 lkm is not complicated and in generally speaking, the uppermost crust is homogeneous in laterally and conspicuously layered in vertically. The shallow cover above the depth of 400m is comparatively loose weathered layer and its P velocity is about 1.65~1.8km/s; The velocity isogram above the depth of 3.0km is nearly horizon, it indicates that the structure is homogeneous in laterally. Though there is a strong velocity gradient at the depth from 400m to 3.0kmt, the velocity isograms are even, i.e. the constant velocity gradient structure appears. Its P wave velocity varies from 1.8km/s to 4.5km/s. It is the second layer, which possible is the bottom boundary of clastic-barbonic acid rock build of well developed marine platform - basin facies of Tarim in Pz; The third layer is located at the depth of about 3.0km~9.0km, the upper part of which shows different velocity structure features from the lower part. The P velocity isograms of upper part are comparatively even and its depth range is from 3.0km to 5.5km, and the P wave velocity gradient obviously decreases. The lower part is a relatively weak gradient layer and its velocity is about 5.1km/s~6.0km/s, and it possible is the light
metamorphic rock in middle and late Proterozoic Group. The fourth layer locats beneath the depth of about 9.5km, which is relatively even and its P wave velocity is about 6.3km/s. This layer may be the middle and deep metamorphic rock in late Archean or early and middle Proterozoic Group.The upper part of third layer is almost even and the average P wave velocity is about 4.8km/s. There is a strong velocity perturbation in this layer, which implyes crustal faults exists at the depths from 4.5kmto 6.0km , and three faults Fl, F2, F3 are inferred. According to the petroleum-geological data, Maigaiti fault passes through the high resolving refraction profile near the post number of 32km and Xiasuhong - Maigaiti fault has two branches, which are located beneath the post numbers of 45km and 58km, respectively. Fl fault maybe just is the Maigaiti fault and its location is slightly close to North. F2 and F3 faults should be the parts of Xiasuhong-Maigaiti fault but there is no any evidence from our results which shows that the fault Fl, F2, F3 pass through the crystalline basement or stretch to the free surface. It is suggested that each of the faults Fl, F2, F3 is not joined with the deep crustal fault in the studied area i.e. the shallow and deep structures of the crust are independent from each other.(6) The Pg data from the high resolving refraction profile which nearly perpendicularly goes through Animaqing suit zone in North Sichuan and South Gansu province were processed in detail by means of Huygens principle wavefront imaging method and the ray hit analysis technique. The fine structures of the upper crust and complex crystalline basement patterns in Animaqing suit zone and its adjacent area are obtained. The Pg wave velocity structures of the studied area show that the velocity structure in the lower part is more complex than in the upper part if roughly take the depth of 2.0km as the boundary. The velocity near surface is lower in Ruoergai basin where it is 3.86km/s~4.4km/s and higher in west Qinling fold zone of South Guansu where it is 4.0km/s~4.9km/s. Above the depth of 2.0km or so, the gradient variation of P velocity is not obvious vertically and the velocity structure is nerally homogeneous laterally in Ruoergai basin, while in the west fold
basin of South Gansu, the velocity gradient variation is obvious vertically and the velocity presents rather strong inhomogeneous laterally. The P velocity structure above the depth of 2.0km is much more complex than the area beneath the depth of 2.0 km, and there is a low velocity zone in the upper crust near North Xianman where Kusehu-Marqing fault goes through. , which possibly is the north boundary between Ganzhi and Songpan block. Wudu-Diebu fault and Zhouqu-Liangdan fault pass through the high resolving refraction profile at 9km, 30km in north Lanmusi, respectively. Beneath the depth of 2.0km, Wudu-Diebu fault shows a low velocity zone with small dimensions in velocity structure figure, and its depth is less than 6.0km. At 30km or so away Lanmusi in north, there is seemly a more large scale a low velocity distribution beneath the depth of 2.0km, and its location corresponds to Zhouqu-Liangdan fault. Because of its closing to the edge of inverse model and lack of the data, we can not give more detail patterns of this low velocity distribution. Generally speaking, Pg wave velocity structure can be divided into three regions, i.e.,south, middle and north regions, if taking both Xiaman and Lanmusi as geological block boundies. The south part is located in Ruoergai basin where the velocity variation is comparatively small and the middle part is located in Animaqin suit zone where its velocity varies greatly, while the north part lies in the West Qinling fold zone and its velocity variation is comparatively large.In the studied area, the differences of crystalline basement structures among different regions are more obvious than that of velocity distributions. Xiaman and Langmusi are the structure variation boundaries. The southern part lies in Zoige basin. Its crystalline basement is deeper and the structure is comparatively complicated. The largest depth of the crystalline basement is about 3.5km there and the basement velocity is about 5.65~5.8km/s; The middle part is located in Alimaqing suit zone and the crystalline basement is more complex than the southern and northern segments, and its depth variation is also large. The deformation of the basement is strongest in the middle segment. There are seemly double layer structure features beneath the depth of 2.5km, one of
which is located at the depth of 3.0km or so and its basement velocity is about 5.65-5.8km/s, and the another lies in the depth of 8.5km and the Pg velocity is about 6.2km/s; The basement structure in northern part is relatively simple and its depth is about 2.2km, and the basement velocity is about 5.5km/s. In Alimaqing suit zone, the crystalline basement is destroyed by Kusehu-Maqing fault zone and characterized by relatively large low velocity zone. Another two fault zones, i.e. Wudu-Diebu fault zone and Zhouqu-Liandan fault zone, are much small comparied with Kusehu-Maqing fault zone in the studied area.(7) Based on Matlab tool packages, the visualized programs for all results had been coded.
复杂地区精细地壳结构的探测与研究主要用地震方法。大量的地震测深结果表明,地壳内部的结构,特别是造山带地区上部地壳的结构,是十分复杂的。一般说来,上部地壳的结构较中下地壳要复杂得多,而且在地表结构复杂的地区更是如此。出于对研究成本和成果产出及探测效率的考虑,宽角反射/折射和高分辨折射己成为精细地壳结构探测的最主要方法之一。上部地壳精细结构的知识则主要来自于宽角反射/折射和高分辨折射人工地震探测Pg波资料的处理结果。一般来说,探测的越细,要求的分辨就越高,所面对的研究目标就越复杂,传统地震数据处理方法的不足和缺陷就逐渐暴露出来。在复杂地壳结构地区,传统的基于程函方程有限差分解的Pg波资料处理方法就存在着以下缺陷:①有限差分反演结果只能给出速度分布,不能给出上部地壳内的可能的界面结构特征,而地壳内的界面是客观存在的,尤其是上部地壳内的折射界面较为发育,其结晶基底具有全球性,近几年在城市活断层探测中大量的高分辨折射地震资料表明,海相沉积盆地的结晶基底一般为一强折射界面;②在程函方程的有限差分解中假设波阵面处处可微,但在一些结构模型特别是在较复杂速度结构介质中波阵面可能自身相交,意味着其梯度并不存在,将导致数值计算上的不稳定性;③计算速度较慢,
The knowledge of earth crust, especially the upper crust, is very important for the progress of geoscience itself and the long term strategy of disaster reduction as well as further exploitation of non-regeneration resources. The learning of the upper crustal structures builds a bridge between the geological results and the deeper crustal tectonics, it is also a first barrier for us to deeply understand the lower crust and even the whole lithosphere. Especially, the upper crust in complex tectonic areas is just like the very thick "cloud" which prevents us to use the "astronomical telescope" of seismic methods to deeply understand the seisgenic background of interior continent earthquakes with great disaster and the crustal evolution history along with the deep geodynamic process. Only after clearly comprehending the upper crustal structure like very thick "cloud" which never disappears in geological history can we possibly further recognize the seisgenic background and its mechanics, make more scientific estimation for the future earthquake disaster and provide more reliable deep seismic evidence for geotectonic theories. For this reason, the first step should be to gain a clear idea about the crustal structures, particularly the fine structure of upper crust.Seismic methods are mainly used to investigate the fine crust structures in complex region. A large number of results from DSS (Deep Seismic Sounding) show that the interior crust structure, particularly the upper crust structure in orogenic area is very complex. Generally speaking, the interior structure of upper crust is more complex than the middle and lower crust, and even more just so in the areas with complex surface structures. Considering the balance among the study cost and achievements as well as investigation efficiency, the wide angle reflection/refraction and high resolution refraction profiles have been main methods for the fine crust structure survey. The learning of the fine upper crust structure mainly come from the results of Pg data of wide angle reflection/refraction and high resolution refraction lines. By and
large, the more fine structure to be investigated, the more high resolution to be required, and the more complex the survey objective to be faced, the inadequacy and defects of traditional processing method for Pg data are gradually exposed. So far, there are following drawbacks in traditional processing method for Pg data, which is based on the finite difference solution of eikonal equation:D The finite difference results only provide velocity distribution and can not supply the possible interface structure features in the upper crust. In fact, the interior crust interfaces do exist, especially in upper crust the refraction interface well developed and the crystalline basement exists globally. A large number of high resolution seismic refraction dada from urban active fault survey in recent several years present that in generally, the crystalline basement of the deposit basin with marine facies is a strong refracting interface;D It is presumed that everywhere on wavefront is derivable in solving the eikonal equation with finite difference method, but the wavefront may be crossed by itself in some structural model, especially in the comparatively complex media model, which means that its gradient does not exist and causes the algorithm to be unstable;□ The computing speed is not fast, which can not well suit the needs of the survey for high resolution and fine complex structures and results in very consuming cost;D The model is parameterized with square meshes, which enable the eikonal equation to be solved conveniently with the finite difference approximation,but it is not reasonable. The complicated level of real crust is different from place to place and even along the same profile, it has comparative difference among the different segments, maybe some part of which is simple and another one is very complex. Furthermore, the data amount along the different part at the same line may be different while denser observation is carried out at the complex structure part. Therefore, whether the model mesh lines are dense or sparse should correspond to both the model complex level and the data amount.In the wake of the quickly progressing of seismic observation
technique and continuously improving of survey system as well as accumulating with a great number of seismic data, especially the needs of high resolution and fine complex structure study in urban active fault investigation and basic researches for some important geoscience problems, the high resolution refraction method will be more and more widely used. For this reason, it is very necessary to improve the drawback in artificial seismic data processing discussed above and develop seismic data processing technique that can be applied to arbitrarily complex crustal structure, particularly the fine structure survey of complex upper crust.In this paper some problems mentioned above in processing seismic data have been settled and the complete process of Pg data interpretation along with some typically applied examples, such as in Jiashi strong earthquake swarm area and Animqing suit belt as well as its adjacent region, are given. Following contents are included in this thesis:(1) In this paper, the histories of seismic refraction investigation and wavefront imaging methods in artificial seismic survey were reviewed. The characteristics of wavefront tracing algorithm in complex crustal structures and its defects in processing the Pg data of artificial seismic investigation were expounded rather completely. Some drawbacks have been improved and a new method for Pg data processing, which had been used to process and to analyze the Pg data obtained from the fine crust structure investigation profiles for urban active fault survey and some key basic geoscience study projects, was proposed. Good results were achieved.( 2 ) The defects during the forward computing for the first arrival Pg wave of the artificial seismic survey in complex structure areas had been removed and the finite different wavefront imaging method based on Huygens principle for velocity structure was proposed. In comparatively complex structures, the ray trace of Pg wave is not just a single refracted path which is presumed in traditional data processing. It presents direct wave near shot points, and at the certain distances away from the shot points it maybe appears as refracted wave, head waved or
diffracted wave which depends on different upper crust and crystalline basement structure features. The improving scheme as follows: based on Hole (1992) finite difference algorithm, using the 5 operators of Lecomte et al.(2000) to simulate the advancing, expanding, evolving and propagating of Pg wave field and to achieve the tracking for the Pg wavefront, and applying backprojectkm of SIRT(Simultaneous Iterative Reconstruction Technique) to reconstruct the P wave velocity structures.(3) The RHA(Ray Hit Analysis) method for determining the interface patterns of upper crust basements was proposed and had been used to process the artificial seismic data obtained from the high resolving refraction profiles carried out in north Sichuan and south Gansu region(NSFC Project) and in Xinjian Jiashi earthquake swarm region. Based on the P velocity model determined by wavefront imaging method, the ray hit numbers for each mesh are computed and the ray hit distribution maps are drawn up, which can be used to determine possible interface patterns. This is because if the crustal interface exists, it is sure that there are denser seismic rays ,with emergent angles near to 90,close to the interface, and the larger velocity contrast between the two sides of the interface, the denser the rays to be distributed. So except the ray dense regions caused by shot point positions, the distribution features of comparatively dense ray hits, with emergent angles near to 90°, should reflect the crustal interface positions.(4) The drawbacks of HRIM(Hagedoorn Refractor Imaging Method) had been removed, and this method was applied in seismic data processing in Fuzhou, Zhangzhou, North Hainan and Jiashi strong earthquake swarm area(Xu Zhaofan, Zhang Xiankang et al.,2002,2005; Jia Shixue,Xu Zhaofan et al., 2005; Urban active fault survey reports of Geophysical Prospecting Center, CEA, 2002,2003, 2004, 2005).Hagedoorn proved that the interface is located at the point where the travel time is just equal to the sum of forward and inverse propagation times. According to this principle, the refractor can be reconstructed. Four conditions must be satisfied using Hagedoorn principle to image the
refractors, i.e. U measured reciprocal travel time data from both forward and reverse layout; D accurate reciprocal time; D known the velocity structure above the refractor; D reconstruction method for seismic wave field. Condition □ and □ may be satisfied during data collection. The velocity structures can be determined by finite difference inversion or wavefront imaging in terms of Huygens principle. An applied example was given by Aldridge et al.(1992) and Lecomte et al.(2000), respectively. In their algorithm, the velocity structures above the refractor need to be provided by other method and, it is difficult in real applications. Based on wavefront expanding in terms of Huygens principle and applying the Lecomte's (2000) five wavefront expanding operator to simulate the seismic wave propagating so as to realize forward computing of seismic wave field, and using Hole's (1992) back-project method for inversion along with visualizing of the results by Matlab algorithms, a complete refractor imaging method in terms of Hagedoorn principle can be completed(Xu Zhao fan & Zhang Xiangkang et al., 2002) and has been used to process the seismic data from high resolving refraction survey for the exploration of fine upper crust structures in urban active fault investigations in some cities and seismic active areas, and good results are achieved.(5) The data of the high resolving refraction profile in Jiashi strong earthquake swarm area were processed in detail with wavefront imaging method in terms of Huygens' principle, synthetical ray hit analysis and Hagedoorn wavefront refractor imaging principle. The fine P wave velocity structures of upper crust and the clear patterns of basement in Jiashi strong earthquake swarm area were obtained, and the possible faults were inferred. Based on the comparatively complete observation system in Jiashi strong earthquake swarm area, the relations among the inverse parameters (such as initial model, iterative times, smoothing scale etc.), wavefront imaging results and ray distribution characteristics were analyzed and discussed in detail. In Jiashi strong earthquake swarm area, the basement with the depth of about 3.0km , where the Pg velocity is about 4.6km/s, appears continuously and completely, which possibly is
the bottom boundary between the complete developed marine platform basin faces and carbonic acid rock build and under which it is light metamorphic rock in middle and late Proterozoic Group ; The interface located at the depth of about 9.0km shows comparatively large variation and its depth gradually deepens from NE to S W, the depthes of which are 9.0km at southwest end and 8.5km at northeast end. This interface is characterized by the crystalline basement where the propagating velocity of Pg is about 6.25km/s. The crystalline basement between the post numbers 38km and 65km is continuous and complete and its depth range is from 8.5km to 9.0km. From the post number 37km to 38km, its depth goes down abruptly and it is about 11.5km between the post numbers 25km and 37km, the location of which is consistent with the fault beneath Jiashi inferred from wide angle reflection/refraction profile JA1 (Zhang Xiankang et al. 2002).From the wavefront imaging results of the velocity structures we can find that the crustal velocity structure above the depth of 1 lkm is not complicated and in generally speaking, the uppermost crust is homogeneous in laterally and conspicuously layered in vertically. The shallow cover above the depth of 400m is comparatively loose weathered layer and its P velocity is about 1.65~1.8km/s; The velocity isogram above the depth of 3.0km is nearly horizon, it indicates that the structure is homogeneous in laterally. Though there is a strong velocity gradient at the depth from 400m to 3.0kmt, the velocity isograms are even, i.e. the constant velocity gradient structure appears. Its P wave velocity varies from 1.8km/s to 4.5km/s. It is the second layer, which possible is the bottom boundary of clastic-barbonic acid rock build of well developed marine platform - basin facies of Tarim in Pz; The third layer is located at the depth of about 3.0km~9.0km, the upper part of which shows different velocity structure features from the lower part. The P velocity isograms of upper part are comparatively even and its depth range is from 3.0km to 5.5km, and the P wave velocity gradient obviously decreases. The lower part is a relatively weak gradient layer and its velocity is about 5.1km/s~6.0km/s, and it possible is the light
metamorphic rock in middle and late Proterozoic Group. The fourth layer locats beneath the depth of about 9.5km, which is relatively even and its P wave velocity is about 6.3km/s. This layer may be the middle and deep metamorphic rock in late Archean or early and middle Proterozoic Group.The upper part of third layer is almost even and the average P wave velocity is about 4.8km/s. There is a strong velocity perturbation in this layer, which implyes crustal faults exists at the depths from 4.5kmto 6.0km , and three faults Fl, F2, F3 are inferred. According to the petroleum-geological data, Maigaiti fault passes through the high resolving refraction profile near the post number of 32km and Xiasuhong - Maigaiti fault has two branches, which are located beneath the post numbers of 45km and 58km, respectively. Fl fault maybe just is the Maigaiti fault and its location is slightly close to North. F2 and F3 faults should be the parts of Xiasuhong-Maigaiti fault but there is no any evidence from our results which shows that the fault Fl, F2, F3 pass through the crystalline basement or stretch to the free surface. It is suggested that each of the faults Fl, F2, F3 is not joined with the deep crustal fault in the studied area i.e. the shallow and deep structures of the crust are independent from each other.(6) The Pg data from the high resolving refraction profile which nearly perpendicularly goes through Animaqing suit zone in North Sichuan and South Gansu province were processed in detail by means of Huygens principle wavefront imaging method and the ray hit analysis technique. The fine structures of the upper crust and complex crystalline basement patterns in Animaqing suit zone and its adjacent area are obtained. The Pg wave velocity structures of the studied area show that the velocity structure in the lower part is more complex than in the upper part if roughly take the depth of 2.0km as the boundary. The velocity near surface is lower in Ruoergai basin where it is 3.86km/s~4.4km/s and higher in west Qinling fold zone of South Guansu where it is 4.0km/s~4.9km/s. Above the depth of 2.0km or so, the gradient variation of P velocity is not obvious vertically and the velocity structure is nerally homogeneous laterally in Ruoergai basin, while in the west fold
basin of South Gansu, the velocity gradient variation is obvious vertically and the velocity presents rather strong inhomogeneous laterally. The P velocity structure above the depth of 2.0km is much more complex than the area beneath the depth of 2.0 km, and there is a low velocity zone in the upper crust near North Xianman where Kusehu-Marqing fault goes through. , which possibly is the north boundary between Ganzhi and Songpan block. Wudu-Diebu fault and Zhouqu-Liangdan fault pass through the high resolving refraction profile at 9km, 30km in north Lanmusi, respectively. Beneath the depth of 2.0km, Wudu-Diebu fault shows a low velocity zone with small dimensions in velocity structure figure, and its depth is less than 6.0km. At 30km or so away Lanmusi in north, there is seemly a more large scale a low velocity distribution beneath the depth of 2.0km, and its location corresponds to Zhouqu-Liangdan fault. Because of its closing to the edge of inverse model and lack of the data, we can not give more detail patterns of this low velocity distribution. Generally speaking, Pg wave velocity structure can be divided into three regions, i.e.,south, middle and north regions, if taking both Xiaman and Lanmusi as geological block boundies. The south part is located in Ruoergai basin where the velocity variation is comparatively small and the middle part is located in Animaqin suit zone where its velocity varies greatly, while the north part lies in the West Qinling fold zone and its velocity variation is comparatively large.In the studied area, the differences of crystalline basement structures among different regions are more obvious than that of velocity distributions. Xiaman and Langmusi are the structure variation boundaries. The southern part lies in Zoige basin. Its crystalline basement is deeper and the structure is comparatively complicated. The largest depth of the crystalline basement is about 3.5km there and the basement velocity is about 5.65~5.8km/s; The middle part is located in Alimaqing suit zone and the crystalline basement is more complex than the southern and northern segments, and its depth variation is also large. The deformation of the basement is strongest in the middle segment. There are seemly double layer structure features beneath the depth of 2.5km, one of
which is located at the depth of 3.0km or so and its basement velocity is about 5.65-5.8km/s, and the another lies in the depth of 8.5km and the Pg velocity is about 6.2km/s; The basement structure in northern part is relatively simple and its depth is about 2.2km, and the basement velocity is about 5.5km/s. In Alimaqing suit zone, the crystalline basement is destroyed by Kusehu-Maqing fault zone and characterized by relatively large low velocity zone. Another two fault zones, i.e. Wudu-Diebu fault zone and Zhouqu-Liandan fault zone, are much small comparied with Kusehu-Maqing fault zone in the studied area.(7) Based on Matlab tool packages, the visualized programs for all results had been coded.
引文
丁林,钟大赉,潘裕生等,1995,东喜马拉雅构造结上新世以来快速抬升的列裂变经迹证据,科学通报,40(16):1497～1500
丁志峰,何正勤,孙为国等,1999b,青藏高原东部及其边缘地区的地壳上地幔三维速度结构,地球物理学报,42(2),197—205。
丁志峰,何正勤,吴建平等,2001,青藏高原地震波三维速度结构的研究,中国地震,17(2),202-209
陈杰,曲国胜,沈军,伽师强震群区及帕米尔东北侧地震构造环境的研究,“九五”国家科技攻关计划,96—913—07,1999。
程裕淇,1994,中国区域地质概论,北京:地质出版社
段永红,张先康,杨卓欣,2002,帕米尔东北侧基底结构研究,地震学报,24(4),p378～384
段永红,张先康,杨卓欣等,2003,长白山天池火山区基地结构研究,地震地质,25(3),p501～508
段永红,张先康等,2002,华北地区上部地壳结构的三维有限差分层析成像,地球物理学报,Vol.45,No.3,p362～369
樊计昌,李松林,赖晓玲,等,新疆伽师地震区三维Q值结构,地震学报,Vol.23,No.6,2001。
嘉世旭,张先康,2005,华北不同构造块体地壳结构及其对比研究,地球物理学报,Vol48,N3,611—620
嘉世旭,张先康,方盛敏,2001,华北裂谷盆地不同块体地壳结构及演化研究,地学前缘,8(2):259～266
刘福田,1984,震源位置和速度结构的联合反演(1)——理论和方法,地球物理学报,27(2):167—175。
刘福田,曲克信,吴华等,1986,华北地区地震层面成像,地球物理学报,29(5):442—449
刘福田,李强,吴华等,1989a,用于速度图像重建的层析成像法,地球物理学报,32(1):46-61
刘福田,曲克信,吴华等,1989b,中国大陆及其邻区的地震层析成像,地球物理学报,3293):281-291
李吉均,1979,青藏高原隆起的时代、幅度和形式探讨,中国科学(B辑),6:608—616
李松林,张先康等,2001,多条人工地震测深剖面资料联合反演首都圈三维地壳结构,地球物理学报,Vol44,No3,p360～368
李松林,张先康,W. D. Mooney,等,伽师地震区地壳细结构及发震断层的初步研究,地球物理学报,Vol.45,No.1,2002。
刘启元,陈九辉,李顺成,等新疆伽师强震群区三维地壳上地幔S波速度结构及其地震成因的探讨,地球物理学报,Vol.43,No.3,2000。
刘志,张先康,周雪松,等,帕米尔东北侧地壳物性结构及其发震环境探讨,地震学报,Vol.25,No.3,2003。
潘裕生,1990,西昆仑山构造特征与演化,3:224—232
曲国胜,陈杰,陈新发,等,伽师中强震群成因及其未来地震趋势初探.见:中国地震局地质研究
孙洪烈,郑度,1998,青藏高原形成演化与发展,广州:广东科学出版社
滕吉文,王国正,刘道洪,等,1974,华北平原中部地区深部构造背景及邢台地震(一),地球物理学报,17(4),255—271
滕吉文,2003,固体地球物理学概论,北京:地震出版社
滕吉文,张中杰,白武明等,2004,岩石圈物理学,北京:科学出版社
王椿镛,二维射线合成地震图的计算及平凉—渭南地震测深剖面的解释,地震研究,8:615—626
王椿镛,1997,中国岩石层结构的回顾与展望,地球物理学报,40:(增刊)82—109
王椿镛,张先康,丁志峰等,1997,大别山造山带上部地壳结构的有限差分层析成像,地球物理学报,40(4),p459～501
王椿镛,张先康,吴庆举等,1994,华北盆地滑脱构造的地震学证据,地球物理学报,37(5),p620～631
于椿镛,张先康,吴庆举等,1994,翼中坳陷内深地震反射剖面揭示的滑脱构造,科学通报,39,625—628
王椿镛,1994,中国深地震反射剖面探测与研究,见:陈永泰主编,中国固体地球物理学进展,北京:海洋出版社,120—127
王椿镛,林中洋,1995,地壳地震学研究进展,见:陈永泰主编,地球与空间科学观测技术进展,北京:地震出版社,272—282
王椿镛,王贵美,林中洋,等,1993,用深地震反射方法研究邢台地震区细结构,地球物理学报,36(4):445-452
王椿镛,吴建平,楼海,等,2003,川西藏东地区的地壳P波速度结构,中国科学(D辑),增刊,181—189
王椿镛,Moony,W D,王溪莉,等,2002,川滇地区的地壳上地幔三维速度结构研究,地震学报,24:1—6
王椿镛,韩渭滨,吴建平,等,2003,松潘—甘孜造山带的地壳结构,地震学报,25(3),229—241
王椿镛,张先康,林中洋,等,1994,束鹿断陷盆地及其邻近地区的地壳结构特征,地震学报,16:472—479
王椿镛,林中洋,陈学波,等,1995,青海门源至福建宁德地学断面综合地球物理解释,地球物理学报,38:590—598
王椿镛,陆原,陈光英,等,1989,地震测深数据处理软件包,地震地磁观测与研究,10(5):50—58
王椿镛,陆原,谢建波,等,1991,地震测深数据文件系统,地球物理学报,34:455—464
王椿镛,楼海,宋亦海,等,1994,人工地震测深数据库系统,地震学报,13:41—52
伍致中,刘东海.1996,塔里木盆地西南坳陷的形成演化[J].新疆石油地质,Vol.17,No.3,
王琪,丁国瑜,乔学军,等,2000,用GPS研究南天山(伽师)地区现今内壳变形,地震学报,Vol.22,No.3,。
徐朝繁,张先康,杨建等,1998,三维壳幔结构中两点射线追踪的分步搜索法,中国地震,Vol.20,No3,p272～278
徐朝繁,张先康,杨键等,1999,非纵观测系统地震测深资料分析处理的人机对话过程,地球物理学进展,Vol 14,No 2,27—32
徐朝繁,张先康,朱金芳等2002,复杂结构中界面的哈格多恩原理折射波前成像方法,地震地质,Vol.24,No.4,p542～548
徐朝繁,张先康,段永红等,2005,新疆伽师强震群区上部地壳细结构的高分辨折射地震探测,地震学报,28(1),待发表
徐朝繁,张先康,胡修奇,2005,射线数分布分析法及其在复杂地壳结构探测中的应用,地震学报,待发表
徐朝繁,张先康,嘉世旭等,2005,琼东北及雷州半岛深地震测深剖面Pg波资料处理结果—有限差分成像和射线数分布分析法,地球物理学报,待发表
徐朝繁,张先康,段永红等,2005,川北甘南地区上部地壳结构的Pg波波前成像和射线数分布分析法结果 已投地球物理学报
胥颐,刘福田,刘建华,等,天山地震带地壳结构与强震构造环境,地球物理学报,Vol43,No2,2000。
肖序常,李廷栋,李光岑等,1988,喜马拉雅岩石圈构造演化总论,北京:地质出版社,1～236
杨欣,高国英,伽师强震群序列特征和震源机曾融生,阚荣举,何传大,1960,柴达木盆地的低频地震工作,科学通报,10:313～316
杨卓欣,赵金仁,张先康,等,伽师强震群区上地壳三维速度层析成像,地震学报,Vol.24,No.2,2002。
张先康,王椿镛等,1996,中国华北地壳上地幔结构研究,第二届海峡两岸地震学术研讨会论文集,p220～233
张先康、赵金仁、张成科,等,帕米尔东北侧地壳结构研究,地球物理学报Vol.45,No.5,2002。
张先康,杨玉春,赵平,等,1994,唐山滦县震区的三维地震透射研究——中、上地壳速度层析成像,地球物理学报,Vol.37,No.6,759-766.
张先康,祝志平,张成科等,1998,张家口—渤海地震带及其两侧地壳上地幔构造与速度结构研究,活动断裂研究,地震出版社,p1～16
张先康,张成科,赵金仁等,1995,唐山地震的深部构造背景,见:《发展中的地震科学研究》,北京:地震出版社
张先康,杨卓欣,杨玉春等,1994,地壳三维结构的层析成像方法——爆炸和地震资料的联合反演,地震学报,vol 17,No 4 422—431
张先康,王椿镛,刘国栋等,1996,延庆—怀来地区地壳细结构——利用深地震反射剖面,地球物理学报,vol 39 No 3 356—364
张先康,祝志平,张成科,等,1998,张家口—渤海地震带及其两侧地壳上地幔构造与速度结构研究,活动断裂研究(6),北京:地震出版社,1—16
张先康,赵金仁,刘国华等,2002,三河—平谷80级大震区震源细结构的深地震反射探测研究,中国地震,Vol 45,No4,326—336
张先康,刘国栋,刘泰升,等,1998,华北地壳结构的三维探测研究,见:寸丹集—庆祝刘光鼎院士工作50周年学术论文集,北京:科学出版社
张先康,张成科,赵金仁等,2002,长白山天池火山区岩浆系统深部结构的深地震测深研究,地震学报,Vol24,No 2,135—143
张先康,李松林,王夫运,等,2003,青藏高原东北缘、鄂尔多斯和唐山的地壳结构差异——深地震测深结果,地震地质,25(1):52—62
周仕勇,许忠淮,韩京,等,主震定位法分析以及1997年新疆伽师强震群高精度定位,地震学报,Vol.21,No.3,1999。
周仕勇,许忠淮,由震源谱推断1997年新疆伽师强震群破裂特性.地震学报.Vol.22,No.2,2000。
漳州城市活断层探测报告,2003,中国地震局地球物理勘探中心
沈阳城市活断层探测报告,2004,中国地震局地球物理勘探中心
西安及周边地区活断层探测报告,2005,中国地震局地球物理勘探中心
赵金仁,张先康,张成科,等,伽师—阿图什震区地壳深部结构特征的探测与研究,中国地 震,Vol.18,No.4,2002。
朱令仁,苏乃秦,杨马陵.1998,年新疆伽师强震群及三次成功的临震预报.中国地震,14(2):
曾融生,阚荣举,何传大,1961,柴达木盆地的低频地震探测的基岩首波和大角度反射波,地球物理学报,10(1):54～66
Abets, G. A. 1994. Three-dimensional inversion of regional P and S arrival times in the East Aleutians and sources of subduction zone gravity highs, J. Geophys. Res., 99, 4395-4412
Achauer, U., 1994. New ideas on the Henya rift based on the inversion of the combined dataset of the 1985 and 1989/90 seismic tomography experiments, Tectonophysics, 236, 305-329.
Aki, K. & Lee, W. H. K., 1976. Determination of the three-dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes l. A homogeneous initial model, J. Geophys. Res., 81, 4381-4399.
Aki, K. & Richards, P. G., 1980. Quantitative seismology: theory and methods, W. H. Freeman, San Francisco.
Aki, K., Christoffersson, A., & Husebye, E. S., 1977. Determination of the three-dimensional seismic structure of the lithosphere, J. Geophys. Res., 82, 277-296.
Aldridge, A. E. & Oldenburg, D. W. 1992. Refractor imaging using an automated wavefront reconstruction method. Geophysics, 57,378～385
Banda, E., Deichman N, Braile L. W. et al., 1982, Amplitude study of the Pg phase. Journal of Geophysics, 51:153～164
Benz, H. M., Zandt, G., & Oppenheimer, D. H., 1992. Lithospheric structure of northern California from teleseismic images of the upper mantle, J. Geophy. Res., 97, 4791-4807.
Blundell, C. A., 1993. Resolution analysis of seismic P—wave velocity estimates using reflection tomography inversion, PhD thesis, Monash University
Bishop, T. P. Bube, K. P., Cutler, R. T. et al., 1985. Tomographic determination of velocity and depth in laterally varying media, Geophysics, 50, 903-923.
Cara, M. & Leveque, J. J., I987. Waveform inversion using secondary observables, Geophy. Res. Lett., 14, 1046-1049 Chiarabba et al, 1997
Carroll, S. & Beresford, G., 1996. Combining reflection tomography with layer replacement for velocity analysis of near-surface reefs, Geophysics, 61,561-569
Cerveny, V., Molotkov, I. A. and Psencik, I. 1977, Ray method in seismology. Praha: Universita karlova
Cerveny, V., 1987. Ray tracing algorithms in three-dimensional laterally varying layered structures, in Seismic tomography: With applications in global seismology and exploration geophysics, edited by G. Nolet, pp. 99-133.
Cerveny, V., 2002. Seismic Ray Theory, Cambridge University Press, Cambridge.
Cheng, N. & House, L., 1996. Minimum traveltime calculations in 3-D graphy theory, Geophysics. 61, 1895-1898.
Chiarabba, C., Amoto, A., and Meghraoui, M., 1997, Tomographic images of the El Asnam fault zone. J. Geophys. Res. 102, 24485-24498
Chiu, S. K. L. , Kanasewich, E. R. , & Phadke, S. , 1986. Thre-dimensionaldetermation of structure and velocity by seismic tomography, geophysics, 51, 1559-1571.
Debayle, E. W. & Kennet, B. L. N. , 2000. The Australian continental upper mantle: Structure and deformation inferred from surface waves, J. Geophys. Res., 105, 25423-25450.
Eberhart-Phillip, D& Reyners, M. , 1997. Continental subduction and three-dimensional crustal structure: The northern South Island, New Zealand, J. Geophys. Res., 102,11848-11861.
Eberhart-Phillips, D.,1986. Three-dimensional velocity structure in northern California Coast ranges from inversion of local earthquake arrivals times, Mull. Seism. Soc. Am., 76,1025-1052.
Eberhart-Phillips, D., 1990. Three-dimensional P and S velocity structure in the Coalinga region, California, J. Geophys. Res., 95,15343-15363.
Eberhart-Phillips, D. & Michael, A. J. , 1993. three-dimensional velocity structure, seismicity and fault structure in the Parkfield region, central California, J. Geophys. Res., 98, 15737-15758.
England, P. C. & Mackenzie, D. P., 1982. A thin viscous sheet model for continental deformation, Geophys. J. R. 70:295-321
England, P.C. & Molnar, P. 1990. Right-lateral shear and rotation as the expiation for strike-slip faultinf in eastern Tibet. Nature, 344:140-142
Farra, E. L. & Madariaga, R.,1988. Non-linear reflection tomography, Geophys. J., 95, 135-147.
Frederiksen, A. W. , Bostock, M. G. & VanDecar, J. C. et al, 1998. Seismic structure of the upper mantle beneath the northern Canadian Cordillera from teleseismic travel-time inversion, Tectonophysics, 294, 43-55.
Fuis,G.S. , Ryberg,T, Lutter, W. J., et. al. , 2001, Seismic mapping, of shallow fault zones in the San Grabriel Mountains from the Los Angeles Region Seismic Experiment, southern California, J. G. R. 106(B4), 6549-6568
Graeber, F. M. & Asch, G. , 1999. Three-dimensional models of P wave velocity and P-to-S velocity ratio in the southern central Andes by simultaneous inversion of local earthquake data, J. Geophys. Res., 104, 20237-20256.
Graeber, F. M., Houseman, G. A., & Greenhalgh, S. A. , 2002. Regional teleseismic tomography of the western Lachlan Orogen and the Newer Volcanic province, southeast Australia, Geophys. J. Int., 149,249-266.
Glahn, A. & Granet, M. , 1993. Southern Rhine Graben: Small-wavelength tomographic study and implications for the dynamic evolution of the graben, Geophys. J. Int., 113,399-418.
Hagedoorm, J. G. 1959. The plus-minus method of interpreting seismic refraction sections [J]. Geophys. Prosp., 7, 158-182
Haslinger,F., Kissling, E. , & Ansorge, J. et al. , 1999. 3-D crustal structure from local earthquake tomography around the Gulf of Arta(Ionian region, NW Greece), Tectonophysics,304,201-218.
Hasegawa, A., Umino, N., & Takagi, A. , 1978. double-planed deep seismic zone and
upper-mantle structure in the northeastern Japan arc, Geophys. J. Royal Astr. Soc., 54, 281-296.
Humphreys, E. & Clayton, R. W. , 1988. Adaption of back projection tomography to seismic travel time problems, J. Geophys. Res. , 93, 1073-1085.
Humphreys, E. & Clayton, R. W. , 1990. Tomography image of the Southern California Mantle, J. Geophys. Res., 95, 19725-19746.
Hildebrand, J.A., Dorman, L. M. , & Hammer, P. T. C. et al. , 1989. Seismic tomography of Jasper Seamount, Geophys. Res. Lett. , 16, 1355-1358.
Hammer, P. T. C. , Dorman, L. M. , Hildebrand, J. H. , et al, 1994. Jasper Seamount structure : Seafloor seismic refraction tomography, J. Geophys. Res., 99, 6731-6752
Hole, J. A. 1992. Nonlinear high-resolution three-dimensional seismic travel time tomography[J], J Geophys Res., 135:6553-6562
Hole, J. A. & Zelt, B.C., 1995. 3-D finite-difference reflection travel times, Geophys. J. Int. , 121, p427~434
Hole, J.A., Clowes, R. M. & Ellis, R. M. 1992. Interface inversion using broadside seismic refraction data and three-dimensional traveltime calculations, J. Geophys. Res. , 97, p3417-3429
Iyer, H. &Hirahara, K. , 1993. Seismic Tomography: Theory and Practice, Chapman & Hall, London.
Julian, B. R. &Gubbins, D. , 1977. Three-dimensional seismic ray tracing, J. Geophys., 43, 95-113.
Kao, H. & Rau, R. J. , 1999. Detailed structures of the subducted Philippine Sea plate beneath northeast Taiwan: A new type of double seismic zone, J. Geophys. Res. , 104, 1015-1033.
Kawakatsu, H. , 1985. Double seismic zone in Tonga, Nature, 316, 53-55.
Kennet, B. L. N., 2002. The Seismic Wavefield, Volume 2; Interpretation of Seismograms on Regional and Global Scales, Cambridge University Press, Cambridge.
Klimes. L. & Kvasnicka, M. , 1994. 3-D network ray tracing, Geophys. J. Int. ,116,7 26-738
Kodaira, S. , Mjelde, R. & Gunnarsson, K. et al. 1998. Structure of the Jan Mayen microcontinent and implicatios for its evolution, Geophys. J. Int., 132, 383-400.
Kosloff , D. , Sherwood,J., & Koren, Z., et al. 1996. Velocity and interface depth determination by tomography of depth migrated gathers,' Geophysics, 61, 1511-1523.
Leveque, J. J. , &masson, F. , 1999., From ACH tomographic models to absolute velocity models, Geophys. J. Int., 137,621-629.
Leveque, J. J. , Rivern, L. , & Wittlinger, G. , 1993. On the use of the checker-board test to assess the resolution of tomographic inversions, Geophys. J. Int., 115, 313-318
Lecomte, I, Gjoystal, H, et al. 2002. Improving modeling and inversion in refraction seismic with a first-order Eikonal solver, Geophysica Prospecting, 48, p437~ 454
Lee,W. H.K.& Pereyra, V. , 1993. Mathematical introduction to seismic tomography, in
Seismic Tomography: theory and pratice, edited by H. M. Iyer & K. Hirahara, pp. 9-12,Chapman & Hall, London
Lutter, W. J. , Nowack, R. L. , & Brail, L., 1990. Seismic imaging of upper crustal structure using travel times from the PASSCAL Ouachita experiment, J. Geophys. Res. , 95, 4621-4631.
Mckenzie.D. 1990. Spinning continents. Nature, 344: 109-110.
McCaughey, M. & Singh, S. C. , 1997. Simultaneous velocity and interface tomography of normal-incidence and wide-aperture seismic traveltime data, Geophys. J. Int., 131, 87-99
Menke, W. , 1989. Geophysical data analysis: Discrete inverse theory, Academic Press, New York.
McQueen, H. W. S. & Lambeck, K. , 1996. Determination of crustal structure in central Australia by inversion of traveltime residuals, Geophys. J. Int., 126,645-662.
Moser, T. J. , 1991. Shortest path calculation of seismic rays, Geophysics, 56,59-67.
Nakanishi, I. & Yamaguchi, K. , 1986. A numerical experiment on nonlinear image reconstruction from first-arrival times for two-dimensional island arc structure, J. Phys. Earth, 34, 195-201
Neele, F., VanDecar, J. , & Snieder, R., 1993. The use of P wave amplitude data in a joint inversion with travel times for upper mantle velocity structure, J. Geophys. Res., 98, 12033-12054.
Nolet,G.,1987. Waveform tomography, in seismic Tomography: With applications in global seismology and exploration geophysics, edited by G. nolet,pp.301-322
Nolet, G. , 1990. Partitioned waveform inversion and two-dimensional structure under the network of autonomously recording seismographs, J. Geophys. Res. , 95, 8499-8512
Oncescu, M. C., Burlacu, V., & Anghel, M. et al, 1984. Three-dimensional P-wave velocity image under the Carpathian Arc, Tectonophysics, 106,305-319.
Pereyra, V. , Lee, W. H. K. , & Keller, H. B. , 1980. Solving two-point seismic-ray tracing problems in a heterogeneous medium, Bull. Seism. Soc. Am. , 70, 79-99.
Parker, R. L. , 1994. Geophysical inverse theory, Princeton University Press, Princeton.
Podvin, P. & Lecomete, L., 1991. Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools, Geophys. J. Int. , 105, 271-284
Qin, F., Luo, Y. , Olsen, Cai, K. B. , et al, 1992. Finite-difference solution of the eikonal equation along expanding wavefronts, Geophysics 57,478-471
Rawlinson , N. , Houseman, G. A., & Collins, C. D. N. , 2001a. Inversion of seismic refraction and wide-angle reflection traveltimes for 3-D layered crustal structure, Geophys. J. Int. , 145, 381-401.
Rawlinson , N. , & Houseman, G. A., 1998. Inversion for interface structure using teleseismic traveltime residuals, Geophys. J. Int., 133, 756-772.
Rawlinson, N. , &Sambridge, M. , 2004. Wave front evolution in strongly heterogeneous layered media using the fast marching method, Geophys. J. Int. , 156, 631-647
Riahi , M. A. , Lund, C. E. , & Pederson, L. B., 1997. Three-dimensional image of the
Moho undulations beneath the Gulf of Bothnia, Tectonophysics,239,149-164.
Ritsema, J. , Nyblade, A. A., & Owens, T. J. et al. , 1998. Upper mantle seismic velocity structure beneath Tanzania, East Africa: Implications for the stability of cratonic lithosphere, J. Geophys. Res., 103, 21201-21213.
Roller, J. C. , and W. H. Jackson, 1966, Seismic Wave Propagation in the Upper Mantle: Lake Superior, Wisconsin to Central. Arizona. Journal of Geophysical Research, 71:5933-5941.
Hill, N. R., 1987, Downward continuation of refracted arrivals to determine shallow structure, Geophyaics, 52, 1188-1198
Saltzer, R. L. & Humphreys, E. D., 1997. Upper mantle P wave velocity structure of the eastern Snake River Plain and its relationship to geodynamic models of the region, J. Geophys. Res., 102, 11829-11841.
Sambridge, M. S. & Kennett, B. L. N. , 1990. Boundary ray tracing in a heterogeneous medium: A simple and versatile algorithm, Geophys. J. Int., 101, 157-168.
Seber, D. , Barazangi, M. , &Tadili, B. A., et al.1996. Three-dimensional upper mantle structure beneath the intraplate Atlas and interplate Rif mountains of Morocco, J. Geophys. Res., 101, 3125-3138.
Sethian, J. A. , 1999. Level Set Methods and Fast Marching Methods, Cambridge University Press, Cambridge.
Sethian, J. A. , & Popovici, A.M., 1999. 3-D traveltime computation using the fast marching method, Geophysics 64, 516-523
Steinhart, J. S., 1964, Lake Superior Seismic Experiment: Shots and Travel times. Journal of Geophysical Research, 69:5335-5352
Sherif, R. E. & Geldart, L. P., 1995. Exploration Seismology (Second Edition), Cambridge University Press.
Scott, J. S. , Masters, T. G. & Vernon, F. L. , 1994. 3-D velocity structure of the San Jacinto fault zone near Anza, California-1. P waves, Geophys. J. Int., 119,611-626.
Smith, W. H. F. & Wessel, P., 1990, Gridding with continuous curvature splines in tension, Geophysics, 55,293-305.
Shalev, E., 1993. Cubic B- splines: Strategies of translating a simple structure to B-spline parameterization, Bull. Seism. Soc. Am., 83, 1617-1627
Steck,L. K., Thurber, C. H. , & Fehler, M., et al. , 1998. Crust and upper mantle P wave velocity structure beneath Valles Caldera, New mexico: results from the Jemez teleseismic tomography experiment, J. Geophys. Res., 103,24301-24320.
Tarantola, A., 1987. Inverse Problem Theory, Elsevier, Amsterdam
Tapponnier, P. &Molnar, P. 1976. Slip-line field theory and large scale continental tectonocs. Nature, 264: 319
Telford, W. M. , Geldart, L. P., and Sheriff, R. E. , 1990, "Applied Geophysics. " Cambridge University Press, Cambridge
Toomey, D. R., Solomon, S. C. , and Purdy, G. M. , 1994, Tomographic imaging of the shallow crustal structure of the East Pacific Rise at 930N, J. Geophys. Res., 99, 24135-24157
Thurber, C. H. , 1983. Earthquake locations and three-dimensional crustal structure
in the Coyote Lake area, central California, J. Geophys. Res., 88, 8226-8236.
Thurber, C. H. , 1993., Local earthquake tomography: Velocities and Vp/Vs-theory, in Seismic Tomography: theory and practice, edited by H. M. Iyer & K. Hirahara, pp. 563-583, Chapman & Hall, London.
Thomson, C. J. &Gubbins,D. 1982., Three-dimensional lithospheric modeling at NORSAR: linearity of the method and amplitude variations from the anomalies, Geophys. J. Royal Astr. Soc. , 71, 1-36.
Um, J. & Thurber, C. , 1987. A fast algorithm for two-point seismic ray tracing, Bull. Seism. Soc. Am., 77,972-986.
VanDecar, J. C. , James, D. E. , &Assumpcao, M. , 1995. Seismic evidence for a fossil mantle plume beneath South America and implications for plate driving forces, Nature, 378, 25-31.
Van Trier, J. & Symes, W. W. , 1991. Upwind finite-difference calculation of traveltimes, Geophysics 56, 812-821
Vidale, J. E. , 1988. Finite-difference calculations of traveltimes, Bull. Seism. Soc. Am., 78, 2062-2076
Vidale, J.E., 1990. Finite-difference calculations of traveltimes in three dimensions, Geophysics 55, 521-526
Wang, Y. & Pratt, R. G., 1997. Sensivities of seismic traveltimes and amplitudes in reflection tomography, Geophys. J. Int., 131,618-642.
Weiland, C.M., Steck, L. K. , & Dawson, P. B. et al, 1995. Nonlinear teleseismic tomography at long Valley Caldera, using three-dimensional minimum travel time ray tracing, J. Geophys. Res., 100, 20379-20390.
White, D. J., 1989. Two-dimensional seismic refraction tomography, Geophys. J. , 97, 223-245.
Wiggins, S. M. , Dorman, L. M., Cornuelle, B. D. , et al, 1996. Hess deep rift valley structure from seismic tomography, J. Geophys. Int., 130, 17-182.
Williamson, P. R. , 1990. Tomographic inversion in reflection seismology, Geophys. J. Int. , 100, 255-274.
Wu Lingzhao, Dave, A. Y., 1987. Injection of Indian crust into Tibetan lower crust: A temperature-dependent viscous model. Tectonics,6(4):505-514.
Yilmaz, O., 1987. Seismic data processing, Society of Exploration Geophysics, Tulsa.
Zelt, C. A. , Ellis, R. M. , 1988, Practical and efficient ray tracing in two-dimensional media for rapid traveltime and amplitude forward modelling, Can. J. expl. Geophys. , 24,16-31
Zelt, C. A. , Ellis, R. M. et al. 2001, Three - dimensional crustal velocity structure beneath the Strait of Georgia British Columbia, Geophys. J. Int., 144, p695~712
Zelt, C. A. , Barton, P.J., 1998, Three-dimensional seismic refraction tomography: A comparison of two methods applied to data from the Faeroe Basin, J. Geophys. Res., 103, p7187~7210
Zelt, C. A. & Smith, R. B. , 1992. Seismic traveltime inversion for 2-D crustal velocity structure, Geophys. J. Int., 108, 16-34
Zelt, C. A. , Hojka ,A. M. 1999. 3-D simultaneous seismic refraction and reflection tomography of wide-angle data from the central Chilean margin, Geophys. Res.
Lett. , 26, p2577~2580
Zelt, C. A. , 1998; Lateral velocity resolution from three-dimensional seismic refraction data, Geo-phys. J. Int., 135, 1101-1112.
Zelt, C. A., 1999. Modelling strategies and model assessment for wide-angle seismic traveltime data, Geophys. J. Int. , 139, 183-204.
Zhang, J. & Toksoz , M. N.,1998. Nonlinear refraction traveltime tomography, Geophysics, 63, 1726-1737.
Zhao, D. & Kanamori, H., 1992. P-wave image of the crust and uppermost mantle in southern Califormia, Geophys. Res. Lett., 23,2329-2332.
Zhao, D. , Hasegawa, A. , & Hoiuchi, S. , 1992. Tomographic imaging of P and S wave velocity structure beneath Northeastern Japan, J. Geophys. Res. , 97, 19909-19928.
Zhao, D. , Hasegawa, A. , & Kanamori, H. , 1994. Deep structure of Japan subduction zone as derived from local, regional, and teleseismic events, J. Geophys. Res. , 99, 22313-22329.
Zhao, D. .Kanamoi, H. , & Humphreys, E. , 1996. Simultaneous inversion of local and teleseismic data for the crust and mantle structure of southern California , phys. Earth planet. Inter.,93, 191-214.
Zhu, H. & Ebel, J. E. , 1994. Tomographic inversion for the seismic velocity structure beneath northern New England using seismic refraction data, J. Geophys. Res. , 99, 15, 331-15, 357.
Zielhuis, A & yan der hilst, R. D. , 1996. Upper-mantle shear velocity beneath easten Australia from inversion of waveforms from SKIPPY portable arrays, Geophys. J. Int. ,127, 1-16
丁志峰,何正勤,孙为国等,1999b,青藏高原东部及其边缘地区的地壳上地幔三维速度结构,地球物理学报,42(2),197—205。
丁志峰,何正勤,吴建平等,2001,青藏高原地震波三维速度结构的研究,中国地震,17(2),202-209
陈杰,曲国胜,沈军,伽师强震群区及帕米尔东北侧地震构造环境的研究,“九五”国家科技攻关计划,96—913—07,1999。
程裕淇,1994,中国区域地质概论,北京:地质出版社
段永红,张先康,杨卓欣,2002,帕米尔东北侧基底结构研究,地震学报,24(4),p378～384
段永红,张先康,杨卓欣等,2003,长白山天池火山区基地结构研究,地震地质,25(3),p501～508
段永红,张先康等,2002,华北地区上部地壳结构的三维有限差分层析成像,地球物理学报,Vol.45,No.3,p362～369
樊计昌,李松林,赖晓玲,等,新疆伽师地震区三维Q值结构,地震学报,Vol.23,No.6,2001。
嘉世旭,张先康,2005,华北不同构造块体地壳结构及其对比研究,地球物理学报,Vol48,N3,611—620
嘉世旭,张先康,方盛敏,2001,华北裂谷盆地不同块体地壳结构及演化研究,地学前缘,8(2):259～266
刘福田,1984,震源位置和速度结构的联合反演(1)——理论和方法,地球物理学报,27(2):167—175。
刘福田,曲克信,吴华等,1986,华北地区地震层面成像,地球物理学报,29(5):442—449
刘福田,李强,吴华等,1989a,用于速度图像重建的层析成像法,地球物理学报,32(1):46-61
刘福田,曲克信,吴华等,1989b,中国大陆及其邻区的地震层析成像,地球物理学报,3293):281-291
李吉均,1979,青藏高原隆起的时代、幅度和形式探讨,中国科学(B辑),6:608—616
李松林,张先康等,2001,多条人工地震测深剖面资料联合反演首都圈三维地壳结构,地球物理学报,Vol44,No3,p360～368
李松林,张先康,W. D. Mooney,等,伽师地震区地壳细结构及发震断层的初步研究,地球物理学报,Vol.45,No.1,2002。
刘启元,陈九辉,李顺成,等新疆伽师强震群区三维地壳上地幔S波速度结构及其地震成因的探讨,地球物理学报,Vol.43,No.3,2000。
刘志,张先康,周雪松,等,帕米尔东北侧地壳物性结构及其发震环境探讨,地震学报,Vol.25,No.3,2003。
潘裕生,1990,西昆仑山构造特征与演化,3:224—232
曲国胜,陈杰,陈新发,等,伽师中强震群成因及其未来地震趋势初探.见:中国地震局地质研究
孙洪烈,郑度,1998,青藏高原形成演化与发展,广州:广东科学出版社
滕吉文,王国正,刘道洪,等,1974,华北平原中部地区深部构造背景及邢台地震(一),地球物理学报,17(4),255—271
滕吉文,2003,固体地球物理学概论,北京:地震出版社
滕吉文,张中杰,白武明等,2004,岩石圈物理学,北京:科学出版社
王椿镛,二维射线合成地震图的计算及平凉—渭南地震测深剖面的解释,地震研究,8:615—626
王椿镛,1997,中国岩石层结构的回顾与展望,地球物理学报,40:(增刊)82—109
王椿镛,张先康,丁志峰等,1997,大别山造山带上部地壳结构的有限差分层析成像,地球物理学报,40(4),p459～501
王椿镛,张先康,吴庆举等,1994,华北盆地滑脱构造的地震学证据,地球物理学报,37(5),p620～631
于椿镛,张先康,吴庆举等,1994,翼中坳陷内深地震反射剖面揭示的滑脱构造,科学通报,39,625—628
王椿镛,1994,中国深地震反射剖面探测与研究,见:陈永泰主编,中国固体地球物理学进展,北京:海洋出版社,120—127
王椿镛,林中洋,1995,地壳地震学研究进展,见:陈永泰主编,地球与空间科学观测技术进展,北京:地震出版社,272—282
王椿镛,王贵美,林中洋,等,1993,用深地震反射方法研究邢台地震区细结构,地球物理学报,36(4):445-452
王椿镛,吴建平,楼海,等,2003,川西藏东地区的地壳P波速度结构,中国科学(D辑),增刊,181—189
王椿镛,Moony,W D,王溪莉,等,2002,川滇地区的地壳上地幔三维速度结构研究,地震学报,24:1—6
王椿镛,韩渭滨,吴建平,等,2003,松潘—甘孜造山带的地壳结构,地震学报,25(3),229—241
王椿镛,张先康,林中洋,等,1994,束鹿断陷盆地及其邻近地区的地壳结构特征,地震学报,16:472—479
王椿镛,林中洋,陈学波,等,1995,青海门源至福建宁德地学断面综合地球物理解释,地球物理学报,38:590—598
王椿镛,陆原,陈光英,等,1989,地震测深数据处理软件包,地震地磁观测与研究,10(5):50—58
王椿镛,陆原,谢建波,等,1991,地震测深数据文件系统,地球物理学报,34:455—464
王椿镛,楼海,宋亦海,等,1994,人工地震测深数据库系统,地震学报,13:41—52
伍致中,刘东海.1996,塔里木盆地西南坳陷的形成演化[J].新疆石油地质,Vol.17,No.3,
王琪,丁国瑜,乔学军,等,2000,用GPS研究南天山(伽师)地区现今内壳变形,地震学报,Vol.22,No.3,。
徐朝繁,张先康,杨建等,1998,三维壳幔结构中两点射线追踪的分步搜索法,中国地震,Vol.20,No3,p272～278
徐朝繁,张先康,杨键等,1999,非纵观测系统地震测深资料分析处理的人机对话过程,地球物理学进展,Vol 14,No 2,27—32
徐朝繁,张先康,朱金芳等2002,复杂结构中界面的哈格多恩原理折射波前成像方法,地震地质,Vol.24,No.4,p542～548
徐朝繁,张先康,段永红等,2005,新疆伽师强震群区上部地壳细结构的高分辨折射地震探测,地震学报,28(1),待发表
徐朝繁,张先康,胡修奇,2005,射线数分布分析法及其在复杂地壳结构探测中的应用,地震学报,待发表
徐朝繁,张先康,嘉世旭等,2005,琼东北及雷州半岛深地震测深剖面Pg波资料处理结果—有限差分成像和射线数分布分析法,地球物理学报,待发表
徐朝繁,张先康,段永红等,2005,川北甘南地区上部地壳结构的Pg波波前成像和射线数分布分析法结果 已投地球物理学报
胥颐,刘福田,刘建华,等,天山地震带地壳结构与强震构造环境,地球物理学报,Vol43,No2,2000。
肖序常,李廷栋,李光岑等,1988,喜马拉雅岩石圈构造演化总论,北京:地质出版社,1～236
杨欣,高国英,伽师强震群序列特征和震源机曾融生,阚荣举,何传大,1960,柴达木盆地的低频地震工作,科学通报,10:313～316
杨卓欣,赵金仁,张先康,等,伽师强震群区上地壳三维速度层析成像,地震学报,Vol.24,No.2,2002。
张先康,王椿镛等,1996,中国华北地壳上地幔结构研究,第二届海峡两岸地震学术研讨会论文集,p220～233
张先康、赵金仁、张成科,等,帕米尔东北侧地壳结构研究,地球物理学报Vol.45,No.5,2002。
张先康,杨玉春,赵平,等,1994,唐山滦县震区的三维地震透射研究——中、上地壳速度层析成像,地球物理学报,Vol.37,No.6,759-766.
张先康,祝志平,张成科等,1998,张家口—渤海地震带及其两侧地壳上地幔构造与速度结构研究,活动断裂研究,地震出版社,p1～16
张先康,张成科,赵金仁等,1995,唐山地震的深部构造背景,见:《发展中的地震科学研究》,北京:地震出版社
张先康,杨卓欣,杨玉春等,1994,地壳三维结构的层析成像方法——爆炸和地震资料的联合反演,地震学报,vol 17,No 4 422—431
张先康,王椿镛,刘国栋等,1996,延庆—怀来地区地壳细结构——利用深地震反射剖面,地球物理学报,vol 39 No 3 356—364
张先康,祝志平,张成科,等,1998,张家口—渤海地震带及其两侧地壳上地幔构造与速度结构研究,活动断裂研究(6),北京:地震出版社,1—16
张先康,赵金仁,刘国华等,2002,三河—平谷80级大震区震源细结构的深地震反射探测研究,中国地震,Vol 45,No4,326—336
张先康,刘国栋,刘泰升,等,1998,华北地壳结构的三维探测研究,见:寸丹集—庆祝刘光鼎院士工作50周年学术论文集,北京:科学出版社
张先康,张成科,赵金仁等,2002,长白山天池火山区岩浆系统深部结构的深地震测深研究,地震学报,Vol24,No 2,135—143
张先康,李松林,王夫运,等,2003,青藏高原东北缘、鄂尔多斯和唐山的地壳结构差异——深地震测深结果,地震地质,25(1):52—62
周仕勇,许忠淮,韩京,等,主震定位法分析以及1997年新疆伽师强震群高精度定位,地震学报,Vol.21,No.3,1999。
周仕勇,许忠淮,由震源谱推断1997年新疆伽师强震群破裂特性.地震学报.Vol.22,No.2,2000。
漳州城市活断层探测报告,2003,中国地震局地球物理勘探中心
沈阳城市活断层探测报告,2004,中国地震局地球物理勘探中心
西安及周边地区活断层探测报告,2005,中国地震局地球物理勘探中心
赵金仁,张先康,张成科,等,伽师—阿图什震区地壳深部结构特征的探测与研究,中国地 震,Vol.18,No.4,2002。
朱令仁,苏乃秦,杨马陵.1998,年新疆伽师强震群及三次成功的临震预报.中国地震,14(2):
曾融生,阚荣举,何传大,1961,柴达木盆地的低频地震探测的基岩首波和大角度反射波,地球物理学报,10(1):54～66
Abets, G. A. 1994. Three-dimensional inversion of regional P and S arrival times in the East Aleutians and sources of subduction zone gravity highs, J. Geophys. Res., 99, 4395-4412
Achauer, U., 1994. New ideas on the Henya rift based on the inversion of the combined dataset of the 1985 and 1989/90 seismic tomography experiments, Tectonophysics, 236, 305-329.
Aki, K. & Lee, W. H. K., 1976. Determination of the three-dimensional velocity anomalies under a seismic array using first P arrival times from local earthquakes l. A homogeneous initial model, J. Geophys. Res., 81, 4381-4399.
Aki, K. & Richards, P. G., 1980. Quantitative seismology: theory and methods, W. H. Freeman, San Francisco.
Aki, K., Christoffersson, A., & Husebye, E. S., 1977. Determination of the three-dimensional seismic structure of the lithosphere, J. Geophys. Res., 82, 277-296.
Aldridge, A. E. & Oldenburg, D. W. 1992. Refractor imaging using an automated wavefront reconstruction method. Geophysics, 57,378～385
Banda, E., Deichman N, Braile L. W. et al., 1982, Amplitude study of the Pg phase. Journal of Geophysics, 51:153～164
Benz, H. M., Zandt, G., & Oppenheimer, D. H., 1992. Lithospheric structure of northern California from teleseismic images of the upper mantle, J. Geophy. Res., 97, 4791-4807.
Blundell, C. A., 1993. Resolution analysis of seismic P—wave velocity estimates using reflection tomography inversion, PhD thesis, Monash University
Bishop, T. P. Bube, K. P., Cutler, R. T. et al., 1985. Tomographic determination of velocity and depth in laterally varying media, Geophysics, 50, 903-923.
Cara, M. & Leveque, J. J., I987. Waveform inversion using secondary observables, Geophy. Res. Lett., 14, 1046-1049 Chiarabba et al, 1997
Carroll, S. & Beresford, G., 1996. Combining reflection tomography with layer replacement for velocity analysis of near-surface reefs, Geophysics, 61,561-569
Cerveny, V., Molotkov, I. A. and Psencik, I. 1977, Ray method in seismology. Praha: Universita karlova
Cerveny, V., 1987. Ray tracing algorithms in three-dimensional laterally varying layered structures, in Seismic tomography: With applications in global seismology and exploration geophysics, edited by G. Nolet, pp. 99-133.
Cerveny, V., 2002. Seismic Ray Theory, Cambridge University Press, Cambridge.
Cheng, N. & House, L., 1996. Minimum traveltime calculations in 3-D graphy theory, Geophysics. 61, 1895-1898.
Chiarabba, C., Amoto, A., and Meghraoui, M., 1997, Tomographic images of the El Asnam fault zone. J. Geophys. Res. 102, 24485-24498
Chiu, S. K. L. , Kanasewich, E. R. , & Phadke, S. , 1986. Thre-dimensionaldetermation of structure and velocity by seismic tomography, geophysics, 51, 1559-1571.
Debayle, E. W. & Kennet, B. L. N. , 2000. The Australian continental upper mantle: Structure and deformation inferred from surface waves, J. Geophys. Res., 105, 25423-25450.
Eberhart-Phillip, D& Reyners, M. , 1997. Continental subduction and three-dimensional crustal structure: The northern South Island, New Zealand, J. Geophys. Res., 102,11848-11861.
Eberhart-Phillips, D.,1986. Three-dimensional velocity structure in northern California Coast ranges from inversion of local earthquake arrivals times, Mull. Seism. Soc. Am., 76,1025-1052.
Eberhart-Phillips, D., 1990. Three-dimensional P and S velocity structure in the Coalinga region, California, J. Geophys. Res., 95,15343-15363.
Eberhart-Phillips, D. & Michael, A. J. , 1993. three-dimensional velocity structure, seismicity and fault structure in the Parkfield region, central California, J. Geophys. Res., 98, 15737-15758.
England, P. C. & Mackenzie, D. P., 1982. A thin viscous sheet model for continental deformation, Geophys. J. R. 70:295-321
England, P.C. & Molnar, P. 1990. Right-lateral shear and rotation as the expiation for strike-slip faultinf in eastern Tibet. Nature, 344:140-142
Farra, E. L. & Madariaga, R.,1988. Non-linear reflection tomography, Geophys. J., 95, 135-147.
Frederiksen, A. W. , Bostock, M. G. & VanDecar, J. C. et al, 1998. Seismic structure of the upper mantle beneath the northern Canadian Cordillera from teleseismic travel-time inversion, Tectonophysics, 294, 43-55.
Fuis,G.S. , Ryberg,T, Lutter, W. J., et. al. , 2001, Seismic mapping, of shallow fault zones in the San Grabriel Mountains from the Los Angeles Region Seismic Experiment, southern California, J. G. R. 106(B4), 6549-6568
Graeber, F. M. & Asch, G. , 1999. Three-dimensional models of P wave velocity and P-to-S velocity ratio in the southern central Andes by simultaneous inversion of local earthquake data, J. Geophys. Res., 104, 20237-20256.
Graeber, F. M., Houseman, G. A., & Greenhalgh, S. A. , 2002. Regional teleseismic tomography of the western Lachlan Orogen and the Newer Volcanic province, southeast Australia, Geophys. J. Int., 149,249-266.
Glahn, A. & Granet, M. , 1993. Southern Rhine Graben: Small-wavelength tomographic study and implications for the dynamic evolution of the graben, Geophys. J. Int., 113,399-418.
Hagedoorm, J. G. 1959. The plus-minus method of interpreting seismic refraction sections [J]. Geophys. Prosp., 7, 158-182
Haslinger,F., Kissling, E. , & Ansorge, J. et al. , 1999. 3-D crustal structure from local earthquake tomography around the Gulf of Arta(Ionian region, NW Greece), Tectonophysics,304,201-218.
Hasegawa, A., Umino, N., & Takagi, A. , 1978. double-planed deep seismic zone and
upper-mantle structure in the northeastern Japan arc, Geophys. J. Royal Astr. Soc., 54, 281-296.
Humphreys, E. & Clayton, R. W. , 1988. Adaption of back projection tomography to seismic travel time problems, J. Geophys. Res. , 93, 1073-1085.
Humphreys, E. & Clayton, R. W. , 1990. Tomography image of the Southern California Mantle, J. Geophys. Res., 95, 19725-19746.
Hildebrand, J.A., Dorman, L. M. , & Hammer, P. T. C. et al. , 1989. Seismic tomography of Jasper Seamount, Geophys. Res. Lett. , 16, 1355-1358.
Hammer, P. T. C. , Dorman, L. M. , Hildebrand, J. H. , et al, 1994. Jasper Seamount structure : Seafloor seismic refraction tomography, J. Geophys. Res., 99, 6731-6752
Hole, J. A. 1992. Nonlinear high-resolution three-dimensional seismic travel time tomography[J], J Geophys Res., 135:6553-6562
Hole, J. A. & Zelt, B.C., 1995. 3-D finite-difference reflection travel times, Geophys. J. Int. , 121, p427~434
Hole, J.A., Clowes, R. M. & Ellis, R. M. 1992. Interface inversion using broadside seismic refraction data and three-dimensional traveltime calculations, J. Geophys. Res. , 97, p3417-3429
Iyer, H. &Hirahara, K. , 1993. Seismic Tomography: Theory and Practice, Chapman & Hall, London.
Julian, B. R. &Gubbins, D. , 1977. Three-dimensional seismic ray tracing, J. Geophys., 43, 95-113.
Kao, H. & Rau, R. J. , 1999. Detailed structures of the subducted Philippine Sea plate beneath northeast Taiwan: A new type of double seismic zone, J. Geophys. Res. , 104, 1015-1033.
Kawakatsu, H. , 1985. Double seismic zone in Tonga, Nature, 316, 53-55.
Kennet, B. L. N., 2002. The Seismic Wavefield, Volume 2; Interpretation of Seismograms on Regional and Global Scales, Cambridge University Press, Cambridge.
Klimes. L. & Kvasnicka, M. , 1994. 3-D network ray tracing, Geophys. J. Int. ,116,7 26-738
Kodaira, S. , Mjelde, R. & Gunnarsson, K. et al. 1998. Structure of the Jan Mayen microcontinent and implicatios for its evolution, Geophys. J. Int., 132, 383-400.
Kosloff , D. , Sherwood,J., & Koren, Z., et al. 1996. Velocity and interface depth determination by tomography of depth migrated gathers,' Geophysics, 61, 1511-1523.
Leveque, J. J. , &masson, F. , 1999., From ACH tomographic models to absolute velocity models, Geophys. J. Int., 137,621-629.
Leveque, J. J. , Rivern, L. , & Wittlinger, G. , 1993. On the use of the checker-board test to assess the resolution of tomographic inversions, Geophys. J. Int., 115, 313-318
Lecomte, I, Gjoystal, H, et al. 2002. Improving modeling and inversion in refraction seismic with a first-order Eikonal solver, Geophysica Prospecting, 48, p437~ 454
Lee,W. H.K.& Pereyra, V. , 1993. Mathematical introduction to seismic tomography, in
Seismic Tomography: theory and pratice, edited by H. M. Iyer & K. Hirahara, pp. 9-12,Chapman & Hall, London
Lutter, W. J. , Nowack, R. L. , & Brail, L., 1990. Seismic imaging of upper crustal structure using travel times from the PASSCAL Ouachita experiment, J. Geophys. Res. , 95, 4621-4631.
Mckenzie.D. 1990. Spinning continents. Nature, 344: 109-110.
McCaughey, M. & Singh, S. C. , 1997. Simultaneous velocity and interface tomography of normal-incidence and wide-aperture seismic traveltime data, Geophys. J. Int., 131, 87-99
Menke, W. , 1989. Geophysical data analysis: Discrete inverse theory, Academic Press, New York.
McQueen, H. W. S. & Lambeck, K. , 1996. Determination of crustal structure in central Australia by inversion of traveltime residuals, Geophys. J. Int., 126,645-662.
Moser, T. J. , 1991. Shortest path calculation of seismic rays, Geophysics, 56,59-67.
Nakanishi, I. & Yamaguchi, K. , 1986. A numerical experiment on nonlinear image reconstruction from first-arrival times for two-dimensional island arc structure, J. Phys. Earth, 34, 195-201
Neele, F., VanDecar, J. , & Snieder, R., 1993. The use of P wave amplitude data in a joint inversion with travel times for upper mantle velocity structure, J. Geophys. Res., 98, 12033-12054.
Nolet,G.,1987. Waveform tomography, in seismic Tomography: With applications in global seismology and exploration geophysics, edited by G. nolet,pp.301-322
Nolet, G. , 1990. Partitioned waveform inversion and two-dimensional structure under the network of autonomously recording seismographs, J. Geophys. Res. , 95, 8499-8512
Oncescu, M. C., Burlacu, V., & Anghel, M. et al, 1984. Three-dimensional P-wave velocity image under the Carpathian Arc, Tectonophysics, 106,305-319.
Pereyra, V. , Lee, W. H. K. , & Keller, H. B. , 1980. Solving two-point seismic-ray tracing problems in a heterogeneous medium, Bull. Seism. Soc. Am. , 70, 79-99.
Parker, R. L. , 1994. Geophysical inverse theory, Princeton University Press, Princeton.
Podvin, P. & Lecomete, L., 1991. Finite difference computation of traveltimes in very contrasted velocity models: a massively parallel approach and its associated tools, Geophys. J. Int. , 105, 271-284
Qin, F., Luo, Y. , Olsen, Cai, K. B. , et al, 1992. Finite-difference solution of the eikonal equation along expanding wavefronts, Geophysics 57,478-471
Rawlinson , N. , Houseman, G. A., & Collins, C. D. N. , 2001a. Inversion of seismic refraction and wide-angle reflection traveltimes for 3-D layered crustal structure, Geophys. J. Int. , 145, 381-401.
Rawlinson , N. , & Houseman, G. A., 1998. Inversion for interface structure using teleseismic traveltime residuals, Geophys. J. Int., 133, 756-772.
Rawlinson, N. , &Sambridge, M. , 2004. Wave front evolution in strongly heterogeneous layered media using the fast marching method, Geophys. J. Int. , 156, 631-647
Riahi , M. A. , Lund, C. E. , & Pederson, L. B., 1997. Three-dimensional image of the
Moho undulations beneath the Gulf of Bothnia, Tectonophysics,239,149-164.
Ritsema, J. , Nyblade, A. A., & Owens, T. J. et al. , 1998. Upper mantle seismic velocity structure beneath Tanzania, East Africa: Implications for the stability of cratonic lithosphere, J. Geophys. Res., 103, 21201-21213.
Roller, J. C. , and W. H. Jackson, 1966, Seismic Wave Propagation in the Upper Mantle: Lake Superior, Wisconsin to Central. Arizona. Journal of Geophysical Research, 71:5933-5941.
Hill, N. R., 1987, Downward continuation of refracted arrivals to determine shallow structure, Geophyaics, 52, 1188-1198
Saltzer, R. L. & Humphreys, E. D., 1997. Upper mantle P wave velocity structure of the eastern Snake River Plain and its relationship to geodynamic models of the region, J. Geophys. Res., 102, 11829-11841.
Sambridge, M. S. & Kennett, B. L. N. , 1990. Boundary ray tracing in a heterogeneous medium: A simple and versatile algorithm, Geophys. J. Int., 101, 157-168.
Seber, D. , Barazangi, M. , &Tadili, B. A., et al.1996. Three-dimensional upper mantle structure beneath the intraplate Atlas and interplate Rif mountains of Morocco, J. Geophys. Res., 101, 3125-3138.
Sethian, J. A. , 1999. Level Set Methods and Fast Marching Methods, Cambridge University Press, Cambridge.
Sethian, J. A. , & Popovici, A.M., 1999. 3-D traveltime computation using the fast marching method, Geophysics 64, 516-523
Steinhart, J. S., 1964, Lake Superior Seismic Experiment: Shots and Travel times. Journal of Geophysical Research, 69:5335-5352
Sherif, R. E. & Geldart, L. P., 1995. Exploration Seismology (Second Edition), Cambridge University Press.
Scott, J. S. , Masters, T. G. & Vernon, F. L. , 1994. 3-D velocity structure of the San Jacinto fault zone near Anza, California-1. P waves, Geophys. J. Int., 119,611-626.
Smith, W. H. F. & Wessel, P., 1990, Gridding with continuous curvature splines in tension, Geophysics, 55,293-305.
Shalev, E., 1993. Cubic B- splines: Strategies of translating a simple structure to B-spline parameterization, Bull. Seism. Soc. Am., 83, 1617-1627
Steck,L. K., Thurber, C. H. , & Fehler, M., et al. , 1998. Crust and upper mantle P wave velocity structure beneath Valles Caldera, New mexico: results from the Jemez teleseismic tomography experiment, J. Geophys. Res., 103,24301-24320.
Tarantola, A., 1987. Inverse Problem Theory, Elsevier, Amsterdam
Tapponnier, P. &Molnar, P. 1976. Slip-line field theory and large scale continental tectonocs. Nature, 264: 319
Telford, W. M. , Geldart, L. P., and Sheriff, R. E. , 1990, "Applied Geophysics. " Cambridge University Press, Cambridge
Toomey, D. R., Solomon, S. C. , and Purdy, G. M. , 1994, Tomographic imaging of the shallow crustal structure of the East Pacific Rise at 930N, J. Geophys. Res., 99, 24135-24157
Thurber, C. H. , 1983. Earthquake locations and three-dimensional crustal structure
in the Coyote Lake area, central California, J. Geophys. Res., 88, 8226-8236.
Thurber, C. H. , 1993., Local earthquake tomography: Velocities and Vp/Vs-theory, in Seismic Tomography: theory and practice, edited by H. M. Iyer & K. Hirahara, pp. 563-583, Chapman & Hall, London.
Thomson, C. J. &Gubbins,D. 1982., Three-dimensional lithospheric modeling at NORSAR: linearity of the method and amplitude variations from the anomalies, Geophys. J. Royal Astr. Soc. , 71, 1-36.
Um, J. & Thurber, C. , 1987. A fast algorithm for two-point seismic ray tracing, Bull. Seism. Soc. Am., 77,972-986.
VanDecar, J. C. , James, D. E. , &Assumpcao, M. , 1995. Seismic evidence for a fossil mantle plume beneath South America and implications for plate driving forces, Nature, 378, 25-31.
Van Trier, J. & Symes, W. W. , 1991. Upwind finite-difference calculation of traveltimes, Geophysics 56, 812-821
Vidale, J. E. , 1988. Finite-difference calculations of traveltimes, Bull. Seism. Soc. Am., 78, 2062-2076
Vidale, J.E., 1990. Finite-difference calculations of traveltimes in three dimensions, Geophysics 55, 521-526
Wang, Y. & Pratt, R. G., 1997. Sensivities of seismic traveltimes and amplitudes in reflection tomography, Geophys. J. Int., 131,618-642.
Weiland, C.M., Steck, L. K. , & Dawson, P. B. et al, 1995. Nonlinear teleseismic tomography at long Valley Caldera, using three-dimensional minimum travel time ray tracing, J. Geophys. Res., 100, 20379-20390.
White, D. J., 1989. Two-dimensional seismic refraction tomography, Geophys. J. , 97, 223-245.
Wiggins, S. M. , Dorman, L. M., Cornuelle, B. D. , et al, 1996. Hess deep rift valley structure from seismic tomography, J. Geophys. Int., 130, 17-182.
Williamson, P. R. , 1990. Tomographic inversion in reflection seismology, Geophys. J. Int. , 100, 255-274.
Wu Lingzhao, Dave, A. Y., 1987. Injection of Indian crust into Tibetan lower crust: A temperature-dependent viscous model. Tectonics,6(4):505-514.
Yilmaz, O., 1987. Seismic data processing, Society of Exploration Geophysics, Tulsa.
Zelt, C. A. , Ellis, R. M. , 1988, Practical and efficient ray tracing in two-dimensional media for rapid traveltime and amplitude forward modelling, Can. J. expl. Geophys. , 24,16-31
Zelt, C. A. , Ellis, R. M. et al. 2001, Three - dimensional crustal velocity structure beneath the Strait of Georgia British Columbia, Geophys. J. Int., 144, p695~712
Zelt, C. A. , Barton, P.J., 1998, Three-dimensional seismic refraction tomography: A comparison of two methods applied to data from the Faeroe Basin, J. Geophys. Res., 103, p7187~7210
Zelt, C. A. & Smith, R. B. , 1992. Seismic traveltime inversion for 2-D crustal velocity structure, Geophys. J. Int., 108, 16-34
Zelt, C. A. , Hojka ,A. M. 1999. 3-D simultaneous seismic refraction and reflection tomography of wide-angle data from the central Chilean margin, Geophys. Res.
Lett. , 26, p2577~2580
Zelt, C. A. , 1998; Lateral velocity resolution from three-dimensional seismic refraction data, Geo-phys. J. Int., 135, 1101-1112.
Zelt, C. A., 1999. Modelling strategies and model assessment for wide-angle seismic traveltime data, Geophys. J. Int. , 139, 183-204.
Zhang, J. & Toksoz , M. N.,1998. Nonlinear refraction traveltime tomography, Geophysics, 63, 1726-1737.
Zhao, D. & Kanamori, H., 1992. P-wave image of the crust and uppermost mantle in southern Califormia, Geophys. Res. Lett., 23,2329-2332.
Zhao, D. , Hasegawa, A. , & Hoiuchi, S. , 1992. Tomographic imaging of P and S wave velocity structure beneath Northeastern Japan, J. Geophys. Res. , 97, 19909-19928.
Zhao, D. , Hasegawa, A. , & Kanamori, H. , 1994. Deep structure of Japan subduction zone as derived from local, regional, and teleseismic events, J. Geophys. Res. , 99, 22313-22329.
Zhao, D. .Kanamoi, H. , & Humphreys, E. , 1996. Simultaneous inversion of local and teleseismic data for the crust and mantle structure of southern California , phys. Earth planet. Inter.,93, 191-214.
Zhu, H. & Ebel, J. E. , 1994. Tomographic inversion for the seismic velocity structure beneath northern New England using seismic refraction data, J. Geophys. Res. , 99, 15, 331-15, 357.
Zielhuis, A & yan der hilst, R. D. , 1996. Upper-mantle shear velocity beneath easten Australia from inversion of waveforms from SKIPPY portable arrays, Geophys. J. Int. ,127, 1-16