[go: up one dir, main page]

CN114595544B - A comprehensive safety evaluation method for buried pipelines in goafs of coal mines - Google Patents

A comprehensive safety evaluation method for buried pipelines in goafs of coal mines Download PDF

Info

Publication number
CN114595544B
CN114595544B CN202210500367.4A CN202210500367A CN114595544B CN 114595544 B CN114595544 B CN 114595544B CN 202210500367 A CN202210500367 A CN 202210500367A CN 114595544 B CN114595544 B CN 114595544B
Authority
CN
China
Prior art keywords
pipeline
early warning
goaf
subsidence
level
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.)
Expired - Fee Related
Application number
CN202210500367.4A
Other languages
Chinese (zh)
Other versions
CN114595544A (en
Inventor
李虎
安兆暾
刘思铭
唐雪梅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ocean University of China
Southwest Petroleum University
Original Assignee
Ocean University of China
Southwest Petroleum University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ocean University of China, Southwest Petroleum University filed Critical Ocean University of China
Priority to CN202210500367.4A priority Critical patent/CN114595544B/en
Publication of CN114595544A publication Critical patent/CN114595544A/en
Application granted granted Critical
Publication of CN114595544B publication Critical patent/CN114595544B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Geometry (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Animal Husbandry (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Agronomy & Crop Science (AREA)
  • Computational Mathematics (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Pipeline Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a comprehensive safety evaluation method for a buried pipeline in a coal mine goaf, which relates to the field of oil and gas transmission and comprises the steps of obtaining a surface subsidence displacement curve of the goaf; dividing a goaf surface displacement early warning interval according to the surface subsidence displacement curve; acquiring a pipeline dynamic stress curve, and dividing a pipeline dynamic early warning interval; and comprehensively evaluating the buried pipeline in the goaf to obtain a comprehensive evaluation grade by combining the displacement early warning interval of the surface of the goaf and the dynamic early warning interval of the pipeline. The invention provides a comprehensive safety evaluation method for a buried pipeline in a coal mine goaf, which aims to solve the problem that the comprehensive safety evaluation of the buried pipeline penetrating through the coal mine goaf is difficult in the prior art, realize the comprehensive evaluation of the dynamic risk of the buried pipeline penetrating through the coal mine goaf, provide a complete model for actual prediction and early warning and provide a guidance basis for engineering evaluation.

Description

一种煤矿采空区埋地管道的综合安全评价方法A comprehensive safety evaluation method for buried pipelines in coal mine goaf

技术领域technical field

本发明涉及油气输送领域,具体涉及一种煤矿采空区埋地管道的综合安全评价方法。The invention relates to the field of oil and gas transportation, in particular to a comprehensive safety evaluation method for buried pipelines in goafs of coal mines.

背景技术Background technique

随着煤矿资源的大量开采,采空区的形成范围逐渐加大。从而对地面诸多建筑物、构筑物的影响不断加大。石油、天然气作为主要能源之一,其长输油气管道作为最大的运输载体素有国民经济动脉之称。油气管道在运输过程中不可避免的会穿越很多复杂地域,进而会遇到许多地质灾害,而煤矿采空沉陷作为比较典型的灾害之一,油气管道在长距离输送过程中难免会穿越煤矿采空区,从而对油气管道造成巨大的破坏。With the massive mining of coal resources, the formation range of gobs has gradually increased. As a result, the impact on many buildings and structures on the ground continues to increase. Oil and natural gas are one of the main energy sources, and their long-distance oil and gas pipelines are known as the national economic artery as the largest transportation carrier. Oil and gas pipelines will inevitably pass through many complex regions during transportation, and will encounter many geological disasters. As one of the more typical disasters, coal mine goaf subsidence, oil and gas pipelines will inevitably pass through coal mine goafs during long-distance transportation. area, causing huge damage to oil and gas pipelines.

以西气东输管线敷设区域为例,西气东输管线一线工程敷设沿线途经9个省,路线途经约76个煤矿开采采空区,其中有8个重大采空区位于山西、宁夏等四个省内,总共受到煤矿采空影响的管段长度约为388千米,油气管道穿越煤矿采空区不仅会影响地下煤矿资源的开采、也会严重影响到油气管道的运输性能。Taking the West-East Gas Pipeline laying area as an example, the first-line project of the West-East Gas Pipeline passes through 9 provinces, and the route passes through about 76 coal mining goafs, of which 8 major goafs are located in Shanxi and Ningxia. In the province, the total length of pipeline sections affected by coal mine goafs is about 388 kilometers. Oil and gas pipelines passing through coal mine goafs will not only affect the mining of underground coal resources, but also seriously affect the transportation performance of oil and gas pipelines.

对于煤矿采空区而言,过量的开采极有可能造成地表开裂甚至坍塌等次生地质灾害。在油气管道穿越采空区时,由于地表的沉陷作用使管道以弯曲受拉为主,若穿越采空形成的迅速塌陷区可能使管道发生悬空甚至断裂。但是现有技术中,对于油气管道穿越煤矿采空区的临界破坏,还没有具体的、综合性的安全评判标准。For gobs in coal mines, excessive mining is very likely to cause secondary geological disasters such as surface cracking or even collapse. When the oil and gas pipeline passes through the goaf, due to the subsidence of the surface, the pipeline is mainly bent and pulled. If it passes through the goaf, the rapid subsidence area may cause the pipeline to hang or even break. However, in the prior art, there is no specific and comprehensive safety evaluation standard for the critical failure of oil and gas pipelines passing through gobs in coal mines.

发明内容Contents of the invention

本发明提供一种煤矿采空区埋地管道的综合安全评价方法,以解决现有技术中难以对穿越煤矿采空区的埋地管道进行综合安全评价的问题,实现综合评价穿越煤矿采空区的埋地管道的动态风险,为实际预测预警提供完整模型、为工程评价提供指导依据的目的。The invention provides a comprehensive safety evaluation method for buried pipelines in goafs in coal mines, to solve the problem in the prior art that it is difficult to conduct comprehensive safety evaluations for buried pipelines passing through goafs in coal mines, and to realize comprehensive evaluation of goafs in coal mines The purpose of providing a complete model for actual prediction and early warning, and providing a guiding basis for engineering evaluation.

本发明通过下述技术方案实现:The present invention realizes through following technical scheme:

一种煤矿采空区埋地管道的综合安全评价方法,包括:A comprehensive safety evaluation method for buried pipelines in goafs of coal mines, including:

获取采空区的地表下沉位移曲线;根据所述地表下沉位移曲线,划分采空区地表位移预警区间;Obtain the surface subsidence displacement curve of the goaf; divide the early warning interval of the goaf surface displacement according to the surface subsidence displacement curve;

获取管道动态应力曲线,划分管道动态预警区间;Obtain the pipeline dynamic stress curve and divide the pipeline dynamic early warning interval;

结合采空区地表位移预警区间、管道动态预警区间,对采空区埋地管道进行综合评价,得到综合评价等级。Combined with the goaf surface displacement early warning interval and the pipeline dynamic early warning interval, the buried pipeline in the goaf is comprehensively evaluated, and the comprehensive evaluation grade is obtained.

针对现有技术中难以对穿越煤矿采空区的埋地管道进行综合安全评价的问题,本发明提出一种煤矿采空区埋地管道的综合安全评价方法。本案发明人在研究过程中发现,地表沉陷变形是由地下煤层逐渐采空而形成的一个自下而上的逐渐影响过程,地表的下沉直接影响到此区域内埋地管道的安全运营,而现有技术完全忽略了这一因素的存在。为此,本方法首先获取采空区的地表下沉位移曲线,然后根据得到的地表下沉位移曲线划分采空区地表位移预警区间;之后,获取管道动态应力曲线,基于应力失效准则来划分管道动态预警区间;最后结合所得到的采空区地表位移预警区间、管道动态预警区间,来进行采空区埋地管道的综合评价。可以看出,本申请结合了采空区地表下沉位移、和穿越采空区的埋地管道的动态应力,建立起了适用于煤矿采空区埋地管道的综合评价方法,填补了现有技术的空白,可综合评价穿越煤矿采空区的埋地管道的动态安全风险,为实际预测预警提供完整模型、为工程评价提供指导依据的效果。Aiming at the problem in the prior art that it is difficult to perform comprehensive safety evaluation on buried pipelines passing through goafs in coal mines, the present invention proposes a comprehensive safety evaluation method for buried pipelines in goafs in coal mines. The inventors of this case found in the research process that the subsidence and deformation of the ground surface is a gradual bottom-up influence process formed by the gradual mining of underground coal seams. The subsidence of the ground surface directly affects the safe operation of buried pipelines in this area, while The prior art completely ignores the existence of this factor. To this end, this method first obtains the surface subsidence displacement curve of the goaf, and then divides the early warning interval of the goaf surface displacement according to the obtained surface subsidence displacement curve; after that, obtains the pipeline dynamic stress curve, and divides the pipeline based on the stress failure criterion Dynamic early-warning interval; finally, combined with the obtained goaf surface displacement early-warning interval and pipeline dynamic early-warning interval, the comprehensive evaluation of goaf-buried pipeline is carried out. It can be seen that this application combines the subsidence displacement of the goaf surface and the dynamic stress of the buried pipeline passing through the goaf, and establishes a comprehensive evaluation method suitable for the buried pipeline in the coal mine goaf, which fills the existing The gap in technology can comprehensively evaluate the dynamic safety risks of buried pipelines passing through coal mine goafs, provide a complete model for actual prediction and early warning, and provide guidance for engineering evaluation.

进一步的,获取采空区的地表下沉位移曲线的方法包括:Further, the method for obtaining the surface subsidence displacement curve of the goaf includes:

设置计算点,计算采空区地表的最大下沉值;Set the calculation point to calculate the maximum subsidence value of the goaf surface;

基于分段Knothe时间函数得到动态下沉方程,根据动态下沉方程绘制地表下沉位移曲线。The dynamic subsidence equation is obtained based on the segmented Knothe time function, and the surface subsidence displacement curve is drawn according to the dynamic subsidence equation.

Knothe时间函数模型符合蠕变理论,认为地表沉陷时刻的下沉速度是地表的最终下沉值,与此时刻动态下沉值之差成一定比例关系。分段Knothe时间函数是对原始Knothe时间函数的改进,其预测得到的地表沉降量更加精确、误差更小。本方案基于分段Knothe时间函数和最大下沉值的结合,得到地表下沉塌陷的方程表达,即可绘制出地表下沉位移曲线。The Knothe time function model conforms to the creep theory, and it is believed that the subsidence velocity at the moment of surface subsidence is the final subsidence value of the surface, which is proportional to the difference between the dynamic subsidence values at this moment. The piecewise Knothe time function is an improvement to the original Knothe time function, and its predicted surface subsidence is more accurate and has smaller errors. This scheme is based on the combination of the segmented Knothe time function and the maximum subsidence value to obtain the equation expression of the surface subsidence and subsidence, and then the surface subsidence displacement curve can be drawn.

此外,本领域技术人员应当理解,Knothe时间函数为本领域通用术语,无标准中文翻译。In addition, those skilled in the art should understand that the Knothe time function is a general term in this field, and there is no standard Chinese translation.

进一步的,所述采空区地表的最大下沉值通过如下公式计算:W m =hqcosψFurther, the maximum subsidence value of the goaf surface is calculated by the following formula: W m = hq cos ψ ;

式中,W m 为最大下沉值;h为煤层厚度;q为下沉系数;ψ为煤层倾角;In the formula, W m is the maximum subsidence value; h is the thickness of the coal seam; q is the subsidence coefficient; ψ is the inclination angle of the coal seam;

所述动态下沉方程为:The dynamic sinking equation is:

Figure 593591DEST_PATH_IMAGE001
Figure 593591DEST_PATH_IMAGE001

式中,W(t)为最大下沉点动态下沉值;t为时间;c为时间系数;τ为地表点出现最大下沉速度的时刻;T为下沉总时长;Ф 1(t)、Ф 2(t)为分段Knothe时间函数;e为自然对数。In the formula, W ( t ) is the dynamic subsidence value of the maximum subsidence point; t is the time; c is the time coefficient; τ is the moment when the maximum subsidence velocity appears at the surface point; , Ф 2 ( t ) is the piecewise Knothe time function; e is the natural logarithm.

进一步的,获取管道动态应力曲线的方法包括:Further, the method for obtaining the dynamic stress curve of the pipeline includes:

设置计算点,计算采空区地表的最大下沉值;Set the calculation point to calculate the maximum subsidence value of the goaf surface;

计算采空区对煤层的影响半径,以所述影响半径的边缘作为管道变形的下沉起始点;Calculate the radius of influence of the goaf on the coal seam, and use the edge of the radius of influence as the sinking starting point of the pipeline deformation;

计算管道在所述下沉起始点的弯矩;Calculate the bending moment of the pipeline at the sinking starting point;

根据所述最大下沉值、影响半径、管道在下沉起始点的弯矩,计算管道动态应力σ 1Calculate the dynamic stress σ 1 of the pipeline according to the maximum subsidence value, the radius of influence, and the bending moment of the pipeline at the starting point of subsidence;

计算管道因内压受到的应力σ 0Calculate the stress σ 0 of the pipeline due to internal pressure;

计算管道的最终应力σσ=σ 1+σ 0Calculate the ultimate stress σ of the pipeline: σ=σ 1 + σ 0 ;

根据最终应力σ绘制管道动态应力曲线。Draw the pipeline dynamic stress curve according to the ultimate stress σ .

管道动态应力曲线可反应埋地管道在穿过煤矿采空区时的应力应变情况,本方案的动态应力曲线以管道动态应力σ 1为基准,在此基础上考虑管道因内压所受的应力情况,以两者之和作为最终应力来绘制管道动态应力曲线。相较于常规的直接以动态应力作为应力分析依据的方式而言,具有更实际的工程应用价值。当然,其中的动态应力可由本领域技术人员根据现有的应力分析手段获得,在此不做赘述。The dynamic stress curve of the pipeline can reflect the stress and strain of the buried pipeline when it passes through the goaf of the coal mine. The dynamic stress curve of this scheme is based on the dynamic stress σ 1 of the pipeline. On this basis, the stress on the pipeline due to internal pressure is considered In this case, the sum of the two is used as the final stress to draw the pipeline dynamic stress curve. Compared with the conventional way of directly using dynamic stress as the basis of stress analysis, it has more practical engineering application value. Certainly, the dynamic stress therein can be obtained by those skilled in the art according to the existing stress analysis means, and will not be repeated here.

进一步的,所述采空区地表的最大下沉值通过如下公式计算:W m =hqcosψ;式中,W m 为最大下沉值;h为煤层厚度;q为下沉系数;ψ为煤层倾角;Further, the maximum subsidence value of the goaf surface is calculated by the following formula: W m = hq cos ψ ; where W m is the maximum subsidence value; h is the thickness of the coal seam; q is the subsidence coefficient; ψ is coal seam inclination;

所述影响半径通过如下公式计算:r=H/tanα;式中,r为影响半径,H为煤层采深,α为开采影响角;The influence radius is calculated by the following formula: r = H /tan α ; where r is the influence radius, H is the coal seam mining depth, and α is the mining influence angle;

若埋地管道在采空区属于管-土协同变形,则管道在下沉起始点的弯矩通过如下公式计算:If the buried pipeline belongs to the pipe-soil collaborative deformation in the goaf, the bending moment of the pipeline at the starting point of subsidence is calculated by the following formula:

Figure 350326DEST_PATH_IMAGE002
Figure 350326DEST_PATH_IMAGE002

式中,M A 为管道在下沉起始点的弯矩,W(l 1,γ)为管道下沉最大点的下沉值,L为下沉起始点到计算点的距离,EI为管道的弯曲刚度,λ为设定系数,且

Figure 692839DEST_PATH_IMAGE003
k为弹性地基系数;In the formula, M A is the bending moment of the pipeline at the sinking starting point, W ( l 1 , γ ) is the sinking value of the maximum sinking point of the pipeline, L is the distance from the sinking starting point to the calculation point, and EI is the bending moment of the pipeline stiffness, λ is the set coefficient, and
Figure 692839DEST_PATH_IMAGE003
, k is the coefficient of elastic foundation;

若埋地管道在采空区属于管-土非协同变形,则管道在下沉起始点的弯矩通过如下公式计算:If the buried pipeline belongs to the pipe-soil non-cooperative deformation in the goaf, the bending moment of the pipeline at the starting point of subsidence is calculated by the following formula:

Figure 717164DEST_PATH_IMAGE004
Figure 717164DEST_PATH_IMAGE004

式中,M A 为管道在下沉起始点的弯矩,EI为管道的弯曲刚度,y D 为暗悬空区域尾端的下沉值,x D 为下沉起始点到暗悬空区域尾端的距离,i为暗悬空区域尾端所对应的地表倾斜度,λ为设定系数,且

Figure 110100DEST_PATH_IMAGE005
k为弹性地基系数;In the formula, M A is the bending moment of the pipeline at the sinking starting point, EI is the bending stiffness of the pipeline, y D is the sinking value of the end of the dark suspended area, x D is the distance from the sinking starting point to the end of the dark suspended area, i is the surface slope corresponding to the end of the dark suspended area, λ is the setting coefficient, and
Figure 110100DEST_PATH_IMAGE005
, k is the coefficient of elastic foundation;

管道因内压受到的应力σ 0通过如下公式计算:σ 0=PD/2S;式中,P为管道内压,D为管道外径,S为管道壁厚。The stress σ 0 suffered by the pipeline due to internal pressure is calculated by the following formula: σ 0 = PD /2 S ; where P is the internal pressure of the pipeline, D is the outer diameter of the pipeline, and S is the wall thickness of the pipeline.

其中,管-土协同变形是指在开采初期,地表变化较小,此时埋地管道在小的地表变化过程中呈现出管道与土体一起变形的情况;Among them, the pipe-soil synergistic deformation refers to the situation that the buried pipeline deforms together with the soil during the small surface change at the initial stage of mining when the ground surface changes little;

管-土非协同变形是指随着开采的不断进行,地面沉降不断加大,当管道下方土体的变形大于管道变形,此时管道和土体的变形表现为在最大下沉点管道与土体发生分离,出现局部的管道下方暗悬空状态,此时称之为管-土非协同变形。在管-土非协同变形阶段,暗悬空区域的管道承受了管道自重以及管道承受到的来自上方的土体重量。而未有暗悬空状态的管-土部分,其管道下方土体的支撑力越靠近采空区边缘,其对管道的支撑力越大。随着开采的不断进行,开采的区域越来越大,形成的暗悬空区域就越来越大。当形成的暗悬空区域够大时,管道上方土体会因为周围土体的下沉而发生破裂。因此,管-土非协同变形是更为复杂的力学模型,但也是工程现场更容易遇到的实际情况。Pipe-soil non-coordinated deformation means that with the continuous mining, the land subsidence increases continuously. When the deformation of the soil below the pipeline is greater than the deformation of the pipeline, the deformation of the pipeline and the soil is manifested at the point of maximum subsidence between the pipeline and the soil. The body is separated, and a local dark suspended state under the pipe appears, which is called pipe-soil non-cooperative deformation. In the pipe-soil non-cooperative deformation stage, the pipe in the dark overhang area bears the weight of the pipe and the weight of the soil from above. For the pipe-soil part without dark overhang, the closer the support force of the soil below the pipe is to the edge of the goaf, the greater the support force for the pipe. As the mining continues, the mining area becomes larger and larger, and the dark overhang area formed becomes larger and larger. When the formed dark overhang area is large enough, the soil above the pipeline will be broken due to the subsidence of the surrounding soil. Therefore, pipe-soil non-cooperative deformation is a more complicated mechanical model, but it is also an actual situation that is more likely to be encountered in engineering sites.

本方案在计算弯矩的过程中,考虑了埋地管道在采空区属于管-土协同变形或管-土非协同变形两种情况,并针对性的分别提出了各自适用的弯矩计算公式,显著提高了后续动态应力的计算准确性和精度。In the process of calculating the bending moment, this scheme considers the two situations that the buried pipeline belongs to the pipe-soil synergistic deformation or the pipe-soil non-synergistic deformation in the goaf, and puts forward the respective applicable bending moment calculation formulas , significantly improving the accuracy and precision of subsequent dynamic stress calculations.

进一步的,采空区地表位移预警区间的划分方法包括:Further, the division method of the early warning interval of surface displacement in the goaf includes:

根据地表下沉位移曲线,绘制位移曲线切线角变化图;According to the surface subsidence displacement curve, draw the tangent angle change diagram of the displacement curve;

在位移曲线切线角变化图中,找出切线角等于45°的两个临界时间点,记为T 1T 2,其中T 1T 2In the tangent angle change diagram of the displacement curve, find two critical time points when the tangent angle is equal to 45°, which are marked as T 1 and T 2 , where T 1 < T 2 ;

划分采空区地表位移预警区间:Divide the early warning interval of surface displacement in the goaf:

若时间t满足t<T 1,则划分为初级预警;If the time t satisfies t< T 1 , it is classified as primary warning;

若时间t满足T 1≤t<T 2,则划分为中级预警;If the time t satisfies T 1 ≤ t < T 2 , it is classified as an intermediate warning;

若时间t满足T 2≤t,则划分为高级预警。If the time t satisfies T 2 ≤t, it is classified as advanced warning.

本方案引入位移曲线切线角对地表的位移变化进行科学合理的预警区间划分。This program introduces the tangent angle of the displacement curve to carry out scientific and reasonable early warning interval division for the displacement change of the surface.

本案发明人在大量研究过程中发现,采空区地表变形分为地表下沉初始期、活跃期与衰退期三个过程,其速度为由慢变快再变慢的过程。因此单位时间内地表变化速率形成的切线角小于45°时,可定义为地表下沉初始期即初始变形阶段;单位时间内地表变化速率形成的切线角大于45°时定义为地表位移活跃期即地表加速变形阶段;经过加速变形阶段后可将单位时间内地表变化速率形成的切线角小于45°时定义为地表位移衰退期即降速变形阶段。The inventors of this case found out during a large amount of research that the surface deformation of the goaf is divided into three processes: the initial period of surface subsidence, the active period and the recession period, and the speed is a process from slow to fast and then slow. Therefore, when the tangent angle formed by the surface change rate per unit time is less than 45°, it can be defined as the initial period of surface subsidence, that is, the initial deformation stage; when the tangent angle formed by the surface change rate per unit time is greater than 45°, it is defined as the active period of surface displacement, namely Surface accelerated deformation stage; after the accelerated deformation stage, when the tangent angle formed by the surface change rate per unit time is less than 45°, it can be defined as the surface displacement decay period, that is, the deceleration deformation stage.

因此本方案采用了位移曲线切线角变化图中切线角等于45°的两个临界时间点T 1T 2作为采空区地表位移预警区间划分的临界点。其中本领域技术人员应当理解,初级预警、中级预警、高级预警的风险等级依次增高。Therefore, this scheme adopts two critical time points T 1 and T 2 when the tangent angle of the displacement curve is equal to 45° in the change diagram of the tangent angle of the displacement curve are used as the critical points for the division of the early warning interval of the surface displacement of the goaf. Among them, those skilled in the art should understand that the risk levels of primary early warning, intermediate early warning, and advanced early warning increase sequentially.

当t<T 1时,地下煤矿开采面积不大,岩层应力重分布区域不大,导致岩层移动缓慢,岩层达到应力重新平衡状态所需的时间较长,使得地表的位移也随之缓慢移动。此阶段由于地表刚开始发生移动变形,移动速率不大,因此对管道的力学行为变化影响不大,采空过程对管道的威胁程度较小,此时管道随地表的缓慢变形而处于初始变形阶段。因此在此阶段管道与采空区地表危险程度都较小,所以将采空区地表位移预警区间划分为初级预警。When t < T 1 , the mining area of the underground coal mine is not large, and the stress redistribution area of the rock formation is not large, resulting in slow movement of the rock formation, and it takes a long time for the rock formation to reach the state of stress rebalance, so that the displacement of the surface also moves slowly. At this stage, the ground surface has just begun to move and deform, and the movement rate is not large, so it has little effect on the change of the mechanical behavior of the pipeline, and the threat of the goaf process to the pipeline is small. At this time, the pipeline is in the initial deformation stage with the slow deformation of the ground surface. . Therefore, at this stage, the pipeline and the goaf surface are less dangerous, so the early warning interval of the goaf surface displacement is divided into primary early warning.

T 1≤t<T 2时,煤矿开采面积达到一定规模,造成大面积的岩层移动,由于采空面积较大,岩层原始应力破坏范围较广,使上覆岩层需要更大面积进行应力重新平衡,此时岩层移动较快。由于此阶段地表变形速度较快,变形位移较大,因此在此阶段管道的力学行为将发生巨大的变化,采空区的逐渐开采对管道具有巨大的威胁。所以将此阶段采空区地表位移预警区间划分为中级预警。When T 1 ≤ t < T 2 , the mining area of the coal mine reaches a certain scale, resulting in a large area of rock movement. Due to the large goaf area, the original stress damage range of the rock layer is wide, and the overlying rock layer needs a larger area for stress restoration. Balanced, the rock formations move faster at this time. Due to the fast surface deformation and large deformation displacement at this stage, the mechanical behavior of the pipeline will change greatly at this stage, and the gradual mining of the goaf will pose a huge threat to the pipeline. Therefore, the early warning interval of surface displacement in the goaf at this stage is divided into intermediate early warning.

T 2≤t时,煤矿开采基本结束,由于在此之前开采的面积经过岩层活跃期的运动,使得在此阶段仅需要少量的岩层位移移动即可达到岩层应力的重新平衡,且此阶段的位移变化速率较慢。但是由于此阶段地表仍有部分位移,且此部分位移在某些时候决定着管道是否被损害,采空区存在塌陷、地质复杂多变,使管道应力应变增加,增加管道危险程度。因此将此阶段采空区地表位移预警区间划分为高级预警。When T 2 ≤ t, the coal mining is basically over. Since the area mined before that has gone through the movement of the active period of the strata, only a small amount of stratum displacement is required to achieve the rebalance of stratum stress at this stage, and the The rate of change of displacement is slow. However, due to the partial displacement of the ground surface at this stage, and this partial displacement determines whether the pipeline is damaged at some point, the goaf has subsidence, and the geology is complex and changeable, which increases the stress and strain of the pipeline and increases the risk of the pipeline. Therefore, the early warning interval of surface displacement in the goaf at this stage is divided into advanced early warning.

进一步的,管道动态预警区间的划分方法包括:Further, the division method of pipeline dynamic early warning interval includes:

若时间t满足t<T 1,则划分为一级预警;If the time t satisfies t< T 1 , it is classified as a first-level early warning;

若时间t满足T 1≤t<T 3,则划分为二级预警;If the time t satisfies T 1 ≤ t < T 3 , it is classified as a second-level early warning;

若时间t满足T 3≤t<T 4,则划分为三级预警;If the time t satisfies T 3 ≤t< T 4 , it is divided into three levels of early warning;

若时间t满足T 4≤t,则划分为四级预警;If the time t satisfies T 4 ≤ t, it is divided into four levels of early warning;

其中,T 3为三级预警临界时间点、T 4均为四级预警临界时间点;所述一级预警、二级预警、三级预警、四级预警的风险等级依次增高。Among them, T3 is the critical time point of the third-level early warning, and T4 is the critical time point of the fourth -level early warning; the risk levels of the first-level early warning, the second-level early warning, the third-level early warning, and the fourth-level early warning increase sequentially.

通过如下公式计算T 4 T 4 is calculated by the following formula:

Figure 478764DEST_PATH_IMAGE006
Figure 478764DEST_PATH_IMAGE006

式中,D为管道外径,M A 为管道在下沉起始点的弯矩,I为管道惯性矩,L为下沉起始点到计算点的距离,λ为设定系数,且

Figure 185820DEST_PATH_IMAGE007
k为弹性地基系数;Ф(T 4)为自变量取T 4时的Knothe时间函数;σ 为管道屈服应力;In the formula, D is the outer diameter of the pipeline, MA is the bending moment of the pipeline at the starting point of subsidence, I is the moment of inertia of the pipeline, L is the distance from the starting point of subsidence to the calculation point, λ is the setting coefficient, and
Figure 185820DEST_PATH_IMAGE007
, k is the elastic foundation coefficient; Ф ( T 4 ) is the Knothe time function when the independent variable is T 4 ; σyield is the yield stress of the pipeline;

通过如下公式计算T 3:T3= T4-T’;式中,T’ 为管道维护所需总时长。Calculate T 3 by the following formula: T 3 = T 4 - T'; where T' is the total time required for pipeline maintenance.

所述管道维护所需总时长通过如下公式计算:T’=t1+ t2+ t3;其中t1为应急动员时长,t2为人员出发到达灾害地点所需时长,t3为维护开工直至完成时长。The total time required for the maintenance of the pipeline is calculated by the following formula: T'=t 1 + t 2 + t 3 ; where t 1 is the length of emergency mobilization, t 2 is the time required for personnel to arrive at the disaster site, and t 3 is the maintenance start until completion time.

穿越采空区的埋地油气管道由于受到煤层逐渐采空的影响,其管道也会随地表的移动而以不同的速度进行变形,划分管道动态预警区间即是为了定量的评估这种变形所带来的安全风险。The buried oil and gas pipeline passing through the goaf is affected by the coal seam's gradual goafing, and its pipeline will also deform at different speeds with the movement of the ground surface. The purpose of dividing the pipeline dynamic warning interval is to quantitatively evaluate the deformation caused coming security risks.

管道正常运行穿越采空区时,管道应力应变以及管道位移都处于缓慢变形阶段,管道动态预警区间划分为一级预警;When the pipeline passes through the goaf in normal operation, the stress and strain of the pipeline and the displacement of the pipeline are in the slow deformation stage, and the pipeline dynamic warning interval is divided into a first-level warning;

T 1≤t<T 3时,管道上的应力应变突然开始加速或者管道上某段上的应力应变发生突变,也就是自T 1时刻之后,管道动态预警区间进入二级预警,此阶段由于管道应力应变变化速度较快,必须实时关注区域内管道内最大应力变化规律,并且需要开始做出应急响应的准备工作。When T 1 ≤ t < T 3 , the stress and strain on the pipeline suddenly start to accelerate or the stress and strain on a certain section of the pipeline suddenly change, that is, after T 1 , the dynamic warning interval of the pipeline enters the second-level early warning. The stress and strain of the pipeline change rapidly, so it is necessary to pay attention to the maximum stress change law in the pipeline in real time, and to start preparations for emergency response.

T 3≤t<T 4时,即是管道动态预警区间在经过二级预警阶段后,采空区持续下沉,使得管道快要接近屈服极限时,采空区地表变形可能仍处于活跃期,即管道变形速率仍然较快,对管道影响程度较大,将预警级别升级为三级预警。此时为了能在管道屈服之前做出及时的应急响应动员、出发、到达、动工、维护等一系列工作,又考虑某些煤矿工作区地质复杂不便于行动等因素,需在管道达到红色预警之前完成维护工作。When T 3 ≤ t < T 4 , it means that the goaf continues to sink after the pipeline dynamic warning interval passes through the second-level warning stage, making the pipeline approaching the yield limit, and the surface deformation of the goaf may still be in an active period. That is, the deformation rate of the pipeline is still relatively fast, and the impact on the pipeline is relatively large, so the warning level is upgraded to a third-level warning. At this time, in order to be able to make timely emergency response mobilization, departure, arrival, construction, maintenance and other series of work before the pipeline yields, and considering factors such as complex geology and inconvenient actions in some coal mine working areas, it is necessary to ensure that the pipeline reaches the red warning. Complete maintenance work.

二级预警和三级预警之间的时间界限为T 3,其具体取值可根据最大最大应力变化速率、应急点与管道灾害位置距离、维抢修队伍的演训技术水平、应急装备技术水平、路况自然条件水平、预期失效严重程度而定。The time limit between the second-level early warning and the third-level early warning is T 3 . It depends on the level of natural road conditions and the expected severity of failure.

T 4≤t时,管道上最大应力即将达到屈服极限,若管道应力继续增加则可能发生管道破坏等严重的安全事故。因此需要在管道达到屈服极限之前完成管道抢修维护工作,确保管道能安全运行。考虑到若管道即将达到屈服极限时煤矿采空区仍然处于活跃期,地表变形仍然处于加速变形状态,即管道应力变化将迅速达到屈服极限。因此,取管道应力达到90%屈服极限为四级预警时间点,此时间点应为应急响应维护工作结束时间点,通过如下公式来计算出T 4When T 4 ≤t, the maximum stress on the pipeline will soon reach the yield limit, and if the pipeline stress continues to increase, serious safety accidents such as pipeline damage may occur. Therefore, it is necessary to complete the pipeline repair and maintenance work before the pipeline reaches the yield limit to ensure the safe operation of the pipeline. Considering that if the pipeline is about to reach the yield limit, the coal mine goaf is still in an active period, and the surface deformation is still in a state of accelerated deformation, that is, the stress change of the pipeline will quickly reach the yield limit. Therefore, take the pipeline stress reaching 90% yield limit as the fourth-level early warning time point, which should be the end time point of emergency response maintenance work, and T 4 is calculated by the following formula:

Figure 457795DEST_PATH_IMAGE008
Figure 457795DEST_PATH_IMAGE008

其中,Ф(T 4)为自变量取T 4时的Knothe时间函数;本领域技术人员应当理解,若采用的是分段Knothe时间函数,则根据T 4属于哪一分段时间内,采用对应公式即可。Among them, Ф ( T 4 ) is the Knothe time function when the independent variable is T 4 ; those skilled in the art should understand that if the piecewise Knothe time function is used, then according to which segment time T4 belongs to, the corresponding The formula will do.

进一步的,对采空区埋地管道进行综合评价的方法包括:Furthermore, the methods for comprehensive evaluation of buried pipelines in gobs include:

若管道动态预警区间属于四级预警,则综合评价等级为灾变级;If the pipeline dynamic early warning interval belongs to the four-level early warning, the comprehensive evaluation level is catastrophic level;

若管道动态预警区间属于三级预警,则综合评价等级为警报级;If the pipeline dynamic early warning interval belongs to the third-level early warning, the comprehensive evaluation level is the warning level;

若管道动态预警区间属于二级预警,则综合评价等级为警示级;If the pipeline dynamic early warning interval belongs to the second-level early warning, the comprehensive evaluation level is the warning level;

若管道动态预警区间属于一级预警、且采空区地表位移预警区间为初级预警,则综合评价等级为注意级;If the pipeline dynamic early warning interval belongs to the first-level early warning, and the surface displacement early warning interval of the goaf is the primary early warning, the comprehensive evaluation level is the attention level;

若管道动态预警区间属于一级预警、且采空区地表位移预警区间为中级预警或高级预警,则综合评价等级为警示级。If the pipeline dynamic early warning interval belongs to the first-level early warning, and the goaf surface displacement early warning interval is the intermediate early warning or high-level early warning, the comprehensive evaluation level is the warning level.

由于现有技术中还没有对穿越煤矿采空区的埋地管道进行综合预警的成熟模型,本申请考虑到在实际工程中,管道的应力应变不会受到理想状态下采空区开采的影响,管道可能会由于局部受剪或其他行为而导致管道力学行为突变,因此,本申请以管道力学行为为基础的管道动态预警区间模型作为主要预警预报指标,以表征地标沉降的采空区地表位移预警区间模型作为辅助预警预报指标,将两者结合,得到能够有效的评价采空区油气管道安全级别的综合模型。Since there is no mature model for comprehensive early warning of buried pipelines passing through coal mine goafs in the prior art, this application considers that in actual engineering, the stress and strain of pipelines will not be affected by the mining of goafs under ideal conditions. The pipeline may cause a sudden change in the mechanical behavior of the pipeline due to local shear or other behaviors. Therefore, this application uses the pipeline dynamic warning interval model based on the pipeline mechanical behavior as the main early warning and forecast index to represent the early warning of the ground displacement of the goaf in the landmark settlement. The interval model is used as an auxiliary early warning and forecast index, and the two are combined to obtain a comprehensive model that can effectively evaluate the safety level of oil and gas pipelines in goaf areas.

本发明与现有技术相比,具有如下的优点和有益效果:Compared with the prior art, the present invention has the following advantages and beneficial effects:

1、本发明一种煤矿采空区埋地管道的综合安全评价方法,结合了采空区地表下沉位移、和穿越采空区的埋地管道的动态应力,建立起了适用于煤矿采空区埋地管道的综合评价方法,填补了现有技术的空白,可综合评价穿越煤矿采空区的埋地管道的动态安全风险,为实际预测预警提供完整模型、为工程评价提供指导依据的效果。1. A comprehensive safety evaluation method for buried pipelines in goafs of coal mines according to the present invention combines the surface subsidence displacement of goafs and the dynamic stress of buried pipelines passing through goafs, and establishes a method suitable for goafs in coal mines. The comprehensive evaluation method of buried pipelines in the area fills the gaps in the existing technology, and can comprehensively evaluate the dynamic safety risks of buried pipelines passing through coal mine goafs, provide a complete model for actual prediction and early warning, and provide guidance for engineering evaluation. .

2、本发明一种煤矿采空区埋地管道的综合安全评价方法,充分考虑了采空区的地表下沉位移对埋地油气管道的作用,具有实际的工程应用价值。2. The present invention provides a comprehensive safety evaluation method for buried pipelines in goafs in coal mines, which fully considers the effect of surface subsidence and displacement of goafs on buried oil and gas pipelines, and has practical engineering application value.

3、本发明一种煤矿采空区埋地管道的综合安全评价方法,考虑管道因内压所受的应力情况,以内压和动态应力的组合作为最终应力来绘制管道动态应力曲线,具有更实际的工程应用价值。3. A comprehensive safety evaluation method for buried pipelines in goafs of coal mines according to the present invention considers the stress conditions of pipelines due to internal pressure, and uses the combination of internal pressure and dynamic stress as the final stress to draw pipeline dynamic stress curves, which is more practical. engineering application value.

4、本发明一种煤矿采空区埋地管道的综合安全评价方法,在计算弯矩的过程中考虑了埋地管道在采空区属于管-土协同变形或管-土非协同变形两种情况,并针对性的分别提出了各自适用的弯矩计算公式,显著提高了动态应力的计算准确性和精度。4. A comprehensive safety evaluation method for buried pipelines in goafs of coal mines according to the present invention. In the process of calculating bending moments, it is considered that buried pipelines in goafs belong to two types of pipe-soil synergistic deformation or pipe-soil non-cooperative deformation. According to the situation, the applicable bending moment calculation formulas are proposed respectively, which significantly improves the calculation accuracy and precision of dynamic stress.

5、本发明一种煤矿采空区埋地管道的综合安全评价方法,采用了位移曲线切线角变化图中切线角等于45°的两个临界时间点作为采空区地表位移预警区间划分的临界点,实现了采空区地表位移预警区间的划分。5. A comprehensive safety evaluation method for buried pipelines in goafs of coal mines according to the present invention adopts two critical time points when the tangent angles of the displacement curves are equal to 45° in the change diagram of the tangent angles of the displacement curves as the critical points for the division of early warning intervals of surface displacements in goafs. Points, the division of the early warning interval of surface displacement in the goaf is realized.

6、本发明一种煤矿采空区埋地管道的综合安全评价方法,通过多个临界时间点的设置和计算,实现了管道动态预警区间的划分。6. The present invention provides a comprehensive safety evaluation method for buried pipelines in goafs in coal mines, through the setting and calculation of multiple critical time points, the division of pipeline dynamic early warning intervals is realized.

7、本发明一种煤矿采空区埋地管道的综合安全评价方法,以管道力学行为为基础的管道动态预警区间模型作为主要预警预报指标,以表征地标沉降的采空区地表位移预警区间模型作为辅助预警预报指标,将两者结合,得到了能够有效评价采空区油气管道安全级别的综合模型。7. The present invention is a comprehensive safety evaluation method for buried pipelines in goafs in coal mines. The pipeline dynamic early warning interval model based on pipeline mechanical behavior is used as the main early warning and forecast index to represent the goaf surface displacement early warning interval model of landmark settlement As an auxiliary early warning and forecast index, combining the two, a comprehensive model that can effectively evaluate the safety level of oil and gas pipelines in goaf areas is obtained.

附图说明Description of drawings

此处所说明的附图用来提供对本发明实施例的进一步理解,构成本申请的一部分,并不构成对本发明实施例的限定。在附图中:The drawings described here are used to provide a further understanding of the embodiments of the present invention, constitute a part of the application, and do not limit the embodiments of the present invention. In the attached picture:

图1为本发明具体实施例的流程示意图;Fig. 1 is the schematic flow chart of the specific embodiment of the present invention;

图2为本发明具体实施例中的地表下沉位移曲线;Fig. 2 is the surface subsidence displacement curve in the specific embodiment of the present invention;

图3为本发明具体实施例中的位移曲线切线角变化图。Fig. 3 is a graph showing the variation of the tangent angle of the displacement curve in a specific embodiment of the present invention.

具体实施方式Detailed ways

为使本发明的目的、技术方案和优点更加清楚明白,下面结合实施例和附图,对本发明作进一步的详细说明,本发明的示意性实施方式及其说明仅用于解释本发明,并不作为对本发明的限定。在本申请的描述中,需要理解的是,术语“前”、“后”、“左”、“右”、“上”、“下”、“竖直”、“水平”、“高”、“低”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请保护范围的限制。In order to make the purpose, technical solutions and advantages of the present invention clearer, the present invention will be further described in detail below in conjunction with the examples and accompanying drawings. As a limitation of the present invention. In the description of this application, it is to be understood that the terms "front", "rear", "left", "right", "upper", "lower", "vertical", "horizontal", "high", The orientation or positional relationship indicated by "low", "inner", "outer" and so on is based on the orientation or positional relationship shown in the drawings, and is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying the referred device Or elements must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be construed as limiting the scope of protection of the present application.

实施例1:Example 1:

如图1所示的一种煤矿采空区埋地管道的综合安全评价方法,包括:As shown in Figure 1, a comprehensive safety evaluation method for buried pipelines in gobs in coal mines includes:

获取采空区的地表下沉位移曲线;根据所述地表下沉位移曲线,划分采空区地表位移预警区间;Obtain the surface subsidence displacement curve of the goaf; divide the early warning interval of the goaf surface displacement according to the surface subsidence displacement curve;

获取管道动态应力曲线,划分管道动态预警区间;Obtain the pipeline dynamic stress curve and divide the pipeline dynamic early warning interval;

结合采空区地表位移预警区间、管道动态预警区间,对采空区埋地管道进行综合评价,得到综合评价等级。Combined with the goaf surface displacement early warning interval and the pipeline dynamic early warning interval, the buried pipeline in the goaf is comprehensively evaluated, and the comprehensive evaluation grade is obtained.

其中,获取采空区的地表下沉位移曲线的方法为:Among them, the method to obtain the surface subsidence displacement curve of the goaf is:

设置计算点,计算采空区地表的最大下沉值;本实施例优选的将计算点设置在管道最大位移处。Set the calculation point to calculate the maximum subsidence value of the goaf surface; in this embodiment, the calculation point is preferably set at the maximum displacement of the pipeline.

基于分段Knothe时间函数得到动态下沉方程,根据动态下沉方程绘制地表下沉位移曲线。The dynamic subsidence equation is obtained based on the segmented Knothe time function, and the surface subsidence displacement curve is drawn according to the dynamic subsidence equation.

本实施例中的分段Knothe时间函数表达为:The segmented Knothe time function in the present embodiment is expressed as:

Figure 21631DEST_PATH_IMAGE009
Figure 21631DEST_PATH_IMAGE009

式中,Ф(t)代表Knothe时间函数;Ф 1(t)、Ф 2(t)为分段Knothe时间函数;t为时间;c为时间系数;τ为地表点出现最大下沉速度的时刻;T为下沉总时长;e为自然对数。where Ф ( t ) represents the Knothe time function; Ф 1 ( t ) and Ф 2 ( t ) are the piecewise Knothe time functions; t is the time; c is the time coefficient; ; T is the total time of sinking; e is the natural logarithm.

其中,采空区地表的最大下沉值通过如下公式计算:W m =hqcosψAmong them, the maximum subsidence value of the goaf surface is calculated by the following formula: W m = hq cos ψ ;

式中,W m 为最大下沉值;h为煤层厚度;q为下沉系数;ψ为煤层倾角;In the formula, W m is the maximum subsidence value; h is the thickness of the coal seam; q is the subsidence coefficient; ψ is the inclination angle of the coal seam;

动态下沉方程表达为:The dynamic sinking equation is expressed as:

Figure 487379DEST_PATH_IMAGE001
Figure 487379DEST_PATH_IMAGE001

式中,W(t)为最大下沉点动态下沉值;t为时间;c为时间系数;τ为地表点出现最大下沉速度的时刻;T为下沉总时长;Ф 1(t)、Ф 2(t)为分段Knothe时间函数;e为自然对数。In the formula, W ( t ) is the dynamic subsidence value of the maximum subsidence point; t is the time; c is the time coefficient; τ is the moment when the maximum subsidence velocity appears at the surface point; , Ф 2 ( t ) is the piecewise Knothe time function; e is the natural logarithm.

其中,获取管道动态应力曲线的方法包括:Among them, the methods for obtaining the dynamic stress curve of the pipeline include:

设置计算点,计算采空区地表的最大下沉值;Set the calculation point to calculate the maximum subsidence value of the goaf surface;

计算采空区对煤层的影响半径,以所述影响半径的边缘作为管道变形的下沉起始点;Calculate the radius of influence of the goaf on the coal seam, and use the edge of the radius of influence as the sinking starting point of the pipeline deformation;

计算管道在所述下沉起始点的弯矩;Calculate the bending moment of the pipeline at the sinking starting point;

根据所述最大下沉值、影响半径、管道在下沉起始点的弯矩,计算管道动态应力σ 1Calculate the dynamic stress σ 1 of the pipeline according to the maximum subsidence value, the radius of influence, and the bending moment of the pipeline at the starting point of subsidence;

计算管道因内压受到的应力σ 0Calculate the stress σ 0 of the pipeline due to internal pressure;

计算管道的最终应力σσ=σ 1+σ 0Calculate the ultimate stress σ of the pipeline: σ=σ 1 + σ 0 ;

最终应力σ绘制管道动态应力曲线。The ultimate stress σ draws the dynamic stress curve of the pipe.

其中,采空区地表的最大下沉值通过如下公式计算:W m =hqcosψAmong them, the maximum subsidence value of the goaf surface is calculated by the following formula: W m = hq cos ψ ;

式中,W m 为最大下沉值;h为煤层厚度;q为下沉系数;ψ为煤层倾角;In the formula, W m is the maximum subsidence value; h is the thickness of the coal seam; q is the subsidence coefficient; ψ is the inclination angle of the coal seam;

所述影响半径通过如下公式计算:r=H/tanα;式中,r为影响半径,H为煤层采深,α为开采影响角;The influence radius is calculated by the following formula: r = H /tan α ; where r is the influence radius, H is the coal seam mining depth, and α is the mining influence angle;

若埋地管道在采空区属于管-土协同变形,则管道在下沉起始点的弯矩通过如下公式计算:If the buried pipeline belongs to the pipe-soil collaborative deformation in the goaf, the bending moment of the pipeline at the starting point of subsidence is calculated by the following formula:

Figure 558978DEST_PATH_IMAGE002
Figure 558978DEST_PATH_IMAGE002

式中,M A 为管道在下沉起始点的弯矩,W(l 1,γ)为管道下沉最大点的下沉值,L为下沉起始点到计算点的距离,EI为管道的弯曲刚度,λ为系数,且

Figure 652836DEST_PATH_IMAGE010
k为弹性地基系数;In the formula, M A is the bending moment of the pipeline at the sinking starting point, W ( l 1 , γ ) is the sinking value of the maximum sinking point of the pipeline, L is the distance from the sinking starting point to the calculation point, and EI is the bending moment of the pipeline stiffness, λ is a coefficient, and
Figure 652836DEST_PATH_IMAGE010
, k is the coefficient of elastic foundation;

若埋地管道在采空区属于管-土非协同变形,则管道在下沉起始点的弯矩通过如下公式计算:If the buried pipeline belongs to the pipe-soil non-cooperative deformation in the goaf, the bending moment of the pipeline at the starting point of subsidence is calculated by the following formula:

Figure 653153DEST_PATH_IMAGE004
Figure 653153DEST_PATH_IMAGE004

式中,M A 为管道在下沉起始点的弯矩,EI为管道的弯曲刚度,y D 为暗悬空区域尾端的下沉值,x D 为下沉起始点到暗悬空区域尾端的距离,i为暗悬空区域尾端所对应的地表倾斜度,λ为设定系数,且

Figure 842082DEST_PATH_IMAGE011
k为弹性地基系数;In the formula, M A is the bending moment of the pipeline at the sinking starting point, EI is the bending stiffness of the pipeline, y D is the sinking value of the end of the dark suspended area, x D is the distance from the sinking starting point to the end of the dark suspended area, i is the surface slope corresponding to the end of the dark suspended area, λ is the setting coefficient, and
Figure 842082DEST_PATH_IMAGE011
, k is the coefficient of elastic foundation;

管道因内压受到的应力σ 0通过如下公式计算:σ 0=PD/2S;式中,P为管道内压,D为管道外径,S为管道壁厚。The stress σ 0 suffered by the pipeline due to internal pressure is calculated by the following formula: σ 0 = PD /2 S ; where P is the internal pressure of the pipeline, D is the outer diameter of the pipeline, and S is the wall thickness of the pipeline.

实施例2:Example 2:

一种煤矿采空区埋地管道的综合安全评价方法,在实施例1的基础上:A comprehensive safety evaluation method for buried pipelines in goafs of coal mines, on the basis of Embodiment 1:

(一)首先划分采空区地表位移预警区间:(1) First, divide the early warning interval of surface displacement in the goaf:

根据地表下沉位移曲线,绘制位移曲线切线角变化图;According to the surface subsidence displacement curve, draw the tangent angle change diagram of the displacement curve;

在位移曲线切线角变化图中,找出切线角等于45°的两个临界时间点,记为T 1T 2,其中T 1T 2In the tangent angle change diagram of the displacement curve, find two critical time points when the tangent angle is equal to 45°, which are marked as T 1 and T 2 , where T 1 < T 2 ;

划分采空区地表位移预警区间:Divide the early warning interval of surface displacement in the goaf:

若时间t满足t<T 1,则划分为初级预警;If the time t satisfies t< T 1 , it is classified as primary warning;

若时间t满足T 1≤t<T 2,则划分为中级预警;If the time t satisfies T 1 ≤ t < T 2 , it is classified as an intermediate warning;

若时间t满足T 2≤t,则划分为高级预警。If the time t satisfies T 2 ≤t, it is classified as advanced warning.

(二)再划分管道动态预警区间:(2) Re-dividing the pipeline dynamic early warning interval:

若时间t满足t<T 1,则划分为一级预警;If the time t satisfies t< T 1 , it is classified as a first-level early warning;

若时间t满足T 1≤t<T 3,则划分为二级预警;If the time t satisfies T 1 ≤ t < T 3 , it is classified as a second-level early warning;

若时间t满足T 3≤t<T 4,则划分为三级预警;If the time t satisfies T 3 ≤t< T 4 , it is divided into three levels of early warning;

若时间t满足T 4≤t,则划分为四级预警;If the time t satisfies T 4 ≤ t, it is divided into four levels of early warning;

其中,T 3为三级预警临界时间点、T 4均为四级预警临界时间点;所述一级预警、二级预警、三级预警、四级预警的风险等级依次增高。Among them, T3 is the critical time point of the third-level early warning, and T4 is the critical time point of the fourth -level early warning; the risk levels of the first-level early warning, the second-level early warning, the third-level early warning, and the fourth-level early warning increase sequentially.

其中,通过如下公式计算T 4Among them, T 4 is calculated by the following formula:

Figure 890941DEST_PATH_IMAGE006
Figure 890941DEST_PATH_IMAGE006

式中,D为管道外径,M A 为管道在下沉起始点的弯矩,I为管道惯性矩,L为下沉起始点到计算点的距离,λ为设定系数,且

Figure 370464DEST_PATH_IMAGE012
k为弹性地基系数;Ф(T 4)为自变量取T 4时的Knothe时间函数;σ 为管道屈服应力;In the formula, D is the outer diameter of the pipeline, MA is the bending moment of the pipeline at the starting point of subsidence, I is the moment of inertia of the pipeline, L is the distance from the starting point of subsidence to the calculation point, λ is the setting coefficient, and
Figure 370464DEST_PATH_IMAGE012
, k is the elastic foundation coefficient; Ф ( T 4 ) is the Knothe time function when the independent variable is T 4 ; σyield is the yield stress of the pipeline;

通过如下公式计算T 3:T3= T4-T’;式中,T’ 为管道维护所需总时长。Calculate T 3 by the following formula: T 3 = T 4 - T'; where T' is the total time required for pipeline maintenance.

管道维护所需总时长通过如下公式计算:T’=t1+ t2+ t3;其中t1为应急动员时长,t2为人员出发到达灾害地点所需时长,t3为维护开工直至完成时长。The total time required for pipeline maintenance is calculated by the following formula: T'=t 1 + t 2 + t 3 ; where t 1 is the length of emergency mobilization, t 2 is the time required for personnel to arrive at the disaster site, and t 3 is the time for maintenance to be completed duration.

(三)最后对采空区埋地管道进行综合评价:(3) Finally, make a comprehensive evaluation of the buried pipeline in the goaf:

若管道动态预警区间属于四级预警,则综合评价等级为灾变级;If the pipeline dynamic early warning interval belongs to the four-level early warning, the comprehensive evaluation level is catastrophic level;

若管道动态预警区间属于三级预警,则综合评价等级为警报级;If the pipeline dynamic early warning interval belongs to the third-level early warning, the comprehensive evaluation level is the warning level;

若管道动态预警区间属于二级预警,则综合评价等级为警示级;If the pipeline dynamic early warning interval belongs to the second-level early warning, the comprehensive evaluation level is the warning level;

若管道动态预警区间属于一级预警、且采空区地表位移预警区间为初级预警,则综合评价等级为注意级;If the pipeline dynamic early warning interval belongs to the first-level early warning, and the surface displacement early warning interval of the goaf is the primary early warning, the comprehensive evaluation level is the attention level;

若管道动态预警区间属于一级预警、且采空区地表位移预警区间为中级预警或高级预警,则综合评价等级为警示级。If the pipeline dynamic early warning interval belongs to the first-level early warning, and the goaf surface displacement early warning interval is the intermediate early warning or high-level early warning, the comprehensive evaluation level is the warning level.

本实施例以某煤矿采空区为例进行计算和评价。In this embodiment, a coal mine goaf is taken as an example for calculation and evaluation.

穿越该煤矿采空区的油气管道为X80管材,外径1024mm,壁厚18mm,设计内压10MPa。煤矿开采的工作面倾角4°-6°,开采走向尺寸571m,宽164m,平均采深260m,平均采厚7.5m。采用综采放顶煤开采,顶板管理方法为全部垮落法。开采沉陷影响传播角为86.2°,走向方向主要影响角正切值为1.9,倾斜方向主要影响角正切值平均为2.1,下沉系数为0.79,水平移动系数为0.35。观测站沿开采走向方向上布置29个点位。The oil and gas pipeline passing through the goaf of the coal mine is made of X80 pipe, with an outer diameter of 1024mm, a wall thickness of 18mm, and a design internal pressure of 10MPa. The inclination angle of the coal mining face is 4°-6°, the mining strike size is 571m, the width is 164m, the average mining depth is 260m, and the average mining thickness is 7.5m. Fully mechanized top-coal caving mining is adopted, and the roof management method is the total caving method. The influence propagation angle of mining subsidence is 86.2°, the tangent value of the main influence angle in the strike direction is 1.9, the average tangent value of the main influence angle in the inclined direction is 2.1, the subsidence coefficient is 0.79, and the horizontal movement coefficient is 0.35. Observation stations are arranged at 29 points along the mining direction.

计算该采空区地表的最大下沉值为-5.925m;下沉起始点的弯矩M A 为-8272.2N·m;c=0.037;下沉总时长T=364d;τ=0.5T=182d;r=136.84m;L=422.34m;k=4500N/m3E=2.1×1011pa;I=0.0037m4λ=0.0347;σ 0 =284.4MPa。The maximum subsidence value of the goaf surface is calculated to be -5.925m; the bending moment M A at the starting point of the subsidence is -8272.2N·m; c =0.037; the total subsidence time T =364d; τ =0.5 T =182d ; r =136.84m; L =422.34m; k =4500N/m 3 ; E= 2.1×10 11 pa; I = 0.0037m 4 ; λ= 0.0347 ;

管道动态应力σ 1为:The dynamic stress σ 1 of the pipeline is:

Figure 243479DEST_PATH_IMAGE013
Figure 243479DEST_PATH_IMAGE013

地表动态下沉方程为:The dynamic subsidence equation of the ground surface is:

Figure 808453DEST_PATH_IMAGE014
Figure 808453DEST_PATH_IMAGE014

根据地表动态下沉方程绘制地表下沉位移曲线,如图2所示。According to the dynamic surface subsidence equation, the surface subsidence displacement curve is drawn, as shown in Figure 2.

通过管道的最终应力公式σ=σ 1+σ 0,绘制管道动态应力曲线。Draw the dynamic stress curve of the pipeline through the final stress formula of the pipeline σ=σ 1 + σ 0 .

根据图2绘制位移曲线切线角变化图,如图3所示,可求得T 1=74d、T 2=301d。According to Fig. 2, draw the tangent angle change diagram of the displacement curve, as shown in Fig. 3, T 1 =74d, T 2 =301d can be obtained.

根据该管道所对应的油气输送公司的抢修动员时间,应急动员时长t1 =1d、人员出发到达灾害地点所需时长t2 =2d ,维护开工直至完成时长t3 =7d,所以管道维护所需总时长为T’=10d。According to the repair mobilization time of the oil and gas transportation company corresponding to the pipeline, the emergency mobilization time is t 1 = 1d, the time required for personnel to arrive at the disaster site is t 2 = 2d, and the time for maintenance to be completed is t 3 = 7d, so the pipeline maintenance requires The total duration is T' = 10d.

穿越采空区的油气管道屈服应力为550MPa,则由管道动态应力曲线可得达到屈服时间点的时间为249d。计算得到T4 =202d,T3 =192d。The yield stress of the oil and gas pipeline passing through the goaf is 550MPa, and the time to reach the yield point can be obtained from the pipeline dynamic stress curve as 249d. Calculated to get T 4 = 202d, T 3 = 192d.

通过上述计算,可得出本实施例的综合评价模型:Through the above calculation, the comprehensive evaluation model of this embodiment can be obtained:

Figure 926582DEST_PATH_IMAGE015
Figure 926582DEST_PATH_IMAGE015

以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The specific embodiments described above have further described the purpose, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above descriptions are only specific embodiments of the present invention and are not intended to limit the scope of the present invention. Protection scope, within the spirit and principles of the present invention, any modification, equivalent replacement, improvement, etc., shall be included in the protection scope of the present invention.

需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其它变体,意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。此外,在本文中使用的术语“连接”在不进行特别说明的情况下,可以是直接相连,也可以是经由其他部件间接相连。It should be noted that in this article, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply that there is a relationship between these entities or operations. any such actual relationship or order exists between them. Furthermore, the term "comprises", "comprising", or any other variation thereof, is intended to cover a non-exclusive inclusion such that a process, method, article, or apparatus comprising a set of elements includes not only those elements, but also includes elements not expressly listed. Other elements mentioned above, or also include elements inherent in such process, method, article or equipment. In addition, the term "connection" used herein may be directly connected or indirectly connected via other components unless otherwise specified.

Claims (5)

1.一种煤矿采空区埋地管道的综合安全评价方法,其特征在于,包括:1. A comprehensive safety evaluation method for buried pipelines in goafs of coal mines, characterized in that it comprises: 获取采空区的地表下沉位移曲线;根据所述地表下沉位移曲线,划分采空区地表位移预警区间;Obtain the surface subsidence displacement curve of the goaf; divide the early warning interval of the goaf surface displacement according to the surface subsidence displacement curve; 获取管道动态应力曲线,划分管道动态预警区间;Obtain the pipeline dynamic stress curve and divide the pipeline dynamic early warning interval; 结合采空区地表位移预警区间、管道动态预警区间,对采空区埋地管道进行综合评价,得到综合评价等级;Combined with the goaf surface displacement early warning interval and the pipeline dynamic early warning interval, the buried pipeline in the goaf is comprehensively evaluated, and the comprehensive evaluation grade is obtained; 采空区地表位移预警区间的划分方法包括:The division method of the early warning interval of surface displacement in the goaf includes: 根据地表下沉位移曲线,绘制位移曲线切线角变化图;According to the surface subsidence displacement curve, draw the tangent angle change diagram of the displacement curve; 在位移曲线切线角变化图中,找出切线角等于45°的两个临界时间点,记为T1、T2,其中T1<T2In the tangent angle change graph of the displacement curve, find two critical time points when the tangent angle is equal to 45°, which are marked as T 1 and T 2 , where T 1 <T 2 ; 划分采空区地表位移预警区间:Divide the early warning interval of surface displacement in the goaf: 若时间t满足t<T1,则划分为初级预警;If the time t satisfies t<T 1 , it is classified as primary warning; 若时间t满足T1≤t<T2,则划分为中级预警;If the time t satisfies T 1 ≤ t < T 2 , it is classified as an intermediate warning; 若时间t满足T2≤t,则划分为高级预警;If the time t satisfies T 2 ≤t, it is classified as advanced warning; 管道动态预警区间的划分方法包括:The division method of pipeline dynamic early warning interval includes: 若时间t满足t<T1,则划分为一级预警;If the time t satisfies t<T 1 , it is classified as a first-level early warning; 若时间t满足T1≤t<T3,则划分为二级预警;If the time t satisfies T 1 ≤ t < T 3 , it is classified as a second-level early warning; 若时间t满足T3≤t<T4,则划分为三级预警;If the time t satisfies T 3 ≤t<T 4 , it is divided into three levels of early warning; 若时间t满足T4≤t,则划分为四级预警;If the time t satisfies T 4 ≤ t, it is divided into four levels of early warning; 其中,T3为三级预警临界时间点、T4均为四级预警临界时间点;所述一级预警、二级预警、三级预警、四级预警的风险等级依次增高;Wherein, T3 is the critical time point of the third -level early warning, and T4 is the critical time point of the fourth -level early warning; the risk levels of the first-level early warning, the second-level early warning, the third-level early warning, and the fourth-level early warning are sequentially increased; 通过如下公式计算T4T 4 is calculated by the following formula:
Figure FDA0003816771170000011
Figure FDA0003816771170000011
式中,D为管道外径,MA为管道在下沉起始点的弯矩,I为管道惯性矩,L为下沉起始点到计算点的距离,λ为设定系数,且
Figure FDA0003816771170000012
k为弹性地基系数;Ф(T4)为自变量取T4时的Knothe时间函数;σ为管道屈服应力;EI为管道的弯曲刚度;
In the formula, D is the outer diameter of the pipeline, M A is the bending moment of the pipeline at the starting point of subsidence, I is the moment of inertia of the pipeline, L is the distance from the starting point of subsidence to the calculation point, λ is the setting coefficient, and
Figure FDA0003816771170000012
k is the elastic foundation coefficient; Ф(T 4 ) is the Knothe time function when the independent variable is T 4 ; σyield is the yield stress of the pipeline; EI is the bending stiffness of the pipeline;
通过如下公式计算T3:T3=T4-T’;式中,T’为管道维护所需总时长;Calculate T 3 by the following formula: T 3 =T 4 -T'; where T' is the total time required for pipeline maintenance; 所述管道维护所需总时长通过如下公式计算:T’=t1+t2+t3;其中t1为应急动员时长,t2为人员出发到达灾害地点所需时长,t3为维护开工直至完成时长;The total time required for the maintenance of the pipeline is calculated by the following formula: T'=t 1 +t 2 +t 3 ; where t 1 is the length of emergency mobilization, t 2 is the time required for personnel to arrive at the disaster site, and t 3 is the maintenance start duration until completion; 对采空区埋地管道进行综合评价的方法包括:Methods for comprehensive evaluation of buried pipelines in gobs include: 若管道动态预警区间属于四级预警,则综合评价等级为灾变级;If the pipeline dynamic early warning interval belongs to the four-level early warning, the comprehensive evaluation level is catastrophic level; 若管道动态预警区间属于三级预警,则综合评价等级为警报级;If the pipeline dynamic early warning interval belongs to the third-level early warning, the comprehensive evaluation level is the warning level; 若管道动态预警区间属于二级预警,则综合评价等级为警示级;If the pipeline dynamic early warning interval belongs to the second-level early warning, the comprehensive evaluation level is the warning level; 若管道动态预警区间属于一级预警、且采空区地表位移预警区间为初级预警,则综合评价等级为注意级;If the pipeline dynamic early warning interval belongs to the first-level early warning, and the surface displacement early warning interval of the goaf is the primary early warning, the comprehensive evaluation level is the attention level; 若管道动态预警区间属于一级预警、且采空区地表位移预警区间为中级预警或高级预警,则综合评价等级为警示级。If the pipeline dynamic early warning interval belongs to the first-level early warning, and the goaf surface displacement early warning interval is the intermediate early warning or high-level early warning, the comprehensive evaluation level is the warning level.
2.根据权利要求1所述的一种煤矿采空区埋地管道的综合安全评价方法,其特征在于,获取采空区的地表下沉位移曲线的方法包括:2. the comprehensive safety evaluation method of a kind of coal mine goaf buried pipeline according to claim 1, is characterized in that, the method for obtaining the surface subsidence displacement curve of goaf comprises: 设置计算点,计算采空区地表的最大下沉值;Set the calculation point to calculate the maximum subsidence value of the goaf surface; 基于分段Knothe时间函数得到动态下沉方程,根据动态下沉方程绘制地表下沉位移曲线。The dynamic subsidence equation is obtained based on the segmented Knothe time function, and the surface subsidence displacement curve is drawn according to the dynamic subsidence equation. 3.根据权利要求2所述的一种煤矿采空区埋地管道的综合安全评价方法,其特征在于,所述采空区地表的最大下沉值通过如下公式计算:3. The comprehensive safety evaluation method of a buried pipeline in a coal mine goaf according to claim 2, wherein the maximum subsidence value of the goaf surface is calculated by the following formula: Wm=hqcosψ;式中,Wm为最大下沉值;h为煤层厚度;q为下沉系数;ψ为煤层倾角;W m = hqcosψ; where, W m is the maximum subsidence value; h is the thickness of the coal seam; q is the subsidence coefficient; ψ is the inclination angle of the coal seam; 所述动态下沉方程为:The dynamic sinking equation is:
Figure FDA0003816771170000021
Figure FDA0003816771170000021
式中,W(t)为最大下沉点动态下沉值;t为时间;c为时间系数;τ为地表点出现最大下沉速度的时刻;T为下沉总时长;Ф1(t)、Ф2(t)为分段Knothe时间函数;e为自然对数。In the formula, W( t ) is the dynamic subsidence value of the maximum subsidence point; t is time; c is the time coefficient; τ is the moment when the surface point has the maximum subsidence velocity; , Ф 2 (t) is the piecewise Knothe time function; e is the natural logarithm.
4.根据权利要求1所述的一种煤矿采空区埋地管道的综合安全评价方法,其特征在于,获取管道动态应力曲线的方法包括:4. the comprehensive safety evaluation method of a kind of coal mine goaf buried pipeline according to claim 1, is characterized in that, the method for obtaining pipeline dynamic stress curve comprises: 设置计算点,计算采空区地表的最大下沉值;Set the calculation point to calculate the maximum subsidence value of the goaf surface; 计算采空区对煤层的影响半径,以所述影响半径的边缘作为管道变形的下沉起始点;Calculate the radius of influence of the goaf on the coal seam, and use the edge of the radius of influence as the sinking starting point of the pipeline deformation; 计算管道在所述下沉起始点的弯矩;Calculate the bending moment of the pipeline at the sinking starting point; 根据所述最大下沉值、影响半径、管道在下沉起始点的弯矩,计算管道动态应力σ1Calculate the dynamic stress σ 1 of the pipeline according to the maximum subsidence value, the radius of influence, and the bending moment of the pipeline at the starting point of subsidence; 计算管道因内压受到的应力σ0Calculate the stress σ 0 of the pipeline due to internal pressure; 计算管道的最终应力σ:σ=σ10Calculate the final stress σ of the pipeline: σ=σ 10 ; 根据最终应力σ绘制管道动态应力曲线。Draw the pipeline dynamic stress curve according to the ultimate stress σ. 5.根据权利要求4所述的一种煤矿采空区埋地管道的综合安全评价方法,其特征在于,5. the comprehensive safety evaluation method of a kind of coal mine goaf buried pipeline according to claim 4, is characterized in that, 所述采空区地表的最大下沉值通过如下公式计算:Wm=hqcosψ;式中,Wm为最大下沉值;h为煤层厚度;q为下沉系数;ψ为煤层倾角;The maximum subsidence value of the goaf surface is calculated by the following formula: W m = hqcosψ; where W m is the maximum subsidence value; h is the thickness of the coal seam; q is the subsidence coefficient; ψ is the inclination angle of the coal seam; 所述影响半径通过如下公式计算:r=H/tanα;式中,r为影响半径,H为煤层采深,α为开采影响角;The influence radius is calculated by the following formula: r=H/tanα; in the formula, r is the influence radius, H is the coal seam mining depth, and α is the mining influence angle; 若埋地管道在采空区属于管-土协同变形,则管道在下沉起始点的弯矩通过如下公式计算:If the buried pipeline belongs to the pipe-soil collaborative deformation in the goaf, the bending moment of the pipeline at the starting point of subsidence is calculated by the following formula:
Figure FDA0003816771170000031
Figure FDA0003816771170000031
式中,MA为管道在下沉起始点的弯矩,W(l1,γ)为管道下沉最大点的下沉值,L为管道下沉起始点到计算点的距离,EI为管道的弯曲刚度,λ为设定系数,且
Figure FDA0003816771170000032
k为弹性地基系数;
In the formula, M A is the bending moment of the pipeline at the starting point of sinking, W(l 1 ,γ) is the sinking value of the maximum sinking point of the pipeline, L is the distance from the starting point of the pipeline sinking to the calculation point, and EI is the Bending stiffness, λ is the setting coefficient, and
Figure FDA0003816771170000032
k is the coefficient of elastic foundation;
若埋地管道在采空区属于管-土非协同变形,则管道在下沉起始点的弯矩通过如下公式计算:If the buried pipeline belongs to the pipe-soil non-cooperative deformation in the goaf, the bending moment of the pipeline at the starting point of subsidence is calculated by the following formula:
Figure FDA0003816771170000033
Figure FDA0003816771170000033
式中,MA为管道在下沉起始点的弯矩,EI为管道的弯曲刚度,yD为暗悬空区域尾端的下沉值,xD为下沉起始点到暗悬空区域尾端的距离,i为暗悬空区域尾端所对应的地表倾斜度,λ为设定系数,且
Figure FDA0003816771170000034
k为弹性地基系数;
In the formula, M A is the bending moment of the pipeline at the sinking starting point, EI is the bending stiffness of the pipeline, y D is the sinking value of the end of the dark suspended area, x D is the distance from the sinking starting point to the end of the dark suspended area, i is the surface slope corresponding to the end of the dark suspended area, λ is the setting coefficient, and
Figure FDA0003816771170000034
k is the coefficient of elastic foundation;
管道因内压受到的应力σ0通过如下公式计算:σ0=PD/2S;式中,P为管道内压,D为管道外径,S为管道壁厚。The stress σ 0 of the pipeline due to internal pressure is calculated by the following formula: σ 0 =PD/2S; where P is the internal pressure of the pipeline, D is the outer diameter of the pipeline, and S is the wall thickness of the pipeline.
CN202210500367.4A 2022-05-10 2022-05-10 A comprehensive safety evaluation method for buried pipelines in goafs of coal mines Expired - Fee Related CN114595544B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210500367.4A CN114595544B (en) 2022-05-10 2022-05-10 A comprehensive safety evaluation method for buried pipelines in goafs of coal mines

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210500367.4A CN114595544B (en) 2022-05-10 2022-05-10 A comprehensive safety evaluation method for buried pipelines in goafs of coal mines

Publications (2)

Publication Number Publication Date
CN114595544A CN114595544A (en) 2022-06-07
CN114595544B true CN114595544B (en) 2022-11-11

Family

ID=81811677

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210500367.4A Expired - Fee Related CN114595544B (en) 2022-05-10 2022-05-10 A comprehensive safety evaluation method for buried pipelines in goafs of coal mines

Country Status (1)

Country Link
CN (1) CN114595544B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109372511A (en) * 2018-09-30 2019-02-22 中国矿业大学(北京) A method of utilizing Yellow River mud dynamic filling coal mining subsided land
CN114218518A (en) * 2021-12-16 2022-03-22 河南省地质矿产勘查开发局第二地质环境调查院 Method for measuring and calculating sedimentation amount of coal mine goaf

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1194199A (en) * 1967-10-03 1970-06-10 Pirelli General Cable Works Warning System for Indicating Dangerous Conditions in the Vicinity.
CN102345796B (en) * 2010-07-28 2013-02-13 中国石油天然气股份有限公司 Method for monitoring oil and gas pipeline in mining subsidence area
CN111156425B (en) * 2020-01-15 2021-04-20 中国石油化工股份有限公司 Pipeline condition monitoring method, device and system
CN113435636B (en) * 2021-06-18 2024-12-03 华北水利水电大学 A dynamic prediction method for the sinking process of surface points in coal mining subsidence
CN113670208B (en) * 2021-08-05 2023-10-27 国家石油天然气管网集团有限公司 A displacement monitoring system for oil and gas pipelines in goaf subsidence areas
CN113792429B (en) * 2021-09-14 2024-05-03 国家石油天然气管网集团有限公司 Method for predicting dynamic stress strain of buried pipeline in goaf by using time function

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109372511A (en) * 2018-09-30 2019-02-22 中国矿业大学(北京) A method of utilizing Yellow River mud dynamic filling coal mining subsided land
CN114218518A (en) * 2021-12-16 2022-03-22 河南省地质矿产勘查开发局第二地质环境调查院 Method for measuring and calculating sedimentation amount of coal mine goaf

Also Published As

Publication number Publication date
CN114595544A (en) 2022-06-07

Similar Documents

Publication Publication Date Title
Zhu et al. Challenges and development prospects of ultra-long and ultra-deep mountain tunnels
CN111058855B (en) Deformation control method and evaluation system for shield underpassing structure
CN107067333B (en) A method for monitoring the stability of high and steep slopes in high cold and high altitude
Chen et al. Blasting excavation induced damage of surrounding rock masses in deep-buried tunnels
CN106767205B (en) Micro-vibration comprehensive control blasting method for subway car section tunnel in bustling urban area
Liu et al. A case study of collapses at the yangshan tunnel of the coal transportation channel from the western inner Mongolia to the Central China
Du et al. Conception and Exploration of Using Data as a Service in Tunnel Construction with the NATM
CN106339798A (en) Soft rock inclined shaft supporting and protecting system stability evaluation method based on fuzzy comprehensive evaluation method
CN106897475B (en) Method and system for determining minimum thickness of covering soil layer of shallow underwater shield tunnel
CN106326636B (en) A kind of rock burst fragment ejection speed prediction technique based on releasable elastic strain energy
CN115809498A (en) A Calculation Method for Longitudinal Deformation of Deep-buried Shield Tunnel Passing Through Existing Lines
CN117151348A (en) A method to identify the type of water inrush from high-level separation layers in mining overlying rock
Jia et al. Identification of goaf instability under blasting disturbance using microseismic monitoring technology
CN110991009B (en) Method for determining stress deformation of buried pipeline by soil body loss below pipeline
Li et al. Study on major construction disasters and controlling technology at the Qingdao Kiaochow Bay subsea tunnel
Su et al. Failure mechanism analysis and treatment of tunnels built in karst fissure strata: A case study
Zhang et al. Fracture pattern of overlying strata in multiple coal seam mining in a physical model vis-à-vis MATLAB analysis and geological radar
CN114595544B (en) A comprehensive safety evaluation method for buried pipelines in goafs of coal mines
Wang et al. Study on rock stability and structural response of shield-driven twin tunnels crossing fault fracture zone based on 3d numerical simulation
CN105911606A (en) Mining overburden rock motion law in-situ observation drill hole layout method
Deng et al. Instability mode analysis of surrounding rocks in tunnel blasting construction with thin bedrock roofs
CN106249287B (en) A kind of determining method of coal mine travelling limit angle
Gou et al. Calculation of unusual stratum load for metro tunnel crossing active ground fissure considering different intersection angles
Dai et al. Water inrush mechanism and safety control in drilling and blasting construction of subsea tunnel
CN117211802A (en) Advanced blasting pressure relief method for tunnel

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20221020

Address after: 266100 Shandong Province, Qingdao city Laoshan District Songling Road No. 238

Applicant after: OCEAN University OF CHINA

Applicant after: SOUTHWEST PETROLEUM University

Address before: 610500, Xindu Avenue, Xindu District, Sichuan, Chengdu, 8

Applicant before: SOUTHWEST PETROLEUM University

GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20221111