模糊數(shù)學(xué)方法在唐莊煤礦構(gòu)造復(fù)雜性評價研究中的應(yīng)用

1、第21卷增刊22009年9月doi：10.3969/j.issn.1674-1803.2009.09.004文章編號：16741803（2009）S2000905中國煤炭地質(zhì)COALGEOLOGYOFCHINAVol21Sup.2Sep2009模糊數(shù)學(xué)方法在唐莊煤礦構(gòu)造復(fù)雜性評價研究中的應(yīng)用閆慶磊，朱炎銘，姚海鵬，曹新款（中國礦業(yè)大學(xué)資源與地球科學(xué)學(xué)院，江蘇徐州221008）摘要：模糊數(shù)學(xué)法中的隸屬函數(shù)、模糊綜合評判、模糊聚類等方法均可以用于構(gòu)造復(fù)雜性評價。本文采用模糊綜合評判方法，以唐莊煤礦的構(gòu)造復(fù)雜性評價研究為例，討論了模糊數(shù)學(xué)理論在構(gòu)造復(fù)雜性評價中的應(yīng)用。將唐莊煤礦18煤層劃分為37個單

2、元，選用適當(dāng)?shù)哪：阕臃謩e對各單元進(jìn)行評價，然后按最大隸屬度原則進(jìn)行歸一化，得出最終結(jié)果，從而確定每個評價單元對應(yīng)的構(gòu)造復(fù)雜程度。關(guān)鍵詞：模糊數(shù)學(xué)方法；模糊綜合評判；構(gòu)造復(fù)雜性；唐莊煤礦中圖分類號：TP399；TD163.1文獻(xiàn)標(biāo)識碼：AApplicationofFuzzyMathematicalMethodonStructuralComplexityAssessmentinTangzhuangCoalmineYanQinglei,ZhuYanming,YaoHaipengandCaoXinkuan(SchoolofResourceandGeosciences,CUMT,Xuzhou,Jian

3、gsu221008)Abstract:Themembershipfunction,fuzzysyntheticdiscriminationandfuzzyclusteringmethodscanbeusedinstructuralcomplexityassessment.Throughfuzzysyntheticdiscrimination,carriedoutstructuralcomplexityassessmentintheTangzhuangcoalmine,divided18coalseamsinthecoalmineinto37units,properfuzzyoperatorsa

4、reselectedtoassesseachunitrespectively,thennormalizingaccordingtothemaximummembershipprinciple.Finalresultis:32unitsbelongtoIIIcategories,accountfor86.5%intotal;5unitsbelongtoIIIIVcategories;inthemass,coalminestructureisinclinedtosimpletype,onlyNos.31and36unitsonthesouthboundarystructurallycomplex.K

5、eywords:fuzzymathematicalmethod;fuzzysyntheticdiscrimination;structuralcomplexity;Tangzhuangcoalmine模糊數(shù)學(xué)法是美國加里福尼亞大學(xué)ZadehLA教授1965年創(chuàng)立的。模糊數(shù)學(xué)的評價方法在20世紀(jì)80年代初始被應(yīng)用于地質(zhì)領(lǐng)域，首先被應(yīng)用于礦業(yè)找礦、巖石分類等方面。在20世紀(jì)90年代初，徐鳳銀首次將模糊綜合評判方法應(yīng)用于礦井構(gòu)造的定量評價和預(yù)測，取得了成功的經(jīng)驗1。經(jīng)過多年來的不斷實踐和完善，模糊數(shù)學(xué)方法已經(jīng)成為礦井構(gòu)造評價領(lǐng)域的一個比較成熟的方法。模糊變換，其數(shù)學(xué)模型為2：假定有兩個論域分別為：因

6、素集U=u1，u2，un，ui為評價因素；評語集V=1，2，m，i為評語等級。對U中全部因素分別進(jìn)行單因素評價，則可獲得從U到V的一個模糊關(guān)系R，RF（U×V）模糊綜合評判的正演方法正在得到日益廣泛的應(yīng)用。這種方法是先由因素集U、評浯集V及單因素評判集R（從U到V的一個模糊映射）構(gòu)成一個綜合評判模型，再根據(jù)各因素的相對重要性給定一個因素權(quán)重集A，經(jīng)A與R的模糊合成，得到一個多因素綜合評判集S。模糊綜合評判的實質(zhì)是一個作者簡介：閆慶磊（1983），男，山東濟寧人，構(gòu)造地質(zhì)專業(yè)，中國礦業(yè)大學(xué)資源學(xué)院碩07。收稿日期：2009-05-26責(zé)任編輯：唐錦秀rn1rn2rnm稱其為單因素評判矩

7、陣，一般0rij1（i=1，2，n；j=1，2，m）。若RF（U×V）一經(jīng)確定，而且存在一個模糊向量AF（U），則可唯一確定一個由U到V的模糊變換B：1基本原理和步驟R=r11r12r1mr21r22r2m，（1）rn1rn2rnm=（b1，b2，bn）B=A×R=（a1，a2，an）×r11r12r1mr21r22r2m（2）10中國煤炭地質(zhì)第21卷其中A=a1，a2，an為評價因素權(quán)重集，可視為因素集上的模糊子集，表達(dá)式為：約3.5km，南北寬約1.2km，面積2.18km2。位于賈汪復(fù)向斜的NW翼，總體上為東西向展布的三角形單斜構(gòu)造，傾向NWW，產(chǎn)狀較為平緩

8、，傾角局部最小為10°，最大為41°，一般為15°20°，走向發(fā)育有次級向斜、背斜各一個，斷裂較為發(fā)育，主要斷層走向集中在NE和NW兩方向（圖2）。A=a+a+a，u1u2un（3）nai為第i個評價因素相應(yīng)的權(quán)重系數(shù)，滿足ai=1。i=1利用模糊合成結(jié)果作為評判基礎(chǔ)，按照最大隸屬度原則，則可得出最后評判結(jié)果。其主要步驟見圖3唐莊礦井田構(gòu)造復(fù)雜程度模糊綜合評價根據(jù)模糊數(shù)學(xué)綜合評價的研究思路和方法，以1。18煤層為研究對象，對唐莊礦井田構(gòu)造的復(fù)雜程度進(jìn)行綜合評價。3.1評價指標(biāo)的選取唐莊煤礦斷層構(gòu)造較發(fā)育，褶皺寬緩，地層傾角一般不大，因此，在選擇評價指標(biāo)時

9、，主要考慮18煤斷層、褶皺的基本參數(shù)指標(biāo)，以實用為原則，選取了煤層傾角變異系數(shù)（R）、斷層密度（M）、煤層絕對傾角（A）及煤厚變異系數(shù)（）四個評價指標(biāo)（表1）。煤層傾角變異系數(shù)R：指煤層底板等高線圖中某一單元內(nèi)傾角變化的標(biāo)準(zhǔn)差與其均值之比。它反映煤層傾角的變化，即小褶皺的發(fā)育程度3。表達(dá)式為：R=圖1模糊綜合評價計算程序框圖Figure1Fuzzycomprehensiveevaluationoftheprogramdiagram1軃）（a-ai軃a×100%（4）式中：ai煤層傾角，（°）；確定評價指標(biāo)，建立“評判因素集”，即：U=u1,u2,u3，un；計算各個指標(biāo)在評

10、價過程中相應(yīng)的權(quán)重系數(shù)，構(gòu)成因素集U上的權(quán)重向量A=a1,a2,a3an；權(quán)重的得出來源于灰色系統(tǒng)關(guān)聯(lián)度的歸一化。數(shù)。軃煤層平均傾角，（°）；a統(tǒng)計點數(shù)。N斷層密度M：指單位面積內(nèi)發(fā)育的斷層條表達(dá)式為：確定各個指標(biāo)的評語論域，形成評語集，即V=v1,v2,v3vm；就確定的評語論域，對每個評語指標(biāo)u1,u2,u3，un相對于每個評語等級作出評價，構(gòu)造模糊相關(guān)矩陣R=（rij）n×m；M=N式中：N斷層條數(shù)，條；（5）統(tǒng)計面積，km2。S煤層傾角絕對值A(chǔ)：指統(tǒng)計單元或塊段內(nèi)煤層傾角的絕對值，表達(dá)式為n選擇恰當(dāng)?shù)哪：铣伤阕樱⒕C合評判的數(shù)學(xué)模型，求出評價的基本結(jié)果，再按最大

11、隸屬度原則，得出最終評判結(jié)果，確定出每個評價單元或塊段對應(yīng)的合理性程度。A=1nai=1i（6）式中：ai煤層傾角，（°）；統(tǒng)計點數(shù)。n煤厚變化系數(shù)：該指標(biāo)是反映單元內(nèi)煤層厚度變化偏離平均厚度程度的參數(shù)。計算公式為：=S×100%（7）2唐莊煤礦地質(zhì)概況唐莊煤礦位于徐州市賈汪區(qū)青山泉鎮(zhèn)境內(nèi)，徐淮及相鄰地區(qū)總體構(gòu)造線方向為近EW向，NNE向次之。礦區(qū)北起F29斷層，南到F5斷層止；西自太原組各煤層露頭起，東到F29、F5兩斷層交匯點。東西長軓M增刊2閆慶磊，等：模糊數(shù)學(xué)方法在唐莊煤礦構(gòu)造復(fù)雜性評價研究中的應(yīng)用表1唐莊礦井田各構(gòu)造單元評價指標(biāo)統(tǒng)計表Table1TablesofT

12、angzhuangminefieldevaluationofallstructuralunits11單元號煤層傾角煤層傾角煤厚變異系數(shù)煤層絕對傾角A斷層密度M單元號煤厚變異系數(shù)煤層絕對傾角A斷層密度M變異系數(shù)R變異系數(shù)R123456789101112131415161718190.02330.03950.00850.0130.08270.13650.10540.13590.22520.02860.07710.11520.12990.16390.04680.06450.34220.15290.1010.02560.01890.01140.03270.02260.01890.09190.16650

13、.09370.01890.01860.02770.02770.01670.1670.11770.03920.02930.018922.22522.57520.02519.919.921.67519.1251922.7521.82528.7519.2519.2515.251615.716.2520.2510.650.05580.09760.03080.04340.02160.01340.02550.01120.0448000.09870.45860.09860.23930.0330.01510.02230.02282021222324252627282930313233343536370.170

14、70.2470.40180.38850.12910.12960.13220.20320.07270.080.33230.13910.05230.30540.10690.19260.06360.05130.02720.10350.01140.01540.04570.13330.03720.02510.04450.00910.07080.12130.01130.02140.02140.03720.07290.0815.511.3758.59.2520.521.2519.12516.7519.77518.7519.2523.7530.1520.47520.7526.529.52519.50.0657

15、0.0220.02120.03730.03340.01270.135500.02970.05480.08450.10010.06330.02220.05080.02180.03340.074圖218煤層底板等高線圖Figure218coalfloorcontourmap其中S=式中：S煤厚標(biāo)準(zhǔn)差；姨軓）2（Mi-Mi=1R×M：母指標(biāo)，為煤層傾角變異系數(shù)與斷層密（8）度的乘積，綜合反映褶皺與斷層的影響程度。3.2權(quán)重的計算權(quán)重是指某一評價因素在決定構(gòu)造復(fù)雜程度時的重要性程度4。為了綜合、科學(xué)的評價構(gòu)造復(fù)雜程度，選擇斷層密度M與煤層傾角變異系數(shù)R的乘積為母指標(biāo)來綜合反映斷層與褶曲的共同

16、影響，計算各評價指標(biāo)的權(quán)重系數(shù)（表2）。軓評定區(qū)內(nèi)的平均厚度；M每個見煤點的實測煤厚；Mi參加評定的見煤點數(shù)；ni=1，2，n。12表2構(gòu)造評價指標(biāo)權(quán)重系數(shù)中國煤炭地質(zhì)第21卷井的各個單元進(jìn)行構(gòu)造復(fù)雜程度綜合評價，得出評價的基本結(jié)果，再按最大隸屬度原則進(jìn)行歸一化，得出最終評價結(jié)果，從而確定每個評價單元對應(yīng)的構(gòu)造復(fù)雜程度?？紤]到Matlab語言對于數(shù)組的處理、運算具有其他編程語言無法比擬的便捷優(yōu)點，在計算模糊綜合評價結(jié)果和歸一化處理時，編程采用Table2Weightfactorsofstructuralevaluation地質(zhì)指標(biāo)傾角變異系數(shù)R煤厚變異系數(shù)煤層傾角A斷層密度MR×M0

17、.22645930.2135840.23500590.32495083.3評價等級劃分在已采區(qū)選一部分典型單元，統(tǒng)計各指標(biāo)的值，將每個指標(biāo)的測算值分別按從小到大的順序排列，然后用最優(yōu)分割法求出每個指標(biāo)四分后對應(yīng)每類構(gòu)造等級的取值范圍5。本次評價以統(tǒng)計結(jié)果為基礎(chǔ)，將唐莊礦已揭露的各指標(biāo)進(jìn)行大量統(tǒng)計，得出各指標(biāo)相對于各等級的取值范圍，作為綜合評價的基礎(chǔ)。唐莊礦18煤采用普通機械化方式開采，根據(jù)礦井地質(zhì)規(guī)程將井田構(gòu)造相對復(fù)雜程度劃分為簡單、較簡單、較復(fù)雜、復(fù)雜四個級別（相應(yīng)的評語集為Matlab語言。根據(jù)每個評價單元構(gòu)造復(fù)雜程度的歸一化結(jié)果，對應(yīng)于IIV構(gòu)造等級的評價單元數(shù)量關(guān)系如圖3所示。V=,）

18、，基本上能滿足要求。因此，采用模糊聚類，利用單元各指標(biāo)的統(tǒng)計結(jié)果，按四類的要求取截值進(jìn)行聚類，求出每類所包括的單元號，分為4組，然后按指標(biāo)排序，劃分出各組對應(yīng)的范圍值，對應(yīng)于IIV構(gòu)造等級，最后得出每類構(gòu)造復(fù)雜程度等級所對應(yīng)的各指標(biāo)取值范圍（表3）。表3唐莊礦井田各指標(biāo)構(gòu)造復(fù)雜程度標(biāo)準(zhǔn)Table3ThestandardindicatorsinthecomplexityconstructionofTangzhuangminefield項目煤層傾角變異系數(shù)R煤厚變異系數(shù)圖3各構(gòu)造等級對應(yīng)評價單元數(shù)目直方圖Figure3Histogramofthestructurallevelevaluationc

19、orrespondtothenumberofunits4結(jié)論根據(jù)以上的綜合評價計算結(jié)果，并將計算結(jié)果標(biāo)注到平面圖上，就得到唐莊礦的礦井構(gòu)造模糊綜合評價的最終成果圖（圖4）。從評判成果圖中可以看出，類單元共計簡單較簡單較復(fù)雜復(fù)雜項目0.00850.080.080.20320.20320.30550.30550.4018煤層絕對傾角A0.0090.02930.02930.04570.04570.09370.09370.167斷層密度M32個，占總數(shù)的86.5%，因此唐莊礦的礦井構(gòu)造總體上偏于簡單，這與該礦生產(chǎn)實踐中的驗證也是相吻合的。但在局部地段，仍存在構(gòu)造相對比較復(fù)雜的單元，其中，礦井構(gòu)造最復(fù)雜

20、的地段在井田構(gòu)造較為發(fā)育，切割關(guān)系明顯的地段31、36單元。井田范圍內(nèi)簡單較簡單較復(fù)雜復(fù)雜8.516.7516.7520.02520.02523.7523.7530.1500.01510.01510.04480.04480.09870.09870.45863.4隸屬函數(shù)的建立根據(jù)構(gòu)造評價指標(biāo)間的相互關(guān)系散點圖可知，指標(biāo)的評語子集符合正態(tài)性分布規(guī)律，因此，取正態(tài)型分布函數(shù)作為隸屬函數(shù)的表達(dá)式：軃）2xUx=exp-（x-22(9)所建立的隸屬函數(shù)包括各評價指標(biāo)相對每個評語等級的隸屬函數(shù)，評價指標(biāo)m個，評語等級n個，則隸屬函數(shù)應(yīng)有m×n個。在前面工作的基礎(chǔ)上，選擇恰當(dāng)?shù)哪：阕?，即可對礦

21、圖4唐莊煤礦18煤礦井構(gòu)造模糊綜合評判成果圖Figure4Tang,18coalminewellconstructedfuzzycomprehensiveevaluationresultsmap增刊2閆慶磊，等：模糊數(shù)學(xué)方法在唐莊煤礦構(gòu)造復(fù)雜性評價研究中的應(yīng)用13斷層的發(fā)育主要受邊界兩條正斷層F29和F5控制，次級褶曲較為不發(fā)育，斷層以EW向，NNE向為主，按照模糊數(shù)學(xué)評判的結(jié)果，在劃分的各個區(qū)塊當(dāng)中，絕大部分受構(gòu)造影響程度較小。參考文獻(xiàn)：1徐鳳銀，魏銘康.礦井構(gòu)造預(yù)測與評價的理論、方法及其應(yīng)用M.徐州：中國礦業(yè)大學(xué)出版社，1993.2楊倫標(biāo)，高英儀模糊數(shù)學(xué)原理及應(yīng)用M.廣州：華南理工大學(xué)出版

22、社，1993:1-161.3劉德民，曹代勇，李永軍，連會青礦井構(gòu)造復(fù)雜程度模糊綜合評判ComGIS的實現(xiàn)J.煤炭科學(xué)技術(shù)，2008，36（2）.4錢光謨，曹代勇，徐志斌，等.煤田構(gòu)造研究方法M.北京:煤炭工業(yè)出版社，1994:217-220.5王生全.模糊綜合評判法在礦井構(gòu)造復(fù)雜程度定量評價中的應(yīng)用J.西北地質(zhì)，1997，1.!（上接第3頁）在廣泛收集區(qū)域地質(zhì)、煤田地質(zhì)、重磁、地震資料及其他相關(guān)的地質(zhì)剖面、鉆探、基礎(chǔ)地質(zhì)研究成果資料的基礎(chǔ)上，通過系統(tǒng)綜合分析研究，對煤炭資源潛力進(jìn)行評價，圈定含煤遠(yuǎn)景區(qū)，進(jìn)一步了解含煤遠(yuǎn)景區(qū)含煤地層、煤層賦存特征并進(jìn)行資源潛力預(yù)測。略性礦產(chǎn)勘查新突破兩個方向上來

23、，把煤炭地質(zhì)調(diào)查資金和技術(shù)力量集中到國家急需的稀有稀缺、優(yōu)質(zhì)煤炭資源煤種和資源總量相對豐富、有利于緩解我國東部資源緊缺狀況或增強市場優(yōu)勢地位，并有望成為國家級大型煤炭遠(yuǎn)景區(qū)的重點勘查區(qū)上來。一是繼續(xù)開展優(yōu)質(zhì)煤炭資源潛力評價和稀有稀缺煤種潛力評價，在西部煤炭資源豐富地區(qū)，評價出一批新的大型煤炭后備基地，為我國煤炭工業(yè)可持續(xù)發(fā)展戰(zhàn)略服務(wù)；二是開展東部煤炭枯竭礦區(qū)深部及外圍煤炭資源調(diào)查，為礦山可持續(xù)發(fā)展發(fā)現(xiàn)一批新的接替資源后備基地；三是開展有條件的南方缺煤省份煤炭資源調(diào)查，發(fā)現(xiàn)一批新的煤炭礦產(chǎn)地，拉動商業(yè)性勘查，提高南方地區(qū)能源保障能力；四是開展全國特殊煤種（如液化用煤、氣化用煤和水煤漿用煤等）的選區(qū)評價工作，圈出特殊煤種保護(hù)性開發(fā)的規(guī)劃選區(qū)，為煤炭資源管理提供技術(shù)支撐；五是開展西北地區(qū)主要含煤盆地煤炭資源整體評價，摸清煤炭資源潛力，查明煤炭資源分布規(guī)律，探索新的煤炭地質(zhì)理論，指導(dǎo)今后煤炭地質(zhì)勘查。參考文獻(xiàn)：1林大揚，程愛國.搞好煤炭及相關(guān)領(lǐng)域的資源大