首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
于水温(18.0±1.0)?C、盐度20、p H 8.1±0.2条件下,以日本囊对虾秋繁同生群仔虾[总长(7.643±0.639)mm]为实验对象,以离水干露为胁迫因子,在确认其离水干露160min内的复水成活率为100%和离水干露240min时的复水存活率略低于5%后,借助显微扫描像素测量技术和多元分析方法定量研究了A、B、C、D、E实验群体(依次为干露处理160—180、180—200、200—220、220—240min时段内的死亡群体和干露处理240min时的存活群体)间形态表型特征的差异。结果表明:(1)在所涉15项形态性状中,各实验群体间均无显著差异(P0.05)的性状多达10项,依次为第一触角柄长、额剑长、眼径、第一腹节长、第二腹节长、第三腹节长、第四腹节长、第五腹节长、尾节长和尾节高;(2)在所涉17项形态比例指标中,各实验群体间均无显著差异(P0.05)的指标仅为3项,依次为头胸甲长/总长、头胸甲高/头胸甲长和尾节高/尾节长,群体间的欧氏距离均达到极显著水平(P0.01),且有随耐干露性能差异的增大而呈显著增大的趋势;(3)经主成分分析,提取到的4个特征值均大于1的主成分,累计贡献率达91.018%,其中第1主成分的贡献率高达57.974%,其载荷绝对值大于0.5的主要变量占形态比例指标总数的64.706%;(4)将A、B、C、D实验群体统归为干露处理淘汰群,E实验群体定义为干露处理选留群。采用逐步判别法,所建的Fisher分类函数方程组可较清晰地区分淘汰群和选留群个体,其中选留群个体的判别准确率P1、P2分别为100%和84.75%,淘汰群个体的判别准确率P1、P2分别为82%和100%,两者综合判别准确率为91%。  相似文献   

2.
以水温16oC、p H 8.1±0.2、盐度20的砂滤海水为实验用水,以体长(7.643±0.639)mm的日本囊对虾(Marsupenaeus japonicus)秋繁同生群仔虾为实验对象,经确认其在水温由16oC瞬降至14oC后并保持96h恒温的最终成活率略低于5%后,以此为水温骤降胁迫处理方式,采取显微扫描像素测量技术和多元分析方法定量研究了A、B、C、D、E实验群体,依次为温降胁迫处理后(0—24)h、(24—48)h、(48—72)h、(72—96)h时段内的死亡群体和温降胁迫处理后96h时的存活群体间形态特征的差异。结果表明:(1)各实验群体在本研究所涉15项形态性状中均无显著差异(P0.05)仅为第一触角柄长、额剑长和眼径,在所涉17项形态比例指标均无显著差异(P0.05)仅为尾节高/尾节长;(2)随耐寒性能差异的增大,各实验群体间的欧氏距离也随之增大,且均达到极显著水平(P0.01);(3)经主成分分析,提取到的4个特征值均大于1的主成分,累计贡献率达87.590%,其中载荷绝对值大于0.80的主要影响变量均仅存在于第1主成分中;(4)将E实验群体定义为耐寒处理选留群,其余实验群体统归为耐寒处理淘汰群,采用逐步导入判别法建立的Fisher分类函数方程组可较清晰地区分淘汰群和选留群个体,其中选留群和淘汰群个体的判别准确率P1分别为99%和88.75%,P2分别为89.80%和98.88%,两者综合判别准确率为93.88%。  相似文献   

3.
于水温(18.0±1.0)°C、盐度20、pH 8.1±0.2条件下,以水流为胁迫因子(时长1min),以日本囊对虾秋繁同生群仔虾[总长(7.643±0.639)mm]为实验对象,以溢水口(实验初始时刻仔虾放置处)为起点,按等距离间距法将自制的水流测定装置(总长1m)依次划分为A(0—25cm)、B(25—50cm)、C(50—75cm)、D(75—100cm)和E(100cm)等五个区段,在确认实验终了时刻分布于A区段内仔虾数量占实验仔虾总数5%的水流速度为0.823cm/s后,以此为实验流速,借助显微扫描像素测量技术和多元分析方法定量研究了A、B、C、D、E实验群体(依次为实验终止时刻分布于A、B、C、D、E区段内的仔虾)个体间形态比例性状间的差异。结果表明:(1)在所涉15项形态测量指标中,各实验群体间均无显著差异(P0.05)的形态性状共计8项,依次为X_4(眼径)、X_5(头胸甲长)、X_6(头胸甲高)、X_7(第一腹节长)、X_8(第二腹节长)、X_9(第三腹节长)、X_~(10)(第四腹节长)和X_14(腹节高);(2)在所涉17项形态比例指标中,实验群体间均无显著差异(P0.05)的形态比例性状共计7项,依次为C_2(额剑长/总长)、C_9(第五腹节长/总长)、C_(10)(尾节长/总长)、C_(11)(尾扇长/总长)、C_(15)(头胸甲高/头胸甲长)、C_(16)(腹节高/第一腹节长)和C_(17)(尾节高/尾节长);(3)经主成分分析,提取到的5个特征值均大于1的主成分,累计贡献率达80.795%,其中第1主成分39.561%,其载荷绝对值大于0.5的主要影响变量占形态比例性状总数的47.059%;(4)将A实验群体定义为水流胁迫处理选留群,B、C、D、E实验群体统归为水流胁迫处理淘汰群。采用逐步判别法,以判别贡献率较大的C_1(第一触角柄长/总长)、C3(眼径/总长)、C7(第三腹节长/总长)、C13(额剑长/头胸甲长)为自变量,所建的Fisher分类函数方程组可较清晰地区分淘汰群和选留群个体,其中选留群的判别准确率P_1、P_2分别为98%和84.85%,淘汰群的判别准确率P_1、P_2分别为82.25%和97.63%,两者综合判别准确率为90.12%。综上可知,借助形态表型分型可实现日本囊对虾仔虾不同抗流性能群体间的筛选。  相似文献   

4.
以体长(1.932±0.204)cm、体质量(1.386±0.055)g的棘胸蛙蝌蚪为实验动物,在水温(24±0.2)℃、DO(7.30±0.01)mg/L、pH 7.30±0.01条件下,采用静水停食法开展了氨氮对蝌蚪的急性攻毒实验,并以此为基础,测定了不同氨氮质量浓度胁迫下棘胸蛙蝌蚪的排氨率、耗氧率、氧氮比及窒息点。结果表明:(1)蝌蚪对氨氮急性攻毒具明显的运动避毒行为,濒死个体的背部皮肤和肝脏均具明显的色变症状;(2)氨氮对蝌蚪的急性致死率具明显的剂量—时间效应,24h、48h、72h、96h LC50依次为177.5、151.7、148.6和146.8mg/L,毒性时段蓄积程度系数呈持续下降趋势,24—48h、48—72h、72—96h的时段MAC值分别为84.04%、10.1%和5.86%;(3)蝌蚪夜均、昼均、日均及时段排氨率随氨氮质量浓度增加均呈阶梯式下降趋势,其中时段排氨率与对照组均无显著差异的为9.80mg/L实验组,夜均、昼均、日均及时段排氨率与对照组均无显著差异的仅为2.45mg/L实验组(P0.05);(4)氨氮对蝌蚪呼吸耗氧具低毒兴奋效应,其夜均、昼均和日均耗氧率随氨氮质量浓度增加均呈先升后降趋势,峰值氨氮质量浓度范围均为4.90—7.35mg/L(P0.05),与对照组均无显著差异的仅为14.70mg/L实验组(P0.05),窒息点含氧量随氨氮质量浓度增加呈先降后升趋势,谷底出现于4.90mg/L实验组(P0.05),与对照组无显著差异的仅为9.80mg/L实验组(P0.05);(5)实验期间,蝌蚪氧氮比(O:N)波动于18.87—25.34之间,昼均和日均氧氮比(O:N)随氨氮质量浓度增加均依次呈"上升—稳定—下降—再稳定"之趋势,两者的峰值氨氮质量浓度范围分别为4.90—14.70mg/L和4.90—9.80mg/L,夜均氧氮比(O:N)呈先升后降趋势,峰值出现于7.35mg/L实验组(P0.05)。  相似文献   

5.
中国五大湖三角帆蚌形态差异分析   总被引:49,自引:8,他引:49       下载免费PDF全文
运用三种多元分析方法 ,采用 1 0个形态比例性状 ,对中国五大湖地理种群与诸暨人工群体三角帆蚌的形态差异进行了比较研究。聚类分析与主成分分析结果表明 ,太湖种群与巢湖种群形态最近 ,而洪泽种群的趋异程度最大。主成分分析构建了 3个主成分 ,其贡献率 :主成分 1为 2 7 87% ,主成分 2为 1 4 2 6% ,主成分 3为 1 1 91 % ,累积贡献率为 5 4 0 4 %。判别分析结果表明 ,5种群及人工群体形态差异显著 (P <0 0 1 )。建立了 5种群及人工群体的判别函数 ,其判别准确率P1 为 73 33%— 86 67% ,P2 为 72 73%— 1 0 0 % ,综合判别率为83 33%。  相似文献   

6.
中国5 个海区不同群体毛蚶形态差异分析   总被引:8,自引:0,他引:8       下载免费PDF全文
运用多变量形态度量学分析方法,采用10个形态性状,对中国5个海区5个野生毛蚶群体间的形态差异进行了比较研究.聚类分析和主成分分析结果表明:天津塘沽群体和山东青岛群体形态最为接近,广西北海群体的趋异程度最大.主成分分析建立了3个主成分--主成分1、主成分2和主成分3,其贡献率分别为:34.70%、19.80%、15.00%,积累贡献率为69.50%. 判别分析结果显示,5个群体间的形态差异显著(P<0.01).建立了5个群体毛蚶的判别函数,其判别准确率P1为45.45%~95.45%、P2为36.36%~95.45%,综合判别率为74.50%.  相似文献   

7.
不同地区毛蚶群体的形态差异的比较研究   总被引:1,自引:0,他引:1  
运用多变量形态度量学分析方法,采用13个形态性状,对山东省4个地区:青岛红岛、威海乳山、东营河口、东营广饶,以及辽宁省葫芦岛和海南省三亚6个不同地区的毛蚶(Scapharca subcrenata)群体间的形态差异进行了比较研究。聚类分析和主成分分析结果表明:河口群体和广饶群体形态最为接近,海南群体的趋异程度最大。主成分分析建立了3个主成分—壳质量、楯面长、前端到腹缘,其贡献率分别为:31.25%,24.46%,12.39%,积累贡献率为68.10%。判别分析结果显示,6个毛蚶群体的形态差异显著(P<0.01)。建立了6个毛蚶群体的判别函数,其判别准确率P1为82.76%~98.63%,P2为80.00%~100%,综合判别率为91.7%。研究结果为毛蚶的地理群体识别、遗传特征研究、种质资源保护与恢复提供重要的依据。  相似文献   

8.
从抗氧化酶活性的变化揭示氨氮胁迫对刀鲚稚幼鱼毒性效应的机制,为其健康养殖提供科学依据。分别以刀鲚稚鱼(体长2.33 cm±0.25 cm)和幼鱼(体长5.12 cm±0.45 cm)为试验对象,进行96 h氨氮急性毒性试验,得出各自的96 h半致死浓度96hLC_(50)以及安全浓度。然后,分别用刀鲚稚幼鱼40%96 h LC_(50)进行24h的氨氮胁迫试验并取样,然后,将胁迫组剩余的刀鲚放入曝气自来水中进行24 h的恢复试验,胁迫及恢复试验均取全鱼测定其总抗氧化能力(T-AOC)、过氧化氢酶(CAT)、总超氧化物歧化酶(T-SOD)、谷胱甘肽过氧化物酶(GSH-Px)非特异性免疫酶的活性。结果显示:在胁迫24 h后稚幼鱼T-AOC活力显著下降(P0.05),GSH-Px活力上升但不显著(P0.05),T-SOD、CAT活力均显著上升(P0.05);恢复24 h后各种酶的活性都在向正常水平恢复,但是最终都没有恢复到对照组水平。研究获得了氨氮对刀鲚稚幼鱼的急性毒性结果和在氨氮胁迫下抗氧化酶活性的变化规律,氨氮胁迫引起了氧化应激,恢复后抗氧化能力增强,缓解了氧化应激,但并未恢复到对照组水平。  相似文献   

9.
采用聚类分析、主成分分析和判别分析3种分析方法,对我国东南沿海5个野生波纹巴非蛤群体10个形态性状进行比较研究。聚类分析和主成分分析结果表明,广东东莞与福建厦门群体、广东湛江与海南海口群体的形态最为接近,广西北海群体的趋异程度最大。主成分分析构建了3个主成分,贡献率:主成分1为23.237%,主成分2为16.109%,主成分3为13.009%,累计贡献率为52.355%。判别分析结果表明,5群体间形态差异显著(P0.01)。建立了5群体的判别函数,判别准确率P1为47.1%—80.0%,P2为45.7%—80.0%,5群体的综合判别率为63.8%。Mantel检验结果表明欧氏距离与地理直线距离无明显相关性(r=0.0822,P=0.3820)。广西北海群体明显地与其它4个地方分开,可能已形成一个新的亚种。本研究可为我国巴非蛤种质资源保护及遗传育种提供依据。  相似文献   

10.
本文研究了氨氮对方斑东风螺(Babylonia areolata)的急性毒性及不同浓度的氨氮溶液(22、47.5、102和220 mg/L)对其体内六种免疫酶(CAT、ACP、AKP、T-SOD、GSH-PX和POD)活性的影响。实验结果表明:氨氮浓度越大,其毒性作用越强,方斑东风螺死亡率越高;500 mg/L氨氮处理组于24 h便出现(20.3±2.1)%死亡率,并且96 h后存活率降为0。而相同条件下,随试验时间的延长,方斑东风螺死亡率也越高。氨氮浓度和处理时间均对方斑东风螺免疫酶活力有显著的影响(P < 0.05)。与对照组相比,处理组CAT和AKP活性均表现出“诱导—抑制”的趋势,ACP和T-SOD活性表现出“抑制—诱导”的趋势,CSH-PX活性呈现“诱导—抑制—诱导”的趋势,而氨氮对处理个体POD活性的影响整体表现为抑制作用。可见,养殖水体中的氨氮会对方斑东风螺免疫酶活性产生较大的影响。本研究从生态现象和生理指标相结合的方式阐述了氨氮对方斑东风螺个体的致毒特性,将对方斑东风螺的养殖具有实用意义,充实了氨氮对贝类毒理实验免疫性能层面的基础数据,可为其他贝类氨氮胁迫实验提供数据基础和科学依据。  相似文献   

11.
采用PCR扩增、克隆获得了岱衢族和闽-粤东族大黄鱼(Pseudosciaena crocea)线粒体基因全长序列,并扩增了各25条样本的高变区域,测序结果可分成3个操作分类单元(Operational Taxonomic Unit,OTU);闽-粤东族含3个类型:D1-N1、D1-N2和D2型;而岱衢族样品均为D1-N1型;设计出能区分不同分型的鉴别引物J1和J2,且优化了鉴别条件。结果表明:根据凝胶电泳的图片,如果样品分型为D1-N2和D2型,则100%为闽-粤东族;如果样品分型为D1-N1型,则该样品为岱衢族的概率为86.2%,为闽-粤东族的概率为13.8%。研究结果为岱衢族大黄鱼的种质资源保护、良种选育和市场销售等方面提供技术支持。  相似文献   

12.
黄海小型水母的分布特征   总被引:3,自引:0,他引:3  
孙松  张芳  李超伦  杨波  吉鹏 《海洋与湖沼》2012,43(3):429-437
根据"973"项目2000年10月—2007年8月浮游生物大网样品,对黄海小型水母类的种类组成及其丰度的时空变化进行了研究,讨论了各优势水母种类在黄海的分布格局。结果表明,黄海测区内小型水母整体丰度较低,且主要分布在50m等深线以浅海域,各小型水母类群如水螅水母类,管水母类和栉水母类以及10种优势种类的季节更替非常明显。与其它海域相比,黄海小型水母丰度水平最低,平均为0.8ind/m3,占整个浮游动物丰度的比例小于0.5%,因此从丰度上看,小型水母在黄海浮游动物功能群中并不占优势。  相似文献   

13.
赵苑  赵丽  张武昌  刘诚刚  魏皓  肖天 《海洋与湖沼》2012,43(6):1030-1038
于2007年3—4月在黄海中部海域采用流式细胞术研究了春季水华过程中聚球藻、微微型真核浮游生物和异养细菌的生物量变化。聚球藻和微微真核型浮游生物的生物量与叶绿素a浓度变化基本呈现相反的趋势,在水华前期较高,水华期迅速下降,直至水华后期又有所升高。异养细菌在整个水华过程中变化较小,生物量在水华期最高,与水柱叶绿素a浓度呈极显著正相关(r=0.319,p<0.01)。水华期这三类微微型浮游生物对浮游植物总碳生物量的贡献很低。纤毛虫和鞭毛虫捕食可能是导致聚球藻和微微型真核浮游生物在水华期生物量降低的主要原因。  相似文献   

14.
有害藻华对海水养殖、人类健康和生态安全构成威胁,已成为全球性的生态灾害.近20年来,米氏凯伦藻(Karenia mikimotoi)、东海原甲藻(Prorocentrum donghaiense)等甲藻形成的有害藻华在东海海域连年暴发,其频率高、危害大,形成及演变机理复杂,对藻华原因种进行监测是开展有害藻华预警和防控的...  相似文献   

15.
采用RACE技术克隆了单细胞绿藻蛋白核小球藻820Rubisco活化酶基因(rca)序列,并用荧光定量PCR方法研究了rca和Rubisco小亚基基因(rbcS)的表达规律。结果获得了1703bp的rcacDNA序列,包括66bp的5’非翻译区、1242bp的开放阅读框和395bp的3’非翻译区。序列比较和分析表明该蛋白核小球藻rca序列与其它绿藻的同源性高达78%—85%;偏好使用以G/C/T结尾的密码子;推测的RCA蛋白等电点和分子量分别为8.44和45.71kDa。荧光定量PCR结果表明随光照时间增加rca与rbcS转录表达逐渐降低;不同盐度处理下rbcS的mRNA量变化不大,而rca在37.5盐度下表达量为25盐度的2.69倍;1.0—3.0mmol/L水杨酸处理rca与rbcS表达均降低。  相似文献   

16.
通过2009-2011年5月下旬-7月下旬辽东湾北部近海10m等深线内的大型水母调查数据,分析了辽东湾北部近海近三年中大型水母资源状况,并探讨了辽东湾大型水母的生态类型。结果显示:辽东湾北部近海大型水母种类主要有海蜇(Rhopilema esculentum)、沙蜇(Nemopilema nomurai)、白色霞水母(Cyanea nozakii)、海月水母(Aurelia sp.1),海蜇和沙蜇是优势种。海蜇幼水母阶段主要集中分布在5m等深线以内的近岸河口水域,随着个体增大有略向深水或密度较小的水域扩散的趋势,仍主要分布在5m等深线两侧水域,属于高温低盐种类。6月份调查海区中发现大量的沙蜇幼水母,随着沙蜇个体增大,7月份调查海区中沙蜇数量大幅度减少。辽东湾海月水母在南部海域出现较多,2010、2011年在北部近海部分海域出现。白色霞水母近几年来辽东湾出现较少,栖息在盐度较高的水域。辽东湾各种大型水母中,沙蜇的生长速度最快。辽东湾海蜇幼水母、沙蜇幼水母的海区出现时间要晚于黄、东海。  相似文献   

17.
基于GOCI的2017年南黄海浒苔演变遥感分析   总被引:3,自引:0,他引:3  
本文应用GOCI高时间分辨率遥感影像,提取浒苔敏感波段并构建合适的评价指标,对NDVI、IGAG、KOSC绿潮主流算法探测能力进行比较,结合目视解译和阈值分割法,对2017年南黄海浒苔信息进行提取并进行演变特征分析。结果表明:NDVI算法的探测能力和稳定性显著优于其他两种,GOCI影像7、8近红外波段的选择对浒苔提取结果存在一定影响,最优为7、5波段组合;2017年浒苔总体经历了“出现-发展-暴发-衰退-消亡”五个阶段,持续时间约为68d,影像中浒苔最早出现在5月13日的盐城外海海域,最后出现在7月12日,最大覆盖面积出现在6月4日,为2363.12km2,在风、流场共同作用下,先向北移动,后在6月下旬沿山东半岛南侧转向东北,前锋位置最终停滞于青岛-烟台-威海一线并逐渐消亡。其发展规律与运移路径与往年相似,但持续时间与暴发时期最大覆盖面积显著少于往年,GOCI影像高时间分辨率的优势使浒苔灾害的逐小时动态监测成为可能。  相似文献   

18.
渔山列岛于2008年成为国家级海洋生态保护区,为了解渔山列岛保护区内潮间带软体动物次级生产力的时空变化。于1982年和2010年的3个季节(春、夏、冬)在渔山列岛进行潮间带软体动物调查,采用多元统计分析软体动物时空差异,以Brey经验公式计算分析软体动物次级生产力(P值,P/B值)。结果表明,年际间物种食性结构发生了明显的改变,1982年潮间带软体动物对次级生产力贡献率最大的物种有单齿螺(Monodonta labio)、锈凹螺(Chlorostoma rustica),占软体动物总次级生产力的42.88%,2010年潮间带软体动物对次级生产力贡献率最大的物种有覆瓦小蛇螺(Serpulorbis imbricata)、条纹隔贻贝(Septifer virgatus),占软体动物总次级生产力的75.96%。1982年,年均栖息密度为2758ind/m2,年均生物量为2100.26g/m2,年均次级生产力为240.04g/(m2·a),年均P/B值为0.63a–1;2010年,年均栖息密度为699ind/m2,年均生物量为1101.85g/m2,年均次级生产力为94.82g/(m2·a),年均P/B值为0.55a–1。多因素方差分析表明软体动物次级生产力年际间差异显著(F=5.761,P0.05),季节间和潮位间差异都不显著(F=0.135,P0.05;F=2.076,P0.05),P/B的值较低,表明该海域软体动物世代更替速度较慢,2010年的P/B值(0.55a–1)低于1982年的P/B值(0.63a–1),表明近30年来群落结构受到影响。  相似文献   

19.
自2009年2月—2010年1月,共采集黄河三角洲流域分布的河蚬498个,采用组织学方法分析其性腺的年周期变化,利用考马斯亮蓝方法进行了蛋白质含量的年周期测量。结果表明,黄河三角洲分布的河蚬雌雄异体,性比为1:1,性腺发育分为6个时期:增殖期、生长期、成熟期、排放期、恢复期和排尽休止期。配子发生期开始于2月份,一年内有两个繁殖高峰期,分别为每年的6月中旬和10月上旬。蛋白质含量无性别差异,闭壳肌和外套膜的蛋白质含量全年都处于较低水平,无明显波动,性腺-内脏团和足的蛋白质含量相对较高,并分别在5月份和5、8月份达到峰值,与河蚬的繁殖盛期显著相关。  相似文献   

20.
鉴于池塘养殖条件下三疣梭子蟹(Portunus trituberculatus)种内相残现象严重,螯足缺失率较高之实际。本研究对三疣梭子蟹的螯足进行了缺失处理,设置了无螯组(C0)、1螯足组(C1)和2螯足组(C2),在实验室内食物资源不足条件下,探讨了螯足数量(螯足缺失)对三疣梭子蟹摄食行为和竞食能力的影响。结果表明:(1)螯足数量对三疣梭子蟹的索食时间无显著影响(P0.05);(2)缺失螯足(C1,C2组)的三疣梭子蟹均可正常摄食,其中C0组蟹可用第一步足代替螯足抱持食物,将食物送入口中,表现出明显的行为补偿策略,但其啃食速度较慢,啃食耗时是C2组蟹的2倍;(3)C0、C1、C2组三疣梭子蟹均有部分个体同时抢占2份食物,其比率与个体具有的螯足数量有关,螯足越多,同时抢占2份食物的比率越高,其中C2组中有20%个体同时抢占了2份食物,而C0组仅有3%个体同时抢占了2份食物;(4)从螯足数量不同的三疣梭子蟹遭同类抢食的次数、比率来看,螯足数量越多的三疣梭子蟹,抢食欲望、抢食能力越强,而且更趋向向双螯均缺失个体发起抢夺性攻击,表现出明显的"强欺弱"现象;(5)螯足缺失直接影响了三疣梭子蟹对食物的保护,其中C0组中50%个体持有的食物可被同类抢走。整体上,缺失单个螯足对三疣梭子蟹摄食和竞食行为的影响不大,而缺失双螯将严重影响了三疣梭子蟹的摄食和竞食能力。  相似文献   

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

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