首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
TAN Hao-yuan  WANG Zhi 《地震地质》2019,41(6):1366-1379
3-D VP and VS images of southern Philippines at the 0~100km depths are generated by inverting a large number of travel-time data from the International Seismological Centre(1960-2017)through seismic tomography method. The results show lateral variation exists in the crust and upper mantle:High VP and VS anomalies emerge in mid-west Mindanao and Bohol Island, which might be caused by the combined action of huge magmatism and ophiolite accretion in the lower crust; low velocity anomalies of the upper mantle in the west of Mindanao are consistent with locations of volcanoes on the surface. It, thus, could be inferred that the low velocity anomaly is closely related to magmatic activity. The dense earthquake distribution along plate margin extending to 100km coincides with the strong activity of the Philippine Sea Plate which is located in the northeast and southeast of Mindanao. Relative weak activity of Sulawesi Sea Basin is presented simultaneously. The subduction of the Philippine Sea Plate is mostly concentrated in the crust and the top of the uppermost mantle. Our tomographic images show that lateral heterogeneities exist in the crust and uppermost mantle of the southern Philippines. Low VP and VS anomalies emerge in Philippine Trench and Cotabato Trench, in contrast, high VP and VS anomalies appear in shallow crust of land area where a large number of earthquakes and magmatic activities develop. This may reflect strong tectonic processes between the Philippine Sea Plate and Philippine Mobile Belt. Low VP and VS anomalies in the crust of eastern Mindanao coinciding with the location of volcanoes on the surface may show partial melting of crust material caused by dehydration of the subducting Philippine Sea Plate. Such a similar phenomenon can be also seen in the south of Negros Island and Cotabato Trench. Thus we infer that active tectonic behaviors are constrained within the crust of the Philippine Sea Plate, Sulu Sea Basin and Sulawesi Basin.Low VP and VS anomalies of the mantle in the mid-west of Mindanao island are associated with magmatic activity which may be caused by a collision between the east and west part of Mindanao at 5Ma. The fracture system in the west of Mindanao provides the possible passage ways of mantle hot material upwelling, coinciding with the model of geothermal distribution in this area. According to the geochemical analysis, ophiolite observed in Sanbaoyan and the western part of Mindanao could indicate material composition from crust to upper mantle on Eurasian continental margin which may show the evidence of rapid expansion environment of mid-ocean ridge. High VP and VS anomalies in the mantle of northeast and southeast of Mindanao coinciding with the distribution of massive earthquake along boundaries show a well agreement with the shape of the Philippine Sea Plate. Dense earthquake distribution in south Mindanao at 100km shows the Philippine Sea Plate has strong activity and stress accumulation in the upper mantle. On the contrary, the seismicity in southwest Mindanao and Cotabato Trench reduces rapidly at the depth from 50km to 100km, revealing weak subduciton and stress release of Sulawesi Basin in the mantle.  相似文献   

2.
An oceanic crustal model has been produced for the Nazca plate south of the Nazca Ridge prior to subduction into the Peru-Chile Trench at 18°S latitude. Consistent delays of thePn arrivals and a discontinuity in the tau-p curve indicate a low-velocity zone at the base of the crust. Observed upper mantle velocities are low; however, the mantle velocity increases with depth, at least to 20 km, to a value of 8.5 km/s. A possible petrological cause for the low-velocity zone is partially serpentinized peridotite; however, no clear refracted shear waves were observed to constrain this interpretation.  相似文献   

3.
渤海及其邻区的地震层析成像   总被引:40,自引:8,他引:40       下载免费PDF全文
由渤海周边地区地震台网记录的P波到时进行层析成像研究,得到了该地区地壳、上地幔速度图像的新信息.结果表明:1.渤海湾内地壳薄,是Moho面上隆的中心.2.渤张断裂由几段相接而成,属壳内断裂.3.该区岩石层厚度变化显著,渤海内为110km,燕山地区大于180km.渤海湾内上地幔盖层速度是研究区内最高的,最大值为8.45km/s,其岩石成分是高密度物质.这是该区呈现高重力异常的缘由.4.渤海内庙岛西部中地壳有一近似圆状的低速体,渤海内强震以及近期地震活动沿其周边展布,那里是高、低速块体的交界部位.研究区内其它地区的强震一般也发生在地壳中高速块体的边侧,或高速块体与低速块体的交界位置.  相似文献   

4.
本文利用30个基准台所记录的238条长周期面波资料,经过适配滤波和分格频散反演,得到中国大陆及邻区147个分格10-105s的纯路径频散,进而反演出青藏高原及邻近地区深至170km的剪切波三维速度结构.研究表明,青藏高原中西部地区和东部地区的地壳平均厚度分别为70±7km和65±7km,地壳平均剪切波速度分别为3.55和3.62km/s,上地幔顶盖平均速度分别为4.63和4.61km/s; 岩石层厚度均为120±10km;东部地区下地壳内30-40km深度处普遍存在低速层;青藏高原及其东侧的上地幔低速层内有横贯东西且明显向上隆起的低速腔.滇西缅北地区的地壳厚45±5km,上地壳及下地壳内都有低速层;上地幔顶盖的速度为4.42km/s,比青藏高原本体及恒河平原都低.恒河平原地壳厚34±2km,速度平均为3.45km/s;上地幔顶盖厚86±10km,速度平均为4.63km/s,顶盖内55-83km深处有一个低速夹层.  相似文献   

5.
谭皓原  王志 《地球物理学报》2018,61(12):4887-4900
菲律宾群岛受到欧亚板块、菲律宾海板块和印度-澳大利亚板块的碰撞作用,地质环境复杂,构造因素多样.尽管近几年来已经有了少数关于该区域层析成像的研究,但这些研究的区域主要集中在马尼拉海沟、吕宋岛及中菲律宾地区,而关于群岛周围其他海沟和南菲律宾地区的讨论相对较少.到目前为止,还没有同时获得过关于菲律宾群岛深部纵、横波速度结构的研究,本次研究通过反演155779条P波震相和59642条S波震相,同时获得了菲律宾群岛从地表至150 km深度的纵、横波速度结构.地震层析成像结果表明该地区的壳幔速度结构具有较强的不均一性,地壳内部存在着广泛的低速异常,而表征俯冲板块的高速异常则沿着群岛周边的海沟展布.南海块体在马尼拉海沟中段的俯冲角度和俯冲活动性比南段小;菲律宾海板块在东吕宋海槽南段微弱的俯冲作用很有可能同本哈姆海台的碰撞有关.菲律宾群岛大部分MW>6.0的强震沿着各个板块的边界发生,体现出菲律宾海板块同欧亚板块之间的强耦合作用,群岛西侧的南海块体在马尼拉海沟16°N-20°N之间呈现出的弱耦合状态可能跟北吕宋地区的拉张应力环境有关,南海块体在16°N以南的地区同上覆块体之间的耦合作用较强;群岛东侧的菲律宾海板块在14°N以北的地区没有强震发生,它与菲律宾群岛之间的耦合程度从北向南逐渐增强,在12°N以南的地区要强于12°N以北的地区;此外苏禄海盆和菲律宾构造带之间也存在着强耦合关系.  相似文献   

6.
中国东部海域岩石圈结构面波层析成像   总被引:11,自引:5,他引:6       下载免费PDF全文
本文通过面波层析成像得到了中国东部海域及邻近地区的地壳上地幔S波速度图像,给出了主要构造单元的区划及其结构特征,并讨论了速度结构与现今构造活动及构造演化历史的关系.研究区内中下地壳的平均速度与地震活动存在比较显著的关系,强震基本都发生在低速区内或高低速过渡区.太行山以东地壳内存在几条北西向低速带,其中张家口—渤海地震带下方的低速带最为显著.东部海域划分成北黄海、南黄海、东海、和冲绳海槽等4个构造块体.北黄海具有较薄较高速的岩石圈,与南华北盆地类似,推测是中生代特提斯洋向北俯冲造成岩石圈减薄的遗迹.北华北地区具有低速的地壳和较厚的岩石圈,岩石圈地幔速度偏低且上下比较均匀,可能反映中生代沿北方缝合带持续碰撞作用的特点.南黄海具有相对较厚的岩石圈,较多地保存了下扬子克拉通的特征.在下扬子与华北地块的拼合过程中,洋壳俯冲可能是北黄海和苏皖地区上地幔低速特征的成因.在125°E以东的朝鲜半岛地区未发现这一拼合过程的遗迹.有可能整个朝鲜半岛都是华北地块的一部分;但也有可能是太平洋俯冲和日本海张开的作用完全改造了朝鲜半岛的岩石圈上地幔,抹去了以往构造运动的痕迹.东海地区的地壳厚度,特别是岩石圈厚度向冲绳海槽方向减小,反映出菲律宾海板块俯冲在弧后广大地区都有影响.冲绳海槽地区可见俯冲的菲律宾海板片以及板片上方显著低速的地壳和上地幔,为冲绳海槽的弧后扩张机制提供了证据.  相似文献   

7.
本文是1986年古雷—石城剖面及嵩口—宜城剖面深地震测深资料的初步研究结果。 对古雷—石城的纵剖面资料,分析了震相特征,共识别出五个波组:P_2、P_3~0、P_4~0、P_5~0及P_n(P_n~0)。通过对波的走时反演,正演拟合和理论地震图方法等计算,得到了该区地壳与上地幔结构模型。 古雷—石城地区地壳具有多层结构,并可划分为上、中、下三层。古雷炮点给出的厚度分别为1.0km、15.7km、12.8km,地壳平均速度为6.29km/s,深度为29.5km,上地幔顶面P_n波速度为7.83km/s。石城炮点给出厚度分别为1.8km、18.3km、12.4km。地壳平均速度为6.29km/3,深度为32.4km,土地幔顶面P_n速度为8.00km/s。 在中地壳下部存在一低速层,其厚度为2.8km,速度为5.85km/s。根据其它研究结果,初步判断低速层介质是半熔融物质组成。 测区内横向变化比较强烈。从东向西有长乐—诏安、政和—海丰和邵武—河源三个大断裂穿过该区,并且都深切至莫霍面;在漳州盆地之下莫霍面隆起约3km,戴云山区下莫霍面凹陷近2km;永安—梅州莫霍面隆起接近3km。莫霍面分布显示出从东南向西北逐渐加深。 宜城—连城—嵩口非纵剖面显示了莫霍面在两处有明显断错,错距约2km邵。表明昭武—河源断裂是切割莫霍面的深大断裂。  相似文献   

8.
Tetsuo  Irifune 《Island Arc》1993,2(2):55-71
Abstract Phase transformations in model mantle compositions and those in subducting slabs have been reviewed to a depth of 800 km on the basis of recent high-pressure experimental data. Seismic velocity and density profiles in these compositions have also been calculated using these and other mineral physics data. The nature of the seismic velocity and density profiles calculated for a pyrolite composition was found to generally agree with those determined by seismic observations (e.g. PREM). The locations of the seismic discontinuities at 400 and 670 km correspond almost exactly to the depths where the transformations of the olivine component to denser phases take place. Moreover, the steep gradients in the seismic velocity/density profiles observed between these depths are qualitatively consistent with those expected from the successive transformations in the complementary pyroxene-garnet component in the pyrolite composition. Further, the calculated seismic velocity and density values agree well with those observed in the upper mantle and mantle transition region within the uncertainties attached to these calculations and observations. Pyrolite or peridotite compositions are thus most likely to represent the composition of the mantle above 670 km depth, although some degrees of chemical heterogeneity may exist in the transition region. The observed sharp discontinuous increases of seismic velocities and density at this depth may be attributed either to the phase transformation to a perovskite-bearing assemblage in pyrolite or to chemical composition changes. Density profiles in subducted slabs have been calculated along adequate geotherms assuming that the slabs are composed of the former oceanic crust underlain by a thicker harzburgitic layer. It is shown that the former oceanic crust is substantially less dense than the surrounding pyrolite mantle at depths below 670 km, while it is denser than pyrolite in the upper mantle and the transition region. The subducted former oceanic crust may be trapped in this region, forming a geochemically enriched layer at the upper mantle-lower mantle boundary. Thick and cool slabs may penetrate into the lower mantle, but the chemically derived buoyancy may result in strong deformation and formation of megalith structures around the 670 km seismic discontinuity. These structures are consistent with those detected by recent seismic tomography studies for subduction zones.  相似文献   

9.
We investigated the physical properties in the upper mantle beneath the Philippine Sea using a theoretical relation derived by Karato [Mapping water content in the upper mantle. Subduction factory, AGU Monograph, in press]. From the attenuation model of Shito and Shibutani [Phys. Earth Planet. Interact., in press] and the velocity model of Widiyantoro et al. [Earth Planet. Sci. Lett. 173 (1999) 91], observed attenuation and velocity anomalies were evaluated to explain the temperature, water content, and chemical heterogeneities in the target area. The results indicate that the observed anomalies in the shallower regions (50–200 km) may be due to chemical composition effects (e.g., concentration of iron), in addition to the temperature and water content anomalies. In contrast, for the deep upper mantle (300–400 km), the observations can be explained by only the effects of high water content (10–50 times higher than the average mantle). These inferred properties of the mantle are consistent with the tectonic history of the Philippine Sea region, which has had a long history of subduction and active magmatism.  相似文献   

10.
We have constrained the shear-wave structure of crust and upper mantle beneath Iceland by analyzing fundamental mode Rayleigh waves recorded at the ICEMELT and HOTSPOT seismic stations in Iceland. The crust varies in thickness from 20 to 28 km in western and northern Iceland and from 26 to 34 km in eastern Iceland. The thickest crust of 34–40 km lies in central Iceland, roughly 100 km west to the current location of the Iceland hotspot. The crust at the hotspot is ∼32 km thick and is underlain by low shear-wave velocities of 4.0–4.1 km/s in the uppermost mantle, indicating that the Moho at the hotspot is probably a weak discontinuity. This low velocity anomaly beneath the hotspot could be associated with partial melting and hot temperature. The lithosphere in Iceland is confined above 60 km and a low velocity zone (LVZ) is imaged at depths of 60 to 120 km. Shear wave velocity in the LVZ is up to 10% lower than a global reference model, indicating the influence of the Mid-Atlantic Ridge and the hotspot in Iceland. The lowest velocities in the LVZ are found beneath the rift zones, suggesting that plume material is channeled along the Mid-Atlantic Ridge. At depths of 100 to 200 km, low velocity anomalies appear at the Tjornes fracture zone to the north of Iceland and beneath the western volcanic zone in southwestern Iceland. Interestingly, a relatively fast anomaly is imaged beneath the hotspot with its center at ∼135 km depth, which could be due to radial anisotropy associated with the strong upwelling within the plume stem or an Mg-enriched mantle residual caused by the extensive extraction of melts.  相似文献   

11.
Summary Measurements of phase velocities along several paths in southeastern Europe and along one path in the Eastern Mediterranean have been made in the period range 15 to 60 sec for the fundamental mode of Rayleigh waves. The data are fitted by layered models of the crust and upper mantle. The crust in the Balkan peninsula is thicker than the normal continental crust by about 10 km while the crust in the northern Aegean Sea has the normal continental thickness. The existence of a low shear velocity layer of thickness equal to 140 km with its top in a depth of 90 km interpretes the data well for the longer periods. The data for the eastern Mediterranean show that this region has an oceanic character. The mean crustal thickness in this area is of the order of 20 km.  相似文献   

12.
Cross-correlation functions of noise are constructed on 119 interstation paths from seismic noise records at stations of Eastern Europe. Dispersion curves of the group velocity of Rayleigh waves obtained from the cross-correlation functions are used for constructing the three-dimensional distribution of the velocity of transverse waves on the East European platform and in adjacent regions by methods of surface-wave tomography. The mean velocity in the crust is minimum in the region of the Caspian depression and Black Sea basin (<3.3 km/s) and maximum in the Baltic shield area (>3.7 km/s). The upper mantle beneath the Baltic and Ukrainian shields is characterized by increased velocity and the absence of the asthenospheric layer. Reduced velocities are noted in the upper mantle of the Black Sea basin. A low-velocity anomaly in the shape of a vertical column is revealed at depths of 200–300 km in the central part of the Dnieper-Donets aulacogen, which confirms the existence of a paleorift in this region.  相似文献   

13.
The paper presents a review and analysis of new seismic data related to the structure of the mantle beneath the East European platform. Analysis of observations of long-range profiles revealed pronounced differences in the structure of the lower lithosphere beneath the Russian plate and the North Caspian coastal depression. The highest P-velocities found at depths around 100 km are in the range 8.4–8.5 km s?1. Deep structure of the Baltic shield is different from the structures of both these regions. No evidence of azimuthal anisotropy in the upper mantle was found. A distribution of P-velocity in the upper mantle and in the transition zone consistent with accurate travel-time data was determined. The model involves several zones of small and large positive velocity gradients in the upper mantle, rapid increases of velocity near 400 and 640 km depths and an almost constant positive velocity gradient between the 400 and 640 km discontinuities. The depth of the 640 km discontinuity was determined from observations of waves converted from P to SV in the mantle.  相似文献   

14.
We investigated the upper mantle anelastic structure beneath the northern Philippine Sea region, including the Izu-Bonin subduction zone and the Shikoku Basin. We used regional waveform data from 69 events in the Pacific and the Philippine Sea slabs, recorded on F-net and J-array network broadband stations in western Japan. Using the S–P phase pair method, we obtained differential attenuation factors, δt*, which represent the relative whole path Q. We conducted a tomographic inversion using 978 δt* values to invert for a fine-scale (50–100 km) three-dimensional anelastic structure.

The results shows two high-Q regions (QP>1000) which are consistent with the locations of the Pacific and the Philippine Sea slabs. Also there is a low-Q (QP110) area extending to the deeper parts (350–400 km) of the model just beneath the old spreading center and the Kinan Seamount Chain in the Shikoku Basin. A small depth dependence of the laterally averaged QP was found, with values of 266 (0–250 km), 301 (250–400 km), and 413 (400–500 km).  相似文献   


15.
Seismic data recorded in the upper mantle triplication distance range between 10° and 30° are generated by wave propagation through complex upper mantle structure. They can be used to place constraints on seismic velocity structures in the upper mantle, key seismic features near the major discontinuities, and anisotropic structure varying with depth. In this paper, we review wave propagation of the upper mantle triplicated phases, how different key seismic features can be studied using upper mantle triplicated data, and the importance of those seismic features to the understanding of mantle temperature and composition. We present two examples of using array triplicated phases to constrain upper mantle velocity structures and detailed features of a certain discontinuity, with one for a shallow event and the other for deep events. For the shallow event, we present examples of how the array triplication data can be used to constrain several key properties of the upper mantle: existence of a lithospheric lid, existence of a low velocity zone beneath the lithospheric lid, and P/S velocity ratio as a function of depth. For deep events, we show examples of how array triplication data can be used to constrain the detailed structures of a certain discontinuity: velocity gradients above and below the discontinuity, velocity jumps across the discontinuity and depth extents of different velocity gradients. We discuss challenges of the upper mantle triplication study, its connection to other approaches, and its potential for further studying some other important features of the mantle: the existence of double 660-km discontinuities, existence of low-velocity channels near major discontinuities and anisotropy varying with depth.  相似文献   

16.
In this study, three receiver function stacking methods are used to study the detailed crust and upper mantle structure beneath south-central Alaska. We used teleseismic waveform data recorded by 36 stations in the Broadband Experiment Across the Alaska Range (BEAAR) and 4 permanent stations in Alaska. H − κ stacking method using P-to-S converted wave and its multiply reflected waves between the Earth's surface and the Moho discontinuity is adopted to estimate the crustal thickness (H) and average crustal VP/VS ratio (κ) in this region. The receiver function results for 24 stations show that the crustal thickness under Alaska ranges from 26.0 to 42.6 km with an average value of 33.8 km, and the VP/VS ratio varies from 1.66 to 1.94 with an average value of 1.81 which corresponds to an average Poisson's ratio of 0.277 with a range from 0.216 to 0.320. High Poisson's ratios under some stations are possibly caused by partial melting in the crust and the uppermost mantle. Common converted point (CCP) stacking results of receiver functions along three lines show clear Moho and slab images under this subduction zone. The depths of the slab from our CCP stacking images are consistent with those estimated from the Wadati–Benioff Zone (WBZ). In the area between two stations DH2 (147.8°W, 63.3°N) and DH3 (147.1°W, 63.0°N), a Moho depth offset of about 10 km is found by both the H − κ and CCP stacking techniques. Common depth point (CDP) stacking of receiver functions shows not only the 410-, 520- and 660-km discontinuities, but also significant variations (−30 to 15 km) in the transition zone thickness under the southwest and southeast parts of the study region. The transition zone becomes thinner by 20–30 km, indicating that the temperature there is 150–200 K higher than that of the normal mantle.  相似文献   

17.
中国大陆及邻近海域的Rayleigh波群速度分布   总被引:11,自引:5,他引:11       下载免费PDF全文
利用Rayleigh波群速度资料反演得到中国大陆及其临近海域的(70°E-145°E,10°N-55°N)15-120s周期的群速度分布图像. 塔里木盆地在15s处清楚地显示为低速,在16-33s左右没有显示,但在36-5s显示为高速,说明塔里木盆地有较深的根. 青藏高原块体是44s至120s图像中最为突出的低速块体,南面与印度板块的分界线以及与北面的塔里木盆地、柴达木盆地的分界清晰,其块体中西部的速度低于东部. 泰国清迈附近存在一尺度为1000km左右的低速带,可能是青藏高原块体的物质向东南方向迁移造成上地幔物质上涌的结果. 南北地震带表现为强烈的速度梯度带,西面为低速,东面为高速. 中国南海的中央、日本海中央、菲律宾海表现为海洋性地壳. 菲律宾海的图像与地形及地震带完全吻合. 环绕菲律宾海及日本海存在400km左右宽的低速带,可能是岩浆活动带.  相似文献   

18.
本文描述随县—安阳剖面深地震测深资料分析解释的初步结果。资料处理结果说明,本地区的莫霍界面埋深约28至36km,其上覆盖层(即地壳)的平均速度约6.24km/s,在新郑和淇县一带该界面埋藏较深。在地表附近沉积层速度约3.50km/s左右,在黄河地区该层厚度最大,约5.5km。而在测线南端的大片地区花岗岩出露地表。该沉积层的下侧是结晶基底,其顶面的首波速度为6.00km/s。在它与莫霍面之间还可分出三层,其相应震相为P2,P03P04,其中P03较为稳定,它是地壳内部的反射波震相,平均速度约5.97km/s,深度为18至23km。本文得出了这几层的深度变化剖面图。并发现在新郑与新乡间黄河流经的地区,地壳内速度的横向变化及梯度层的分布情况有其明显特点。从对其他震相的分析,还获得新乡以南沉积层内的速度分布及舞阳附近的断层位置,并表明在本区上地幔内存在高速夹层的可能性。   相似文献   

19.
利用我国地震台站记录的瑞利波观测资料,通过适配滤波频时分析技术进行数据处理,获得了穿越我国东南及陆缘地区的瑞利波频散.使用随机反演理论取得了东南大陆及陆缘地带4°×4°网格的纯路径频散数据.在网格反演的基础上使用Harkrider的面波反演程序求得了该区剪切波的三维速度结构.结果表明:1.华南大陆Moho界面埋深为30-40km,并由西向东逐渐减薄,在陆缘与浅海地域为25-28km,具有明显的分区特征.2.上地幔低速层埋深为60-0km,变化幅度较大,这与深部断裂分布及深层过程有关,但NS向剖面上各界面的起伏变化均比EW向剖面平缓.3.东南陆缘是东亚大陆的海陆过渡带,在深部表现为Moho界面埋深和地壳平均速度降低的地带,地幔深部界面的起伏形态充分表明,深浅介质结构和物质耦合的不均匀性.  相似文献   

20.
香港地区海陆地震联测及深部地壳结构研究   总被引:7,自引:0,他引:7       下载免费PDF全文
为了探明南海北部海陆过渡带的深部地壳结构,我们在香港外海域进行了一次海陆地震联洲的实验,利用固定地震台网远距离接收海上气枪信号,接收距离远达200多km,并利用此次实验的测线1剖面模拟得到了海陆过渡带的深部地壳速度结构.速度结构模型表明:研究区海陆过渡带的地壳结构非均匀性较明显,由陆至海沉积层有一个突然增厚的特点;莫霍面深度约为26~29 km,上地壳P波速度约为5.5~6.4 km/s,下地壳P波速度为6.5~6.9 km/s.在担杆列岛往海方向有一个低速破碎带,其上地壳P波速度为5.2~6.1 km/s,下地壳P波速度为6.2~6.4 km/s,结合野外地质调查的结果,推测它可能为滨海断裂带.在担杆列岛往陆方向香港和深圳之间的研究区域,莫霍面有较大起伏,可能与此处发育的海丰断裂有关.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号