首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 312 毫秒
1.
冯驰  金琦  王艳楠  赵丽娜  吕恒  李云梅 《环境科学》2015,36(5):1557-1564
叶绿素a作为水质参数之一,常用来作为衡量水体富营养化程度的指示标准.利用从太湖及洞庭湖获取的326个实测数据,基于实测遥感反射率对水体光谱进行光学分类,结果表明所采集的样点可分为3种水体类型.结合GOCI的波段设置,建立了不同类型水体的叶绿素a浓度反演模型.水体类型一可以利用490 nm(3波段)和555 nm(4波段)来反演,水体类型二可利用660 nm(5波段)和443 nm(2波段),水体类型三利用745 nm(7波段)和680 nm(6波段).精度分析表明,分类后的平均相对误差明显下降,类型一为38.91%、类型二为24.19%、类型三为22.90%;类型一均方根误差为4.87μg·L-1、类型二为8.13μg·L-1、类型三为11.66μg·L-1;分类前后的总体平均相对误差由49.78%降低到29.59%,总体均方根误差由14.10μg·L-1降低到9.29μg·L-1,分类后反演精度得到了显著提高.利用2013年5月13日8景GOCI影像反演了太湖的叶绿素a浓度,结果表明,2013年5月13日太湖叶绿素a浓度日变化显著,高值区主要集中在竺山湾、梅梁湾、贡湖湾,低值区主要集中在湖心区以及南部区域,10:00以后太湖西南部沿岸的叶绿素a浓度显著降低.这种先分类后反演的方法对于二类水体的模型反演精度的提高具有重要作用.  相似文献   

2.
基于GF-1 WFV影像和BP神经网络的太湖叶绿素a反演   总被引:7,自引:0,他引:7  
叶绿素a浓度是可直接遥感反演的重要水质参数之一,常用来评价湖泊水体的富营养化程度.太湖是典型的二类水体,光学性质复杂,应用一类水体线性反演模式拟合较为片面且难以找到最佳拟合模型.BP神经网络模型具有模拟复杂非线性问题的功能.为研究高分一号卫星16m多光谱相机WFV4结合BP神经网络进行太湖叶绿素a浓度监测的可行性,实验利用GF-1 WFV4影像和实时的地面采样数据,建立了BP神经网络模型,同时采用波段比值经验模型进行对比.经精度检验,BP神经网络模型预测值与实测值之间的可决系数R2高达0.9680,而波段比值模型的R2为0.9541,且均方根误差RMSE由波段比值模型的18.7915降低为BP神经网络模型的7.6068,平均相对误差e也由波段比值模型的19.16%降低为BP神经网络模型的6.75%.结果证明,GF-1 WFV4影像应用BP神经网络模型反演太湖叶绿素a浓度较波段比值模型精度有所提高.将经过水体掩膜的GF-1 WFV4影像用于训练好的BP神经网络反演太湖叶绿素a浓度分布,结果显示,叶绿素a高浓度区集中分布在湖心区北部、竺山湾、梅梁湾区域,与之前的研究一致.本文研究结果验证了采用BP神经网络模型对GF-1 WFV4影像进行太湖叶绿素a浓度反演的可行性.  相似文献   

3.
基于稀疏表达的水体遥感反射率高光谱重构及其应用   总被引:2,自引:2,他引:0  
高光谱重构技术可以有效地突破多光谱卫星传感器波段设置的限制,获得更多更有效的地物光谱信息.本研究基于稀疏表达方法提出了一种针对水体遥感反射率的高光谱重构算法,以太湖、杭州湾的原位水体光谱数据为数据源,在5种常用水色传感器(Sentinel-2A MSI、MERIS、MODIS Aqua、GOCI以及ⅦRS)上进行了高光谱重构实验,最后将该算法应用于GOCI数据,进行了算法适用性验证.结果表明:(1)基于稀疏表达的高光谱重构算法可以在不利用实测光谱数据的情况下实现高光谱重构,光谱重构精度高于多元回归光谱重构算法;(2)基于稀疏表达的高光谱重构算法在5种水色传感器上都取得了较好的效果,平均相对误差均在10%以下,均方根误差均在0.005 sr-1以下;(3)相比于原始GOCI多光谱数据,经稀疏表达高光谱重构后的GOCI数据在叶绿素a浓度和总悬浮物浓度估算精度上有不同程度提升.其中对叶绿素a浓度估算而言,平均相对误差从80.6%减少至51.5%,均方根误差从12.175μg·L~(-1)减少至7.125μg·L~(-1);对悬浮物浓度估算而言,平均相对误差从19.1%减少至18.8%,均方根误差从29.048 mg·L~(-1)减少至28.596 mg·L~(-1).  相似文献   

4.
利用Hyperion高光谱数据的三波段法反演太湖叶绿素a浓度   总被引:6,自引:3,他引:3  
杜聪  王世新  周艺  阎福礼 《环境科学》2009,30(10):2904-2910
以2004-08-19太湖野外试验所获取的水质数据(叶绿素a浓度7.8~154.3μg.L-1,总悬浮物浓度65.0~190.2 mg.L-1,N=38)和同步的Hyperion星载高光谱数据为研究对象,利用三波段算法反演太湖水体的叶绿素a浓度.通过分析太湖固有光学量的特点,提出适用于太湖的3个特征波段的选择依据,并对波段进行优化计算,在此基础上建立了三波段统计模型,最后对模型的反演精度进行分析与评价.结果表明,Hyperion的B34(691.37 nm)、B37(721.90 nm)和B50(854.18 nm)组成三波段模型变量与叶绿素a浓度具有最高的相关系数(r=0.934),模型的决定系数(R2)和均方根误差(RMSE)分别为0.872和13.93μg.L-1,其反演精度优于传统经验统计模型,如比值模型(R2=0.844,RMSE=15.41μg.L-1)和一阶微分模型(R2=0.831,RMSE=16.00μg.L-1).研究结果证实了三波段法适用于内陆富营养化浑浊水体和Hyperion高光谱数据,为今后更精确地反演内陆水体的叶绿素a浓度提供了参考依据.  相似文献   

5.
高光谱重构技术可以有效的突破多光谱卫星传感器波段设置的限制,获得更多更有效的地物光谱信息。本研究基于稀疏表达方法提出了一种针对水体遥感反射率的高光谱重构算法,以太湖、杭州湾的原位水体光谱数据为数据源,在5种常用水色传感器(Sentinel-2A MSI、MERIS、MODIS Aqua、GOCI以及VIIRS)上进行了高光谱重构实验,最后将该算法应用于GOCI数据,进行了算法适用性验证。结果表明:(1)基于稀疏表达的高光谱重构算法可以在不利用实测光谱数据的情况下实现高光谱重构,光谱重构精度高于多元回归光谱重构算法;(2)基于稀疏表达的高光谱重构算法在5种水色传感器上都取得了较好的效果,平均相对误差均在10%以下,均方根误差均在0.005sr-1以下;(3)相比于原始GOCI多光谱数据,经稀疏表达高光谱重构后的GOCI数据在叶绿素a浓度和总悬浮物浓度估算精度上有不同程度提升。其中对叶绿素a浓度估算而言,平均相对误差从80.6%减少至51.5%,均方根误差从12.175μg·L-1减少至7.125μg·L-1;对悬浮物浓度估算而言,平均相对误差从19.1%减少至18.8%,均方根误差从29.048 mg·L-1减少至28.596 mg·L-1。  相似文献   

6.
不同方法估算太湖叶绿素a浓度对比研究   总被引:10,自引:2,他引:8  
基于2006-01-07~2006-01-09和2006-07-29~2006-08-01太湖地面实测高光谱数据以及同步水质参数数据,对比分析了三波段模型、两波段模型、反射峰位置法、一阶微分法4种方法用于估算太湖叶绿素a浓度的精度,并讨论其应用于遥感影像中估算叶绿素a浓度的可行性. 2次采样3类水色参数总悬浮物、叶绿素a浓度和有色可溶性有机物在440 nm处吸收系数的变化范围分别为12.24~285.20 mg·L-1、 4.83~155.11 μg·L-1和0.27~2.36 m-1.前述4种方法在反演太湖水体的叶绿素a浓度时都取得较高的精度;决定系数分别为:0.813、 0.838、 0.872、 0.819,均方根误差分别为:13.04、 12.12、 13.41、 12.13 μg·L-1;相对误差分别为:35.5%、 34.9%、 24.6%、 41.8%.反射峰位置法估算精度最高,但应用到叶绿素a浓度遥感影像估算比较困难.三波段模型和两波段模型的反演结果优于传统的一阶微分法,且在卫星遥感反演中具有良好的应用前景.根据模拟MERIS数据,分别得到最优三波段模型[R-1(665)-R-1(709)]×R(754)和两波段模型R(709)/R(681),其决定系数、均方根误差、相对误差分别为0.788、 13.87 μg·L-1、 37.3%和0.815、 12.96 μg·L-1、 34.8%,反映了MERIS数据能非常好地应用于太湖这类浑浊二类水体叶绿素a浓度的精确估算.  相似文献   

7.
刘忠华  李云梅  檀静  郭宇龙  周莉  刘阁 《环境科学》2012,33(9):3000-3008
总悬浮物浓度是水体重要的水质参数.本研究利用太湖春季、秋季和巢湖夏季多期野外实测数据,通过对生物光学模型进行合理的简化构建适用于太湖、巢湖水体总悬浮物浓度反演的半分析模型,并将该模型应用于MERIS和环境一号卫星高光谱卫星影像上以验证该方法的适用性.结果表明:①针对太湖和巢湖水体,总悬浮物浓度最优反演波段范围为730~832nm(氧气吸收带除外);②针对MERIS数据,波段10(中心波长754 nm)和波段12(中心波长779 nm)均适用于太湖总悬浮物浓度反演,而波段11(中心波长761 nm)由于氧气吸收带的影响不适用于总悬浮物浓度反演;③针对太湖MERIS数据,模型反演结果的相对误差基本上呈现出随距离卫星过境时间增大而逐渐增加的趋势,在卫星过境时间正负3 h内测量的样点,模型反演结果的相对误差均在50%以内,而时间差超过3 h,相对误差则逐渐增大到50%以上;④环境一号卫星高光谱数据17个波段(B83~B99)均能够对巢湖总悬浮物浓度进行较好的反演,其反演效果要好于MERIS数据在太湖的反演结果.  相似文献   

8.
环境一号卫星在大型水体水环境监测与评价中具有独特的优势。为探求遥感影像在水体叶绿素a浓度反演中的应用,基于环境一号卫星CCD数据和同步实测叶绿素a浓度值,通过影像辐射定标、大气校正和几何精校正等预处理获取水体反射率,分别将单波段和不同特征波段组合的反射率与实测叶绿素a值进行皮尔逊相关分析,选取R20.8的波段组合进行建模,通过对3种波段组合反演结果对比和精度验证,发现基于CCD数据第4波段与第3波段反射率比值的二次模型具有良好的反演效果,模型预测值与实测值的最小相对误差为0.76%,平均相对误差10.99%,均方根误差为0.007 6 mg/L,明显低于实测叶绿素a浓度的平均值;最后基于该模型实现了太湖叶绿素a浓度反演,并对叶绿素a的时空分布进行了初步分析。  相似文献   

9.
太湖水域叶绿素a浓度的遥感反演研究   总被引:5,自引:0,他引:5  
利用太湖水域MODIS遥感数据的各波段反射率组合计算值,与实测的叶绿素a浓度进行相关性分析,找到相关性最好的反射率组合,建立反演太湖叶绿素a浓度的遥感模型.结果表明,利用MODIS数据可以较好地实现对太湖水域叶绿素a浓度的定量反演计算,并以MODIS数据第3、第17波段的反射率组合作为遥感指数建立了反演叶绿素a浓度的模型.第3、第17波段的波长范围分别为459nm~479nm、890nm~920nm,这一波段选择与以往使用TM数据得到的结论有所不同.  相似文献   

10.
利用高光谱反演模型评估太湖水体叶绿素a浓度分布   总被引:3,自引:1,他引:2  
叶绿素a浓度是评价水体富营养化和初级生产力的一个重要参数,高光谱遥感是获取叶绿素a浓度的有效手段.为建立太湖水域叶绿素a的最佳高光谱估算模型,选取2015年5—7月共计60组同步实测高光谱数据和叶绿素a浓度数据,在地面光谱反射率和叶绿素a浓度相关性分析的基础上,使用2∶1的数据样本进行太湖水域叶绿素a的最佳高光谱估算模型的建立和验证,筛选模型分别为波段比值、三波段、荧光峰位置、峰谷距离、一阶微分、NDCI(Normalized Difference Chlorophyll Index)、峰面积、荧光峰高度、WCI(Water Chlorophyll-a Index)和四波段模型.结果表明,建模得到的四波段模型决定系数最高,峰面积模型的决定系数相对最低;四波段模型的反演精度最高,均方根误差(RMSE)为0.00376 mg·L~(-1),平均绝对误差(MAPE)为27.86%,而WCI模型的反演精度相对最低,RMSE为0.01231 mg·L~(-1),MAPE为45.11%.将反演精度最高的四波段模型应用于2015年8月3日的两景HSI(Hyperspectral Imaging Radiometer)高光谱影像数据,也得到较高精度,利用同步实测叶绿素a浓度验证的决定系数为0.7643,RMSE为0.00433 mg·L~(-1),MAPE为45.62%.在春、夏季叶绿素对水体光学特性占主导作用且叶绿素分布均匀的情景下,本研究可为太湖水域叶绿素a的高光谱反演和水环境监测提供有价值的参考,其它季节水体光谱特点的研究尚待进一步开展.  相似文献   

11.
The potential harm of heavy metals is a primary concern in application of sludge to the agricultural land. A pot experiment was conducted to evaluate the effect of two sludges on fractionation of Zn and Cu in soil and their phytotoxicity to pakchoi. The loamy soil was mixed with 0%, 20%, 40%, 60% and 80% (by weight) of digested sewage sludge (SS) and composted sludge (SC). The additions of both sludges caused a significant raise in all fractions, resulting in that exchangeable (EXCH) and organic bound (OM) became predominance of Zn and organic bound Cu occupied the largest portion. There was more available amount of Zn and Cu in SS treatments than SC treatments. During the pot experiment, the concentration of Zn in EXCH, carbonate (CAR) and OM and Cu in EXCH and OM fractions decreased in all treatments, so their bioavailability reduced. Germination rate and plant biomass decreased when the addition rate was high and the best yield appeared in 20% mixtures at the harvest of pakchoi. The two sludges increased tissue contents of Zn and Cu especially in the SS treatments. Zn in pakchoi was not only in relationship to ΔEXCH and ΔCAR forms but also in ΔOM forms in the sludge-soil mixtures. Tissue content of Cu in pakchoi grown on SC-soils could not be predicted by ΔEXCH. These correlation rates between Zn and Cu accumulation in pakchoi and variation of different fractions increased with time, which might indicate that sludges represented stronger impacts on the plant in long-term land application.  相似文献   

12.
A hydroponic experiment was carried out to study intraspecific differences in the effects of different concentrations of cadmium (Cd)(0-10 mg/L) and arsenate (As(V)) (0-8 mg/L) on the growth parameters and accumulation of Cd and As in six wheat varieties Jing-9428, Duokang-1, Jingdong-11, Jing-411, Jingdong-8 and Zhongmai-8. The endpoints of wheat seedlings, including seed germination,biomass, root length and shoot height, decreased with increasing the Cd and As concentrations. Significant differences in seed germination, biomass, root length, shoot height and the accumulation of Cd and As were observed between the treatments and among the varieties (p < 0.05). The lethal dosage 50% were about 20, 80, 60, 60, 80 and 20 mg As/L for Jing-9428, Duokang-1, Jingdong-11,Jing-411, Jingdong-8 and Zhongmai-8, respectively, and the corresponding values for Cd were about 30, 80, 20, 40, 60 and 10 mg Cd/L, respectively. Among the six varieties, Duokang-1 was found to be the most resistant to Cd and As toxicity, and Zhongmai-8 was the most sensitive to Cd and As co-contamination. The resistance of the six varieties was found dependant on the seedling uptake of Cd and As. Duokang-1 was the most suitable for cultivation in Cd and As co-contaminated soils.  相似文献   

13.
Toxic effects of two agrochemicals on nifH gene in agricultural black soil were investigated using denaturing gradient gel electrophoresis (DGGE) and sequencing approaches in a microcosm experiment. Changes of soil nifH gene diversity and composition were examined following the application of acetochlor, methamidophos and their combination. Acetochlor reduced the nifH gene diversity (both in gene richness and diversity index values) and caused changes in the nifH gene composition. The effects of acetochlor on nifH gene were strengthened as the concentration of acetochlor increased. Cluster analysis of DGGE banding patterns showed that nifH gene composition which had been affected by low concentration of acetochlor (50 mg/kg) recovered firstly. Methamidophos reduced nifH gene richness that except at 4 weeks. The medium concentration of methamidophos (150 mg/kg) caused the most apparent changes in nifH gene diversity at the first week while the high concentration of methamidophos (250 mg/kg) produced prominent effects on nifH gene diversity in the following weeks. Cluster analysis showed that minimal changes of nifH gene composition were found at 1 week and maximal changes at 4 weeks. Toxic effects of acetochlor and methamidophos combination on nifH gene were also apparent. Different nifH genes (bands) responded differently to the impact of agrochemicals: four individual bands were eliminated by the application of the agrochemicals, five bands became predominant by the stimulation of the agrochemicals, and four bands showed strong resistance to the influence of the agrochemicals. Fifteen prominent bands were partially sequenced, yielding 15 different nifH sequences, which were used for phylogenetic reconstructions. All sequences were affiliated with the alpha- and beta-proteobacteria, showing higher similarity to eight different diazotrophic genera.  相似文献   

14.
The effects of arbuscular mycorrhizal (AM) fungus (Glomus mosseae) and phosphorus (P) addition (100 mg/kg soil) on arsenic (As) uptake by maize plants (Zea mays L.) from an As-contaminated soil were examined in a glasshouse experiment.Non-mycorrhizal and zero-P addition controls were included.Plant biomass and concentrations and uptake of As,P,and other nutrients,AM colonization,root lengths,and hyphal length densities were determined.The results indicated that addition of P significantly inhibited root colonization and development of extraradical mycelium.Root length and dry weight both increased markedly with mycorrhizal colonization under the zero-P treatments,but shoot and root biomass of AM plants was depressed by P application.AM fungal inoculation decreased shoot As concentrations when no P was added,and shoot and root As concentrations of AM plants increased 2.6 and 1.4 times with P addition,respectively.Shoot and root uptake of P,Mn,Cu,and Zn increased,but shoot Fe uptake decreased by 44.6%,with inoculation, when P was added.P addition reduced shoot P,Fe,Mn,Cu,and Zn uptake of AM plants,but increased root Fe and Mn uptake of the nonmycorrhizal ones.AM colonization therefore appeared to enhance plant tolerance to As in low P soil,and have some potential for the phytostabilization of As-contaminated soil,however,P application may introduce additional environmental risk by increasing soil As mobility.  相似文献   

15.
Several main metabolites of benzo[a]pyrene (BaP) formed by Penicillium chrysogenum, Benzo[a]pyrene-1,6-quinone (BP 1,6- quinone), trans-7,8-dihydroxy-7,8-dihydrobenzo[a]pyrene (BP 7,8-diol), 3-hydroxybenzo[a]pyrene (3-OHBP), were identified by high-performance liquid chromatography (HPLC). The three metabolites were liable to be accumulated and were hardly further metabolized because of their toxicity to microorganisms. However, their further degradation was essential for the complete degradation of BaP. To enhance their degradation, two methods, degradation by coupling Penicillium chrysogenum with KMnO4 and degradation only by Penicillium chrysogenum, were compared; Meanwhile, the parameters of degradation in the superior method were optimized. The results showed that (1) the method of coupling Penicillium chrysogenum with KMnO4 was better and was the first method to be used in the degradation of BaP and its metabolites; (2) the metabolite, BP 1,6-quinone was the most liable to be accumulated in pure cultures; (3) the effect of degradation was the best when the concentration of KMnO4 in the cultures was 0.01% (w/v), concentration of the three compounds was 5 mg/L and pH was 6.2. Based on the experimental results, a novel concept with regard to the bioremediation of BaP-contaminated environment was discussed, considering the influence on environmental toxicity of the accumulated metabolites.  相似文献   

16.
Sorption of chlorotoluron in ammonium sulfate, urea and atrazine multi-solutes system was investigated by batch experiments. The results showed application of nitrogen fertilizers to the soil could affect the behavior of chlorotoluron. At the same concentration of N, sorption of chlorotoluron decreased as the concentration of atrazine increased on the day 0 and 6 in soil, respectively. The sorption of chlorotoluron increased from 0 to 6 d when soils were preincubated with deionized water, ammonium sulfate and urea solution for 6 d. That indicated incubation time was one of the most important factors for the sorption of chlorotoluron in nitrogen fertilizers treatments. The individual sorption isotherms of chlorotoluron in rubbery polymer and silica were strictly linear in single solute system, but there were competition sorption between pesticides or between pesticides and nitrogen fertilizers. That indicated the sorption taken place by concurrent solid-phase dissolution mechanism and sorption on the interface of water-organic matter or water-mineral matter.  相似文献   

17.
Laogang landfill near Shanghai is the largest landfill in China, and receives about 10000 t of daily garbage per day, Samples of topsoil and plants were analyzed to evaluate mercury pollution from the landfill. For topsoil samples, there were significant correlations among total mercury (HgT), combinative mercury (Hgc) and gaseous mercury (HgG), and content of total organic carbon (TOC), but, no significantly relationship was found between Hg content and filling time. Hg content changes in vertical profiles with time showed that the average Hgv of profiles 1992, 1996, and 2000 was similar, but their average HgG was quite different. HgT was significantly correlated with Hgc in profile 1992 and 2000, and Hgv was significantly correlated with Hg6 in profile 1996. HgG/Hgv ratio in profile samples decreased in the order of (HgG,/HgT)1992〉(HgG/HgT)1996〉〉(HgG/HgT)2000. A simple outline of Hg release in landfill could be drawn: with increasing of filling time, degradation undergoes different biodegradation, accordingly, gaseous mercury goes through small, more, and small proportion to total mercury. Distribution of Hg in plants was inhomogeneous, following the order of leaf〉root〉stem. The highest value of leaf may be associated with higher atmospheric Hg from landfill. Ligneous plants (e.g. Phyllostachys glanca, Prunus salicina and Ligustrum lucidum) are capable of enriching more Hg than herbaceous plants.  相似文献   

18.
Phytoremediation is a potential cleanup technology for the removal of heavy metals from contaminated soils.Bidens maximowicziana is a new Pb hyperaccumulator,which not only has remarkable tolerance to Pb but also extraordinary accumulation capacity for Pb.The maximum Pb concentration was 1509.3 mg/kg in roots and 2164.7 mg/kg in overground tissues.The Pb distribution order in the B. maximowicziana was:leaf>stem>root.The effect of amendments on phytoremediation was also studied.The mobility of soil Pb and the Pb concentrations in plants were both increased by EDTA application.Compared with CK(control check),EDTA application promoted translocation of Pb to overground parts of the plant.The Pb concentrations in overground parts of plants was increased from 24.23-680.56 mg/kg to 29.07-1905.57 mg/kg.This research demonstrated that B.maximowicziana appeared to be suitable for phytoremediation of Pb contaminated soil,especially,combination with EDTA.  相似文献   

19.
Decomposition of alachlor by ozonation and its mechanism   总被引:1,自引:0,他引:1  
Decomposition and corresponding mechanism of alachlor, an endocrine disruptor in water by ozonation were investigated. Results showed that alachlor could not be completely mineralized by ozone alone. Many intermediates and final products were formed during the process, including aromatic compounds, aliphatic carboxylic acids, and inorganic ions. In evoluting these products, some of them with weak polarity were qualitatively identified by GC-MS. The information of inorganic ions suggested that the dechlorination was the first and the fastest step in the ozonation of alachlor.  相似文献   

20.
The influence of the nonionic surfactant Tween 80 on pentachlorophenol (PCP) oxidation catalyzed by horseradish peroxidase was studied. The surfactant was tested at concentrations below and above its critical micelle concentration (CMC). Enhancement of PCP removal was observed at sub-CMCs. The presence of Tween 80 in the reaction mixture reduced enzyme inactivation which occurred through a combination of free radical attack and sorption by precipitated products. A simple first-order model was able to simulate time profiles for enzyme inactivation in the presence or absence of Tween 80. At supra-CMCs, the surfactant caused noticeable reductions in PCP removal, presumably through micelle partitioning of PCP which precluded the hydrophobic PCP molecule from interacting with the enzyme.  相似文献   

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

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