CN106407569A - 一种厚松散层薄基岩条件下地表下沉值计算方法 - Google Patents
一种厚松散层薄基岩条件下地表下沉值计算方法 Download PDFInfo
- Publication number
- CN106407569A CN106407569A CN201610841258.3A CN201610841258A CN106407569A CN 106407569 A CN106407569 A CN 106407569A CN 201610841258 A CN201610841258 A CN 201610841258A CN 106407569 A CN106407569 A CN 106407569A
- Authority
- CN
- China
- Prior art keywords
- soil body
- value
- layer soil
- unconsolidated formation
- compression settlement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 33
- 238000004364 calculation method Methods 0.000 title claims abstract description 7
- 239000002689 soil Substances 0.000 claims abstract description 91
- 238000007906 compression Methods 0.000 claims abstract description 44
- 230000006835 compression Effects 0.000 claims abstract description 41
- 238000005065 mining Methods 0.000 claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 21
- 239000011435 rock Substances 0.000 claims description 13
- 230000010354 integration Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 235000008331 Pinus X rigitaeda Nutrition 0.000 claims description 2
- 235000011613 Pinus brutia Nutrition 0.000 claims description 2
- 241000018646 Pinus brutia Species 0.000 claims description 2
- 238000005303 weighing Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims 1
- 238000004519 manufacturing process Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 10
- 239000003245 coal Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000000205 computational method Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明提供了一种厚松散层薄基岩条件下地表下沉值计算方法,包括松散层土体原始数据采集,土体压缩沉降的判断,判断工作面采动程度,计算松散层土体压缩沉降值及数值修正匹配等五步。本发明计算方法简单,具有较强的理论基础,计算准确性高,且方便计算机程序语言的实现,一方面有效的提高了计算厚松散层薄基岩矿区地表下沉值的准确性,另一方面也厚松散层薄基岩矿区安全尘产提供了有力的技术支持。
Description
技术领域
本发明涉及一种厚松散层薄基岩条件下地表下沉值计算方法,属矿山开采沉陷与防护技术领域。
背景技术
煤炭是我国的主体能源,在我国的能源结构中大约占到百分之七十的比例,其在我国的国民经济与国计民生中占有举足轻重且不可或缺的地位。在我国华北、华中、华东等矿区赋存有厚松散层的地层,如鹤壁矿区地表的第四系松散层厚度为50~170m,永夏矿区松散层厚度达300多米,淮北矿区地表上部松散层厚度为270m,淮南矿区松散层厚度最大达467m,开滦矿区松散层厚度最大为300m。随着煤炭开采的持续发展,在厚松散层薄基岩矿区地表移动变形规律出现了很多异常的现象,例如,下沉系数大于1.0,地表移动范围较大等等。这些规律和参数的差异,为厚松散层薄基岩矿区地表移动变形预测、村庄保护煤柱留设及各类建构筑物保护带来了不确定性。因此针对这一现状,迫切需要开发一种能够准确计算厚松散层薄基岩矿区地表下沉值的方法,以满足实际使用的需要。
发明内容
本发明目的就在于克服上述不足,提供一种基于物理力学中的卡诺定理及土体力学相关理论的厚松散层薄基岩条件下地表下沉值计算方法。
为实现上述目的,本发明是通过以下技术方案来实现:
一种厚松散层薄基岩条件下地表下沉值计算方法,包括以下步骤:
(1)松散层土体原始数据采集,通过现场勘查及查阅原始地质勘查数据,获取待监测松散层土体的地质结构、各层土体厚度、密度、压缩模量及基岩面向土体传递空间的高度数据及常规概率积分法计算得到的地表下沉值,并存档备用;
(2)土体压缩沉降的判断,根据现场采集数据,计算基岩与松散层交界面形成的移动盆地的下沉系数q1、计算松散层土体表面形成的地表移动盆地的下沉系数q2,判断对比计算得到的q1和q2数值大小关系,且当q2值大于q1值时,则松散层土体存在压缩沉降,并进行下一步计算作业,q2值小于等于q1值时,则松散层土体未发生压缩沉降,并终止计算作业;
(3)判断工作面采动程度,根据现场采集数据,工作面采动程度判断计算函数模型为:L/H,
其中:
L为采掘工作面长度;
H采掘工作面深度;
且当采掘工作面长度与采掘工作面宽度都达到或超过1.2H0~1.4H0时,地表达到充分采动,其中H0为平均采深;
(4)计算松散层土体压缩沉降值,首先构建松散层土体压缩沉降计算函数模型,然后将步骤(1)中采集到的待监测松散层土体的地质结构、各层土体厚度、密度、压缩模量及基岩面向土体传递空间的高度数据带入到松散层土体压缩沉降计算函数模型中进行数据运算,并得出松散层土体压缩沉降值;
(5)数值修正匹配,将步骤(4)计算得到的松散层土体压缩沉降值与常规概率积分法计算得到的地表下沉值求和,并得到最终地表下沉值。
进一步的,步骤(2)中,还可根据上覆基岩面的最大下沉值与地表的最大下沉值比较土体压缩沉降的判断,当上覆基岩面的最大下沉值小于地表的最大下沉值时,松散层土体存在压缩沉降。
进一步的,步骤(3)中,在厚表土层薄基岩矿区,使用基岩厚度作为采深来衡量采动程度。
进一步的,步骤(4)中,松散层土体压缩沉降计算函数模型为:
其中:
W为松散层土体压缩沉降值;
Wi为松散层土体第i层的压缩下沉量;
n为松散层土体最大层数,n≥1;
其中Wi的计算函数模型为:
li为第i层土体的厚度,ρi为第i层土体的密度,Esi为第i层土体的压缩模量,g为重力加速度,Di为基岩面向第i层土体传递空间的高度,
其中,Di=Wi-1+di;
Wi-1为第i-1层土体的压缩沉降值;
di为第i-1层土体的厚度。
本发明计算方法简单,具有较强的理论基础,计算准确性高,且方便计算机程序语言的实现,一方面有效的提高了计算厚松散层薄基岩矿区地表下沉值的准确性,另一方面也厚松散层薄基岩矿区安全生产提供了有力的技术支持。
附图说明
图1为空间传递示意图;
图2为实测下沉值曲线与计算下沉值曲线;
图3为本发明实施流程图。
具体实施方式
以下结合具体实施例对本发明作具体的介绍。
本实施例以山东省某煤矿实际采掘作业情况中厚表土薄基岩地区开采地表移动变形规律、参数与兖州矿区在正常开采条件下获得的地表移动规律、参数存在较大的差异为例进行具体说明,其中工作面采厚2.9m,而地表最大下沉达到3.932m,基岩与松散层交界面形成的移动盆地的下沉系数达到0.7,松散层土体表面形成的地表移动盆地的下沉系数达到1.3。
如图1—3所示,一种厚松散层薄基岩条件下地表下沉值计算方法,包括以下步骤:
(1)松散层土体原始数据采集,通过现场勘查及查阅原始地质勘查数据,获取待监测松散层土体的地质结构、各层土体厚度、密度、压缩模量及基岩面向土体传递空间的高度数据及常规概率积分法计算得到的地表下沉值,并存档备用;
(2)土体压缩沉降的判断,根据现场采集数据,计算基岩与松散层交界面形成的移动盆地的下沉系数q1、计算松散层土体表面形成的地表移动盆地的下沉系数q2,判断对比计算得到的q1和q2数值大小关系,其中,q1为0.7,q2为1.3,q1严重小于q2,说明该工作面上方松散层土体存在压缩沉降;
(3)判断工作面采动程度,根据现场采集数据,以为采掘工作面长度与采掘工作面深度的比值可判断该工作面为充分采动程度。
(4)计算松散层土体压缩沉降值,首先构建松散层土体压缩沉降计算函数模型,然后将步骤(1)中采集到的待监测松散层土体的地质结构、各层土体厚度、密度、压缩模量及基岩面向土体传递空间的高度数据带入到松散层土体压缩沉降计算函数模型中进行数据运算,并得出松散层土体压缩沉降值;
(5)数值修正匹配,将步骤(4)计算得到的松散层土体压缩沉降值与常规概率积分法计算得到的地表下沉值求和,并得到最终地表下沉值。
本实施例中,步骤(2)中,还可根据上覆基岩面的最大下沉值与地表的最大下沉值比较土体压缩沉降的判断,当上覆基岩面的最大下沉值小于地表的最大下沉值时,松散层土体存在压缩沉降。
本实施例中,步骤(3)中,在厚表土层薄基岩矿区,使用基岩厚度作为采深来衡量采动程度。
本实施例中,步骤(4)中,松散层土体压缩沉降计算函数模型为:
其中:
W为松散层土体压缩沉降值;
Wi为松散层土体第i层的压缩下沉量;
n为松散层土体最大层数,n≥1;
其中Wi的计算函数模型为:
li为第i层土体的厚度,ρi为第i层土体的密度,Esi为第i层土体的压缩模量,g为重力加速度,Di为基岩面向第i层土体传递空间的高度,
其中,Di=Wi-1+di;
Wi-1为第i-1层土体的压缩沉降值;
di为第i-1层土体的厚度。
本实施例中,所述的松散层土体压缩沉降计算函数模型为基于卡诺定理基础推导获得。
从图2可以看出,通过该方法的计算结果符合实测下沉曲线,从而验证了该方法的正确性与可行性。
本发明计算方法简单,具有较强的理论基础,计算准确性高,且方便计算机程序语言的实现,一方面有效的提高了计算厚松散层薄基岩矿区地表下沉值的准确性,另一方面也厚松散层薄基岩矿区安全生产提供了有力的技术支持。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (4)
1.一种厚松散层薄基岩条件下地表下沉值计算方法,其特征在于:包括以下步骤:
(1)松散层土体原始数据采集,通过现场勘查及查阅原始地质勘查数据,获取待监测松散层土体的地质结构、各层土体厚度、密度、压缩模量及基岩面向土体传递空间的高度数据及常规概率积分法计算得到的地表下沉值,并存档备用;
(2)土体压缩沉降的判断,根据现场采集数据,计算基岩与松散层交界面形成的移动盆地的下沉系数q1、计算松散层土体表面形成的地表移动盆地的下沉系数q2,判断对比计算得到的q1和q2数值大小关系,且当q2值大于q1值时,则松散层土体存在压缩沉降,并进行下一步计算作业,q2值小于等于q1值时,则松散层土体未发生压缩沉降,并终止计算作业;
(3)判断工作面采动程度,根据现场采集数据,工作面采动程度判断计算函数模型为:L/H,
其中:
L为采掘工作面长度;
H采掘工作面深度;
且当采掘工作面长度与采掘工作面宽度都达到或超过1.2H0~1.4H0时,地表达到充分采动,其中H0为平均采深;
(4)计算松散层土体压缩沉降值,首先构建松散层土体压缩沉降计算函数模型,然后将步骤(1)中采集到的待监测松散层土体的地质结构、各层土体厚度、密度、压缩模量及基岩面向土体传递空间的高度数据带入到松散层土体压缩沉降计算函数模型中进行数据运算,并得出松散层土体压缩沉降值;
(5)数值修正匹配,将步骤(4)计算得到的松散层土体压缩沉降值与常规概率积分法计算得到的地表下沉值求和,并得到最终地表下沉值。
2.根据权利要求1所述的一种厚松散层薄基岩条件下地表下沉值计算方法,其特征在于,步骤(2)中,还可根据上覆基岩面的最大下沉值与地表的最大下沉值比较土体压缩沉降的判断,当上覆基岩面的最大下沉值小于地表的最大下沉值时,松散层土体存在压缩沉降。
3.根据权利要求1所述的一种厚松散层薄基岩条件下地表下沉值计算方法,其特征在于,步骤(3)中,在厚表土层薄基岩矿区,使用基岩厚度作为采深来衡量采动程度。
4.根据权利要求1所述的一种厚松散层薄基岩条件下地表下沉值计算方法,其特征在于,步骤(4)中,松散层土体压缩沉降计算函数模型为:
其中:
W为松散层土体压缩沉降值;
Wi为松散层土体第i层的压缩下沉量;
n为松散层土体最大层数,n≥1;
其中Wi的计算函数模型为:
li为第i层土体的厚度,ρi为第i层土体的密度,Esi为第i层土体的压缩模量,g为重力加速度,Di为基岩面向第i层土体传递空间的高度,
其中,Di=Wi-1+di;
Wi-1为第i-1层土体的压缩沉降值;
di为第i-1层土体的厚度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610841258.3A CN106407569B (zh) | 2016-09-22 | 2016-09-22 | 一种厚松散层薄基岩条件下地表下沉值计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610841258.3A CN106407569B (zh) | 2016-09-22 | 2016-09-22 | 一种厚松散层薄基岩条件下地表下沉值计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106407569A true CN106407569A (zh) | 2017-02-15 |
CN106407569B CN106407569B (zh) | 2019-10-08 |
Family
ID=57997340
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610841258.3A Expired - Fee Related CN106407569B (zh) | 2016-09-22 | 2016-09-22 | 一种厚松散层薄基岩条件下地表下沉值计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106407569B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110398227A (zh) * | 2019-07-29 | 2019-11-01 | 中国电建集团铁路建设有限公司 | 一种可溶岩场地近接地铁车站异常差异沉降监测方法 |
CN110555559A (zh) * | 2019-09-05 | 2019-12-10 | 首都师范大学 | 一种厚松散层非充分采动条件下地面沉降预计方法 |
CN112949106A (zh) * | 2020-11-26 | 2021-06-11 | 中南大学 | 一种岩土工程地质地表移动变形状态的检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609782A (zh) * | 2012-01-10 | 2012-07-25 | 中国矿业大学(北京) | 基于采煤沉陷模拟和复耕率的边采边复时机的优选方法 |
CN105926569A (zh) * | 2016-04-28 | 2016-09-07 | 河北地质大学 | 一种基于沉降监测数据的煤矿老采空区场地稳定性定量评价方法 |
-
2016
- 2016-09-22 CN CN201610841258.3A patent/CN106407569B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609782A (zh) * | 2012-01-10 | 2012-07-25 | 中国矿业大学(北京) | 基于采煤沉陷模拟和复耕率的边采边复时机的优选方法 |
CN105926569A (zh) * | 2016-04-28 | 2016-09-07 | 河北地质大学 | 一种基于沉降监测数据的煤矿老采空区场地稳定性定量评价方法 |
Non-Patent Citations (4)
Title |
---|
HONGRI ZHANG 等: "Mechanism and prediction research of the surface mining subsidence of thick loose layers in mining area", 《APPLIED MECHANICS AND MATERIALS》 * |
何荣 等: "大采深条带开采宽度确定方法研究", 《河南理工大学学报(自然科学版)》 * |
吴新辉: "司马矿厚黄土层薄基岩下开采地表移动变形规律研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
蔡来良 等: "急倾斜煤层开采最大下沉值对数预计法", 《煤矿安全》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110398227A (zh) * | 2019-07-29 | 2019-11-01 | 中国电建集团铁路建设有限公司 | 一种可溶岩场地近接地铁车站异常差异沉降监测方法 |
CN110555559A (zh) * | 2019-09-05 | 2019-12-10 | 首都师范大学 | 一种厚松散层非充分采动条件下地面沉降预计方法 |
CN110555559B (zh) * | 2019-09-05 | 2022-06-14 | 首都师范大学 | 一种厚松散层非充分采动条件下地面沉降预计方法 |
CN112949106A (zh) * | 2020-11-26 | 2021-06-11 | 中南大学 | 一种岩土工程地质地表移动变形状态的检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106407569B (zh) | 2019-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103995947B (zh) | 改进的煤层底板突水脆弱性评价方法 | |
CN105184864B (zh) | 一种用于天然地基换填工程量计算的场地地层三维地质结构模型生成方法 | |
CN110362905B (zh) | 一种盾构隧道开挖面失稳概率确定方法 | |
CN103953392B (zh) | 深埋隧道断面上岩爆风险位置判别方法 | |
CN103744128B (zh) | 一种用于地下洞室岩爆烈度等级的综合预报方法 | |
CN110991048B (zh) | 一种关闭井工矿地表沉陷预测方法 | |
CN104361211A (zh) | 岩石统计损伤本构模型的构建和应用方法 | |
CN106593524A (zh) | 一种固体充填开采地表沉陷动态预计方法 | |
CN106372297A (zh) | 砂土岩溶地层中盾构与溶洞安全垂直距离的确定方法 | |
Wu et al. | Upper limit analysis of stability of the water-resistant rock mass of a Karst tunnel face considering the seepage force | |
CN106407569A (zh) | 一种厚松散层薄基岩条件下地表下沉值计算方法 | |
CN108241772A (zh) | 考虑多重因素的裂隙承压含水层中隧道涌水量的确定方法 | |
CN103778480A (zh) | 一种基于敏感度分析的裂隙带高度预测方法 | |
CN108665146B (zh) | 大型露天矿山边坡稳定性评价精度确定方法 | |
CN107391778B (zh) | 一种圆形隧洞渗流量的解析计算方法 | |
CN101845815A (zh) | 一种软土地层含承压水基坑突涌塑性破坏计算处理方法 | |
CN106251053A (zh) | 一种采煤沉陷区地表建设高压输电线塔的适宜性评估方法 | |
CN104564069B (zh) | 一种基于方格网法的地面动态沉陷预测与复垦方法 | |
CN105631155A (zh) | 一种水库诱发地震概率计算方法 | |
CN103541376A (zh) | 采煤沉陷区地基在重复开采条件下的基础变形预测方法 | |
CN103046526A (zh) | 深基坑底面隆起远程智能监测三维数字预警方法与系统 | |
CN205428142U (zh) | 水气压力传感器监测预警预报系统 | |
Li et al. | Prediction of mining subsidence in shallow coal seam | |
CN111912953B (zh) | 一种基于挖方量监测的深凹矿山开采边坡稳定性测定方法 | |
CN106372295A (zh) | 砂土岩溶地层中盾构与溶洞安全水平距离的确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191008 |
|
CF01 | Termination of patent right due to non-payment of annual fee |