朱廣軼,鄭仰發(fā),張華興

(1.沈陽大學建筑工程學院,遼寧沈陽 110044;2.天地科技股份有限公司開采設計事業(yè)部,北京 100013)

礦區(qū)地質(zhì)動力區(qū)劃趨勢面分析的可視化軟件系統(tǒng)

朱廣軼1,鄭仰發(fā)1,張華興2

(1.沈陽大學建筑工程學院,遼寧沈陽 110044;2.天地科技股份有限公司開采設計事業(yè)部,北京 100013)

地質(zhì)動力區(qū)劃是礦區(qū)防治多種地質(zhì)災害的新方法?；谮厔菝娣治隼碚摵?3次樣條插值函數(shù)及VB調(diào)用 CAD二次開發(fā)成圖技術(shù),構(gòu)建了地質(zhì)趨勢面分析的可視化軟件系統(tǒng)。綜合礦區(qū)地質(zhì)地形圖、震中分布圖、地下勘察及地應力測量結(jié)果,軟件系統(tǒng)可有效區(qū)劃出活動斷裂,從而為煤礦開采設計與瓦斯動力事故等災害預防提供科學理論依據(jù)。

地質(zhì)動力區(qū)劃;地質(zhì)趨勢面;可視化系統(tǒng);活動斷裂

地質(zhì)動力區(qū)劃首創(chuàng)于俄羅斯地質(zhì)力學與礦山測量研究院[1],是近十幾年發(fā)展起來的集地質(zhì)動力學、測量學、板塊學說、地球物理學、巖體力學等多學科的研究成果,旨在探討研究區(qū)域內(nèi)各種動力過程,以及參與該過程的地球能量場。通過對地形地貌、構(gòu)造物理、河流水系的考察分析,同地質(zhì)構(gòu)造、地震分布、地球物理、航片、衛(wèi)星照片等圖紙資料進行對比,確定研究區(qū)域內(nèi)各級 (活動)斷塊圖及斷塊間的相互作用。利用地應力測量、數(shù)值分析、GIS等技術(shù)手段,以斷塊圖為模型,進行原巖應力計算,確定應力場的特性,為工程活動的防災減災提供科學理論依據(jù)。

在采礦領(lǐng)域,地質(zhì)動力區(qū)劃方法對礦區(qū)地質(zhì)災害預測分析具有時間短、工程評價費用低、效果顯著等特點,并且已在國內(nèi)多個礦區(qū)得到了卓有成效的應用[2-3]。地質(zhì)趨勢面分析法是地質(zhì)動力區(qū)劃進行構(gòu)造斷塊劃分最行之有效的方法,屬于常用的數(shù)學地質(zhì)方法[4],適用于對平面和空間中分布的各種觀測資料進行統(tǒng)計與擬合分析。

目前,我國還沒有自主產(chǎn)權(quán)的地質(zhì)趨勢面分析專用軟件。本文構(gòu)建了地質(zhì)趨勢面分析的可視化軟件系統(tǒng),這不僅彌補了國內(nèi)此領(lǐng)域的空白,而且開創(chuàng)了地質(zhì)動力區(qū)劃可視化運用的新思路,對逐步完善地質(zhì)動力區(qū)劃方法具有一定意義。

1 趨勢面分析理論

1.1 理論模型

在理論上,原始觀測面一般可以分解為 3部分:區(qū)域趨勢,局部異常,隨機干擾[4]。所謂區(qū)域趨勢是指遍及全區(qū)的,較大規(guī)模地質(zhì)作用的反映;局部異常是由規(guī)模比研究區(qū)小的地質(zhì)作用所產(chǎn)生的,其影響在區(qū)域內(nèi)并非處處可見的隨機分量;隨機干擾,一般認為是由抽樣誤差和觀測誤差組成,不包括系統(tǒng)誤差。統(tǒng)計學理論指出,隨機成分的分離是有條件的,要求有重復抽樣的觀測數(shù)據(jù),這在地質(zhì)工作中往往難以做到[4]。因此,在實際工作中,只要求分離其中的兩種成分,從而理論模型在實際應用時就可表述為:

上述模型中,局部異常分量必然包含隨機干擾分量,只不過局部異常分量處于主導地位,而隨機分量所占的比重很小。分析思路是分別將分離開的分量作圖,然后作地質(zhì)解釋,一般只作等值線圖,有時也作傾向剖面投影圖。

1.2 數(shù)學模型

理論模型 (1)在數(shù)學上用函數(shù)式表述:

式中,Hi為第 i點的實測值,ui為第 i點處東西方向的坐標值,vi為第 i點處南北向的坐標值,ei為第 i點處的擬合剩余值,p為 u坐標值的最高冪次, q為 v坐標值的最高冪次,ajk是 u為 j次冪,v為 k次冪的多項式系數(shù),n為觀測點數(shù)。

要建立趨勢面擬合函數(shù),重要的是求出待定系數(shù) ajk,一般采用最小二乘準則。根據(jù)微積分原理,可以得到多元線性方程組的矩陣方程:

其中,m為擬合多項式函數(shù)的最高次數(shù),r與m間的關(guān)系式:

將 n組非規(guī)則分布觀測點的三維坐標代入(3)式,解算出系數(shù) a0,a1,…,ar后回代入式(2),就可求得趨勢值與剩余值。

1.3 斷裂構(gòu)造的地質(zhì)特征

斷裂構(gòu)造在觀測面上往往表現(xiàn)為一條直線或曲線,在殘差圖上可表觀為細長的帶狀正異常 (往往對應于隱伏的逆斷層)或負異常 (往往對應于隱伏的正斷層),也可以是由大剩余值構(gòu)成的,并列的正負剩余帶或梯度陡變帶?？傊?只要在殘差圖中有帶狀分布的異常出現(xiàn),就有理由區(qū)劃該異常為斷裂的反映。這種斷裂識別法是一種基礎判別的方法,能大概地確定斷裂的位置和走向。

2 可視化軟件系統(tǒng)

該趨勢面分析計算的可視化系統(tǒng)根據(jù)地質(zhì)測量所得三維數(shù)據(jù)為基礎,其核心程序包括:利用順序、選擇、循環(huán)結(jié)構(gòu)編制矩陣式 (3)的趨勢面擬合函數(shù)常系數(shù)的求解程序,及計算觀測點的趨勢值和剩余值,采用3次樣條插值函數(shù)與VB調(diào)用 CAD二次開發(fā)成圖技術(shù),將所得文件數(shù)據(jù)繪制成趨勢圖、殘差圖。系統(tǒng)具體的運行過程如圖 1示。

圖 1 趨勢面分析可視化系統(tǒng)技術(shù)路線

3 工程應用

大興井田區(qū)域內(nèi)地質(zhì)構(gòu)造發(fā)育[3]。大興礦地處北北東向展布的向斜盆狀構(gòu)造區(qū),西南側(cè)為井田盆狀構(gòu)造中心,也是沉降中心。北翼平緩,為 5°左右,南翼陡,為 20～30°,西翼陡,約 40°,東翼平緩,為 5°左右。盆軸南北兩端呈抬起形態(tài)。

3.1 Ⅱ級斷裂趨勢面分析

在大興井田包括邊界西部的較大區(qū)域內(nèi),依據(jù)地質(zhì)測量所得三維數(shù)據(jù),選用實用性較強的二次擬合函數(shù),利用可視化程序進行趨勢面分析計算,可得到擬合趨勢面函數(shù)相關(guān)參數(shù)值 (表 1)和大興井田Ⅱ級斷裂的二次殘差圖 (圖 2)。

表1 大興井田擬合趨勢面參數(shù)值

圖2 大興井田二次擬合殘差分布

分析圖 2作地質(zhì)解釋:在殘差圖上可明顯發(fā)現(xiàn)較大等高距的梯度陡變帶、并列的正負剩余帶。據(jù)斷裂構(gòu)造的地質(zhì)特征,可初步判定該殘差圖中 0線帶為曲線斷裂。該斷裂延伸區(qū)域范圍較大,可劃為Ⅱ級斷裂。將該Ⅱ級斷裂與相應地形圖,地震震中分布圖對照,其位置方向均吻合較好,因此,可判定該斷裂為活動斷裂。

3.2 大興井田活動斷塊區(qū)劃

對大興井田較大區(qū)域進行地質(zhì)動力區(qū)劃,在不同比例尺的地形圖上,利用趨勢面分析可視化程序系統(tǒng)劃分不同等級的斷裂。殘差圖結(jié)合各類圖紙資料與測量結(jié)果,準確劃分出各級活動斷裂。

Ⅰ,Ⅱ級活動斷裂主要依據(jù)殘差圖與震中分布圖綜合對照進行確定。同 1∶100×104遼寧及渤海周圍地震震中分布圖相比較,Ⅰ,Ⅱ級斷裂吻合的較好,Ⅰ,Ⅱ級斷裂基本上在震中分布帶上,在Ⅰ,Ⅱ級斷裂交界處附近,是地震分布的最大密度位置。以上表明:Ⅰ,Ⅱ級斷裂在現(xiàn)代地殼運動中起主導控制作用。

在劃分Ⅲ,Ⅳ,Ⅴ級活動斷裂時,利用本可視化系統(tǒng)進行趨勢面擬合分析后,考察主要溝谷、坡腳線、河流的膝狀彎曲、鞍部、等高線輪廓等特征,結(jié)合地下地質(zhì)勘查 (節(jié)理面、擦痕等構(gòu)造形跡)、地應力測量結(jié)果等多因素綜合分析來確定。大興井田區(qū)劃具體結(jié)果如圖 3示。

圖3 大興井田活動斷塊區(qū)劃結(jié)果

4 結(jié) 論

(1)地質(zhì)動力區(qū)劃為防治礦區(qū)多種地質(zhì)災害開創(chuàng)了新途徑,趨勢面分析法是地質(zhì)動力區(qū)劃最行之有效的方法。本文基于趨勢面分析理論與可視化編程軟件 VB調(diào)用 CAD二次開發(fā)成圖技術(shù),構(gòu)建了地質(zhì)動力區(qū)劃趨勢面分析計算的可視化系統(tǒng)。通過分析該可視化系統(tǒng)繪制的殘差圖,能夠有效確定斷裂位置和方向。

(2)在大興井田范圍內(nèi)進行趨勢面分析,結(jié)合礦區(qū)地質(zhì)地形圖、震中分布圖、地下勘察及地應力測量結(jié)果等資料,準確區(qū)劃出該井田區(qū)域內(nèi)的活動斷塊?；顒訑嗔褬?gòu)造附近是礦區(qū)地質(zhì)災害的高發(fā)區(qū),活動斷塊的準確劃分對煤礦開采設計、瓦斯動力事故預防等,具有重要指導作用和理論意義。

[1]И . М.佩圖霍夫,И . М.巴杜金娜.地下地質(zhì)動力學 [M].王 麗,陳學華譯.北京:煤炭工業(yè)出版社,2006.

[2]朱廣軼,高智廣,郝 哲,等.地質(zhì)動力區(qū)劃在鶴崗城區(qū)的應用 [J].遼寧工程技術(shù)大學學報,2005,24(5):646-648.

[3]朱廣軼,李保東,張 忠.大興井田地質(zhì)動力區(qū)劃研究 [J].煤礦安全,2005,36(3):1-3.

[4]張啟銳 .地質(zhì)趨勢面分析 [M].北京:科學出版社,1990.

[責任編輯:周景林]

Visualization Software System of Trend Face Analysis ofM ine GeologicalDynam ic D ivision

ZHU Guang-yi1,ZHENG Yang-fa1,ZHANG Hua-xing2

(1.Architecture Engineering School,ShenyangUniversity,Shenyang 110044,China; 2.CoalMining&DesigningDepartment,Tiandi Science&Technology Co.,Ltd,Beijing 100013,China)

Geological dynamic division was a new method of preventing geological disasters.This paper constructed a visualized software system of trend face analysis on the basis of trend face analysis theory,thrice spline interpolation function and CAD second developmentwith VB.This software could effectively find active fault by combining geological figure,epicenter contribution figure,underground exploration and geo-stress survey result,which provided scientific basis formining design and preventingmethane dynamic disaster.

geological dynamic division;geological trend face;visualization system;active fault

TD163

A

1006-6225(2010)06-0031-03

2010-07-15

國家重點科技攻關(guān)項目 (96-223-01-03-04)

朱廣軼 (1962-),男,遼寧喀左人,教授,研究生導師,主要從事礦區(qū)巖層移動控制與防災減災方面的研究。