首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
地球介质的衰减特性可以从地震波中的频谱信息中获取,并且能够提供地球内部结构非均匀性的重要信息.与大量的速度结构研究相比,对衰减结构的研究却较少,且大多仅对日本群岛局部地区做研究.本文运用ML振幅层析成像技术反演获得了整个日本群岛上地壳的品质因子Q0值(1 Hz时的Q值).计算中共利用了日本Hi-net台网971个台站记录到的5559个地震事件,超过60000 条Sg波最大振幅-频率数据,并且满足震中距小于2°、震源深度小于10 km.计算结果表明,低Q0值区域主要存在于日本群岛的中央火山带,与日本新生代以来的火山分布非常一致;同时,位于从神户到京都东北地区、以及四国岛东部等地的非火山区,断层密集,对应的Q0值也很低;另外,西太平洋俯冲海沟与日本群岛东海岸之间的区域也是低Q0值区,这说明高温、岩石强烈破碎、深厚的沉积物是造成地震波强烈衰减的主要原因.高Q0值主要分布在中央火山带前缘与日本东海岸之间.  相似文献   

2.

青藏高原东南缘的龙门山断裂两侧具有陡峭的地形特征,在约50~100 km的水平距离内,地形高程从2000 m增加到4000 m,该区强烈的壳幔变形特征及地球动力学模式一直是研究的热点问题.本文从四川地区49个固定台站记录的远震资料提取了P波接收函数,获得了四川盆地及周边的地壳厚度和泊松比,并以此构建反演的初始模型.在线性反演的基础上,引入了分别拟合低频和高频接收函数的两步反演技术,用以反演台站下方的地壳S波速度结构.数字试验表明,该方法可以有效抑制接收函数反演的不唯一性,为了得到最优解,最后用Bootstrap重采样技术估计解的不确定性.结果表明,四川盆地的地壳厚度在40~46 km,松潘-甘孜块体北部的地壳厚度为46~52 km,而南部增厚到50~60 km.从四川盆地向西跨过龙门山断裂,地壳厚度增加了10~15 km.在四川盆地及周边地区,地壳泊松比在0.26~0.32之间,呈块体分布特征,高泊松比(0.28~0.32)主要沿龙门山断裂以及安宁河-小江断裂分布.地壳S波速度结构表明,来自青藏高原中部的中下地壳低速层可能受到了坚硬的四川盆地阻挡,改变原来的运动方向并沿龙门山断裂展布,由于低速层的囤积导致该区地形陡峭和下地壳增厚.

  相似文献   

3.
Aso Volcano experienced a huge pyroclastic eruption 90 thousand years ago, and formed a large caldera (18 km × 25 km). In order to test the hypothesis of a magma body in the mid and lower crust that has been suggested geophysically and geochemically, we investigated seismic velocity discontinuities and velocity structure beneath Aso Caldera using receiver functions and a genetic algorithm inversion. We confirm the existence of the Moho at depths between 30 km and 35 km and a large velocity anomaly should exist in the deep portion of the crust beneath Aso Caldera, from imaging of receiver functions observed only at stations outside the caldera. As a result of a more detailed examination with GA inversion, a low velocity layer is detected at depths between 10 km and 24 km beneath the western part of the caldera. S-wave velocity of the layer is estimated to be 2.0–2.4 km/s. We estimate that the low velocity layer contains at most 15% melt or 30% aqueous fluid. The layer exists near the Conrad and at the same depths as the swarm of the low frequency earthquakes and a compressional and dilatational deformation source which are expected to be caused by fluid movement beneath the middle-eastern part of the caldera. Fluid contained in the layer might be related with huge pyroclastic eruptions of Aso Volcano.  相似文献   

4.
The Tohoku megathrust earthquake, which occurred on March 11, 2011 and had an epicenter that was 70 km east of Tohoku, Japan, resulted in an estimated ten's of billions of dollars in damage and a death toll of more than 15 thousand lives, yet few studies have documented key spatio-temporal seismogenic characteristics. Specifically, the temporal decay of aftershock activity, the number of strong aftershocks (with magnitudes greater than or equal to 7.0), the magnitude of the greatest aftershock, and area of possible aftershocks. Forecasted results from this study are based on Gutenberg-Richter's relation, Bath's law, Omori's law, and Well's relation of rupture scale utilizing the magnitude and statistical parameters of earthquakes in USA and China (Landers, Northridge, Hector Mine, San Simeon and Wenchuan earthquakes). The number of strong aftershocks, the parameters of Gutenberg-Richter's relation, and the modified form of Omori's law are confirmed based on the aftershock sequence data from the MW9.0 Tohoku earthquake. Moreover, for a large earthquake, the seismogenic structure could be a fault, a fault system, or an intersection of several faults. The seismogenic structure of the earthquake suggests that the event occurred on a thrust fault near the Japan trench within the overriding plate that subsequently triggered three or more active faults producing large aftershocks.  相似文献   

5.
This study seeks to construct a hazard function for earthquake probabilities based on potential foreshocks. Earthquakes of magnitude 6.5 and larger that occurred between 1976 and 2000 in an offshore area of the Tohoku region of northeast Japan were selected as events for estimating probabilities. Later occurrences of multiple events and aftershocks were omitted from targets. As a result, a total of 14 earthquakes were employed in the assessment of models. The study volume spans 300 km (East-West) × 660 km (North-South) × 60 km in depth. The probability of a target earthquake occurring at a certain point in time-space depends on the number of small earthquakes that occurred per unit volume in that vicinity. In this study, we assume that the hazard function increases geometrically with the number of potential foreshocks within a constrained space-time window. The parameters for defining potential foreshocks are magnitude, spatial extent and lead time to the point of assessment. The time parameter is studied in ranges of 1 to 5 days (1-day steps), and spatial parameters in 20 to 100 km (20-km steps). The model parameters of the hazard function are determined by the maximum likelihood method. The most effective hazard function examined was the following case: When an earthquake of magnitude 4.5 to 6.5 occurs, the hazard for a large event is increased significantly for one day within a 20 km radius surrounding the earthquake. If two or more such earthquakes are observed, the model expects a 20,000 times greater probability of an earthquake of magnitude 6.5 or greater than in the absence of such events.  相似文献   

6.

板内火山活动是板内下方深部物质组成的重要视窗,而中国东北是研究板内火山的热点场所.中国东北火山活动主要沿大兴安岭及松辽盆地周缘分布,隐匿在大兴安岭内的各火山群落鲜有地球物理观测.本研究聚焦于大兴安岭北部的诺敏河火山群,利用围绕诺敏河火山布设的43个流动台站记录到的17个远震事件数据,通过时间域波形匹配法计算了300条直达P波震相的Δt*,并进一步采用贝叶斯蒙特卡洛方法反演得到了研究区的二维Δt*模型,为约束诺敏河火山地区的深部结构提供了新的观测证据.结果显示研究区高衰减区域与火山活动在空间位置上有很好的一致性.诺敏河火山与研究区北端临近小古里河火山区域均观测到高Δt*值,可能由区域下方热地幔物质上涌导致.诺敏河火山东部、科洛河火山西部区域观测到低Δt*值,与高衰减区域最大Δt*差值达约0.3±0.05 s,可能由两部分原因导致:该区域下方残留未被侵蚀的岩石圈,以及局部区域或存在的部分熔融导致的残留地幔矿物脱水.

  相似文献   

7.
2011年3月11日日本东北近海9.0级地震及伴发海啸造成了重大的人员伤亡和经济损失.本文介绍了这次灾害的特点和震后应急救灾行动,总结了此次应急的经验和不足,并对我国防震减灾工作提出了几点启示.  相似文献   

8.
李艳娥  陈学忠 《地震》2017,37(4):10-21
2011年3月11日, 一个 MW9.1地震袭击了日本本州地区, 为了分析这次地震前后主震破裂区内应力时空变化, 我们选取1996年1月~2016年6月期间发生在破裂区内的563个5.0≤MS≤6.9地震, 研究了视应力随时间的变化和空间分布。 日本MW9.1地震前从2002年中起视应力开始呈趋势性上升变化, 到2009年初以0.18 MPa/a的速率从0.6 MPa上升到1.76 MPa, 相差约3倍, 直到地震发生前夕一直保持在1.5 MPa之上。 地震发生之后, 直到2016年6月在破裂区内视应力呈缓慢下降变化, 但仍保持在1.5 MPa之上较高水平。 视应力在地面上和断层面上的分布显示, 1996—2005年间破裂区仅存在个别视应力高值, 从2006年到2011年2月, 破裂区大面积出现视应力高值。 在日本MW9.1地震发生之后的近3个月内, 破裂区视应力整体处于高值水平, 之后在较高的水平上缓慢减弱。 视应力是地震断层面上平均应力的下限, 视应力的高低在一定程度上反映的是震源断层面上平均应力的高低。 在日本MW9.1地震前, 发生在破裂区内的地震, 其断层面上的平均应力经历了大约8.5年的趋势上升变化过程。 这次大地震前破裂区所在的地壳应力逐渐增加, 最后导致断层面错动发生日本MW9.1地震。  相似文献   

9.
张冰  符力耘  魏伟  管西竹 《地球物理学报》2014,57(10):3373-3388
澳大利亚西北大陆架卡拉汶盆地含有丰富的油气资源, 但其地质条件复杂, 普遍存在异常地层高压分布,特别是该区区域盖层-中深层巨厚的Muderong页岩层内存在局部异常高压带,给石油钻井带来巨大困难.本文根据卡拉汶盆地北部地区10余口井的实测地层压力资料、声波测井数据和联井地震剖面, 研究该页岩层内异常压力千赫兹尺度声波响应的衰减特征,及其与地层有效压力之间的关联计算模型.通过联井地震数据的井-震相关性分析进行井-震过渡,研究井旁地震数据在异常高压带的衰减特征,及其与地层有效压力之间的关联计算模型,从而实现井震资料联合进行地层压力地震外推反演.研究表明,Muderong页岩层地层有效压力与井中声波速度/阻抗品质因子和井旁地震阻抗品质因子具有良好的正向相关性,根据二者交汇数据拟合经验公式反算得到的各井有效压力结果基本反映出实测有效压力曲线的变化特征.卡拉汶盆地北部地区地震数据具有良好的井-震相关性,地震资料品质相对可靠,确保了声-压相关分析的井-震过渡.最后,通过地震阻抗品质因子与地层有效压力之间的关联计算模型,进行该区异常地层压力地震外推反演,预测结果与钻井证实的异常压力分布特征基本符合.本研究为卡拉汶盆地中深部区域盖层局部异常高压带的地层压力外推反演提供了一种有效途径.  相似文献   

10.
利用波形互相关方法识别乳山震群中的重复地震,挑选一组在时间上跨越2015年5月22日ML 5.0地震的重复地震序列(repeating earthquake sequences,RES),利用测算单台直达P波第一个周期信号频散衰减特性的方法,估算乳山震群震源区介质品质因子,结果显示,不同路径下Qm值随时间的变化形态基本一致,即此次ML 5.0地震前Qm值升高,临震及震后Qm值下降并恢复至震前水平。同时,通过对震源区近场及远场应力状态的相关性分析,认为此次ML 5.0地震发生在远场作用下区域构造应力状态调整过程中,与区域构造的介质状态、断层分布、速度结构密切相关。  相似文献   

11.

收集华北克拉通地区188个宽频带流动台站观测资料进行处理.通过背景噪声面波数据和接收函数双重资料约束联合反演,得到了研究区沉积层厚度、地壳厚度及地壳S波速度结构.结果显示:(1)沉积盖层厚度与地质构造相对应,盆地区与隆起区分界明显.(2)研究区地壳厚度变化范围约29~46 km,自西向东逐渐变薄.(3)中、上地壳华北盆地S波速度偏高,可能与新生代以来多次沉降所造成的相对高的岩石强度有关;(4)下地壳S波速度显示研究区主要存在三个低速区,分别是唐山—天津周边、张北及太行山造山带地区;华北盆地存在显著高速异常,推测可能是由于华北盆地经历下地壳拆沉后,大规模的伸展作用相伴随的幔源基性铁镁质岩浆底侵至下地壳结晶所造成的.(5)多个发生过强震的区域表现出沉积层下方存在较大范围的(约10 km)高速体,并且高速体又被其下低S波速度包裹,壳内岩石强度的差异为应力积累及地震发生提供条件.

  相似文献   

12.
利用西沙琛航岛流动地震台站和永兴岛固定地震台站的资料,提取了远震P波接收函数,结合正演和反演方法模拟了台站下方的地壳结构.模拟结果显示:西沙群岛地壳顶部存在2 km厚的新生代低速沉积层,横波速度只有2.0~2.2 km/s;上地壳为一速度梯度带,横波速度由2 km处的3.4 km/s逐渐增加到12 km深度时的3.8 km/s;下地壳存在明显低速层,厚度达到12 km,平均横波速度3.5 km/s;莫霍面埋深26~28 km,也表现为一速度梯度带,横波速度从3.8 km/s变化到4.6 km/s左右,并保持稳定;该地区的地壳泊松比值大于0.3,推测西沙群岛的壳内低速层和异常泊松比值与地幔热活动引起的韧性流变构造和岩石矿物的各向异性排列有关.  相似文献   

13.

地震激发的地球自由振荡中, 环型自由振荡(特别是低频环型自由振荡)的精确探测相比于球型自由振荡更加困难. 环型自由振荡对SH波速度结构更加敏感, 能更好地监测走滑断层运动, 因此精确估算环型自由振荡的谱峰参数可为地球内部大尺度SH波速度结构和大地震的震源机制提供更强约束. 本研究利用2011年日本Tohoku地震后四分量钻孔应变仪的观测资料, 提取了该地震激发的1.3 mHz以下低频环型和球型自由振荡信号, 与STS-1宽频带地震仪的结果相比较, 发现四分量钻孔应变仪对1 mHz以下的环型自由振荡具有更强的探测能力. 本研究还采用频率域自回归方法精确估计了0T20T6的简并频率, 结果与PREM模型理论值相差非常小, 平均相对误差仅为2.52×10-4.此外还用钻孔应变仪观测资料探测到了0T2谱峰分裂现象. 研究结果表明, 钻孔应变仪能有效探测到大地震激发的低频地球环型自由振荡信号, 可为地幔大尺度各向异性结构研究提供高质量的基础数据.

  相似文献   

14.
In order to investigate crustal structure beneath the eastern Marmara region, a seismic refraction survey was conducted across the North Anatolian Fault (NAF) zone in north west Turkey. Two reversed profiles across two strands of the NAF zone were recorded in the Armutlu Highland where a tectonically active region was formed by different continents. We used land explosions in boreholes and quarry blasts as seismic sources. A reliable crustal velocity and depth model is obtained from the inversion of first arrival travel times. The velocity-depth model will improve the positioning of the earthquake activities in this active portion of the NAF. A high velocity anomaly (5.6–5.8 km s−1) in the central highland of Armutlu block and the low velocity (4.90 km s−1) pattern north of Iznik Lake are the two dominant features. The crustal thickness is about 26 ± 2 km in the north and increases to about 32 ± 2 km beneath the central Armutlu block in the south. P-wave velocities are about 3.95 km s−1 to 4.70 km s−1 for the depth range between about 1 km and 5 km in the upper crust. The eastern Marmara region has different units of upper crust with velocities varying with depth to almost 8 km. The high upper crust velocities are associated with Armutlu metamorphic rocks, while the low velocity anomalies are due to unconsolidated sedimentary sequences. The western side of Armutlu block has complex tectonics and is well known for geothermal sources. If these sources are continuous throughout the portions of the crust, it may be associated with a granitic intrusion and deformation along the NAF zone. That is, the geothermal sources associated with the low velocity may be due to the occurrence of widespread shear heating, even shear melting. The presence of shear melting may indicate the presence of crustal fluid imposed by two blocks of the NAF system.  相似文献   

15.
利用2002~2003年中国地震局地质研究所台阵实验室以唐山大震区为中心布设的40个流动宽频带地震台站和首都圈数字台网的33个宽频带台站的远震数据,采用接收函数非线性反演方法得到其中72个宽频带台站下方60 km深度范围内的S波速度结构.根据得到的各台站下方地壳上地幔的S波速度结构,并综合刘启元等(1997)用接收函数非线性反演方法得到的延怀盆地15个宽频带流动台站下方的地壳上地幔S波速度结构模型,给出了39°N~41°N,114°E~119.5°E区域内沿不同走向、不同深度S波速度分布.由于综合了利用首都圈数字地震台网的宽频带台站以及流动地震台阵的观测数据,本文给出了较前人同类研究空间分辨率更好的结果.结果表明: (1)研究区的速度结构,特别是怀来以东的速度结构十分复杂.在10~20 km深度范围内,研究区地壳具有高速和低速异常块体的交错结构.研究区中上地壳速度结构主要被与张渤地震带大体重合的NW向高速条带和穿越唐山大震区的NE向高速条带所控制,而其中下地壳的速度结构主要为延怀—三河—唐山地区上地幔隆起所控制.(2)研究区内存在若干壳内S波低速体,它们主要分布在唐山,三河及延怀盆地等地区.在这些地区,壳内低速体伴随着壳幔界面的隆起和上地幔顶部速度结构的横向变化.(3)地表断层分布与地壳速度结构分区有较好的相关性,表明断层对不同块体有明显的控制作用.其中,宝坻断裂,香河断裂和唐山断裂均为超壳断裂.(4)首都圈内大地震的分布与壳内低速体及上地幔顶部的速度结构有密切关系.对于唐山大地震的成因,仅考虑板块作用引起的水平应力场是不够的,有必要充分重视由于上地幔变形引起的地壳垂直变形和上地幔物质侵入造成的热效应.  相似文献   

16.
用转换函数方法研究喜马拉雅地区速度结构   总被引:4,自引:2,他引:2       下载免费PDF全文
利用流动数字地震台网提供的三分量地震波形记录,应用转换函数及快速模拟退火算法对喜马拉雅山脉地区46个地震站下的地壳横波速度结构进行了反演,为进一步揭示青藏高原喜马拉雅山脉地区的动力学演化过程提供了可靠的地球物理证据.根据本文结果可清晰看到,喜马拉雅山脉地区作为当今地壳活动最活跃的地区,物质交换非常活跃,地下结构远远未达到平衡,地壳速度有很大差异,在板块边界处莫霍界面速度间断不是非常明显,自喜马拉雅南坡向高原腹地,地壳厚度大致从55 km增加到80 km;沿经度方向,莫霍面也有一定的起伏.通过研究得到另外一个证据是,在喜马拉雅的主中央逆冲断裂,由大陆碰撞产生的主要构造,其深度可能要大于80 km.  相似文献   

17.
本文利用宽频流动台阵记录的远震波形资料和接收函数波动方程叠后偏移方法,获得了华北克拉通东北部边界及其邻近地区的地壳和地幔转换带的间断面结构图像.结果显示研究区域的地壳厚度存在显著的横向变化:以南北重力梯度带为界,西北部的兴蒙造山带地壳较厚(~40 km),东南部的燕山带、松辽盆地和辽东台隆地壳明显较薄(30~35 km).这有可能反映,研究区南北重力梯度带两侧地壳在中-新生代区域构造伸展过程中经历了不同程度的改造和减薄.地幔转换带成像结果显示,研究区410 km和660 km间断面结构存在横向差异.经度121°E-122°E之间,上地幔底部出现双重间断面,深度分别为660 km和690 km.经度122.5°E以东(北黄海地区),410 km间断面有5~20 km幅度的下沉,660 km间断面有5~15 km幅度的抬升;该地区地幔转换带厚度相对全球平均偏薄10~20 km,指示着该地区较热的上地幔底部温度环境.我们认为太平洋俯冲板块可能停滞在研究区119°E-122°E经度范围的地幔转换带中,但未延伸至118°E以西;而俯冲板块在124°E以东可能局部穿透了上地幔底部而进入下地幔,同时引起小尺度的地幔对流,导致北黄海地区下地幔物质的上涌.  相似文献   

18.

通过对南北地震带北段区域所布设的676个流动地震台站观测资料进行处理,联合反演面波频散与接收函数数据,获得了研究区内地壳厚度、沉积层厚度的分布情况以及地壳上地幔高分辨率S波速度结构成像结果.反演结果显示研究区地壳厚度从青藏高原东北缘向外总体逐渐变薄,秦岭造山带地壳厚度较同属青藏高原东北缘的北祁连块体明显减薄;鄂尔多斯盆地及河套盆地分布有非常厚的沉积层,阿拉善块体部分区域也有一定沉积层分布,沉积层与研究区内盆地位置较为一致;松潘—甘孜块体、北祁连造山带等青藏高原东北缘总体表现为S波低速异常;在中下地壳,松潘—甘孜块体下方的低速体比北祁连造山带下方的低速体S波速度值更小、分布深度更浅,更有可能对应于部分熔融的地壳;鄂尔多斯盆地在中下地壳以及上地幔内有着较大范围的高速异常一直延伸到120 km以下,而河套盆地地幔只在80 km以上部分有着高速异常的分布,此深度可能代表了河套盆地的岩石圈厚度,来自深部地幔的热物质上涌造成了该区域的岩石圈减薄;阿拉善块体在地壳和上地幔都表现出高低速共存的分布特征,暗示阿拉善块体西部岩石圈可能受青藏高原东北缘的挤压作用发生改造.

  相似文献   

19.

地震波衰减研究是了解岩石圈构造特征的有效方法.本文基于OBS2016-2测线的海底地震仪(OBS)数据, 首次对南海东北部洋陆过渡区域地震波衰减特征进行研究, 通过正演模拟获得该区域二维纵波衰减(QP)结构.结果表明, 下陆坡的上地壳存在一个宽约40 km, 厚度约为4~5 km的高衰减区, 其特征为低纵波和横波速度(VP为5.5~6.3 km·s-1VS为3.1~3.6 km·s-1)以及较低的波速比(VP/VS为1.72~1.80), 对应较低的QP(280~410), 推测与断裂发育有关, 且受到火山活动的影响.洋陆过渡及洋壳区域的上地壳高衰减区具有低QP(300~400)和高VP/VS(1.90~1.96)特征, 可能对应较多的断裂发育及流体运移.洋陆过渡区域下地壳高速异常体表现为相对低的QP(550~600), 对应较高的VP(7.0~7.8 km·s-1)和VS(3.5~3.8 km·s-1)以及较高的VP/VS(1.85~1.96), 推测与蛇纹石化作用有关.蛇纹石化可能进一步增加岩石的孔隙度并导致更多的流体运移, 使得洋陆过渡及洋壳区域存在较高的地震波衰减.QP结构有助于我们分析南海大陆边缘的地震波衰减特征, 结合VPVS以及VP/VS, 可以更好地了解该区域地质结构和岩石属性, 对进一步挖掘OBS数据信息有重要参考价值.

  相似文献   

20.

利用郯庐断裂带鲁苏皖段及邻区的261个宽频带地震台站5年记录的763个远震波形数据,计算并筛选得到了10846条远震P波接收函数.采用P波接收函数H-κ法得到了该区的地壳厚度和壳内平均泊松比值,并采用共转换点叠加法进一步揭示台站下方Moho界面的起伏形态.研究发现:⑴研究区Moho界面埋深在27~40 km范围内变化,平均深度在~34 km,总体上以郯庐断裂带为界呈现出东薄西厚的特征.地壳厚度在不同块体之间或者是块体内部存在着明显差异,表明不同的地质构造单元经历了不同的构造演化过程.⑵研究区地壳泊松比在0.15~0.32之间变化,平均泊松比为0.24,略低于全球陆壳和中国陆壳平均泊松比值;然而,较大的泊松比浮动范围却意味着研究区内地壳物质具有强烈的横向非均匀性及物质组成的复杂性.沿郯庐断裂带展布着一条NNE-SSW方向的泊松比高值异常带,推测是镁铁质基性岩沿郯庐断裂带上涌至地壳所致,亦或是高温高压的幔源热物质底侵至下地壳所致.⑶研究区的地壳厚度和壳内平均泊松比存在着反相关的关系;地壳厚度和地表地貌特征呈镜向关系,即造山隆起区Moho界面埋藏较深,而平原盆地区Moho界面埋藏较浅.Moho埋深等值线分布图和研究区布格重力异常特征对应良好.⑷共转换点(Common Conversion Point,CCP)叠加法对Moho界面的刻画与H-κ叠加法求得的地壳厚度结果具有较好的一致性.CCP剖面表明郯庐断裂带不仅是扬子断块区和华北块体的分界断裂,更是一条切割Moho面、深抵上地幔的深大超壳断裂带,错距达4~7 km.⑸研究区内部分台站下方存在壳内分界面,仍能通过改变壳内速度后采用H-κ法获得其埋藏深度.

  相似文献   

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

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