APP下载

岩层真厚度计算公式推导与编程

2021-05-20范茂煌

中国金属通报 2021年4期

范茂煌

(云南省有色地质局三〇六队,云南 昆明 650217)

在实测地质剖面数据处理过程中,岩层真厚度计算是核心问题。查阅相关文献发现,被引用较多的岩层真厚度计算公式为:

(1)岩层真厚度计算行业标准式[1]:D=L×(sinαcosβsinγ±sinβcosα)。

(2)岩层真厚度计算列昂托夫斯基式[2]:①当岩层倾向与地面坡向相反时:D=L(sinαcosβsinγ+cosαsinβ);②当岩层倾向与地面坡向相同且α>β时:D=L(sinαcosβsinγ-cosαsinβ);③当岩层倾向与地面坡向相同且α<β时:D=L(cosαsinβ-sinαcosβsinγ)。

(3)岩层真厚度计算徐开礼式[3]:D=L×|sinαcosβsinγ±sinβcosα|。以上式中:D为岩层真厚度,L为导线斜距,α为岩层真倾角,β为导线坡度角(非负数),γ为剖面导线与地层走向线的锐夹角,坡向与岩层倾向相反取“+”号,坡向与岩层倾向相同取“-”号。郭福生和林银山[2]提出当岩层倾向与地面坡向相同且α>β时,列昂托夫斯基式公式②有矛盾,例如:α=25°,β=23°,γ=60°,带入公式②为负值,与客观不符合,并指出列昂托夫斯基式②式、③式运用条件分别为α′>β、β>α′,其中α′为岩层在导线方向上的似倾角,tanα′=tanαsinγ。徐开礼式对取“-”号时取绝对值,未说明原由,且γ的具体计算未提及。刘纯瑶、吴盼泉、殷鉴等人在列昂托夫斯基岩层真厚度公式的推导[4]中,公式“±”号的取值推导借助岩层橫剖面,并在岩层橫剖面上分析:测线所在地形坡度与岩层相反、测线所在地形坡度与岩层相同且倾斜更缓、测线所在地形坡度与岩层相同且倾斜更陡三种情况,但未对其他情况分析。本文通过地质实体与几何模型相融合,采用五级地质实体参数作为讨论划分要素,确定了列昂托夫斯基②式、③式运用条件分别为α>β′、β′>α,其中,β′∈[β,90],且β′=arctan(tanβ/sinγ),详细介绍了剖面导线与地层走向线的锐夹角γ的计算公式,最终推导出各种岩层和测线空间组合下岩层真厚度计算的综合公式,并进行了编程。……

登录APP查看全文