鄒　波

巖層真厚度的自動化計算

鄒波

（四川鑫順礦業(yè)股份有限公司，成都 610091）

目前地質(zhì)工作中巖層真厚度的計算一般都是采用列昂托夫斯基的傳統(tǒng)公式（中國地質(zhì)調(diào)查局地質(zhì)調(diào)查技術(shù)標(biāo)準(zhǔn)（DD2006-01），但公式中有一個正負(fù)號的判斷經(jīng)常困擾著廣大的地質(zhì)技術(shù)人員。文章里介紹了通過EXCEL編輯公式來實(shí)現(xiàn)巖層真厚度的自動化計算的方法，以方便廣大的地質(zhì)技術(shù)人員進(jìn)行巖層真厚度批量自動計算，提高工作效率，減少運(yùn)算錯誤。

；巖層；真厚度；計算

在野外地質(zhì)工作中，計算巖層真厚度是每個地質(zhì)技術(shù)人員經(jīng)常面臨的問題。實(shí)測剖面地層厚度計算受諸多因素影響（如坡向、地形坡度角、巖層產(chǎn)狀、導(dǎo)線方位等），在運(yùn)用列氏厚度計算公式H=L（sin α×cosβ×sinγ±cosα×sinβ）時，存在“+”、“-”號選取及正負(fù)地形坡度角等問題，稍微不慎將符號判斷錯誤，將導(dǎo)致計算結(jié)果錯誤。通過研究、圖形驗(yàn)證、結(jié)合實(shí)際工作，運(yùn)用EXCEL公式編輯，做成方便實(shí)用的程序，利于廣大地質(zhì)技術(shù)人員快速、準(zhǔn)確的計算巖層真厚度。

1　列氏公式概述

列氏H=L（sinα×cosβ×sinγ±cosα×sinβ）：式中：H-巖層真厚度（m）；L-導(dǎo)線斜距（m）；α -巖層傾角（°）；β-導(dǎo)線坡度角（ ±°）； γ -剖面導(dǎo)線與地層走向線的銳夾角（°）。當(dāng)坡向與巖層傾向相反時，公式中用加號計算；當(dāng)坡向與巖層傾向相同時，公式中用減號計算。

圖1　

2　應(yīng)用列氏公式計算巖層真厚度實(shí)例

用列式公式H=L（sinα×cosβ×sinγ± cosα×sinβ）計算時，導(dǎo)線坡度角是仰角（+）還是俯角（-），坡度角β的取值，在中國地質(zhì)調(diào)查局地質(zhì)調(diào)查技術(shù)標(biāo)準(zhǔn)（DD2006-01）中并未做說明，經(jīng)查閱資料，實(shí)際運(yùn)用證實(shí)計算時應(yīng)該對導(dǎo)線坡度角β取絕對值進(jìn)行計算才對。

示例：當(dāng)巖層測量的導(dǎo)線方位為90°，巖層產(chǎn)狀為270°∠70°，導(dǎo)線坡角有+30°和-30°的兩種情況。如果用正負(fù)號，當(dāng)導(dǎo)線坡角為負(fù)時，公式D=L（sinα×cosβ×sinγ±cosα×sinβ）的最后一項(xiàng)sinβ為負(fù)值，由于導(dǎo)線坡角為負(fù)時，坡向與巖層傾向相反，公式中的±號應(yīng)取“+”號，相當(dāng)于公式里括號中的兩項(xiàng)為相減；而當(dāng)導(dǎo)線坡角為正時，sinβ雖為正值，但坡向與巖層傾向相同，公式中的±號應(yīng)取“-”號，公式里括號中的兩項(xiàng)仍為相減。假設(shè)導(dǎo)線斜距均為1m時，兩種情況計算出來的真厚度均為0.64m。導(dǎo)線坡角使用了正負(fù)號導(dǎo)致兩種截然不同的情況計算出來的結(jié)果一樣，這樣，導(dǎo)線坡角使用了正負(fù)號就不正確了。同樣，而導(dǎo)線坡角取絕對值后，假設(shè)導(dǎo)線斜距均為1m，坡角為正時，計算的真厚度為0.64m，為負(fù)時，計算的真厚度為0.98m。

為了驗(yàn)證取導(dǎo)線坡度角β絕對值計算出來的結(jié)果的正確性，筆者在AutoCAD2004軟件中，將長度、坡度按實(shí)際比例尺畫圖，通過測量標(biāo)注的方式進(jìn)行驗(yàn)證，經(jīng)實(shí)際測量出的真厚度證實(shí)：無論導(dǎo)線坡度角是仰角（+）還是俯角（-），導(dǎo)線坡度角β必須取絕對值才是正確的。驗(yàn)證圖形如圖1。

從數(shù)值上看，導(dǎo)線坡度角β始終取正值即可使計算的結(jié)果正確，但為了實(shí)現(xiàn)使用EXCEL自動化計算的目的以及真實(shí)反映導(dǎo)線的坡度情況，導(dǎo)線坡度角β仍要有正負(fù)之分，以判斷坡向與巖層傾向，如后面3.2節(jié)所述。上面的圖形同時證實(shí)了列式公式中提到的剖面導(dǎo)線與地層走向線的夾角，除了銳夾角外，實(shí)際上要包含90度直角的情況。當(dāng)剖面導(dǎo)線與巖層走向正好垂直時，其夾角為90度，該公式同樣能適用。

圖2　

3　使用EXCEL自動化計算巖層真厚度

使用EXCEL進(jìn)行公式編輯，能較方便的實(shí)現(xiàn)巖層真厚度的自動化計算。實(shí)現(xiàn)該功能的難點(diǎn)在于剖面導(dǎo)線與地層走向線的銳夾角的計算以及列氏公式中正負(fù)號的判斷。下面對各功能的實(shí)現(xiàn)進(jìn)行詳細(xì)介紹。3.1 剖面導(dǎo)線與地層走向線的銳夾角的計算

1）巖層產(chǎn)狀的傾向小于270°的，先將巖地產(chǎn)狀加上90°，得到地層的走向線方位角；如果巖層產(chǎn)狀的傾向大于等于270°，則將傾向減去90°，得到地層的走向線方位角。

2）將剖面導(dǎo)線的方位角與地層走向線相減，再取絕對值，即使用較大的方位角度數(shù)減去較小的方位角度數(shù)，即得剖面導(dǎo)線與地層走向線的夾角。

3）分析上面第二步得到的剖面導(dǎo)線與地層走向線的夾角度數(shù)，如果這個度數(shù)大于等于90°，則減去90°，直到得到的夾角度數(shù)小于或等于90°為止，即得到剖面導(dǎo)線與地層走向線的銳夾角。

3.2列氏公式中正負(fù)號的判斷

1）坡向可能有兩個方向，當(dāng)?shù)匦纹露冉菫檎龝r，坡向應(yīng)與導(dǎo)線方位相反，即坡向的方位角為導(dǎo)線方位角加上180°；當(dāng)?shù)匦纹露冉菫樨?fù)時，坡向應(yīng)與導(dǎo)線方位相同，即坡向的方位角等于導(dǎo)線方位角。

2）將導(dǎo)線坡向的方位角與巖層傾向的方位角相減，再取絕對值，（即使用較大的方位角度數(shù)減去較小的方位角度數(shù)）即得坡向與巖層傾向的夾角。

3）分析上面第二步得到的坡向與巖層傾向的夾角，這個夾角可分為3種情況，即夾角小于90°，夾角大于90°且小于270°和夾角大于270 °。如果這個度數(shù)大于90°，坡向與巖層傾向相反，使用列氏公式時取正號；反之，如果這個度數(shù)小于90°，則坡向與巖層傾向相同，使用列氏公式時取負(fù)號。實(shí)際工作時，還會遇到坡向與巖層傾向分別位于第一、四象限的情況，有時兩者相減會得到一個大于270°的值，這種情況夾角實(shí)際也是小于90°，如坡向45度與傾向340 °之間的夾角實(shí)際是65°，而不是295°，畫方位草圖即可很方便的分析出兩者之間的關(guān)系是相同還是相反，如圖2。實(shí)測地質(zhì)剖面一般垂直巖層走向進(jìn)行，會有一定的夾角，小夾角的情況是實(shí)測地質(zhì)剖面、布置工程等工作要避免的，因此就不會出現(xiàn)坡向與巖層傾向的夾角等于90°的情況，即順層測量巖層或順層采樣。

3.3其它

在使用EXCEL公式運(yùn)算三角函數(shù)時，需要使用弧度的角度單位。條件判斷在EXCEL中使用IF條件函數(shù)實(shí)現(xiàn)，中間過程保存在臨時單元格中備用。

圖3　

4　制作方法

我們首先在EXCEL中做一個基本的巖層真厚度計算表格，表頭及輔助列如圖3所示，然后在G到S列編輯相應(yīng)的公式，最終在S列得到真厚度值。

從上圖中可以看到，G、H兩列為輔助列，用于保證I列得到不大于90°的夾角；N列也是輔助列，用于計算坡向的方位角，O列則是根據(jù)N列的結(jié)果來判斷當(dāng)前剖面測量的地層傾向與坡向是否相同或相反，以此判斷真厚度計算時應(yīng)該使用列氏公式的“加”法式還是“減”法式，用“+”和“-”來作為判斷標(biāo)志，為S列計算真厚度做準(zhǔn)備；J、K、L三列和P、Q、R三列則是分別計算各參數(shù)的正弦、余弦值；M、R兩列則分別求出列氏公式“±”前后兩部分的乘積。S列則按照列氏公式，根據(jù)B、M、R三列的值計算出巖層的真厚度。

為了便于美觀和使用方便，最后將輔助列G、H、N列全部隱藏。為了防止誤操作損壞程序，可將左側(cè)需要輸入原始編錄數(shù)據(jù)的區(qū)域設(shè)置為可編輯區(qū)域，用于輸入原始數(shù)據(jù)，其它區(qū)域設(shè)置為保護(hù)區(qū)域，或者將該電子表格文件設(shè)置為只讀格式。

5　總結(jié)

樣品真厚度的計算也是礦區(qū)工作中經(jīng)常遇到的。因樣品的真厚度與地質(zhì)剖面測量中某個巖層的真厚度實(shí)質(zhì)是一樣的，樣槽的方位角相當(dāng)于地質(zhì)剖面測量中導(dǎo)線的方位角；樣槽的坡角相當(dāng)于地質(zhì)剖面測量中的坡度角；樣長相當(dāng)于地質(zhì)剖面測量中的某一巖層的導(dǎo)線長，其余類推。因此，樣品真厚度的計算同樣可以使用這個巖層真厚度自動計算表來計算。

需要說明的是：實(shí)測地質(zhì)剖面或采樣時，由于有可能在某段地層的往復(fù)測量或采樣，如果真厚度值總是一個正值，這樣會造成計算的累計真厚度偏大。本計算程序計算的地層真厚度值有正有負(fù)，當(dāng)導(dǎo)線坡向與巖層傾向一致且導(dǎo)線坡度角絕對值大于巖層真傾角時，從巖層頂面向巖層底面的方向產(chǎn)生的厚度為正厚度，反之為負(fù)厚度。這樣就消除了連續(xù)往復(fù)測量地層時，地層真厚度累計偏大的錯誤。如不需要考慮值的正負(fù)，則取真厚度計算結(jié)果的絕對值即可。

［1］ DD2006-01，固體礦產(chǎn)勘查原始地質(zhì)編錄規(guī)程［S］.

［2］郭福生等.巖層真厚度計算公式述評［J］.中國煤田地質(zhì)，1999（2）.

［3］劉書才等.實(shí)測剖面地層厚度計算方法探討［J］.山東國土資源，2006（9）.

Automatic Calculation of the Actual Thickness of Rock Bed

ZOU Bo
（Sichuan Xinshun Mining Co. Ltd.， Chengdu 610091）

At present the actual thickness of rock bed is calculated by the three traditional formulas（CGS Standard :DD2006-01）， but plus or minus judgment in formulas is often very difficult. This paper gives an EXCEL editor formula to calculate automatically the actual thickness of rock bed.

rock bed； actual thickness； calculation

P622.22

A

1006-0995（2015）04-0632-0

10.3969/j.issn.1006-0995.2015.04.037

2015-02-25

鄒波（1977—），男，四川成都人，助理工程師，研究方向：礦山地質(zhì)