CN110162900A - A method of simulation coal-bed flooding flow event - Google Patents

A method of simulation coal-bed flooding flow event Download PDF

Info

Publication number
CN110162900A
CN110162900A CN201910449651.1A CN201910449651A CN110162900A CN 110162900 A CN110162900 A CN 110162900A CN 201910449651 A CN201910449651 A CN 201910449651A CN 110162900 A CN110162900 A CN 110162900A
Authority
CN
China
Prior art keywords
coal
flow event
bed flooding
simulation
coal body
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.)
Pending
Application number
CN201910449651.1A
Other languages
Chinese (zh)
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.)
Liaoning Technical University
Original Assignee
Liaoning Technical 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 Liaoning Technical University filed Critical Liaoning Technical University
Priority to CN201910449651.1A priority Critical patent/CN110162900A/en
Publication of CN110162900A publication Critical patent/CN110162900A/en
Pending legal-status Critical Current

Links

Classifications

    • 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
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10081Computed x-ray tomography [CT]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20024Filtering details
    • G06T2207/20032Median filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

The present invention provides a kind of method for simulating coal-bed flooding flow event, is related to dust prevention and control technical field.The present invention obtains coal body original-gray image by CT scan, median filtering noise reduction process is carried out to coal body original gradation figure obtained by MATLAB software, obtain txt formatted file, it is then introduced into COMSOL software, coal body model is reconstructed with interpolating function, assigns model correlated variables parameter, model initial boundary conditions are set, network is divided, and numerical simulation is carried out to injection water coal sample flow event, obtains flowing law of the water under the true hole of coal body.Method provided by the invention is simple to operate, effectively shortens experimental period and cost, provides theoretical foundation and technical support to flowing law of the accurate description water in coal body, provides directive significance for coal-bed flooding depositing dust technique.

Description

A method of simulation coal-bed flooding flow event
Technical field
The present invention relates to dust prevention and control technical field more particularly to a kind of methods for simulating coal-bed flooding flow event.
Background technique
With the needs that mine coal production is exploited, the efficient mining process cutting tooth cutting shattering process of underground mine is always Along with the generation of a large amount of coal dusts, the coal miner that is exposed on working face under the conditions of high concentrate dust for a long time and easily suffer from serious Pneumoconiosis (CWP), while high concentrate dust is also serious restriction underground coal safety the main reason for inducing dust explosion disaster Production.Coal dust water filling is used as the dust-proof protrusion-dispelling of underground high-efficient, drive row gas main method, is increasingly used in producing Cheng Zhong.In order to study seepage flow flow process of the water in coal seam in the injecting process with pressure, need to be unfolded to coal body internal holes crazing The specific research of gap.However since coal body internal void fissured structure is complicated, the flow event of water can not be monitored effectively, so that coal The research of Moisture migration is difficult to move even one step in body.Currently, being directed to the study limitation in coal body internal void crack in numerical simulation Emulation, and the existing simulation software model basic condition setting based on finite element method design is relatively-stationary, Bu Nengmo Intend the process of the water filling seepage flow under true hole, i.e., the crack in model needs default and can not change, this is practical not with scene Symbol.So the method for coal and rock hole fissured structure can be constructed nondestructively really to probe into fortune of the water in coal body by being badly in need of one kind Dynamic process.
Summary of the invention
The technical problem to be solved by the present invention is in view of the above shortcomings of the prior art, provide a kind of simulation coal-bed flooding infiltration The method of stream process, can accurate description coal-bed flooding when coal sample seepage process, to hole change accurate description, analyze and research water In the mechanism and the characteristics of motion of coal body flowing, directive significance is provided for coal-bed flooding depositing dust technique.
In order to solve the above technical problems, the technical solution used in the present invention is:
A method of simulation coal-bed flooding flow event, comprising the following steps:
Step 1 carries out CT scan to the coal sample after test tape pressure injection water, obtains original-gray image;
Using Industrial X-ray CT detection system, by cone-beam scan and DR real time imagery Through Several Survey Measure, to coal sample into Row scanning imagery nondestructively scans the hole fissured structure inside coal body, obtains original-gray image;
Step 2 carries out median filtering noise reduction process to original-gray image with MATLAB, obtains txt formatted file;
The pore structural information in the original-gray image of CT scan is carried out using the I=imread function in MATLAB Identification reads one layer of picture of coal, and noise spot in picture is carried out noise reduction process using median filtering function, and use Fid=fopen function in MATLAB extracts each information put on image, exports as txt formatted file;
Step 3 imports obtained txt formatted file in COMSOL software, with interpolating function by coal body model weight Structure;
The size for corresponding to picture to txt formatted file obtained models, and picture size is 2.88cm × 3.01cm, Inputted in txt formatted file under %grid and input two groups of numbers, one group be 1-301, another group be 1-288, obtained %data is inputted before data, selects interpolating function to be imported in COMSOL global definition;
Step 4 assigns model correlated variables parameter;
Keying in density under global definition parameter in COMSOL software is 1000kg/m3, dynamic viscosity 0.001Pa, material Expect Attributions selection CoAI;
Step 5, setting model initial boundary conditions;
Physical field selects coefficient form boundary partial differential equation, and dependent variable physical quantity selects Head (m), the choosing of source item physical quantity Head is selected to change over time;Boundary setting is carried out to obtained geometrical model, the coboundary of coal body is set as the head of 20mm, Lower boundary is set as the head of 0mm, generates the head difference of 20mm, and driving fluid flows, and right boundary is impermeable;
Step 6 divides network, and carries out numerical simulation to injection water coal sample flow event;
According to the specification of selected materials, the unit size in crack is refined, the unit size of setting crack position is 57.2mm is arranged in 5mm, other positions;Setting be modeled as steady-flow, by simulation process obtain water in coal body motion process, The accurate description of wetted perimeter.
The beneficial effects of adopting the technical scheme are that simulation coal-bed flooding flow event provided by the invention Method image is handled by MATLAB median filtering using CT scan coal body image, and then obtain each of image Point gray value provides modeling foundation for the accurate microcosmic hole fissured structure of coal body that presents.With COMSOL interpolating function by coal body into Hole crack distribution characteristics on coal body is finely characterized, avoids and directly preset to coal body model by row model reconstruction Crack condition more can accurately express the hole crack inside coal body, more can accurately be consistent with scene.At MATLAB Image is managed, is finely characterized Coal Pore Structure with COMSOL simulation software (all suitable for the image after each CT scan With), and water filling simulation is carried out, it effectively optimizes and is lauched the rule flowed in coal body in micro-scale, using laboratory experiment The method combined with computer simulation further improves the accuracy of model and the reliability of experiment, is that coal-bed flooding is real It tests and provides reliably modeling method and calculation method, provide directive significance for coal-bed flooding depositing dust technique.
Detailed description of the invention
Fig. 1 is the method flow diagram of simulation coal-bed flooding flow event provided in an embodiment of the present invention;
Fig. 2 is CT scan equipment schematic provided in an embodiment of the present invention;
Fig. 3 is CT scan original graph provided in an embodiment of the present invention;
Fig. 4 is that median filtering provided in an embodiment of the present invention filters out noise figure;
Fig. 5 is coal body micropore structure illustraton of model provided in an embodiment of the present invention;
Fig. 6 is that prototype network provided in an embodiment of the present invention divides figure;
Fig. 7 is coal-bed flooding seepage simulation result figure provided in an embodiment of the present invention.
In figure: 1, radiographic source;2, cone-beam x-ray;3, coal sample;4, imaging plane;5, planar array detector.
Specific embodiment
With reference to the accompanying drawings and examples, specific embodiments of the present invention will be described in further detail.Implement below Example is not intended to limit the scope of the invention for illustrating the present invention.
Coal sample is cut into 4cm × 4cm × 10cm cuboid in use for laboratory certified reference coal cutting machine first, to coal body Water filling with pressure is carried out, water flow is made to penetrate into coal body.Coal-bed flooding flow event is simulated using method provided in this embodiment, is such as schemed Shown in 1, it is described that specific step is as follows.
Step 1 carries out CT scan to the coal sample after test tape pressure injection water, obtains original-gray image;
For CT scan equipment principle as shown in Fig. 2, coal sample 3 is lain in a horizontal plane on objective table, opening radiographic source 1 makes cone-beam Ray 2 penetrates coal sample 3, can be imaged on planar array detector 5, and complete image is chosen on planar array detector 5;
Step 2 carries out median filtering noise reduction process to original-gray image with MATLAB, obtains txt formatted file:
CT scan original graph is as shown in figure 3, noise spot in picture is handled using median filtering function, and utilizes The value of pixel each present on image is exported as txt formatted file by MATLAB.That is: a digital signal sequences are carried out When filtering processing, the long window that a length is odd number L, L=2n+1 are defined first, and window signal sample is (x1,x2, x3......xn), n size of data is arranged, x1<x2<x3<......<xn, then have:
Wherein, y is series of windows (x1,x2,x3......xn) intermediate value.
Window is slided in data, by window hit exactly pair pixel value replaced with the intermediate value of each pixel of window, it may be assumed that
Y (i)=med [x (i-N) ..., x (i) ..., x (i+N)]
Original-gray image eliminates the noise point of image, is effectively retained after the median filter process of MATLAB The pore character and boundary information of coal body, it is as shown in Figure 4 that median filtering filters out noise figure.
Txt formatted file is imported COMSOL by step 3, rebuilds coal body model with interpolating function.In group Selection parameter curved surface in geometry is built, the first parameter is s1, minimum value 0, maximum value 301.Second parameter is s2, minimum Value is 0, maximum value 288.Setting relative tolerance is 1e-1, and maximum knot number is 50.Coal body micropore structure illustraton of model such as Fig. 5 It is shown.
Step 4 assigns model correlated variables parameter:
Keying in density under global definition parameter in COMSOL is 1000kg/m3, dynamic viscosity 0.001Pa, material category Property selection CoAI.
Step 5, setting model initial boundary conditions:
The coboundary that coal body is arranged is the head of 20mm, and in the head of lower boundary setting 0mm, which produces the water of 20mm Head is poor, driving fluid flowing, and right boundary is impermeable.
According to but be not limited to cubic law and describe flowing law of the water in coal body, i.e., withAs The mathematical model of seepage simulation, wherein ρ is fluid density, and g is acceleration of gravity, and μ is the dynamic viscosity of fluid, and a is crack Porosity, KsFor hydraulic conductivity.Flux in crack is provided by speedThis depends on hydraulic conductivity Ks, water Head H=H (x, y) and hydraulic gradientAccording to hydraulic conductivity KsWith coefficient of transmissibility TsPer unit crack can also be calculated The displacement of widthAnd ideal conditions there are (for example, rough surface) when deviation, can be used thick Rough coefficient f adjusts cubic law, therefore coefficient of transmissibility isWherein describe what fluid in rock fracture flowed Potential-flow model uses Reynolds equationIt, can be into if the other Fluid Control Equations of researcher's addition in need Row addition manually.
Step 6 divides network, and carries out numerical simulation to injection water coal sample flow event;
According to the specification of selected materials, in order to accurately reach the precision of solution, it is desirable that the unit size in crack is refined. The unit size that crack position is arranged is 5mm, and 57.2mm is arranged in other positions.Setting is modeled as steady-flow.After network divides As shown in Figure 6.
Coal-bed flooding seepage simulation result figure is as shown in fig. 7, arrow direction as shown in the figure is seepage water in coal body Flow path.Water can be obtained during by coal body, it is obvious in main crack flow process, it is shunted in the middle part of secondary interstice It is dynamic, it is not flowed substantially in base portion.For water in coal body flow process, speed overall trend is to be gradually reduced.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;Although Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used To modify to technical solution documented by previous embodiment, or some or all of the technical features are equal Replacement;And these are modified or replaceed, model defined by the claims in the present invention that it does not separate the essence of the corresponding technical solution It encloses.

Claims (7)

1. a kind of method for simulating coal-bed flooding flow event, it is characterised in that: the following steps are included:
Step 1 carries out CT scan to the coal sample after test tape pressure injection water, obtains original image grayscale image;
Step 2 carries out median filtering noise reduction process to original image grayscale image with MATLAB, obtains txt formatted file;
Step 3 imports obtained txt formatted file in COMSOL software, with interpolating function by coal body model reconstruction;
Step 4 assigns model correlated variables parameter;
Step 5, setting model initial boundary conditions;
Step 6 divides network, and carries out numerical simulation to injection water coal sample flow event, obtains water in coal body by simulation process The accurate description of middle motion process, wetted perimeter.
2. the method for simulation coal-bed flooding flow event according to claim 1, it is characterised in that: in the step 1, adopt Coal sample is scanned by cone-beam scan and DR real time imagery Through Several Survey Measure with Industrial X-ray CT detection system Picture nondestructively scans the hole fissured structure inside coal body, obtains original image grayscale image.
3. the method for simulation coal-bed flooding flow event according to claim 1, it is characterised in that: in the step 2, adopt The pore structural information in the original image grayscale image of CT scan is identified with the I=imread function in MATLAB, is read One layer of picture of coal is taken, and noise spot in picture is subjected to noise reduction process using median filtering function, and in MATLAB Fid=fopen function extracts each information put on image, exports as txt formatted file.
4. the method for simulation coal-bed flooding flow event according to claim 1, it is characterised in that: right in the step 3 The size that txt formatted file obtained corresponds to picture is modeled, and picture size is 2.88cm × 3.01cm, in txt format Inputted in file under %grid and input two groups of numbers, one group be 1-301, another group be 1-288, it is defeated before obtained data Enter %data, selects interpolating function to be imported in COMSOL global definition.
5. the method for simulation coal-bed flooding flow event according to claim 1, it is characterised in that: in the step 4, Keying in density in COMSOL software under global definition parameter is 1000kg/m3, dynamic viscosity 0.001Pa, material properties selection CoAI。
6. the method for simulation coal-bed flooding flow event according to claim 1, it is characterised in that: in the step 5, object Manage field select coefficient form boundary partial differential equation, dependent variable physical quantity select Head (m), source item physical quantity selection head with Time change;Boundary setting is carried out to obtained geometrical model, the coboundary of coal body is set as the head of 20mm, lower boundary setting For the head of 0mm, the head difference of 20mm is generated, driving fluid flows, and right boundary is impermeable.
7. the method for simulation coal-bed flooding flow event according to claim 1, it is characterised in that: in the step 6, root According to the specification of selected materials, the unit size in crack is refined, the unit size of setting crack position is 5mm, other positions Install 57.2mm;Setting is modeled as steady-flow.
CN201910449651.1A 2019-05-28 2019-05-28 A method of simulation coal-bed flooding flow event Pending CN110162900A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910449651.1A CN110162900A (en) 2019-05-28 2019-05-28 A method of simulation coal-bed flooding flow event

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910449651.1A CN110162900A (en) 2019-05-28 2019-05-28 A method of simulation coal-bed flooding flow event

Publications (1)

Publication Number Publication Date
CN110162900A true CN110162900A (en) 2019-08-23

Family

ID=67629268

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910449651.1A Pending CN110162900A (en) 2019-05-28 2019-05-28 A method of simulation coal-bed flooding flow event

Country Status (1)

Country Link
CN (1) CN110162900A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110702570A (en) * 2019-09-24 2020-01-17 山东科技大学 Method for realizing visualization of coal body pore fracture dynamic seepage process

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101268935A (en) * 2008-04-25 2008-09-24 清华大学 In vivo fluorescence numerator imaging modelling approach capable of calling multiple imaging algorithms
CN105160050A (en) * 2015-07-01 2015-12-16 中国矿业大学(北京) Loaded coal and rock damage constitutive equation building method based on CT scanning
CN106960070A (en) * 2016-12-28 2017-07-18 山东科技大学 A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT
CN109211666A (en) * 2018-08-31 2019-01-15 山东科技大学 The method of coal body permeability under predicted stresses loading environment based on CT scan

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101268935A (en) * 2008-04-25 2008-09-24 清华大学 In vivo fluorescence numerator imaging modelling approach capable of calling multiple imaging algorithms
CN105160050A (en) * 2015-07-01 2015-12-16 中国矿业大学(北京) Loaded coal and rock damage constitutive equation building method based on CT scanning
CN106960070A (en) * 2016-12-28 2017-07-18 山东科技大学 A kind of seepage simulation method that coal body is reconstructed based on finite element discretization member CT
CN109211666A (en) * 2018-08-31 2019-01-15 山东科技大学 The method of coal body permeability under predicted stresses loading environment based on CT scan

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
白若男: ""基于显微CT技术的煤体微观孔隙结构构建及低压水渗流数值模拟", 《中国优秀硕士学位论文全文数据库工程科技I辑》 *
陈俊国: "煤层气储层孔裂隙多尺度渗透率预测和流固藕合模型", 《中国博士学位论文全文数据库工程科技I辑》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110702570A (en) * 2019-09-24 2020-01-17 山东科技大学 Method for realizing visualization of coal body pore fracture dynamic seepage process

Similar Documents

Publication Publication Date Title
US9031820B2 (en) Method for creating a 3D model of a hydrocarbon reservoir, and method for comparative testing of hydrocarbon recovery techniques
Ishutov et al. 3D printing sandstone porosity models
RU2743121C1 (en) Method of numerical modeling of filtration/damage/stress communication during water injection into coal-bearing mass during regionalization
Zhang et al. An extended grain‐based model accounting for microstructures in rock deformation
US9046509B2 (en) Method and system for estimating rock properties from rock samples using digital rock physics imaging
Ji et al. Characterization of pore structure and strain localization in Majella limestone by X-ray computed tomography and digital image correlation
CN107449707A (en) Three-dimensional characterization determination method and device for quantification of pores with different scales in shale reservoir
CN109211666B (en) The method of coal body permeability under predicted stresses loading environment based on CT scan
NO20121394A1 (en) Procedure for obtaining consistent and integrated physical properties of porous media
Chen et al. A new model of pore structure typing based on fractal geometry
Voltolini et al. Visualization and prediction of supercritical CO2 distribution in sandstones during drainage: An in situ synchrotron X-ray micro-computed tomography study
Islam et al. Structural characterization and numerical simulations of flow properties of standard and reservoir carbonate rocks using micro-tomography
Han et al. Numerical computation of elastic properties for porous rocks based on CT-scanned images using direct mapping method
Payton et al. Pore-scale assessment of subsurface carbon storage potential: implications for the UK Geoenergy Observatories project
CN110162900A (en) A method of simulation coal-bed flooding flow event
Kohanpur et al. Using direct numerical simulation of pore-level events to improve pore-network models for prediction of residual trapping of CO2
Mohammadmoradi et al. Pore network and morphological characterization of pore-Level structures
Izgec Reactive flow in vuggy carbonates: methods and models applied to matrix acidizing of carbonates
Mielniczuk et al. Capillary water in 2‐D drying‐Cracking soil sub‐grain models: Morphology and kinematics of evaporation and Haines jumps
Gharbi et al. A relationship between the fractal dimension and scaling groups of unstable miscible displacements
CN103077556B (en) The Three-dimension Numerical Model method for designing of sand production
Ishutov 3D printing porous proxies as a new tool for laboratory and numerical analyses of sedimentary rocks
Kameda et al. To see a rock in a grain of sand
CN114818236A (en) Seepage simulation method for reconstructing coal body pore crack network structure
Khan et al. Capillary rise in vuggy media

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190823