GENESIS OF NODULAR LIMESTONE IN THE UPPER CAMBRIAN BAOSHAN FORMATION IN WESTERN YUNNAN PROVINCE
-
摘要: 从瘤体与基质的矿物组成、岩石的结构构造出发,深入探讨了滇西施甸地区保山组瘤状灰岩的沉积环境及其形成机理。研究表明,滇西施甸地区上寒武统保山组瘤状灰岩存在3种成因:一种是呈互层状的泥质和灰质沉积物在半固结或弱固结状态下,经差异压实-压溶作用而形成瘤状构造,岩石中各组分含量及后期压溶作用强度的差异导致了瘤状灰岩不同的存在形式;另一种是泥质和灰质沉积物在半固结或弱固结状态下经斜坡滑塌作用形成串珠状、砾屑状灰岩;第三种则是泥灰质沉积层间歇性地暴露干裂,移位磨蚀,并被后期沉积物充填胶结形成瘤状灰岩。Abstract: In this paper, the depositional environment and formation mechanism of the nodular limestone from the Under Cambrian Baoshan Formation in Shidian region, western Yunnan Province was discussed on the basis of mineral components, texture and structure of the node and matrix. Researches show that there are three origins of the nodular limestone in the study area. The first is that the semi-consolidate or weak-consolidate mud and carbonate sediment is altered into the nodular structure by the differential compaction and pressure solution, and the occurrence mode of the nodular limestone is determined by the content of compositions of rock and the pressure solution intensity. The second is that the semi-consolidate or weak-consolidate mud and carbonate sediment is altered into the shoestring and rudaceous nodular limestone by the slope slumping. The third is that the argillo calcareous sediment with exposure and cracking is cemented into the nodular limestone by sediment filling.
-
0. 引言
遥感技术在基础地质工作中应用越来越广泛,但是由于植被等因素的干扰,难以将植被覆盖下的信息从图像中分离出来,植被是遥感对地表探查的一道天然屏障。在基岩被植被完全覆盖的地区,遥感所获得的信息主要是植被信息,给遥感地质勘查工作带来极大不便。因此,找到一种可行的植被剔除方法,还原植被覆盖下的基岩信息,可更好地发挥遥感技术在基础地质工作中的应用能力。
过去研究人员主要基于像元的植被指数对植被等干扰信息进行提取,但这种基于像元的模式不能有效地提高提取的精度。因为遥感图像中每个像元一般是多个地物的混合体,一个像元里记录了多类不同性质的地面目标,图像的光谱特征也是多个地物光谱特征的混合反映。混合像元的存在,产生了同物异谱、同谱异物现象,从而给地物的分类造成许多困难。因此,为了提高干扰信息提取的精度,就必须找到一种有效的进行混合像元分解的方法。线性光谱混合模型(Linear Spectrum Mixture Model,LSMM)是目前应用较多的混合像元分解方法[1~2],但LSMM因为要对整幅图像的每个像元都进行处理,所以处理速度较慢。
蚁群算法是由意大利学者Dorigo等人,通过模拟自然界蚂蚁寻径的行为提出的一种全新仿生进化算法[3~4],是具有离散性、并行性、鲁棒性、正反馈性等特点的一种随机搜索方法。由于其概念简明、实现方便,迅速得到认可,并在优化问题求解、电力系统、计算机、冶金自动化等领域都有成功的应用[5]。但目前很少有人将蚁群算法引入遥感地质领域,由于遥感地质领域的信息提取也可以看作是一个组合优化问题,蚁群算法的全局性、离散性和基于概率选择路径等特点对于遥感图像非常适用。
本文为了剔除植被干扰信息,综合考虑LSMM处理速度慢而蚁群算法识别目标速度快的特点,结合蚁群算法和线性光谱混合模型,建立基于蚁群搜索的光谱分解模型,剔除植被干扰信息,重构不含有植被信息的新的多波段图像,以期为后续基础地质工作提供基础影像。
1. 研究方法
1.1 蚁群算法基本原理
对于采用整数编码求解组合优化问题的蚁群系统,设求解问题因素有N个,蚁群中共有M只蚂蚁,τij(t)表示在t时刻i和j组合之间信息素的数量。蚂蚁m在运动过程中根据各个路径上信息素的数量决定下一步的路径。用pijm(t)表示在t时刻蚂蚁m由城市i转移到城市j的概率,则:
pmij(t)={ταij(t)⋅ηβij(t)∑r∈Tmταir(t)⋅ηβir(t)j∈Tm0otherwise (1) 其中:Tm表示蚂蚁m下尝试组合情况的集合,该集合随蚂蚁m的行进过程而动态改变。
信息量τij(t)随时间的推移会逐步衰减,用1-ρ表示它的衰减程度。经过n个时刻,要根据下式对各路径上的信息量作更新:
τij(t+1)=ρ⋅τij(t)+Δτij (2) Δτij=M∑s=1Δτkij (3) Δτijk表示蚂蚁m在本次循环中在组合i和j之间留下的信息量,其计算方法根据蚁群系统的计算模型而定。
1.2 蚂蚁移动规则
正如Chialvo[6]和Millonas[7]所描述的,单只蚂蚁的状态可以用位置参数r和方向参数θ来表示。蚂蚁的状态转移规则可用一个加权函数表示:
w(σ)=[1+σ1+δσ]β (4) 这个方程描述了移动到信息素浓度为σ(r)的像素r处的相对概率。参数β表示一种随机度。β大,则w(σ)值较大,蚂蚁以较大的权系数跟随外激素浓度大的路径;反之,对蚂蚁路径选择影响不大。1/δ表示了蚂蚁感知外激素的能力。
蚂蚁从像元k运动到像元i的归一化转移概率定义为:
pik=w(σi)w(Δi)∑j/kw(σj)w(Δj) (5) w(Δ)是一加权因子,Δj表示蚂蚁在t-1时刻运动时的方向改变量,它的取值为8个离散的w值,一般蚂蚁的来向w定为1/20,原路径方向定为1,邻近的8个像元顺时针依次为1/20,1/12,1/4,1/2,1,1/2,1/4,1/12。
蚂蚁爬行时,先计算8邻域每个像元的转移概率pik,然后根据轮赌算法选择移动到下一个像元。为了有利于全局搜索最优解,避免陷入局部最优解,设定了蚂蚁的兴奋阈值,当某支蚂蚁达到该值时,采用随机算法从8邻域中选择一个像元进行移动。
1.3 线性光谱混合模型
线性光谱混合模型(Linear Spectrum Mixture Model,LSMM)主要的目的是绘制特定像元内地面物质的相对丰度图。因此,它必须要满足几个假设:
① 在一个像元里端元的混合必须是线性的。这种线性关系下光子在从太阳到传感器的路径上仅仅与一种端元成分交互,这要求端元组分分布在足够大的面积上。
② 在图像上所有的地物类型必须要有足够多的对照光谱,以便对它们进行分离。
③ 在图像里存在的端元组分必须在最少一个像元里保持纯状态,这个像元必须能正确识别且可空间定位。
LSMM假定传感器测量的光谱是像元内所有成分光谱的线性组合[1]。LSMM模型描述如下:
Ri=n∑k=1fkRik+εi (6) 在这里i为光谱波段的数目;k是端元的数目;Ri是像元波段i的光谱反射率;fk是像元内端元k的面积比;Rik是波段i端元k的光谱反射率;εi是波段i的模拟误差。
fk受下式约束:
n∑k=1fk=1且0≤fk≤1 (7) 残差RMS用下式计算:
RMS=√(m∑i=1ε2i)/m (8) RMS残差要对所有图像像元都计算。RMS越大,模型的匹配就越差。所以残差图可以用来估计选择的端元是否适当,选择的端元数目是否足够。
因为在LSMM模型中,端元反射率和像元反射率都是已知的,只有端元在像元里的面积比未知,只要波段的数目大于端元数目,就可以求出端元的面积比,从而得到各分量图像和残差图像,进而得到地物的分类图像。
LSMM的优点在于它相对简单,并且是一种有物理意义的丰度测量方法;理论上有较好的科学性,对于解决像元内的混合现象有一定的效果。
2. 基于蚁群算法的光谱分解算法
本文中基于蚁群算法的光谱分解算法采用的是混合像元分解蚁群算法,输入经过反射率转换预处理后的影像,输出不含有植被信息的新的多波段图像。
【算法开始】
① 初始化
1) 蚁群及相关参数;
2) 信息素矩阵;
3) 访问标志矩阵P;
4) 残差矩阵RMS;
5) 丰度矩阵F;
② 迭代,直到达到停止条件
1) 评价每只蚂蚁,对每只蚂蚁做
如果没有访问过当前像元则完成下述工作
(1) 建立端元系数矩阵x;(矩阵x主要存储像元内各子端元的面积)
(2) 建立当前像元波段矩阵y;(矩阵y主要存储当前像元各波段的光谱反射率值)
(3) 调用最小二乘法计算各端元的比例系数矩阵a;(重新计算各子端元的面积)
(4) 对矩阵a进行负数漂移处理;
(5) 对矩阵a进行归一化处理;
(6) 计算该点残差存放到RMS矩阵;
(7) 记录丰度信息到F;
(8) 根据矩阵a调整该像元的各波段值;
(9) 设置当前像元访问标志;
(10) 如果剔除端元所占比例超过某阈值,调整(增加)信息素;
2) 每只蚂蚁爬行,做
如果 蚁群迭代到某兴奋阈值 则
蚂蚁采用随机算法从8邻域中选择一个像元进行移动
否则
(1) 根据公式(4)计算8邻域每个像元的w(σ)
(2) 根据进入方向分配8邻域每个像元的权值
(3) 根据公式(5)计算8邻域每个像元的转移概率pik
(4) 按照轮赌算法选择移动到下一个像元
③ 按照BSQ格式保存图像数据;
④ 按照BSQ格式保存残差矩阵;
⑤ 按照BSQ格式保存剔除专题的丰度矩阵;
【算法结束】
3. 试验及结果分析
选取青海黄南州吉地地区为试验区,该区为植被高度覆盖区。试验区遥感数据为Landsa7号卫星的ETM,景号为132-36,2000年7月18日获取的影像,太阳高度角为64.2°。试验区图像大小为408×402像素,经过反射率转换预处理后的影像作为输入图像。基于ENVI3.6及VC++6.0开发的蚁群光谱分解软件进行试验。
通过求取植被指数,可以发现该区植被覆盖度为85.7157%(见图 1,绿色部分为植被),高植被覆盖率严重影响了矿化弱信息的提取。选用ETM1,2,3,4,5,7等6个波段进行主成分变换,其中PC1和PC2占95.3%。通过主成分分析,选用PC1和PC2制作二维散点图(见图 2)。通过散点图分析,确定图像上主要有植被及3种未知类别共4个端元。运行蚁群光谱分解软件,最后得到剔除植被后影像(见图 3)及残差图(见图 4)。
从残差图上可以看出,最大残差为0.103,绝大部分残差在0.05以下。说明光谱分解的效果较好。对比图 2和图 3可以发现,进行植被剔除后,原来绿色部分的植被覆盖区被还原为该区本来的面目。
4. 结论
基于蚁群算法的光谱分解方法剔除植被信息,首次将蚁群这种全新的算法引入到遥感地质领域,综合考虑传统的光谱分解植被剔除方法处理速度慢和蚁群算法识别目标速度快的特点,通过残差图分析以及原图与剔除植被后影像对比分析,初步验证了基于蚁群算法的光谱分解方法来剔除植被信息的可行性。
由于蚁群算法的理论体系还有待进一步完善,所以用这种全新的算法来处理遥感地质领域的一些复杂问题,还有许多工作要做。如蚁群的初始参数怎样优化,才能取得理想的效果,有待于进一步深入研究。
-
图 3 瘤状灰岩成因模式
(A)差异压实—压溶成因:1—半固结或弱固结状态下的薄层灰质层与薄层泥质层互层沉积物;2—灰质层经压实作用被拉长变薄,发生细颈化;3—开阔台地相或滩后澙湖相接受进一步压实,细颈化的灰质层被压断,边缘局部溶解、圆化,形成扁豆状、疙瘩状瘤体,部分瘤体受差异压实—压溶作用而发生短距离剪切位错,形成断续状、杂乱状瘤状灰岩,台地前缘斜坡相则发生滑塌作用,形成串珠状瘤状灰岩
(B)暴露成因:4—半固结的含泥碳酸盐沉积物;5—地壳隆升,海水退出,沉积物间歇性暴露地表形成网脉状、不规则状裂隙;6—黏土物质经水流作用充填于裂隙中,胶结先期形成的灰岩碎块形成瘤状构造,后期经受风化溶蚀而成蜂巢状Figure 3. Genetic models of nodular limestone
-
[1] 董兆雄, 朱晓惠, 侯方浩, 等.一种特殊瘤状灰岩的成因研究[J].沉积学报, 2002, 20(1):20~24. http://www.cnki.com.cn/Article/CJFDTOTAL-CJXB200201003.htmDONG Zhao-xiong, ZHU Xiao-hui, HOU Fang-hao, et al. Study on origin of special nodular limestone[J]. Acta Sedimentologica Sinica, 2002, 20(1):20~24. http://www.cnki.com.cn/Article/CJFDTOTAL-CJXB200201003.htm [2] 夏丹, 汪凯明, 罗顺社, 等.燕山地区高于庄组张家峪亚组瘤状灰岩成因研究[J].石油地质与工程, 2009, 23(1):4~7. http://www.cnki.com.cn/Article/CJFDTOTAL-SYHN200901005.htmXIA Dan, WANG Kai-ming, LUO Shun-she, et al. Study on genesis of knollenkalk at Zhangjiayu subgroup of Gaoyuzhuang Group at Yanshan area[J]. Petroleum Geology and Engineering, 2009, 23(1):4~7. http://www.cnki.com.cn/Article/CJFDTOTAL-SYHN200901005.htm [3] 张霞, 林春明, 凌洪飞, 等.浙西地区奥陶系砚瓦山组瘤状灰岩及其成因探讨[J].古地理学报, 2009, 11(5):481~490. http://www.cnki.com.cn/Article/CJFDTOTAL-GDLX200905004.htmZHANG Xia, LIN Chun-ming, LING Hong-fei, et al. Nodular limestone and its genesis from the Ordovician Yanwashan Formation in western Zhejiang Province[J]. Journal of Palaeogeography, 2009, 11(5):481~490. http://www.cnki.com.cn/Article/CJFDTOTAL-GDLX200905004.htm [4] 云南省地质矿产局.云南省区域地质志[M].北京:地质出版社, 1990.Bureau of Geology and Mineral Exploration and Development of Yunnan Province. Regional geology of Yunnan Province[M]. Beijing:Geological Publishing House, 1990. [5] 时言, 龚大明.滇西保山施甸地区泥盆系研究新知[J].成都地质学院学报, 1992, 19(3):21~34. http://www.cnki.com.cn/Article/CJFDTOTAL-CDLG199203002.htmSHI Yan, GONG Da-ming. Advances of Devonian research in Baoshan-Shidian area, western Yunnan[J]. Journal of Chengdu College of Geology, 1992, 19(3):21~34. http://www.cnki.com.cn/Article/CJFDTOTAL-CDLG199203002.htm [6] LIU Ben-pei, FENG Qing-lai, FANG Nian-qiao. Tectonic evolution of the Palaeo-Tethys in Changning-Menglian Belt and adjacent Regions, Western Yunnan[J]. Journal of China University of Geosciences, 1991, 2(1):18~28. [7] Mullins H T, Wilber R J, Bardmen M R. Nodular carbonate sediment on Bahamian slopes[J]. Sediment Petrol, 1980, 50(1):117~138. http://archives.datapages.com/data/sepm/journals/v47-50/data/050/050001/0117.htm?q=%2BauthorStrip%3Aboardman+authorStrip%3Awillhite [8] 孙庆峰.结核状碳酸盐岩研究综述[J].上海地质, 2005, (4):62~67. http://www.cnki.com.cn/Article/CJFDTOTAL-SHAD200504020.htmSUN Qing-feng. On nodular limestones[J]. Shanghai Geology, 2005, (4):62~67. http://www.cnki.com.cn/Article/CJFDTOTAL-SHAD200504020.htm [9] 韩树芬.安微沿长江地区下三叠统瘤状灰岩成因研究[J].地质科学, 1983, (3):232~238. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKX198303003.htmHAN Shu-fen.[J]. Study on the genesis of the Lower Triassic nodular limestone in the area along Yangtze River in Anhui Province[J]. Scientia Geologica Sinica, 1983, (3):232~238. http://www.cnki.com.cn/Article/CJFDTOTAL-DZKX198303003.htm [10] Mullins H T, Neumann A C, Wilber R J, et al. Nodular carbonate sediment on Bahamian slopes:Possible precursors to nodular limestones[J]. Journal of Sedimentary Petrology, 1980, 50(1):117~131. http://archives.datapages.com/data/sepm/journals/v47-50/data/050/050001/0117.htm?q=%2BauthorStrip%3Aboardman+authorStrip%3Awillhite [11] Rais P, Louis Schmid B, Bernasconi S M, et al. Palaeocean ographic and palaeoclimatic reorganization around the Middle Late Jurassic transition[J]. Palaeogeography, Palaeoclimat ology, Palaeoecology, 2007, 251:527~546. doi: 10.1016/j.palaeo.2007.05.008 [12] 蓝光志, 张廷山, 高卫东.川西北地区早志留世瘤状灰岩的类型、成因及意义[J].西南石油学院学报, 1994, 16(3):1~5. http://www.cnki.com.cn/Article/CJFDTOTAL-XNSY403.000.htmLAN Guang-zhi, ZHANG Ting-shan, GAO Wei-dong. The types, genesis and significance of the Early Silurian nodular limestone in the area of northwestern Sichuan[J]. Journal of Southwest Petroleum Institute, 1994, 16(3):1~5. http://www.cnki.com.cn/Article/CJFDTOTAL-XNSY403.000.htm [13] 朱洪发, 王恕一.苏南、皖南三叠纪瘤状灰岩、蠕虫状灰岩的成因[J].石油实验地质, 1992, 14(4):454~460. doi: 10.11781/sysydz199204454ZHU Hong-fa, WANG Shu-yi. The origins of the Triassic nodular and vermicular limestones in south Jiangsu-south Anhui provinces[J]. Experimental Petroleum Geology, 1992, 14(4):454~460. doi: 10.11781/sysydz199204454 [14] 郭福生, 梁鼎新.浙江江山砚瓦山组瘤状灰岩的成因[J].矿物岩石, 1993, 13(3):74~80. http://www.cnki.com.cn/Article/CJFDTOTAL-KWYS199303011.htmGUO Fu-sheng, LIANG Ding-xin. Genesis of the nodular limestone in Washan Formation in Jiangshanyan, Zhejiang Province[J]. Journal of Mineralogy and Petrology, 1993, 13(3):74~80. http://www.cnki.com.cn/Article/CJFDTOTAL-KWYS199303011.htm [15] 冯增昭, 王英华, 李尚武, 等.下扬子地区中下三叠统青龙群岩相古地理研究[M].昆明:云南科技出版社, 1988:1~69.FENG Zeng-zhao, WANG Ying-hua, LI Shang-wu, et al. Research on lithofacies palaeogeography of the Lower-Middle Triassic Qinglong Group in the lower Yangtze region[M]. Kunming:Yunnan Science and Technology Press, 1988:1~69. [16] 郭福生.下扬子地区三叠系下统瘤状灰岩成因研究[J].华东地质学院学报, 1989, 12(4):17~22.GUO Fu-sheng. Study on the genesis of the Lower Triassic nodular limestone in the lower Yangtze region[J]. Journal of East China College of Geology, 1989, 12(4):17~22. [17] 金振奎, 冯增昭.云贵地区二叠系瘤石灰岩的成因[J].岩石矿物学杂志, 1994, 13(2):133~137. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKW402.004.htmJIN Zhen-kui, FENG Zeng-zhao. The origin of the Permian nodular limestones in Yunnan-Guizhou region[J]. Acta Petrologica et Mineralogica, 1994, 13(2):133~137. http://www.cnki.com.cn/Article/CJFDTOTAL-YSKW402.004.htm [18] 黄勇, 邓贵标, 彭成龙, 等.滇西保山南部早中奥陶世沉积缺失的发现及意义[J].贵州地质, 2009, 26(1):1~6. http://www.cnki.com.cn/Article/CJFDTOTAL-GZDZ200901003.htmHUANG Yong, DENG Gui-biao, Peng Cheng-long, et al. The discovery and significance of absence in Early-Middle Ordovician in southern Baoshan, western Yunnan[J]. Guizhou Geology, 2009, 26(1):1~6. http://www.cnki.com.cn/Article/CJFDTOTAL-GZDZ200901003.htm [19] 陆松年.初论"泛华夏造山作用"与加里东和泛非造山作用的对比[J].地质通报, 2004, 23(9):952~958. http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD2004Z2018.htmLU Song-nian. Comparison of the Pan-Cathaysian orogeny with the Caledonian and Pan-African orogenies[J]. Geological Bulletin of China, 2004, 23(9):952~958. http://www.cnki.com.cn/Article/CJFDTOTAL-ZQYD2004Z2018.htm [20] Wopfner H. Gondwana origin of the Baoshan and Tengchong terrenes of west Yunnan[C]//Hall R, Blundel D. Tectonic evolution of Suntheast Asia. London:Geological Society, 1996:539~547. [21] Yin A, Harrison T M. Geologic evolution of the Himalayan-Tibetan orogeny[J]. Annual Review of Earth and Planetary Sciences, 2000, 28:211~280. doi: 10.1146/annurev.earth.28.1.211 [22] 钟大赉.滇川西部古特提斯造山带[M].北京:科学出版社, 1998.ZHONG Da-ben. Palaeo-Tethysides in the west of Yunnan and Sichuan provinces[M]. Beijing:Science Press, 1998. [23] 陈福坤, 李秋立, 王秀丽, 等.云南特提斯带保山-腾冲地块早古生代岩浆岩[J].地球学报, 2005, 26(增刊):93. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDJ200509001031.htmCHEN Fu-kun, LI Qiu-li, WANG Xiu-li, et al. Early Paleozoic magmatism in Baoshan-Tengchong block of the Tethyan belt, Yunnan Province[J]. Acta Geoscientica Sinica, 2005, 26(Supp.):93. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGDJ200509001031.htm -