首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 734 毫秒
1.
  相似文献   

2.
Characteristics of the seismicity in depth ranges 0–33 and 34–70 km before ten large and great (M w = 7.0−9.0) earthquakes of 2000–2008 in the Sumatra region are studied, as are those in the seismic gap zones where no large earthquakes have occurred since at least 1935. Ring seismicity structures are revealed in both depth ranges. It is shown that the epicenters of the main seismic events lie, as a rule, close to regions of overlap or in close proximity to “shallow” and “deep” rings. Correlation dependences of ring sizes and threshold earthquakes magnitudes on energy of the main seismic event in the ring seismicity regions are obtained. Identification of ring structures in the seismic gap zones (in the regions of Central and South Sumatra) suggests active processes of large earthquake preparation proceed in the region. The probable magnitudes of imminent seismic events are estimated from the data on the seismicity ring sizes.  相似文献   

3.
The cause for prolific seismicity in the Koyna region is a geological enigma. Attempts have been made to link occurrence of these earthquakes with tectonic strain as well as the nearby reservoirs. With a view to providing reliable seismological database for studying the earth structure and the earthquake process in the Koyna region, a state of the art digital seismic network was deployed for twenty months during 1996–97. We present preliminary results from this experiment covering an area of 60 × 80 km2 with twenty seismic stations. Hypocentral locations of more than 400 earthquakes confined to 11×25 km2 reveal fragmentation in the seismicity pattern — a NE — SW segment has a dip towards NW at approximately 45°, whilst the other two segments show a near vertical trend. These seismic segments have a close linkage with the Western Ghat escarpment and the Warna fault. Ninety per cent of the seismicity is confined within the depth range of 3–10 km. The depth distribution of earthquakes delimits the seismogenic zone with its base at 10 km indicating a transition from an unstable to stable frictional sliding regime. The lack of shallow seismicity between 0 and 3 km indicates a mature fault system with well-developed gouge zones, which inhibit shallow earthquake nucleation. Local earthquake travel time inversion for P- and S-waves show ≈ 2% higher velocity in the seismogenic crust (0–10 km) beneath the epicentral tract relative to a lower velocity (2–3%) in the adjoining region. The high P- and S-wave velocity in the seismogenic crust argues against the presence of high pressure fluid zones and suggests its possible linkage with denser lithology. The zone of high velocity has been traced to deeper depths (≈ 70 km) through teleseismic tomography. The results reveal segmented and matured seismogenic fault systems in the Koyna region where seismicity is possibly controlled by strain build up due to competent lithology in the seismic zone with a deep crustal root.  相似文献   

4.
We analyse the seismicity pattern including b-value in the north Sumatra-Great Nicobar region from 1976 to 2004. The analysis suggests that there were a number of significant, intermediate and short-term precursors before the magnitude 7.6 earthquake of 2 November 2002. However, they were not found to be so prominent prior to the magnitude 9.0 earthquake of 26 December 2004 though downward migration of activity and a 50-day short-term quiescence was observed before the event. The various precursors identified include post-seismic and intermediate-term quiescence of 13 and 10 years respectively, between the 1976 (magnitude 6.3) and 2002 earthquakes with two years (1990–1991) of increase in background seismicity; renewed seismicity, downward migration of seismic activity and foreshocks in 2002, just before the mainshock. Spatial variation in b-value with time indicates precursory changes in the form of high b-value zone near the epicenter preceding the mainshocks of 2004 and 2002 and temporal rise in b-value in the epicentral area before the 2002 earthquake.  相似文献   

5.
基于集集强震群序列地震特征的地震追踪预测   总被引：1，自引：0，他引：1  
郑魁香  赵汝仁 《地学前缘》2002,9(2):493-498
分析集集强震群前余震序列的 7年 (1993/ 0 9/ 2 1— 2 0 0 0 / 0 9/ 2 0 )中震级规模在M =3 0以上的地震目录 ,可以找到前震类型、孕震空区特征、孕震条带特征、前震丛集性活动与信号震特征、主震前平静以及余震序列的二次余震等至少 6项清楚的地震序列特征。利用已发展出的年度强震趋势分析步骤的经验 ,佐以依据地震序列特征进一步加以追踪的观念 ,以集集地震序列分析为例 ,试图将地震趋势分析由年的时间尺度 ,追踪到更短的月的时间范围 ;并尝试建立台湾地区西部地震带浅源强震的追踪分析步骤 ,并为以测震学为基础的地震预测提供逼近短临时间尺度的分析方法。  相似文献   

6.

The seismicity of South Australia over the period 1980–92 is presented as a follow‐up to earlier studies. The South Australian seismic network has undergone a significant expansion in the last decade, and with it an increase in the number and precision of located earthquakes. The distribution of recent seismic activity is similar to the historical pattern of earthquakes and the previous instrumental seismicity maps, all of which show the three main areas as being the Flinders‐Mt Lofty Ranges, Eyre Peninsula, and the southeast. The one notable exception in the recent study is the presence of earthquake activity in the Musgrave Block, a previously aseismic region. Intensity characteristics are reported for earthquakes that were sufficiently widely felt. Fault plane solutions for three Flinders Ranges earthquakes (previously unpublished) are also presented; the focal mechanisms are consistent with predominant northeast‐southwest compression. The seismic moment method was used to estimate the seismic risk for the major population centres in terms of probability of exceedance of seismic intensity within a given period. These estimates are based on the recurrence parameters and intensity attenuation function for the region. The results place Adelaide close to the AS2121 ‐ 1979 Earthquake Code Zone I/Zone 2 boundary.  相似文献   

7.
We studied the temporal behavior of the background shallow seismicity rate in 700 circular areas across inland Japan. To search for and test the significance of the possible rate changes in background seismicity, we developed an efficient computational method that applies the space–time ETAS model proposed by Ogata in 1998 to the areas. Also, we conducted Monte Carlo tests using a simulated catalog to validate the model we applied. Our first finding was that the activation anomalies were found so frequently that the constant background seismicity hypothesis may not be appropriate and/or the triggered event model with constraints on the parameters may not adequately describe the observed seismicity. However, quiescence occasionally occurs merely by chance. Another outcome of our study was that we could automatically find several anomalous background seismicity rate changes associated with the occurrence of large earthquakes. Very significant seismic activation was found before the M6.1 Mt. Iwate earthquake of 1998. Also, possible seismic quiescence was found in an area 150 km southwest of the focal region of the M7.3 Western Tottori earthquake of 2000. The seismicity rate in the area recovered after the mainshock.  相似文献   

8.
The Burmese Arc seismic activity is not uniform for its ∼ 1100 km length; only the Northern Burmese Arc (NBA) is intensely active. Six large earthquakes in the magnitude range 6.1–7.4 have originated from the NBA Benioff zone between 1954–2011, within an area of 200 × 300 km2 where the Indian plate subducts eastward to depths beyond 200 km below the Burma plate. An analysis on seismogenesis of this interplate region suggests that while the subducting lithosphere is characterized by profuse seismicity, seismicity in the overriding plate is rather few. Large earthquakes occurring in the overriding plate are associated with the backarc Shan-Sagaing Fault (SSF) further east. The forecasting performance of the Benioff zone earthquakes in NBA as forerunner is analysed here by: (i) spatial earthquake clustering, (ii) seismic cycles and their temporal quiescence and (iii) the characteristic temporal b-value changes. Three such clusters (C1–C3) are identified from NBA Benioff Zones I & II that are capable of generating earthquakes in the magnitude ranges of 7.38 to 7.93. Seismic cycles evidenced for the Zone I displayed distinct quiescence (Q1, Q2 and Q3) prior to the 6th August 1988 (M 6.6) earthquake. Similar cycles were used to forecast an earthquake (Dasgupta et al. 2010) to come from the Zone I (cluster C1); which, actually struck on 4 February 2011 (M 6.3). The preparatory activity for an event has already been set in the Zone II and we speculate its occurrence as a large event (M > 6.0) possibly within the year 2012, somewhere close to cluster C3. Temporal analysis of b-value indicates a rise before an ensuing large earthquake.  相似文献   

9.
The earthquake of 6 October 1987 (M = 6.6), which occurred near the Shipunsky Cape, Kamchatka, was the largest crustal event in the vicinity of the main city of Kamchatka — Petropavlovsk-Kamchatsky — during the last three decades. It was followed by numerous aftershocks. This earthquake allowed us to test the effectiveness of the seismic hazard monitoring in Kamchatka, including the seismological, geodetic and hydrogeochemical surveys. The seismic survey provided the location and source nature of the main shock and aftershocks and the seismic environment of the main shock. The geodetic and hydrogeochemical surveys have yielded data on the response to earthquakes of the Earth's surface deformations, water level, and chemical elements concentration in the underground water. As a result, the following data were obtained:

u

  • The earthquake of 6 October had a seismic moment 4–10 E18 Nm, thrust type of faulting and the source volume of 20 × 20 × 10 km3. The maximum intensity was VI–VII (MSK-64 scale) and maximum acceleration 88 cm/s2.
  • Before this event, a relative increase in the number of the upper mantle (depth more than 100 km) moderate magnitude earthquakes during 5 years and a one-year period of seismic quiescence for small shallow earthquakes, were recognized. Significant anomalies in HCO3 and H3BO3 concentrations in the underground waters were observed in the wells a week before the main shock.
    •   相似文献   

    10.
    Seismic hazard analysis of the northwest Himalayan belt was carried out by using extreme value theory (EVT). The rate of seismicity (a value) and recurrence intervals with the given earthquake magnitude (b value) was calculated from the observed data using Gutenberg–Richter Law. The statistical evaluation of 12,125 events from 1902 to 2017 shows the increasing trend in their inter-arrival times. The frequency–magnitude relation exhibits a linear downslope trend with negative slope of 0.8277 and positive intercept of 4.6977. The empirical results showed that the annual risk probability of high magnitude earthquake M?≥?7.7 in 50 years is 88% with recurrence period of 47 years, probability of M?≤?7.5 in 50 years is 97% with recurrence period of 27 years, and probability of M?≤?6.5 in 50 years is 100% with recurrence period of 4 years. Kashmir valley, located in the NW Himalaya, encompasses a peculiar tectonic and structural setup. The patterns of the present and historical seismicity records of the valley suggest a long-term strain accumulation along NNW and SSE extensions with the decline in the seismic gap, posing a potential threat of earthquakes in the future. The Kashmir valley is characterized by the typical lithological, tectono-geomorphic, geotechnical, hydrogeological and socioeconomic settings that augment the earthquake vulnerability associated with the seismicity of the region. The cumulative impact of the various influencing parameters therefore exacerbates the seismic hazard risk of the valley to future earthquake events.  相似文献   

    11.
    Seismotectonics of the Nepal Himalaya from a local seismic network   总被引：3，自引：0，他引：3  
    The National Seismological Network of Nepal consists of 17 short period seismic stations operated since 1994. It provides an exceptional view of the microseismic activity over nearly one third of the Himalayan arc, including the only segment, between longitudes 78°E and 85°E, that has not produced any M>8 earthquakes over the last century. It shows a belt of seismicity that follows approximately the front of the Higher Himalaya with most of the seismic moment being released at depths between 10 and 20 km. This belt of seismicity is interpreted to reflect interseismic stress accumulation in the upper crust associated with creep in the lower crust beneath the Higher Himalaya. The seismic activity is more intense around 82°E in Far-Western Nepal and around 87°E in Eastern Nepal. Western Nepal, between 82.5 and 85°E, is characterized by a particularly low level of seismic activity. We propose that these lateral variations are related to segmentation of the Main Himalayan Thrust Fault. The major junctions between the different segments would thus lie at about 87°E and 82°E with possibly an intermediate one at about 85°E. These junctions seem to coincide with some of the active normal faults in Southern Tibet. Lateral variation of seismic activity is also found to correlate with lateral variations of geological structures suggesting that segmentation is a long-lived feature. We infer four 250–400 km long segments that could produce earthquakes comparable to the M=8.4 Bihar–Nepal earthquake that struck eastern Nepal in 1934. Assuming the model of the characteristic earthquake, the recurrence interval between two such earthquakes on a given segment is between 130 and 260 years.  相似文献   

    12.
    In estimating the likelihood of an earthquake hazard for a seismically active region, information on the geometry of the potential source is important in quantifying the seismic hazard. The damage from an earthquake varies spatially and is governed by the fault geometry and lithology. As earthquake damage is amplified by guided seismic waves along fault zones, it is important to delineate the disposition of the fault zones by precisely determined hypocentral parameters. We used the double difference (DD) algorithm to relocate earthquakes in the Koyna-Warna seismic zone (KWSZ) region, with the P- and S-wave catalog data from relative arrival time pairs constituting the input. A significant improvement in the hypocentral estimates was achieved, with the epicentral errors <30 m and focal depth errors <75 m i.e. errors have been significantly reduced by an order of magnitude from the parameters determined by HYPO71. The earthquake activity defines three different fault segments. The seismogenic volume is shallower in the south by 3 km, with seismicity in the north extending to a depth of 11 km while in the south the deepest seismicity observed is at a depth of 8 km. By resolving the structure of seismicity in greater detail, we address the salient issues related to the seismotectonics of this region.  相似文献   

    13.
    地震迁移的类型、特征及机制讨论   总被引：2，自引：2，他引：0  
    地震活动的迁移是指地震沿着某一方向循序地发生,是地震活动总体无序中局部出现的有序结构。通过系统分析区域上典型的地震迁移现象可以发现,迁移可划分为沿断裂走向的纵向迁移、垂直断裂走向的横向迁移与岩石圈尺度的深源迁移三大类。结合具体的实例分析,可初步归纳出纵向迁移（包括单向、双向、反复和跳跃式迁移等常见形式）、横向迁移和深源地震迁移的主要特征,并初步估算出了不同类型迁移的速度值,其中沿全球板块边界纵向迁移平均速度约为V=569 km/a,沿亚板块边界的平均速度约为V=120 km/a,沿大陆内部断裂带平均速度约为V=50 km/a。横向迁移相对比较复杂,初步发现在东亚存在两种速度分别为约20 km/a、50 km/a的迁移现象。而深源地震迁移速度的全球平均值大约为360 km/a。地震的纵向与横向迁移都存在不同层次和级别,也存在多种不同频率、能量与速度的形变波与迁移现象,这很可能是区域上大地震丛集活动过程中断层相互作用、地震应力触发和岩石圈尺度的形变波传播等因素共同作用的结果,而这种大空间尺度上的地震迁移现象的存在及其所具有的规律性特征,显然可为开展区域地震危险性分析提供新的思路和方法参考。   相似文献   

    14.
    Seismic potential of Southern Italy   总被引：1，自引：2，他引：1  
    To improve estimates of the long-term average seismic potential of the slowly straining South Central Mediterranean plate boundary zone, we integrate constraints on tectonic style and deformation rates from geodetic and geologic data with the traditional constraints from seismicity catalogs. We express seismic potential (long-term average earthquake recurrence rates as a function of magnitude) in the form of truncated Gutenberg–Richter distributions for seven seismotectonic source zones. Seismic coupling seems to be large or even complete in most zones. An exception is the southern Tyrrhenian thrust zone, where most of the African–European convergence is accommodated. Here aseismic deformation is estimated to range from at least 25% along the western part to almost 100% aseismic slip around the Aeolian Islands. Even so, seismic potential of this zone has previously been significantly underestimated, due to the low levels of recorded past seismicity. By contrast, the series of 19 M6–7 earthquakes that hit Calabria in the 18th and 19th century released tectonic strain rates accumulated over time spans up to several times the catalog duration, and seismic potential is revised downward. The southern Tyrrhenian thrust zone and the extensional Calabrian faults, as well as the northeastern Sicilian transtensional zone between them (which includes the Messina Straits, where a destructive M7 event occurred in 1908), all have a similar seismic potential with minimum recurrence times of M ≥ 6.5 of 150–220 years. This potential is lower than that of the Southern Apennines (M ≥ 6.5 recurring every 60 to 140 years), but higher than that of southeastern Sicily (minimum M ≥ 6.5 recurrence times of 400 years). The high seismicity levels recorded in southeastern Sicily indicate some clustering and are most compatible with a tectonic scenario where the Ionian deforms internally, and motions at the Calabrian Trench are small. The estimated seismic potential for the Calabrian Trench and Central and Western Sicily are the lowest (minimum M ≥ 6.5 recurrence times of 550–800 years). Most zones are probably capable of generating earthquakes up to magnitudes 7–7.5, with the exception of Central and Western Sicily where maximum events sizes most likely do not exceed 7.  相似文献   

    15.
    Ten moderate to large (magnitude 6–7) earthquakes have occurred in southwestern British Columbia and northwestern Washington in the last 130 years. A future large earthquake close to Vancouver, Victoria, or Seattle would cause tens of billions of dollars damage and would seriously impact the economies of Canada and the United States. An improved understanding of seismic hazards and risk in the region has been gained in recent years by using geologic data to extend the short period of instrumented seismicity. Geologic studies have demonstrated that historically unprecedented, magnitude 8 to 9 earthquakes have struck the coastal Pacific Northwest on average once every 500 years over the last several thousand years; another earthquake of this size can be expected in the future. Geologic data also provide insights into the likely damaging effects of future large earthquakes in the region. Much of the earthquake damage will result directly from ground shaking, but damage can also be expected from secondary phenomena, including liquefaction, landslides, and tsunamis. Vancouver is at great risk from earthquakes because important infrastructure, including energy and transportation lifelines, probably would be damaged or destroyed by landslides and liquefaction-induced ground failure.  相似文献   

    16.
    Volcanic earthquakes on Kamchatka can be divided into two large groups: earthquakes with depths of 0–40 km generated by stresses which arise during magma migration in the Earth's crust under volcanos (the first group), and the earthquakes directly connected with the eruptions (volcanic tremor, explosive earthquakes, etc.—the second group). This paper presents a review of some energetic, spectral and spatio-temporal characteristics of the Kamchatkan volcanic earthquakes of the first group and their relationship with volcanic phenomena.

    Seismicity related to volcanic activity has the following specific features: a local and predominantly swarm-like pattern of earthquake origination; iteration of earthquake swarms in the same seismically active zones; many shallow and relatively small events; a small magnitude limit (up to 5.5–6); the existence of longer-period variations of volcanic earthquake foci as compared to the tectonic one; and a comparatively high value of the slope of the earthquake recurrence plot. At the same time, similarity in behaviour of some parameters of the seismic regime during the preparation and development of eruptions and prior to large earthquakes, as well as the destruction of samples, are noted.  相似文献   


    17.
    东南亚地区是“21世纪海上丝绸之路”（以下简称“海洋丝路”）的重要组成部分,该区历史上曾发生十余次巨大地震,地震及其次生地质灾害是威胁东南亚地区经济社会发展和国际合作的主要自然灾害。系统梳理该区地震活动的时空分布特征及评估未来灾害风险格局,对于推进“一带一路”倡议实施及区域经济社会可持续发展具有重要意义。文章基于东南亚地区1900年以来M≥5地震的时空分布统计分析和地震b值计算,揭示出该区的地震活动在时间上表现出活跃期与平静期交替变化的特征;空间上表现出明显的聚集效应,成丛性强且主要集中在5个地震统计区内,其中印尼—马来多岛弧盆系地震区和菲律宾群岛地震区的地震活动最为活跃。总体而言,东南亚5个地震区的b值偏低,在0.42~0.91之间。该区内的地震b值也存在时空差异,受大地震事件、俯冲带年龄、活动断裂带和震源深度等众多因素影响,但主控因素在不同区域有所不同。地震b值时空变化特征对区域地震活动预测具有启示作用。上述认识为推进“海洋丝路”工程建设和“一带一路”防灾减灾对策提供了科学支撑。  相似文献   

    18.
    The Japan Trench subduction zone, located east of NE Japan, has regional variation in seismicity. Many large earthquakes occurred in the northern part of Japan Trench, but few in the southern part. Off Miyagi region is in the middle of the Japan Trench, where the large earthquakes (M > 7) with thrust mechanisms have occurred at an interval of about 40 years in two parts: inner trench slope and near land. A seismic experiment using 36 ocean bottom seismographs (OBS) and a 12,000 cu. in. airgun array was conducted to determine a detailed, 2D velocity structure in the forearc region off Miyagi. The depth to the Moho is 21 km, at 115 km from the trench axis, and becomes progressively deeper landward. The P-wave velocity of the mantle wedge is 7.9–8.1 km/s, which is typical velocity for uppermost mantle without large serpentinization. The dip angle of oceanic crust is increased from 5–6° near the trench axis to 23° 150 km landward from the trench axis. The P-wave velocity of the oceanic uppermost mantle is as small as 7.7 km/s. This low-velocity oceanic mantle seems to be caused by not a lateral anisotropy but some subduction process. By comparison with the seismicity off Miyagi, the subduction zone can be divided into four parts: 1) Seaward of the trench axis, the seismicity is low and normal fault-type earthquakes occur associated with the destruction of oceanic lithosphere. 2) Beneath the deformed zone landward of the trench axis, the plate boundary is characterized as a stable sliding fault plain. In case of earthquakes, this zone may be tsunamigenic. 3) Below forearc crust where P-wave velocity is almost 6 km/s and larger: this zone is the seismogenic zone below inner trench slope, which is a plate boundary between the forearc and oceanic crusts. 4) Below mantle wedge: the rupture zones of thrust large earthquakes near land (e.g. 1978 off Miyagi earthquake) are located beneath the mantle wedge. The depth of the rupture zones is 30–50 km below sea level. From the comparison, the rupture zones of large earthquakes off Miyagi are limited in two parts: plate boundary between the forearc and oceanic crusts and below mantle wedge. This limitation is a rare case for subduction zone. Although the seismogenic process beneath the mantle wedge is not fully clarified, our observation suggests the two possibilities: earthquake generation at the plate boundary overridden by the mantle wedge without serpentinization or that in the subducting slab.  相似文献   

    19.
    The Himalayas has experienced varying rates of earthquake occurrence in the past in its seismo-tectonically distinguished segments which may be attributed to different physical processes of accumulation of stress and its release, and due diligence is required for its inclusion for working out the seismic hazard. The present paper intends to revisit the various earthquake occurrence models applied to Himalayas and examines it in the light of recent damaging earthquakes in Himalayan belt. Due to discordant seismicity of Himalayas, three types of regions have been considered to estimate larger return period events. The regions selected are (1) the North-West Himalayan Fold and Thrust Belt which is seismically very active, (2) the Garhwal Himalaya which has never experienced large earthquake although sufficient stress exists and (3) the Nepal region which is very seismically active region due to unlocked rupture and frequently experienced large earthquake events. The seismicity parameters have been revisited using two earthquake recurrence models namely constant seismicity and constant moment release. For constant moment release model, the strain rates have been derived from global strain rate model and are converted into seismic moment of earthquake events considering the geometry of the finite source and the rates being consumed fully by the contemporary seismicity. Probability of earthquake occurrence with time has been estimated for each region using both models and compared assuming Poissonian distribution. The results show that seismicity for North-West region is observed to be relatively less when estimated using constant seismicity model which implies that either the occupied accumulated stress is not being unconfined in the form of earthquakes or the compiled earthquake catalogue is insufficient. Similar trend has been observed for seismic gap area but with lesser difference reported from both methods. However, for the Nepal region, the estimated seismicity by the two methods has been found to be relatively less when estimated using constant moment release model which implies that in the Nepal region, accumulated strain is releasing in the form of large earthquake occurrence event. The partial release in second event of May 2015 of similar size shows that the physical process is trying to release the energy with large earthquake event. If it would have been in other regions like that of seismic gap region, the fault may not have released the energy and may be inviting even bigger event in future. It is, therefore, necessary to look into the seismicity from strain rates also for its due interpretation in terms of predicting the seismic hazard in various segments of Himalayas.  相似文献   

    20.
    针对2015年4月25日尼泊尔M8.1地震后喜马拉雅造山带的未来地震危险性问题,通过对喜马拉雅带历史大地震应变能释放和在尼泊尔地震发震前后的区域地震活动图像进行了分析研究。结果发现喜马拉雅带很可能已进入新-轮的地震活跃期。此次尼泊尔大地震不足以将喜马拉雅带中段的地壳应变能全部释放,喜马拉雅带中段的地震活动和藏南裂谷带地震活动具有密切的关联,在喜马拉雅带中段和藏南裂谷带还将有大地震活动。同时研究结果还显示现今在喜马拉雅带的东段存在阿萨姆围空区和不丹围空区,在喜马拉雅的西段出现噶尔围空区,喜马拉雅西段新德里和西藏接壤地区以及喀喇昆仑断裂上噶尔县地区地震危险性很高,喜马拉雅东段林芝山南地区以南的阿萨姆和不丹地区危险性很高,应引起重视。  相似文献   

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

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