DESIGN AND IMPLEMENTATION OF THE OPERATION MANAGEMENT AND SECURITY SAFEGUARD SYSTEM FOR FIELD GEOLOGICAL SURVEY
-
摘要: 基于北斗导航定位、遥感卫星以及基础地理空间信息承载与综合显示平台而建立的野外地质调查作业管理与安全保障系统,采用组件化、平台化的系统体系结构,具备强大且灵活的可扩展性和可集成性。该系统可实现野外地质调查区域基础遥感影像、GIS和DEM数据等的分层显示、北斗通讯和导航信息的实时显示、通讯历史记录和定位信息的查询、野外地质调查人员的通讯与定位等;同时可实现野外地质调查人员与各级管理部门的互联互通,管理部门对野外作业人员的作业态势、作业进度可进行综合查询,对野外作业人员的外勤安全、遇险救援提供决策支撑。该系统在部分野外地质调查示范单位运行效果良好。Abstract: Based on the Beidou satellite navigation service, remote sensing satellites, as well as the platform of basic geospatial information bearing with integrated displaying, the system of operation management and security safeguard for field geological survey is established. That uses a structure of components and platformwith a powerful and flexible scalability and integration. The hierarchical displays of remote sensing images, GIS and DEM data for field geological survey area, also the query of communication and position, are implemented. At the same time, the interoperability of the personnel about field geological survey and the management departments, the integrated query of the work situation and the progress of the operation, the decision support provided for security and relief, are implemented. Currently, the system has deployed several demonstrational units of field geological survey and gained recognition.
-
0. 前言
断裂作用可导致两种不同形式应变:连续应变和非连续应变。由断裂位移来估算的应变是非连续应变[1~3] (图 1)。由于应变概念只能应用于连续变形, 所以Jamison (1989) [1]引进了断裂应变这一词, 用来描述由一系列断裂位移导致的非连续变形。实际上, 断裂应变是非连续应变还是连续变形取决于所研究的范围和断裂位移的相对大小[1, 4]。如果所研究地区的尺度与断裂位移相比大很多, 断裂应变就可以视为连续变形。到目前为止, 断裂应变在构造地质界尤其在中国构造地质界没有得到足够重视。这主要是由于计算断裂应变必须获得每一条断裂位移大小和方向。这在野外工作中是比较困难的。因为往往缺乏合适的和足够的被错动标志体。但是对于露头条件很好的地区, 尤其是现代断裂活动区, 比较而言容易观察到断裂运动标志体。这些标志物很少经过后期的破坏和沉积物覆盖。
对于断裂应变, 有两种计算方法。一种是计算单个断裂引起的应变量, 也就是垂直断裂方向的剖面上求水平方向的应变。另一种通过研究区内所有断裂数据求断裂应变椭球体, 这样就可以知道主应变的方向和大小。本文只介绍第一种方法。
1. 断裂作用过程的连续应变
变形实验表明, 在断裂发生前, 有一定的弹性应变积累。在静岩应力条件下(σv=σH=σh=ρgz), 其表达式为:
(1) 这里εh为静岩应力条件下的水平应变, E是扬氏模量, ν为泊松比, ρgz就是垂直应力。如果岩体变形已经达到断裂阶段, 这种弹性应变已不可恢复。但是它与断裂引起的塑性应变相比通常可以忽略不计。
断裂作用过程中, 并不总是断块的刚性运动, 而是可以引起地块一定的塑性变形。这种塑性应变与断裂最大位移量和断裂长度有关。Schultz和Fossen (2002) [5]给出了一个计算公式:
(2) 这里D为断裂最大位移, L为断裂长度, δ为断裂倾角, 而σ为断裂面上有效应力的校正值, 其具体的计算表达式见Schultz和Fossen (2002) [5]。
据Schultz和Fossen (2002), 对于正断层, 计算出的伸展应变为2 %~ 3 %; 对于逆断层, 计算出的压缩应变为4 %~ 5 %。如图 2所示, 这种应变与断裂最大位移(D)与断裂长度(L)的比值成正相关。即比值越大, 应变量越大。
2. 断裂作用过程的非连续应变
2.1 断裂类型
断裂应变与三个因素有关:断裂几何形态、断裂的旋转性、断裂的位移大小。这三种因素的不同组合, 给出的计算方法不一样。断裂有很多的分类方法, 在这里不一一列举。我们只介绍Wernicke (1982) [6]的分类方法。他是根据断裂形态和断裂旋转性来进行分类的(表 1)。按照这个分类, 只有平面状断裂不会引起地层旋转情况。铲状断裂可以引起地层旋转但本身不一定旋转。如果同时有地层和断裂旋转, 可以是平面状断裂也可以铲状断裂引起。
2.2 断裂旋转机制
到目前为止, 已提出了三种断裂旋转机制。最早认为, 断裂旋转是刚性的, 断块内部没有任何变形[7]。这种机制存在的问题是断裂旋转产生的空隙, 这些空隙没有用其他的机制加以完满解释。而且这种空隙与旋转幅度成正比(图 3)。由于上述原因产生了垂直简单剪切和斜向简单剪切模型。垂直剪切机制提出断裂上盘由于断裂作用, 发生垂直方向的简单剪切作用。越靠近断裂, 剪切作用越强, 同时认为这种剪切机制也是地层旋转的原因[8~10]。然而有作者认为, 在铲状断裂的上盘, 剪切方向不一定是垂直的, 它可以在任何方向进行, 这取决于断裂的几何形态和远场应力作用状态。这就是所谓的斜向简单剪切机制[11] (图 4)。
图 3 断裂刚性旋转示意图(a)表示断裂还没有位移时的状态; (b)表示断裂发生位移同时发生旋转, 断裂倾角变小。在问号处留下的空隙没有得到很好的解释。断裂旋转的角度等于地层的倾角, 也就是θ =δ0-δFigure 3. Diagram of the rigid -body mechanism(a)The initial state in which the faults are with no displacement; (b)The fault dips decrease with the rotation of faults.The spaces with interrogation marks are not well explained. The rotated angle of the bed is equal to that of the faults, that is to say, θ =δ0-δ2.3 断裂应变的计算
当断裂和地层都不发生旋转时(δ =δ0), 断裂作用前后地层的长度不变(图 5)。断裂引起的伸长CB (对正断层而言)或缩短(对逆断层而言)等于断裂水平断距或平错。因为CB等于dcos (δ), 所以应变量应为:
(3) 从(3)式可以看到, 应变的大小与断裂倾角成反相关关系; 也就是说, 对相同的总位移, 断裂倾角越大, 应变越小。同时也能看到, 应变与位移大小成正相关关系; 也就是说, 相同的断裂倾角, 总位移越大, 应变越大。
对于刚性旋转机制断块(图 6), 由于在断块内部没有变形, 断裂旋转的角度等于地层的倾角, 即θ=δ0-δ。断裂的平错等于L0cot(δ)sin(θ), 断裂作用后的长度为D′C′=D′B+BC′=L0[cot(δ)sin(θ)+cos(θ)], 所以应变大小为:
(4) 由(4)式可以看到, 应变大小与断裂现在的倾角δ和地层的倾角θ有着很复杂的关系, 而不是我们想象的那么简单。特别地, 我们由断裂现在的倾角δ可以计算断裂形成时的倾角δ0, 其表达式是:
(5) 对于垂直剪切机制, 断块内部发生了垂直方向简单剪切。这样一来, 断块旋转以后的地层长度比原来的要长, 但在水平方向的投影与原始的长度一致。如图 7所示, 断裂平错为h =Lbcot(δ)sin(θ), 变形前的长为L0=Lbcos(θ), 所以应变大小为:
(6) 由(6)式可以推断, 应变大小与断裂现在的倾角成反相关关系, 而与地层的倾角成正相关关系。有兴趣的读者可以计算一下, 对于相同的断裂旋转角, 垂直剪切引起的水平应变要大于刚性旋转引起的应变量[10]。
同样地, 由断裂现在的倾角δ可以计算断裂形成时的倾角δ0, 其表达式是:
(7) 可以看到, 对于垂直剪切机制, 变形前后的断裂倾角不是简单的代数关系。比较式(7)和式(5), 对于相同的地层倾角和断裂现在的倾角, 由垂直剪切机制得到的断裂原始倾角相对较小。
斜向剪切机制与垂直剪切机制有些类似(图 8)。但是二者相比, 斜向剪切机制引起的拉伸量比垂直剪切机制引起的拉伸量大, 用等式表示为:
(8) 其应变量为:
(9) 从这个方程可以推断, 当α等于零时, 式(9)等同于式(6)。式(9)可以进一步变为:
(10) 该式表明, 只要tan(α)tan(θ)小于1, 斜向剪切机制引起的拉伸量比垂直剪切机制引起的拉伸量大。据White等(1986)[11], α的大小一般小于45度, θ的大小也小于45度, 因此, tan(α)tan(θ)小于1。
2.4 研究实例
我们考察来自于墨西哥中央桌子山San Miguelito地区渐新世火山岩区的断裂。该区发育有一系列的多米诺式的正断层(图 9)。通过研究认为, 这些断裂经历了垂直剪切作用。
断裂参数都是通过野外质地测量所得(表 2)。测量标志体为Cantera未熔结凝灰岩。地层倾角为断块内的平均值。也就是通过测量一系列的倾角值, 然后求得平均值。在一个断块内, 可以得到5~10个地层数据。我们看到, 每一断块的应变量各不相同。各断块的地层的拉伸应变的最大值15. 5 % (断块7), 而水平应变的最大值达31. 2 % (断块8)。对于整个剖面, 求得的应变量为20 %。
3. 结论
(1) 断裂作用可导致连续应变即塑性变形与非连续应变。他们分别有不同的计算公式。
(2) 计算断裂的非连续应变应考虑断裂几何形态、断裂的旋转与否、断裂的位移大小等三个因素。
(3) 刚性旋转时, 没有断块内变形。它引起的水平非连续应变最小。垂直剪切作用使断块内地层变形, 即有地层的连续性应变。在相同条件下, 它引起的非连续应变量比刚性旋转机制引起的非连续应变量大。斜向剪切也使断块内地层变形, 也有地层的连续性应变。在相同条件下, 它引起的非连续应变比垂直剪切机制引起的非连续应变应变大。
-
表 1 2011年度野外地质调查作业管理与安全保障系统部署情况统计
Table 1. Deployment statistics of operation management and security safeguard system for field geological survey in 2011
-
[1] 毕思文.全球变化与地球系统科学统一研究的最佳天然实验室——青藏高原[J].系统工程理论与实践, 1997, (5):72~77. http://www.cnki.com.cn/Article/CJFDTOTAL-XTLL705.011.htmBI Si-wen. A best laboratory of the universal research for the earth's global change and earth system science:The Qinghai-Tibet Plateau[J]. Systems Engineering Theory & Practice, 1997, (5):72~77. http://www.cnki.com.cn/Article/CJFDTOTAL-XTLL705.011.htm [2] 吴中海, 赵希涛, 吴珍汉, 等.西藏纳木错地区约120 ka BP以来的古植被、古气候与湖面变化[J].地质学报, 2004, 78(2):242~252. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200402013.htmWU Zhong-hai, ZHAO Xi-tao, WU Zhen-han, et al. Palaeovegetation, palaeeclimate and lake-level change since 120 ka BP in Nam Co, central Xizang[J]. Acta Geologica Sinica, 2004, 78(2):242~252. http://www.cnki.com.cn/Article/CJFDTOTAL-DZXE200402013.htm [3] 吴绍洪, 尹云鹤, 郑度, 等.青藏高原近30年气候变化趋势[J].地理学报, 2005, 60(1):3~11. doi: 10.11821/xb200501001WU Shao-hong, YIN Yun-he, ZHENG Du, et al. Climate changes in the Tibetan Plateau during the last three decades[J]. Acta Geographica Sinica, 2005, 60(1):3~11. doi: 10.11821/xb200501001 [4] 倪晋宇, 胡道功, 周春景.东昆仑造山带纳赤台群形成的大地构造环境探讨[J].地质力学学报, 2010, 16(1):11~20. http://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?flag=1&file_no=20100102&journal_id=dzlxxbNI Jin-yu, HU Dao-gong, ZHOU Chun-jing. Discussion on tectonic environment of Naij Tal Group, East Kunlun orogenic belt[J]. Journal of Geomechanics, 2010, 16(1):11~20. http://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?flag=1&file_no=20100102&journal_id=dzlxxb [5] 邓吉秋, 鲍光淑.基于GIS的湖南省国土资源遥感综合调查信息系统分析与设计[J].中南工业大学学报:自然科学版, 2003, 34(2):192~195. http://www.cnki.com.cn/Article/CJFDTOTAL-JSGG200405064.htmDENG Ji-qiu, BAO Guang-shu. Analysis and design of information system of remote sensing integrated survey about land resources in Hunan Province based on GIS[J]. Jourmal of Central South University of Technology:Natural Science, 2003, 34(2):192~195. http://www.cnki.com.cn/Article/CJFDTOTAL-JSGG200405064.htm [6] Fortier P J. 数据库技术大全[M]. 林瑶, 范建华, 赵刚, 译. 北京: 电子工业出版社, 1999.Fortier P J.Database Systems Handbook [M]. LIN Yao, FAN Jian-hua, ZHAO Gang, Translated. Beijing:Publishing House of Electronics Industry, 1999. [7] Walters R E, Coles M, Ferracchiati F, 等. 深入SQL Server 2008[M]. 任斌, 刘芳芳, 译. 北京: 人民邮电出版社, 2011.Walters R E, Coles M, Ferracchiati F, et al. Accelerated SQL Server 2008[M]. JIA Hong-feng, Translated. Beijing:Tsinghua University Press, 2011. [8] Sack J. SQL Server 2008实战[M]. 金迎春, 译. 北京: 人民邮电出版社, 2010.Sack J. SQL Server 2008transact-SQL recipes[M]. JIN Ying-chun, Translated. Beijing:Posts & Telecom Press, 2010. [9] Lobel L, Brust A J, Forte S. 精通SQL Server 2008程序设计[M]. 贾洪峰, 译. 北京: 清华大学出版社, 2010.Lobel L, Brust A J, Forte S. Mastering SQL Server 2008 program design[M]. JIA Hong-feng, Translated. Beijing:Tsinghua University Press, 2010. [10] Skeet J. 深入理解C#(第2版)[M]. 周靖, 朱永光, 姚琪琳, 译. 北京: 人民邮电出版社, 2012: 18~19.Skeet J. C# in depth (Second Edition) [M]. ZHOU Jing, ZHU Yong-guang, YAO qi-lin, Translated. Beijing:Posts & Telecom Press, 2012:18~19. -