首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 109 毫秒
1.
在巴楚地区露头高频层序研究及塔中地区钻测井三级层序划分的基础上,对塔里木盆地塔中地区塔中162井和塔中43井下奥陶统鹰山组层序OSQ2的自然伽马能谱测井ln(Th/K)曲线进行了滤波处理、快速傅里叶变换处理、小波时频分析以及调谐处理,识别出具有米兰科维奇旋回特征的高频旋回。结果表明:斜度(黄赤交角)旋回最为明显,周期为37.0 ka,是以向上变浅为主的六级米级旋回的主要控制因素,其主旋回平均厚度在台地边缘相带的塔中162井及塔中43井中分别为4.55 m及3.97 m,在巴楚地区半局限-开阔台地相带中为2~4 m;95 ka的短偏心率周期是形成五级准层序的主要控制因素,而大致代表四级准层序组形成时限的413 ka长偏心率周期在地层记录中表现并不明显。据高频层序叠置关系分析及平均主旋回个数初步估算,塔中-巴楚地区下奥陶统鹰山组层序OSQ2的形成时限大约为4.92 Ma,塔中43井鹰山组层序OSQ2受后期构造隆升而被剥蚀约192  m。结合四级层序(准层序组)发育时限(413 ka)及其叠置关系,以及与其调谐的小波变换曲线特征,初步将塔中-巴楚地区鹰山组层序OSQ2划分出12个准层序组。自然伽马能谱测井ln(Th/K)曲线包含沉积古水深相对变化的信息,是反映气候变化生成的米兰科维奇旋回层序的良好指标。  相似文献   

2.
四川盆地中二叠统茅口组形成于中二叠世中晚期。通过分析四川盆地中二叠统茅口组露头剖面岩性特征、对比广元上寺长江沟露头剖面碳氧同位素变化趋势与沉积层序发育特征,将茅口组划分为2个可全盆地追踪的Ⅱ型三级层序,并分析了茅口组三级层序的主控因素。在此基础上,对钻穿中二叠统茅口组典型井的自然伽马能谱测井曲线开展频谱分析,识别出多种具有米兰科维奇旋回特征的高频旋回,并计算出茅口组沉积的平均速率及沉积时限等相关参数。最后,采用数字滤波消除掉其他次要旋回因素的影响而仅保留与主控因素相关的旋回信息,建立茅口组的高频层序划分方案。结果表明,茅口组三级层序在形成过程中主要受控于构造升降及全球海平面变化;茅口组米兰科维奇旋回特征明显,其中长偏心率旋回(413.0 ka)和短偏心率旋回(123.0 ka)分别是形成四级层序(准层序组)和五级层序(准层序)的主控因素,与之对应的平均旋回厚度在龙17井区分别为13.44m和4.31m,在安平1井区分别为16.03m和4.68m;茅口组大约发育15个四级层序,其发育时限大约为6.11Ma;根据构造背景曲线和长偏心率旋回曲线的叠加曲线划分高频层序,其高频层序界面更加接近实际地层的发育情况。  相似文献   

3.
旋回地层学的研究对象由海相地层逐渐转移到陆相深水地层.为了探究通过旋回地层学研究方法在陆相浅水湖盆建立高频层序格架的可能,采用频谱分析、天文调谐、“米氏”沉积速率拾取等手段对东海陆架盆地西湖凹陷黄岩区渐新统坳陷湖盆的岩相敏感曲线-GR数据开展研究.频谱分析结果显示黄岩区花港组受405 ka长偏心率周期(E)、121 ka和97 ka短偏心率周期(e)、28 ka斜率周期(O)以及约22 ka岁差周期(P)的调控.在浅水湖泊环境的富泥区,以La2004天文解决方案给出的65°N平均日照序列为参照进行天文调谐,建立了浮动天文年代标尺;在浅水三角洲环境的富砂区,通过滑动窗口频谱分析得到“米氏”沉积速率曲线,结合锆石铀铅测年资料的校验和岩石组合类型联合约束,完成高频层序格架划分.研究表明:黄岩区花港组持续时间约为10.9 Ma,可识别出27个405 ka长偏心率周期,根据偏心率滤波曲线和日照量各级次包络面的对应关系共划分11个四级层序,27个五级层序,分别对应着0.8~1.6 Ma和0.4 Ma的基准面旋回,每个旋回在沉积速率曲线上表现为高?低?高的特征.本次研究提供了一种利用“米氏”沉积速率进行高频层序格架搭建及层序?时间单元拾取的新方法,拓宽了旋回地层学的应用范围.   相似文献   

4.
Sedimentary cyclic sequences deposited during the Late Palaeozoic Ice Age are widespread. Glacio‐eustatic control of the cyclic patterns is commonly accepted, and the durations of the cyclothems generally match the short‐ and long‐eccentricity Milankovitch orbital parameters. Nevertheless, geochemical fingerprints of orbital parameters are poorly documented in deep‐time sedimentary records. Here, we report on well‐exposed Bashkirian cyclothems of c. 123 ka and c. 400 ka durations from the Valdorria platform. The shorter‐term cyclothems can be grouped into longer‐term composite sequences that are consistent with generally accepted durations of c. 125 ka and c. 400 ka for Milankovitch eccentricity cycles. The stratigraphic pattern is mirrored by the isotope geochemical signals, which show distinct recurring trends. These trends are confirmed by statistical tests. Whereas intrinsic factors and/or subaerial exposure related to sea‐level lowstands might have truncated cycle patterns in tectonically stable basins, rapid subsidence of the Valdorria platform's foreland basin appears to have contributed to a faithful recording of cyclothems of different orders. The patterns and biostratigraphic constraint revealed in this study demonstrate the power of orbital forcing in imprinting sedimentary and geochemical signals in the rock record.  相似文献   

5.
显生宙沉积旋回会受到地球轨道参数偏心率、地轴斜率和岁差的周期性变化(米兰科维奇旋回)的控制,但目前对前寒武系的相关研究较为薄弱。天津蓟县剖面中元古界洪水庄组-铁岭组为一套潟湖-潮坪相沉积,主要呈砂岩-页岩互层的韵律性产出。为探索这种岩性韵律是否与轨道旋回有关,对其进行了高分辨率的岩性刻画,并以磁化率和伽马能谱数据作为古气候-古环境替代性指标,通过频谱分析等方法进行旋回地层学分析。结果表明各指标均记录了完整的米兰科维奇旋回。铁岭组一段下部由短偏心率、斜率和岁差控制的地层旋回厚度分别为1.2~1.5 m、0.4 m和0.17~0.19 m,其中0.1 m的半岁差周期也较明显;洪水庄组二段顶部由长偏心率、短偏心率、斜率和岁差控制的地层旋回厚度分别为1.1~1.8 m、0.34~0.54 m、0.14~0.16 m和0.08~0.09 m。以识别出的稳定长偏心率周期405 ka旋回对洪水庄组进行了天文年代校准,并估计出当时的短偏心率、斜率以及岁差周期分别为100 ka、22~24 ka和15 ka。在洪水庄组中记录的偏心率振幅变化周期为~2 Ma,比现今~2.4 Ma的周期略短。  相似文献   

6.
北黄海东部坳陷始新统米兰科维奇旋回特征   总被引:1,自引:0,他引:1       下载免费PDF全文
吴淑玉  刘俊 《地球科学》2015,40(11):1933-1944
米兰科维奇旋回是记录在沉积地层中的表现形式, 其代表的时间涵义是进行高分辨率地层划分和对比的有效手段, 从地层中揭示的米兰科维奇旋回, 可以完善地层学尤其是旋回地层学理论.以北黄海东部坳陷为研究对象, 根据J. Laskar的解决方案计算出该区始新统的米兰科维奇旋回周期为: 125 ka和99 ka的偏心率周期, 51 ka和39 ka的地轴斜率周期, 23 ka和19 ka的岁差周期.通过对3口井的GR(自然伽马)和SP(自然电位)测井曲线进行频谱分析, 发现其频谱峰值与天文周期存在着良好的对应关系, 因此可以判定该区域地层中保存着完好的米兰科维奇沉积旋回.地层旋回厚度存在13.03~15.89 m的长周期、3.70~5.21 m的中周期和2.17~2.94 m的短周期, 并由此计算地层的沉积速率为121.20~127.12 m/Ma.从隆起沿着斜坡往湖盆中心, 沉积持续时间越长且沉积厚度也越厚, 但沉积速率相对稳定.通过连续小波变换对始新统地层进行小层划分, 划分出6期沉积体, 以每个沉积体为独立窗口进行频谱分析计算出沉积的持续时间和速率, 从气候变化的影响分析每个阶段的沉积环境.从下往上层序地层E6期为低水位体系域, E5和E4为湖泊扩张体系域, E3为高位体系域, E2和E1为湖泊收缩体系域.以上方法可以证明米兰科维奇进行沉积旋回分析是一种有效的方法.   相似文献   

7.
上二叠统记录了地质历史时期最大规模的生物灭绝事件和最深刻的环境变化。对上二叠统的层序地层格架进行精细 描绘,建立高分辨率的地层序列,是深入了解此次事件及其演化的基础和关键。基于钻测井、岩心观测及地球化学分析测 试结果等资料综合分析,并运用小波分析技术,对黔西盘县上二叠统煤系进行了米氏旋回的识别和划分,结果表明,研究 区上二叠统煤系记录了稳定的米兰科维奇旋回,天文轨道周期对其沉积过程具有明显影响,由长、短偏心率、地轴斜率和 岁差周期引起的地层旋回厚度分别为16.06~17.24 m、5.39~5.70 m、2.11~2.15 m、1.12~1.21 m,长偏心率周期对地层中沉积 旋回的控制和影响最强。对长、短偏心率周期进行滤波分析后,建立了上二叠统煤系“浮动”天文年代标尺,为约束同沉 积火山事件层(Tonstein) 的形成及其持续时限提供了年代学依据。以区域等时对比效果明显的长偏心率旋回为标尺,并结 合旋回沉积序列、旋回界面特征等,将上二叠统煤系划分为4个三级层序,并进一步划分为16个四级层序(对应于中期旋 回),建立了研究区高频层序地层格架。  相似文献   

8.
上二叠统记录了地质历史时期最大规模的生物灭绝事件和最深刻的环境变化。对上二叠统的层序地层格架进行精细 描绘,建立高分辨率的地层序列,是深入了解此次事件及其演化的基础和关键。基于钻测井、岩心观测及地球化学分析测 试结果等资料综合分析,并运用小波分析技术,对黔西盘县上二叠统煤系进行了米氏旋回的识别和划分,结果表明,研究 区上二叠统煤系记录了稳定的米兰科维奇旋回,天文轨道周期对其沉积过程具有明显影响,由长、短偏心率、地轴斜率和 岁差周期引起的地层旋回厚度分别为16.06~17.24 m、5.39~5.70 m、2.11~2.15 m、1.12~1.21 m,长偏心率周期对地层中沉积 旋回的控制和影响最强。对长、短偏心率周期进行滤波分析后,建立了上二叠统煤系“浮动”天文年代标尺,为约束同沉 积火山事件层(Tonstein) 的形成及其持续时限提供了年代学依据。以区域等时对比效果明显的长偏心率旋回为标尺,并结 合旋回沉积序列、旋回界面特征等,将上二叠统煤系划分为4个三级层序,并进一步划分为16个四级层序(对应于中期旋 回),建立了研究区高频层序地层格架。  相似文献   

9.
将天文轨道周期与不同级别的旋回联系起来,旨在使米兰科维奇周期这一高精度时间标尺纳入高频层序地层划分中,实现高精度旋回的识别与划分。以哈萨克斯坦南图尔盖盆地Aryskum地堑Ary301井为例,基于不同沉积特征,分别对卡拉甘塞组I~IV段自然伽马测井数据进行频谱分析和连续小波变换,结果显示沉积地层中保存完好的米兰科维奇旋回,天文轨道周期对Aryskum地堑沉积过程具有明显影响,并将31.9~39.5 m旋回厚度解释为受400 kyr长偏心率周期控制,11.9~14.2 m,6.7~8.8 m旋回厚度分别受125 kyr和95 kyr短偏心率周期控制。对长、短偏心率周期进行滤波分析后,与天文模型理论周期曲线进行对比,建立卡拉甘塞组的浮动天文年代标尺,分别以400 kyr、125 kyr偏心率周期滤波曲线作为中期和短期旋回划分的参考曲线,共识别出11.5个中期旋回和47个短期旋回,为高频旋回划分提供了一种不受人为因素影响的天然标准,保证了研究区旋回划分的科学性和统一性。  相似文献   

10.
曹海洋  王华  赵睿 《地球科学》2017,42(12):2299-2311
在旋回地层学研究中,高分辨率地球化学数据作为一重要的替代指标,其一般获取方式(ICP-MS或WD-XRF)较耗费时间和经费,而手持X射线衍射仪(ED-XRF)能够高效的解决这一问题,目前将ED-XRF运用在沉积地层研究中,尤其是旋回地层学研究领域中较少.以酒泉盆地青西凹陷早白垩世下沟组深湖相泥岩、白云质泥岩为研究对象,将连续采集的岩心样品通过ED-XRF测试,该测试方法能够保证在不破坏样品的前提下,并在误差允许的范围内的条件下,短时间内(5 min)完成对单一样品的测试.通过测试精确得出24种主量及微量元素的含量,将得出的地球化学特征作为旋回地层学研究的替代指标,并针对K/Al,Ca/Al,Si/Al,Fe/Al,Ti/Al,Mn/Al,Ni/Al,Zn/Al和Pb/Al的比值进行频谱分析.数据表明W2井下沟组一段(SQK1g1)地层旋回性表现出与天文旋回明显的对应关系.沉积地层中的旋回厚度比值为20.25:4.75:1.75:1.00,与米兰科维奇旋回中长偏心率、短偏心率、斜率和岁差的周期比值:400 ka:95 ka:37 ka:20 ka相对应.该地区的沉积记录受天文轨道周期的影响得到了证明,与此同时,为手持X射线衍射仪在旋回地层学中的应用开辟了新领域.   相似文献   

11.
选择黔西地区盘1井、威宁1井、纳雍1井和织金1井,对长兴组自然伽马测井数据进行频谱分析,识别出其中的米兰科维奇周期,并采用滤波手段提取出其中的稳定周期,进行连井对比。结果表明频谱峰值与米兰科维奇周期存在良好的对应关系,研究区长兴组地层中保存着123. 0 ka的偏心率周期、45. 0 ka和35. 6 ka的斜率周期、21. 2 ka和17. 7 ka的岁差周期,各级周期的旋回厚度分别为4. 59~5. 73 m、1. 37~2. 15 m和0. 66~0. 98 m;各井的短偏心率旋回曲线显示研究区长兴组地层记录了约18个短偏心率旋回,据此估计出长兴组沉积时限约为2. 21 Ma;通过短偏心率旋回层和所计算沉积速率的连井对比,发现地层中海相标志岩层段等时性对比效果明显,反映了区域快速海侵的过程,而沉积速率的变化可能受到了海平面变化和物源供应的共同影响。  相似文献   

12.
The lowermost Cretaceous (Berriasian) Sierra del Pozo Formation is divisible into metre‐scale cycles that are bundled in a four‐tiered hierarchy of cycles and cyclic sequences. At the smallest scale, sixth‐order cycles, thought to be the product of precessionally forced sea‐level fluctuations (c. 20 ka), average less than a metre in thickness, shallow upward and are bounded by surfaces where deeper facies abruptly overlie shallower facies. Bundles of sixth‐order cycles, called fifth‐order sequences, are recognized by two types of asymmetric patterns in facies distribution. First, more pronounced facies changes occur at sixth‐order cycle boundaries lower in a sequence, whereas smaller facies changes occur at cycle boundaries higher in a sequence. Secondly, subtidal marine limestones (i.e. biomicrites and biosparites) are the dominant facies lower in a cyclic sequence, whereas more restricted or shallower water facies (i.e. mud‐cracked microbial laminites, dolomites, shale and clay) are more predominant higher in a cyclic sequence. The bundling of sixth‐order cycles is explained as a product of periodic change (100 ka) in the eccentricity of the earth's orbit. The degree of orbital eccentricity modulates the magnitude of precessional sea‐level fluctuations, which in turn determine the relative facies contrasts at sixth‐order cycle boundaries. Larger scale fourth‐ and third‐order sequences are defined by similar patterns in facies contrast at successive sequence boundaries and by a change in predominant facies type from bottom to top. These patterns are explained as the product of 400 ka and 2 Ma periodic variations in the eccentricity of the earth's orbit. In summary, the strength of the precessional signal varies in consort with periodic changes in orbital eccentricity at three time scales producing a stacked hierarchy of cycles and sequences of cycles. Change in obliquity of the earth's axis may modulate the effects of the precessional signal and thereby modify patterns of preserved cyclic structure. This interpretation of cyclicity in the Sierra del Pozo section is markedly different from that of Jiménez de Cisneros & Vera (1993), who attributed all rock cycles to a single process (obliquity) even though the cycles they described varied in thickness from less than 1 m to more than 4 m. The cycles they described are shown here to be either precessional cycles or composite sets of precessional cycles (100 ka sequences). In this new hierarchical interpretation, ninety‐six 100 ka and twenty‐four 400 ka sequences are recognized, extending the time of deposition of the Sierra del Pozo section to more than 9·6 Ma.  相似文献   

13.
地球围绕太阳的轨道要素周期性变化引起气候旋回变化,进而影响到地质历史中沉积物及生物面貌的旋回变化。牛38井始新统沙河街组三段旋回地层学研究表明,沉积物和古生物对404.8ka的天文偏心率旋回具有明显的响应特征。在一个偏心率旋回之内,随着地球围绕太阳公转轨道位置的变化,气候由热变冷再变热,古生物属种数量、沉积物粒度出现相应的的旋回变化。若以404.8ka偏心率旋回为单位来划分高频层序,可以保证所划分的地层具有异旋回性和严格的等时性。  相似文献   

14.
ABSTRACT Sedimentary cycles recorded in young sediments are often attributed to fluctuations of the Earth's climate on a 104−106-year scale which in turn is governed by periodic variations in solar insolation linked to orbital (Milankovitch) parameters. A spectacular example of cyclic stratal patterns in ancient deposits is the Middle Triassic Latemar carbonate platform (W Dolomites, N Italy). Based on spectral analyses from previous studies, a superimposition of precession (∼20 ka) and eccentricity (∼100 ka) controlled sea-level fluctuations has been suggested to account for the stacking hierarchy at Latemar, with ∼20 ka being assigned to each highest-order depositional cycle. Zircon U–Pb isotopic ages from volcanic-ash layers within the cyclic succession, corroborated by biostratigraphic constraints, suggest that the average time interval for every individual cycle is significantly smaller than the shortest Milankovitch period and therefore challenge previously published interpretations relating distinct spectral peaks to the above mentioned hierarchy. However, our new spectral data indicate that cyclicities resembling Milankovitch characteristics might exist, but on an entirely different scale. Our findings show that frequency spectra should only be interpreted in combination with robust age control. They also encourage the search for complementary mechanisms controlling carbonate deposition.  相似文献   

15.
针对应用常规方法进行高频层序划分和对比容易受人为因素影响,所建立的高级别层序地层格架具有多解性等问题,本文引入天文地层学中气候旋回受天文周期驱动的理论,选择涠西南凹陷WZ11-4N油田4口钻井中的流一段地层为研究对象,进行米兰科维奇旋回特征分析。对自然伽马数据进行频谱分析的结果显示,在不同钻井中,流一段地层主要受偏心率周期(401ka、125 ka和95 ka)控制。对该周期滤波分析后,选择控制流一段发育的主要米兰科维奇周期曲线,建立了该区的高分辨率天文年代标尺。最后以偏心率125 ka周期曲线作为6级层序划分的参考曲线,对流一段层序进行高频旋回地层划分与对比。在此基础上,最终实现高频地层格架下沉积相的精细对比。  相似文献   

16.
松辽盆地东岭地区泉头组高频沉积旋回的叠加型式分析   总被引:1,自引:1,他引:0  
根据自然伽马测井曲线计算的旋回厚度,采用Fischer图解识别高频旋回的叠加型式,试图揭示湖相沉积中GR测井曲线反映的泥质含量波动与气候旋回的关系。选择松辽盆地东岭地区顶底清楚,地层连续的SN109井泉头组测井曲线中作为实验数据,根据极值点法共识别出1166个高频沉积旋回,旋回厚度变化于0.25~2.75 m之间,旋回平均厚度为0.62 m;根据Fischer图解,估算高频旋回的平均周期为17.2 ka,中频旋回周期在86.1~103.3 ka之间,低频旋回周期在305~416 ka之间,分别与米兰科维奇旋回岁差周期19~24 ka,轴斜率周期85~140 ka,偏心率周期350~400 ka存在着对应关系,说明SN109井泉头组湖相沉积旋回的形成受米兰科维奇气候旋回机制的控制。  相似文献   

17.
赵军  赵凯  张金宇 《沉积学报》2022,40(3):801-812
基于米兰科维奇天文旋回理论在三角洲前缘开展砂岩组和单层级别的油层对比研究,频谱分析和小波变换揭示萨尔图油层Ⅲ油层组存在米兰科维奇天文旋回,可识别出比值接近5∶2∶1的3个旋回厚度,分别对应于偏心率100 ka、斜率40 ka和岁差20 ka周期,推测地层在沉积时受到米兰科维奇天文轨道周期性变化的影响。根据岁差20 ka周期带通滤波(BPF)和堆积速率变化特征,将萨尔图油层Ⅲ油层组划分为9个砂岩组和20个单层。研究表明,岁差20 ka周期调控气候波动的机制与超短期基准面升降变化之间存在着一定的协调作用,利用米兰科维奇天文旋回划分出的单层具有严格的等时性。20个单层砂体的展布特征表明,萨尔图油层Ⅲ油层组经历了多期水进—水退事件,与砂岩组的划分具有较好的对应关系。  相似文献   

18.
中国西南地区晚二叠世泥炭地净初级生产力及其控制因素   总被引:2,自引:0,他引:2  
邵龙义  汪浩  LargeDJ 《古地理学报》2011,13(5):473-480
米兰科维奇旋回理论是古环境研究中重要的时间“度量”工具.文中以贵州普安糯东17号煤层和云南富源天佑10号煤层为例,对晚二叠世煤层的地球物理测井信号进行频谱分析,以获得其中的米兰科维奇轨道周期参数.研究发现,测井信号所反映出的煤层灰分含量变化受泥炭地发育时期的米兰科维奇轨道周期(123 ka(偏心率):35.6 ka(斜...  相似文献   

19.
Based on outcrops, drilling, cores, and logging data, the Liangshan Formation and Qixia Formation in the northwestern Sichuan region are analyzed using the theory and method of sequence stratigraphy. Three third-order sequence boundaries are identified, including the basal boundary of Liangshan Formation (type Ⅰ), the boundary between Qi 1 Member and Qi 2 Member (type Ⅱ), and the boundary between the Qixia Formation and the Maokou Formation (type Ⅰ). The Middle Permian Qixia Stage can be subdivided into two third-order sequences (from base to top: SQ1 and SQ2). The SQ1 sequence corresponds to the Liangshan Formation and the Qi 1 Member, and the SQ2 sequence corresponds to the Qi 2 Member. Based on the sequence stratigraphic framework reconstruction of the Qixia Stage, the sedimentary filling pattern of the “onlapping at the base and truncation at the top” is discovered. The contour maps of SQ1 and SQ2 sequence stratigraphic thickness are calculated and plotted. Combined with the stratigraphic filling pattern, the tectonic-palaeogeographical framework of the early SQ1 and SQ2 phases of the Qixia Stage, their geological significance is discussed. At the Qixia Stage, the differentiation between uplift and depression occurs in the northwest and northeast trends. The Hannan uplift, the northern Sichuan uplift, the northern margin uplift and the Guangyuan-Wangcang sag have controlled the sedimentary pattern in the study area. It is indicated that prototype of the Guangyuan-Wangcang Trough began in the Middle Permian Qixia Stage. In the study area, the karst cave-type dolostone reservoir and the karst-type limestone reservoir around the Guangyuan-Wangcang sag rim are favorable exploration areas.  相似文献   

20.
古气候是影响陆相层序发育的一个重要因素,其周期性变化是源于天文旋回的驱动力(即米兰柯维奇假说)。米兰柯维奇天文旋回包括偏心率旋回、倾斜率旋回和岁差旋回,这几个轨道参数所驱动的古气候变迁,分别形成了周期为10万年(或40万年)、4万年及2万年的高频层序,即小层序组、小层序及小层单元,研究结果表明,泌阳断陷湖盆的高频层序是由地球轨道旋回所驱动的古气候周期性变迁所形成的。  相似文献   

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

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