CN116541904A - Gas extraction drilling design method - Google Patents

Gas extraction drilling design method Download PDF

Info

Publication number
CN116541904A
CN116541904A CN202310314588.7A CN202310314588A CN116541904A CN 116541904 A CN116541904 A CN 116541904A CN 202310314588 A CN202310314588 A CN 202310314588A CN 116541904 A CN116541904 A CN 116541904A
Authority
CN
China
Prior art keywords
drilling
hole
borehole
holes
gas extraction
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
CN202310314588.7A
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.)
Guoneng Wuhai Energy Huangbaici Mining Co ltd
CHN Energy Wuhai Energy Co Ltd
CCTEG Chongqing Research Institute Co Ltd
Original Assignee
Guoneng Wuhai Energy Huangbaici Mining Co ltd
CHN Energy Wuhai Energy Co Ltd
CCTEG Chongqing Research Institute Co Ltd
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 Guoneng Wuhai Energy Huangbaici Mining Co ltd, CHN Energy Wuhai Energy Co Ltd, CCTEG Chongqing Research Institute Co Ltd filed Critical Guoneng Wuhai Energy Huangbaici Mining Co ltd
Priority to CN202310314588.7A priority Critical patent/CN116541904A/en
Publication of CN116541904A publication Critical patent/CN116541904A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21FSAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F7/00Methods or devices for drawing- off gases with or without subsequent use of the gas for any purpose
    • 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/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/043Optimisation of two dimensional placement, e.g. cutting of clothes or wood
    • 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
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/10Geothermal energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Geometry (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Agronomy & Crop Science (AREA)
  • Game Theory and Decision Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Development Economics (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Computational Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a gas extraction drilling design method. The gas extraction drilling design method comprises the following steps: establishing a bedding drilling model: designing the drill holes of the sections of each coal seam and forming a drill hole group; establishing a layer-penetrating drilling model: arranging drilling holes on the plane and/or the section, determining two-dimensional coordinates of the drilling holes, and acquiring at least one of an inclination angle, an azimuth angle and a hole depth of the drilling holes; and establishing a drilling analysis effect model, wherein the drilling analysis effect model performs drilling arrangement effect analysis through a drilling arrangement effect judging method based on space linear density analysis and/or a drilling arrangement effect judging method based on space nuclear density analysis. The invention effectively solves the problem of poor universality of gas extraction drilling design in the prior art.

Description

Gas extraction drilling design method
Technical Field
The invention relates to the technical field of gas extraction, in particular to a gas extraction drilling design method.
Background
At present, gas disasters are always one of main disasters of coal mines, wherein one of main technical means and methods for gas management is gas extraction. Specifically, the design of the extraction drilling holes of the coal mines in China is mainly designed by extraction department designers through experience or by using prepared design templates.
However, the above-mentioned extraction borehole design is only suitable for the occurrence of stable near-horizontal coal seam with small occurrence variation and the extraction borehole with low requirement on borehole design, but the adaptation is poor for the occurrence of unstable coal seam with large inclination variation and the extraction borehole design of outstanding coal seam with strict requirement on borehole control range, so that the borehole design which strictly accords with regulations and has no extraction blank zone is difficult to realize.
Disclosure of Invention
The invention mainly aims to provide a gas extraction drilling design method, which aims to solve the problem of poor universality of gas extraction drilling design in the prior art.
In order to achieve the above object, the present invention provides a gas extraction borehole design method, including: establishing a bedding drilling model: designing the drill holes of the sections of each coal seam and forming a drill hole group; establishing a layer-penetrating drilling model: arranging drilling holes on the plane and/or the section, determining two-dimensional coordinates of the drilling holes, and acquiring at least one of an inclination angle, an azimuth angle and a hole depth of the drilling holes; and establishing a drilling analysis effect model, wherein the drilling analysis effect model performs drilling arrangement effect analysis through a drilling arrangement effect judging method based on space linear density analysis and/or a drilling arrangement effect judging method based on space nuclear density analysis.
Further, the method for designing the drill holes of the sections of each coal seam comprises a divergent drill hole design, wherein the divergent drill hole design comprises a fixed drill hole number design method, and the fixed drill hole number design method comprises the following steps: step S11: performing drilling hole arrangement, and obtaining hole coordinates K (x, y, z) of each drilling hole; step S12: determining the drilling interval, and acquiring the final Kong Zuobiao D (X, Y, Z) of each drilling hole; step S13: at least one of the inclination angle, azimuth angle and hole depth of the borehole is obtained from the hole coordinates K (X, Y, Z) and the final Kong Zuobiao D (X, Y, Z).
Further, in step S11, the hole distribution area is formed after the two sides of the area to be drilled are removed, and the method for opening the hole distribution area in the horizontal direction includes: reserving a region for placing a drilling machine at the two sides, and uniformly drilling holes in the hole distribution region.
Further, in step S12, the method of determining the borehole spacing includes: step S51: optionally selecting a first preset value as a drilling interval; step S52: hole distribution is carried out according to a vertical distance type drilling interval calculation mode; step S53: judging the drilling condition, if the drilling can not completely cover the control area when the number of the drilling holes reaches the set number, selecting a second preset value which is larger than the first preset value as the drilling hole interval, and repeating the step S52; if the control area is completely covered when the set number of the drill holes is not reached, selecting a third preset value smaller than the first preset value as the drill hole interval, and repeating the step S52; the method for acquiring the final hole coordinates of each drilled hole comprises the following steps: after the drill hole distance is determined, the coordinates of the boundary intersection point of each drill hole and the control area are final hole coordinates.
Further, in step S13, the method for obtaining at least one of the inclination angle θ, the azimuth angle α and the hole depth L of the drill hole according to the hole coordinates K (X, Y, Z) and the final Kong Zuobiao D (X, Y, Z) includes: the calculation formula of the hole depth L is as follows:
the calculation formula of the inclination angle theta of the drilling hole is as follows:
the calculation formula of the azimuth angle alpha of the drilling hole is as follows: />
Further, the divergent borehole design also includes a fixed borehole spacing design method, the fixed borehole spacing design method including: step S21: drilling arrangement is carried out at preset drilling intervals; step S22: and (5) hole distribution is carried out according to a hole distribution mode with the fixed number of holes to be drilled, and each drilling parameter is obtained.
Further, the method for designing the drilling holes of the sections of each coal seam further comprises parallel drilling hole design.
Further, before the drilling arrangement is performed on the plane and/or the section plane and the two-dimensional coordinates of the drilling are determined, establishing the layer-penetrating drilling model further comprises: and determining a drilling hole bottom and a coal point boundary corresponding to the drilling control range boundary, and taking smaller values in the drilling hole bottom and the coal point boundary as drilling control ranges.
Further, a method of arranging a borehole in a plane and/or a section and determining two-dimensional coordinates of the borehole comprises: the vertical direction drilling arrangement method comprises the following steps: acquiring coordinate values K (0, X, Z) of the opening point of each drilling hole in the front-back direction and the vertical direction, and coordinate values D (X, X, Z) of the coal outlet point in the front-back direction and the vertical direction; the horizontal drilling arrangement method comprises the following steps: acquiring coordinate values K '(/ Y) of the opening point of each drilling hole in the left-right direction and coordinate values D' (/ Y) of the coal outlet point in the left-right direction; the coordinates of the hole opening points of each drilling hole in the front-back direction, the vertical direction and the left-right direction are integrated to obtain a three-dimensional coordinate K "(0, Y, Z), and the coordinates of the coal outlet points of each drilling hole in the front-back direction, the vertical direction and the left-right direction are integrated to obtain a three-dimensional coordinate D" (X, Y, Z).
Further, the method for analyzing the borehole placement effect based on the spatial linear density analysis comprises the following steps: step S31: dividing the grid cells; step S32: determining a search radius r; step S33: taking the searching radius r as a preset circle, and calculating the length of each line segment falling into the preset circle to obtain the importance of the line segment; step S34: calculating the drilling line density of the grid unit; step S35: and displaying the judging result in a cloud picture form, and judging the blank zone of the drilling hole.
Further, the search radius R and the effective extraction radius R of the drill hole satisfy the following conditions: r is more than or equal to 2 and less than or equal to 3R; in step S34, the calculation method of the grid cell borehole line density includes: multiplying the length of each line segment by the importance of the drilling hole, summing the multiplied values, and dividing the summed values by the area of a preset circle; the importance of the drill hole can be calculated according to the effective extraction radius R of the drill hole.
Further, the method for analyzing the borehole placement effect based on the method for evaluating the borehole placement effect based on the spatial kernel density analysis comprises the following steps: step S41: determining a kernel function; step S42: determining the influence range of the line segment; step S43: acquiring the importance degree of the line segment falling in the influence range of the line segment; step S44: calculating the density of the drilling nuclei of the grid cells; step S45: and displaying the judging result in a cloud picture form, and judging the blank bands.
Further, in step S44, the method for calculating the density of the mesh unit drilling nuclei includes: the distance weights of all line segments falling within the range of influence of the line segment are summed.
When the technical scheme of the invention is applied, when the gas extraction drilling hole design is needed, the drilling holes on the sections of each coal seam are designed and a drilling group is formed to establish a bedding drilling model, then the drilling holes are arranged on the plane and/or the section and the two-dimensional coordinates of the drilling holes are determined, at least one of the inclination angle, the azimuth angle and the hole depth of the drilling holes is obtained to establish a bedding drilling hole model, finally the drilling hole arrangement effect analysis is carried out by a drilling hole arrangement effect judging method based on space linear density analysis and/or a drilling hole arrangement effect judging method based on space nuclear density analysis, so that an auxiliary means is provided for the quality supervision of gas extraction drilling hole engineering, the accurate grasp of the mine gas extraction drilling hole design and execution conditions is realized, the problem that the universality of the gas extraction drilling hole design in the prior art is poor is solved, and the intelligent design of the gas extraction drilling hole is realized. Meanwhile, the intelligent design of the gas extraction drilling hole is carried out by utilizing a computer, the intelligent design method has the characteristics of wide application range and rapid and convenient design, the designed extraction drilling hole can meet the requirements of relevant control ranges and construction positions such as anti-outburst rules, coal mine safety regulations and the like, the generation of a drilling blank zone is effectively avoided, in addition, the designed drilling hole control range, uniformity and the like can be evaluated, the optimal adjustment of the drilling hole design by a designer according to the evaluation result is facilitated, and finally, the optimal extraction drilling hole design which meets the requirements of relevant regulations and enables the construction engineering quantity to be minimum is obtained.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention. In the drawings:
FIG. 1 is a diagram showing a comparison of two opening patterns, namely a fixed two-wall distance and a fixed opening pitch, according to a first embodiment of a gas extraction borehole design method of the present invention;
FIG. 2 shows a schematic borehole diagram of an trial-and-error layout in the gas drainage borehole design method of FIG. 1;
FIG. 3 is a schematic diagram showing the coordinates of the calculation of drilling parameters in the gas extraction drilling design method of FIG. 1;
FIG. 4 is a schematic diagram of a hole bottom and a coal point boundary corresponding to a boundary of a drilling control range in the gas extraction drilling design method in FIG. 1;
FIG. 5 shows a schematic diagram of a linear density analysis method of the gas drainage borehole design method of FIG. 1;
FIG. 6 is a graph showing the relationship between the gas content of the extracted coal and the borehole distance in the gas extraction borehole design method in FIG. 1;
FIG. 7 is a schematic diagram of a nuclear density analysis method in a second embodiment of a gas extraction borehole design method according to the present invention;
FIG. 8 shows a graph of the gas content of the extracted coal body versus the borehole distance in the gas extraction borehole design method of FIG. 7.
Detailed Description
It should be noted that, in the case of no conflict, the embodiments and features in the embodiments may be combined with each other. The invention will be described in detail below with reference to the drawings in connection with embodiments.
It is noted that all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs unless otherwise indicated.
In the present invention, unless otherwise indicated, terms of orientation such as "upper" and "lower" are used generally with respect to the orientation shown in the drawings or to the vertical, vertical or gravitational orientation; also, for ease of understanding and description, "left, right" is generally directed to the left, right as shown in the drawings; "inner and outer" refer to inner and outer relative to the outline of the components themselves, but the above-described orientation terms are not intended to limit the present invention.
In order to solve the problem of poor universality of gas extraction drilling design in the prior art, the application provides a gas extraction drilling design method.
Example 1
As shown in fig. 1 to 6, the gas extraction borehole design method includes:
establishing a bedding drilling model: designing the drill holes of the sections of each coal seam and forming a drill hole group;
establishing a layer-penetrating drilling model: arranging drilling holes on the plane and/or the section, determining two-dimensional coordinates of the drilling holes, and acquiring at least one of an inclination angle, an azimuth angle and a hole depth of the drilling holes;
and establishing a drilling analysis effect model, wherein the drilling analysis effect model performs drilling arrangement effect analysis through a drilling arrangement effect judging method based on space linear density analysis.
By applying the technical scheme of the embodiment, when the gas extraction drilling hole design is needed, the drilling holes on the sections of each coal seam are designed and a drilling group is formed to establish a bedding drilling model, then drilling holes are arranged on a plane and/or a section and two-dimensional coordinates of the drilling holes are determined, at least one of the inclination angle, the azimuth angle and the hole depth of the drilling holes is acquired to establish a bedding drilling hole model, finally the drilling hole arrangement effect analysis is carried out by a drilling hole arrangement effect judging method based on space linear density analysis and/or a drilling hole arrangement effect judging method based on space nuclear density analysis, so that an auxiliary means is provided for quality supervision of gas extraction drilling holes, the accurate grasp of mine gas extraction measure design and execution conditions is realized, the problem that the universality of the gas extraction drilling hole design in the prior art is poor is solved, and the intelligent design of the gas extraction drilling holes is realized. Meanwhile, the intelligent design of the gas extraction drilling hole is carried out by utilizing a computer, the intelligent design method has the characteristics of wide application range and rapid and convenient design, the designed extraction drilling hole can meet the requirements of relevant control ranges and construction positions such as anti-outburst rules, coal mine safety regulations and the like, the generation of a drilling blank zone is effectively avoided, in addition, the designed drilling hole control range, uniformity and the like can be evaluated, the optimal adjustment of the drilling hole design by a designer according to the evaluation result is facilitated, and finally, the optimal extraction drilling hole design which meets the requirements of relevant regulations and enables the construction engineering quantity to be minimum is obtained.
In the process of establishing the bedding drilling model, firstly, drilling holes of all sections are respectively designed, then, drilling holes of a plurality of sections are combined according to requirements, and finally, the actual design meeting the field requirements is formed. Specifically, as shown in fig. 1, the double-roadway tunneling bedding drilling design can be formed by combining five units of left-side drill site divergent drilling, left-side roadway head-on parallel drilling, transverse parallel drilling, right-side roadway head-on parallel drilling and right-side drill site divergent drilling. Wherein the borehole designs include divergent borehole designs and parallel borehole designs; the divergent drilling design method comprises a fixed drilling number design method and a fixed drilling interval design method.
In the embodiment, the method for judging the drilling arrangement effect based on the space linear density analysis expands the search area, also considers the difference of the line segment importance in the search area, and can realize the automatic evaluation of the uniformity degree of the gas extraction drilling arrangement. Therefore, a user can control and evaluate the design drilling holes, the control range of construction drilling holes, the number of drilling holes in each area and the uniformity degree of the large-range working surface by using the space linear density judging method, so that the occurrence of blank zones of the drilling holes can be effectively avoided, and the drilling holes in the extraction range are uniformly distributed in the coal seam. In addition, whether the construction drilling hole meets the design requirement or not can be timely found, and whether the constructed drilling hole meets the related regulation requirement or not.
In this embodiment, the method for designing the drill holes of the section of each coal seam includes a divergent drill hole design, the divergent drill hole design includes a fixed drill hole number design method, and the fixed drill hole number design method includes:
step S11: performing drilling hole arrangement, and obtaining hole coordinates K (x, y, z) of each drilling hole;
step S12: determining the drilling interval, and acquiring the final Kong Zuobiao D (X, Y, Z) of each drilling hole;
step S13: at least one of the inclination angle, azimuth angle and hole depth of the borehole is obtained from the hole coordinates K (X, Y, Z) and the final Kong Zuobiao D (X, Y, Z).
Specifically, the design method of the fixed drilling intervals can cover the required control range with the minimum drilling quantity, ensure the minimum drilling engineering quantity, meet the requirements of the control ranges of anti-outburst rules, coal mine safety regulations and the like, and can also enable the drilling construction quantity to be minimum. Under the condition of meeting the design requirement, the construction amount of underground drilling can be greatly saved, the coal exploitation cost is effectively reduced, and the coal mine safety and economic benefit are increased.
In this embodiment, in step S11, a hole distribution area is formed after removing the two sides of the area to be drilled, and the method for opening the hole distribution area in the horizontal direction includes:
reserving a region for placing a drilling machine at the two sides, and uniformly drilling holes in the hole distribution region.
In this embodiment, in step S12, the method for determining the borehole pitch includes:
step S51: optionally selecting a first preset value as a drilling interval;
step S52: hole distribution is carried out according to a vertical distance type drilling interval calculation mode;
step S53: judging the drilling condition, if the drilling can not completely cover the control area when the number of the drilling holes reaches the set number, selecting a second preset value which is larger than the first preset value as the drilling hole interval, and repeating the step S52; if the control area is completely covered when the set number of the drill holes is not reached, selecting a third preset value smaller than the first preset value as the drill hole interval, and repeating the step S52;
specifically, the specific method flow of the intelligent design of the fixed drilling number mode bedding drilling is as follows:
first, a borehole aperture arrangement is performed, and aperture coordinates K (x, y, z) of the respective boreholes are calculated. The horizontal (row) direction openings comprise two modes of fixed two-side distance and fixed opening spacing. When drilling holes are usually arranged, a certain space is reserved in the horizontal direction at two sides, so that a drilling machine can be placed in the drilling process, and the drilling holes can be conveniently operated, and the other middle area allows the drilling holes to be arranged, which is called a hole distribution area temporarily. Under the mode of fixing the distance between two sides and hole distribution, the distance from the boundary drilling holes to the two sides of the roadway is fixed and unchanged, and is equal to the minimum distance reserved for placing the drilling machine, the drilling holes are evenly distributed in the hole distribution area, and the maximization of the drilling hole spacing is ensured. Under the mode of hole distribution with fixed hole spacing, the drill holes are uniformly distributed at the middle position of the hole distribution area according to the given hole spacing, so that the maximization of the distance from the boundary drill hole to the roadway wall is ensured. Two hole pattern pairs are shown in fig. 2. The vertical (column) holes are evenly distributed in the height range of the holes, so that the maximization of the hole spacing is ensured.
The reasonable borehole spacing is then determined and the final Kong Zuobiao D (X, Y, Z) for each borehole is calculated. The drilling interval is measured in a vertical distance mode, the distance between adjacent drilling end hole points is not fixed in a rectangular drilling control area, and the coordinates of the drilling end hole points cannot be directly calculated, so that a hole distribution method based on an attempt approximation principle is designed, and the method specifically comprises the following steps: firstly, taking any value as a drilling interval, and carrying out hole distribution according to a vertical distance type drilling interval calculation mode. When the given fixed number of holes is reached, if the holes cannot fully cover the control area, the hole spacing is too small, and the hole spacing should be increased for rearrangement (shown in fig. 2 (a)); conversely, it is indicated that the hole pitch is too large, and the hole pitch should be reduced for rearrangement (as shown in fig. 2 (b)). According to this method, the drill hole spacing is adjusted and the holes are rearranged until a fixed number of drill holes just cover the control area, at this time, the drill hole spacing is a reasonable drill hole spacing, and the coordinates of the intersection point of each corresponding drill hole and the boundary of the control range are the final Kong Zuobiao of the drill holes (as shown in fig. 2 (c)).
In this embodiment, in step S12, the method for obtaining the final hole coordinates of each borehole includes:
after the drill hole distance is determined, the coordinates of the boundary intersection point of each drill hole and the control area are final hole coordinates.
In this embodiment, in step S13, the method for obtaining at least one of the inclination angle θ, the azimuth angle α, and the hole depth L of the drill hole according to the hole coordinates K (X, Y, Z) and the final Kong Zuobiao D (X, Y, Z) includes:
the calculation formula of the hole depth L is as follows:
the calculation formula of the inclination angle theta of the drilling hole is as follows:
the calculation formula of the azimuth angle alpha of the drilling hole is as follows:
specifically, the specific method flow of the intelligent design of the fixed-borehole interval mode bedding borehole is as follows:
firstly, arranging holes at a hole spacing given by a user until a control range of the designed holes completely covers a range of expected required control, wherein the number of the holes is the minimum, then, hole arrangement is carried out according to a hole arrangement mode with the fixed number of the holes, and finally, parameters of all the holes are obtained.
In this embodiment, the divergent borehole design further includes a fixed borehole spacing design method, the fixed borehole spacing design method including:
step S21: drilling arrangement is carried out at preset drilling intervals;
step S22: and (5) hole distribution is carried out according to a hole distribution mode with the fixed number of holes to be drilled, and each drilling parameter is obtained.
In this embodiment, the method for designing the drill holes of the sections of each coal seam further includes parallel drill hole design.
In this embodiment, before the drilling arrangement is performed on the plane and/or the section plane and the two-dimensional coordinates of the drilling are determined, the establishing of the through-layer drilling model further includes:
and determining a drilling hole bottom and a coal point boundary corresponding to the drilling control range boundary, and taking smaller values in the drilling hole bottom and the coal point boundary as drilling control ranges.
Specifically, the control range of the through-layer drilling is provided with two calculation modes of ' hole bottom ' and ' coal point. For a certain set of holes, the control range of the hole bottom is different from that of the coal point, and even the control range of the coal point is quite different, which is mainly related to the thickness of the coal layer, the included angle between the holes and the coal layer, and the like. In practical application of coal mines, the smaller of the two is usually taken as the control range of the group of drilling holes, but the two are required to be converted in the computer intelligent design of the drilling holes, and the two are uniformly calculated according to one mode. For different spatial position relations of coal seam-tunnel, a calculation model of the hole bottoms and coal-seeing point boundaries of the front, the rear, the left side and the right side corresponding to the boundary of the drilling control range is established, as shown in fig. 4.
In this embodiment, a method of arranging a borehole in a plane and/or a section and determining two-dimensional coordinates of the borehole includes:
the vertical direction drilling arrangement method comprises the following steps: acquiring coordinate values K (0, X, Z) of the opening point of each drilling hole in the front-back direction and the vertical direction, and coordinate values D (X, X, Z) of the coal outlet point in the front-back direction and the vertical direction;
the horizontal drilling arrangement method comprises the following steps: acquiring coordinate values K '(/ Y) of the opening point of each drilling hole in the left-right direction and coordinate values D' (/ Y) of the coal outlet point in the left-right direction;
the coordinates of the hole opening points of each drilling hole in the front-back direction, the vertical direction and the left-right direction are integrated to obtain a three-dimensional coordinate K "(0, Y, Z), and the coordinates of the coal outlet points of each drilling hole in the front-back direction, the vertical direction and the left-right direction are integrated to obtain a three-dimensional coordinate D" (X, Y, Z).
Specifically, in the vertical drilling arrangement, because the measurement of the drilling interval in the sectional view adopts a vertical distance mode, the drilling arrangement is carried out by applying the same trial approach hole distribution method as the sequential drilling, so that the control range of the front-back direction is covered by the minimum drilling column number under the condition that the projection interval of the drilling section is not larger than the given drilling interval. From the cross-sectional arrangement of the drilled holes, coordinate values of the front-rear and vertical directions and the like of each row of drilled hole opening points, coal outlet points and the like, namely, K (0, X, Z) and D (X, Z) are obtained.
In a horizontal drilling arrangement, the holes are distributed in the left-right direction using a conventional hole bottom spacing metering mode. Similar to the sectional view, the drilling arrangement simultaneously satisfies the conditions of full coverage of the left-right direction control range, no more than a given drilling pitch at the bottom of the holes, and a minimum of 3 drilling columns. From the floor plan, coordinate values K '(×y,) and D' (×y,) of the opening point, the coal outlet point, and the like in the right-left direction can be obtained. The coordinates obtained by hole distribution in the horizontal direction and the vertical direction are integrated, so that three-dimensional coordinates of a hole opening point and a coal outlet point of each drilling hole, namely K "(0, Y, Z) and D" (X, Y, Z) can be obtained.
Specifically, the process of calculating the inclination angle theta, the azimuth angle alpha and the hole depth L of the drilling hole based on the space geometry algorithm is as follows:
the calculation method of the inclination angle theta and the azimuth angle alpha of the drilling hole is the same as that of the bedding drilling hole, and the deflection angle, the inclination angle and the hole depth of the drilling hole are rapidly calculated by using a mature space geometrical algorithm according to the relative coordinates K (X, Y, Z) of the drilling hole opening points and the coordinates D (X, Y, Z) of the drilling hole coal outlet points. It should be noted that, the actual through-layer drilling hole still continues to construct a section of length after passing through the coal seam, on one hand, the inspection drilling hole has indeed penetrated the full thickness of the coal seam, on the other hand, a section of drilling hole is reserved downwards to the hole for storing deposited slag (water) so as to prevent the coal hole section from being blocked and the extraction effect from being affected. Therefore, the calculation of the hole depth L requires adding the length of the drilled hole after passing through the coal seam on the basis of the space distance between the hole opening and the coal outlet point.
In this embodiment, the method for performing the borehole placement effect analysis based on the borehole placement effect evaluation method of the spatial linear density analysis includes:
step S31: dividing the grid cells;
step S32: determining a search radius r;
step S33: taking the searching radius r as a preset circle, and calculating the length of each line segment falling into the preset circle to obtain the importance of the line segment;
step S34: calculating the drilling line density of the grid unit;
step S35: and displaying the judging result in a cloud picture form, and judging the blank zone of the drilling hole.
Specifically, the borehole arrangement effect judging method based on space linear density analysis comprises the following specific implementation steps: the search radius r is used for drawing a preset circle by taking the center of each grid pixel as the circle center, the length of each line segment falling into the preset circle is calculated, if the importance degree given by a user to each line segment is considered, the length of each line segment is multiplied by the importance degree and summed, and then the obtained sum is divided by the area of the preset circle (search circle) to obtain the linear density of the grid unit, as shown in fig. 5. In addition, the search area may be triangular, rectangular, hexagonal, or the like.
In this embodiment, the specific method for evaluating the borehole placement effect based on the linear density analysis is as follows:
(1) dividing the grid cells; the size of the grid cells determines the final patterning accuracy, and the larger the grid cells are, the lower the patterning accuracy is; the smaller the grid cell, the higher the patterning accuracy, but the larger the corresponding calculation amount. The size of the grid cells takes two modes, absolute size and relative size. Absolute sizes such as 0.5m×0.5m, 1m×1m, etc., relative sizes such as 1/5000, 1/10000, etc. of the judgment range;
(2) determining a search radius r; the searching radius r has an important influence on the judging result, if the searching range is too large, the difference of drilling density calculation results among the areas can be weakened, and the comparison is not obvious; if the search range is too small, the calculation result shows no regional difference. Determining the effective extraction radius of the drill hole with the reasonable search radius r being 2-3 times through experiments;
(3) determining the importance of the line segment; the influence of different types of holes on the gas extraction of the coal body in the search area is reflected in the gas extraction holes, for example, the gas extraction effects of the directional holes with equal length and the common holes are obviously different, and the gas extraction effects of the holes with different diameters are also different, so that the holes with different types and different diameters should be given with different importance in the evaluation process. The importance of the drilling hole can be calculated according to the effective extraction radius of the drilling hole;
(4) calculating the drilling line density of the grid cells; calculating the linear density of each grid unit according to a linear density calculation method;
(5) and displaying the judgment result and judging the blank zone, wherein the judgment result is finally displayed in a cloud picture form. In order to show the existence of the blank bands of the drill holes in the patterns, and simultaneously facilitate the comparison between different patterns, a unified blank band judgment standard is set.
In this embodiment, the search radius R and the effective extraction radius R of the borehole satisfy: r is more than or equal to 2 and less than or equal to 3R; in step S34, the calculation method of the grid cell borehole line density includes:
multiplying the length of each line segment by the importance of the drilling hole, summing the multiplied values, and dividing the summed values by the area of a preset circle;
the importance of the drill hole can be calculated according to the effective extraction radius R of the drill hole.
Example two
The difference between the gas extraction drilling design method in the second embodiment and the first embodiment is that: the manner in which the borehole analysis effect model is built is different.
In this embodiment, the borehole placement effect evaluation method based on the spatial kernel density analysis performs the borehole placement effect analysis.
In other embodiments not shown in the drawings, the borehole analysis effect model performs the borehole placement effect analysis by a borehole placement effect evaluation method based on spatial line density analysis and a borehole placement effect evaluation method based on spatial kernel density analysis.
In this embodiment, the method for analyzing the borehole placement effect based on the method for evaluating the borehole placement effect of the spatial kernel density analysis includes:
step S41: determining a kernel function;
step S42: determining the influence range of the line segment;
step S43: acquiring the importance degree of the line segment falling in the influence range of the line segment;
step S44: calculating the density of the drilling nuclei of the grid cells;
step S45: and displaying the judging result in a cloud picture form, and judging the blank bands.
Specifically, the borehole placement effect evaluation method based on the spatial kernel density analysis is an improvement of the line density analysis method, and considers the influence of the distance between the grid unit and the line segment on the basis of the line density analysis method, and changes the length of the line segment falling into the search area into the distance weight, and the kernel density of the grid unit is the sum of the distance weights of all the line segments contributing to the same, as shown in fig. 6.
In this embodiment, the specific method for evaluating the borehole placement effect based on the nuclear density analysis is as follows:
(1) determining a kernel function; the kernel function is used to describe the change in distance weight coefficient with distance size. With the kernel density function, the distance weight of 1 total line segments can be dispersed from the center position to the outside, wherein the highest value is given at the center, and gradually decreases to zero in the outside influence range. Common kernel functions include a uniform kernel, a triangle kernel, a quadratic kernel, a quartic kernel, a sextuse kernel, a gaussian kernel, a cosine kernel, and the like, and the kernel functions must be reasonably selected. The kernel function represents the influence of the distance from the grid unit coal body to the drilling hole on the extraction effect in the gas extraction judgment. According to the relation between the gas reduction amount of the gas extraction coal body and the drilling distance (as shown in fig. 7 and 8), a fourth-order kernel function is selected in the embodiment;
(2) determining the line segment influence range; the line section influence range is embodied in the gas extraction judgment and is the influence range of the extraction drilling. The drilling influence range is mainly related to the type of drilling and the diameter of the drilling under the condition that factors such as the negative extraction pressure and the extraction time are not considered. Therefore, in the judging process, different types of gas extraction drilling holes with different diameters are required to be endowed with different influence ranges, namely the effective influence ranges of various drilling holes;
(3) determining the importance of the line segment; in the gas extraction drilling holes, the importance of the line segment reflects the influence of different types of drilling holes on the gas extraction of the coal body in the search area, and the value is taken according to the importance of the line segment in the linear density analysis;
(4) calculating the nuclear density of the grid cells; the density corresponding to each grid unit is equal to the sum of the products of all nuclear surface values superimposed on the grid unit and the importance of the corresponding line segments;
(5) displaying a judgment result and judging a blank zone; and finally, displaying the judgment result in a cloud picture form. According to the nuclear density analysis method, the nuclear density calculation result of the drilling blank zone is 0, and the calculation result of the drilling control area is greater than 0.
In this embodiment, in step S44, the method for calculating the density of the mesh unit drilling nuclei includes:
the distance weights of all line segments falling within the range of influence of the line segment are summed.
From the above description, it can be seen that the above embodiments of the present invention achieve the following technical effects:
when the gas extraction drilling hole design is needed, firstly, drilling holes on the sections of each coal seam are designed and a drilling group is formed to establish a bedding drilling hole model, then drilling holes are arranged on a plane and/or a section and two-dimensional coordinates of the drilling holes are determined, at least one of the inclination angle, the azimuth angle and the hole depth of the drilling holes is obtained to establish a bedding drilling hole model, finally, the drilling hole arrangement effect analysis is carried out through a drilling hole arrangement effect judging method based on space linear density analysis and/or a drilling hole arrangement effect judging method based on space nuclear density analysis, so that an auxiliary means is provided for gas extraction drilling hole engineering quality supervision, the problem that the mine gas extraction measure design and the execution condition are poor in universality is solved, and the intelligent design of the gas extraction drilling holes is realized. Meanwhile, the intelligent design of the gas extraction drilling hole is carried out by utilizing a computer, the intelligent design method has the characteristics of wide application range and rapid and convenient design, the designed extraction drilling hole can meet the requirements of relevant control ranges and construction positions such as anti-outburst rules, coal mine safety regulations and the like, the generation of a drilling blank zone is effectively avoided, in addition, the designed drilling hole control range, uniformity and the like can be evaluated, the optimal adjustment of the drilling hole design by a designer according to the evaluation result is facilitated, and finally, the optimal extraction drilling hole design which meets the requirements of relevant regulations and enables the construction engineering quantity to be minimum is obtained.
It will be apparent that the embodiments described above are merely some, but not all, embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present invention without making any inventive effort, shall fall within the scope of the present invention.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of example embodiments in accordance with the present application. As used herein, the singular is also intended to include the plural unless the context clearly indicates otherwise, and furthermore, it is to be understood that the terms "comprises" and/or "comprising" when used in this specification are taken to specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
It should be noted that the terms "first," "second," and the like in the description and claims of the present application and the above figures are used for distinguishing between similar objects and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used may be interchanged where appropriate such that embodiments of the present application described herein may be implemented in sequences other than those illustrated or described herein.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, but various modifications and variations can be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (13)

1. The gas extraction drilling design method is characterized by comprising the following steps of:
establishing a bedding drilling model: designing the drill holes of the sections of each coal seam and forming a drill hole group;
establishing a layer-penetrating drilling model: arranging drilling holes on a plane and/or a section, determining two-dimensional coordinates of the drilling holes, and acquiring at least one of an inclination angle, an azimuth angle and a hole depth of the drilling holes;
and establishing a drilling analysis effect model, wherein the drilling analysis effect model performs drilling arrangement effect analysis through a drilling arrangement effect judging method based on space linear density analysis and/or a drilling arrangement effect judging method based on space nuclear density analysis.
2. The gas drainage borehole design method according to claim 1, wherein the method of designing the boreholes of each coal seam section comprises a divergent borehole design, the divergent borehole design comprising a fixed borehole number design method, the fixed borehole number design method comprising:
step S11: performing drilling hole arrangement, and obtaining hole coordinates K (x, y, z) of each drilling hole;
step S12: determining the drilling interval, and acquiring the final Kong Zuobiao D (X, Y, Z) of each drilling hole;
step S13: at least one of the inclination angle, azimuth angle, and hole depth of the borehole is obtained from the hole coordinates K (X, Y, Z) and the terminal Kong Zuobiao D (X, Y, Z).
3. The gas extraction drilling design method according to claim 2, wherein in the step S11, a hole distribution area is formed after removing two sides of the area to be drilled, and the hole distribution area is perforated in a horizontal direction, which includes:
reserving a region for placing a drilling machine at the two sides, and uniformly drilling holes in the hole distribution region.
4. The gas extraction borehole design method according to claim 2, wherein in the step S12, the method of determining the borehole spacing comprises:
step S51: optionally selecting a first preset value as a drilling interval;
step S52: hole distribution is carried out according to a vertical distance type drilling interval calculation mode;
step S53: judging the drilling condition, if the drilling can not completely cover the control area when the number of the drilling holes reaches the set number, selecting a second preset value which is larger than the first preset value as a drilling interval, and repeating the step S52; if the number of the holes does not reach the set number of the holes, the control area is completely covered, a third preset value smaller than the first preset value is selected as a hole interval, and the step S52 is repeated;
the method for acquiring the final hole coordinates of each drilled hole comprises the following steps:
and after the drill hole distance is determined, the coordinates of the boundary intersection point of each drill hole and the control area are the final hole coordinates.
5. The gas extraction borehole design method according to claim 2, wherein in the step S13, the method of acquiring at least one of the inclination angle θ, the azimuth angle α, and the hole depth L of the borehole from the hole coordinates K (X, Y, Z) and the terminal Kong Zuobiao D (X, Y, Z) comprises:
the calculation formula of the hole depth L is as follows:
the calculation formula of the inclination angle theta of the drilling hole is as follows:
the calculation formula of the azimuth angle alpha of the drilling hole is as follows:
6. the gas extraction borehole design method of claim 2, wherein the divergent borehole design further comprises a fixed borehole spacing design method comprising:
step S21: drilling arrangement is carried out at preset drilling intervals;
step S22: and (5) hole distribution is carried out according to a hole distribution mode with the fixed number of holes to be drilled, and each drilling parameter is obtained.
7. The gas drainage borehole design method according to claim 1, wherein the method of designing the borehole of each coal seam section further comprises parallel borehole design.
8. The gas extraction borehole design method of claim 1, wherein prior to borehole placement in the plane and/or section and determination of the two-dimensional coordinates of the borehole, the establishing a through-the-horizon borehole model further comprises:
and determining a drilling hole bottom and a coal point boundary corresponding to the drilling control range boundary, and taking smaller values in the drilling hole bottom and the coal point boundary as drilling control ranges.
9. The gas extraction borehole design method according to claim 1, wherein the method of arranging the borehole in a plane and/or a section and determining the two-dimensional coordinates of the borehole comprises:
the vertical direction drilling arrangement method comprises the following steps: acquiring coordinate values K (0, X, Z) of the opening point of each drilling hole in the front-back direction and the vertical direction, and coordinate values D (X, X, Z) of the coal outlet point in the front-back direction and the vertical direction;
the horizontal drilling arrangement method comprises the following steps: acquiring coordinate values K '(/ Y) of the opening point of each drilling hole in the left-right direction and coordinate values D' (/ Y) of the coal outlet point in the left-right direction;
the coordinates of the hole opening points of each drilling hole in the front-back direction, the vertical direction and the left-right direction are integrated to obtain a three-dimensional coordinate K "(0, Y, Z), and the coordinates of the coal outlet points of each drilling hole in the front-back direction, the vertical direction and the left-right direction are integrated to obtain a three-dimensional coordinate D" (X, Y, Z).
10. The gas extraction borehole design method according to claim 1, wherein the borehole layout effect analysis method based on the spatial linear density analysis comprises the steps of:
step S31: dividing the grid cells;
step S32: determining a search radius r;
step S33: taking the searching radius r as a preset circle, and calculating the length of each line segment falling into the preset circle to obtain the importance of the line segment;
step S34: calculating the drilling line density of the grid unit;
step S35: and displaying the judging result in a cloud picture form, and judging the blank zone of the drilling hole.
11. The gas extraction borehole design method according to claim 10, wherein the search radius R and the borehole effective extraction radius R satisfy: r is more than or equal to 2 and less than or equal to 3R; in the step S34, the method for calculating the line density of the grid cell drill holes includes:
multiplying the length of each line segment by the importance of the drilling hole, summing the multiplied values, and dividing the summed values by the area of the preset circle;
the importance of the drilling hole can be calculated according to the effective extraction radius R of the drilling hole.
12. The gas extraction borehole design method according to claim 1, wherein the method for borehole placement effect analysis based on the borehole placement effect evaluation method for spatial nuclear density analysis comprises:
step S41: determining a kernel function;
step S42: determining the influence range of the line segment;
step S43: acquiring the importance degree of the line segment falling in the influence range of the line segment;
step S44: calculating the density of the drilling nuclei of the grid cells;
step S45: and displaying the judging result in a cloud picture form, and judging the blank bands.
13. The gas extraction borehole design method according to claim 12, wherein in the step S44, the calculation method of the grid cell borehole core density includes:
the distance weights of all line segments falling within the line segment influence range are summed.
CN202310314588.7A 2023-03-28 2023-03-28 Gas extraction drilling design method Pending CN116541904A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310314588.7A CN116541904A (en) 2023-03-28 2023-03-28 Gas extraction drilling design method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310314588.7A CN116541904A (en) 2023-03-28 2023-03-28 Gas extraction drilling design method

Publications (1)

Publication Number Publication Date
CN116541904A true CN116541904A (en) 2023-08-04

Family

ID=87456719

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310314588.7A Pending CN116541904A (en) 2023-03-28 2023-03-28 Gas extraction drilling design method

Country Status (1)

Country Link
CN (1) CN116541904A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117408085A (en) * 2023-12-12 2024-01-16 天津矿智科技有限公司 Method, system, equipment and medium for designing coal mine gas extraction hole position
CN117473792A (en) * 2023-12-22 2024-01-30 天津矿智科技有限公司 Method, system, equipment and medium for optimally designing gas extraction hole

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117408085A (en) * 2023-12-12 2024-01-16 天津矿智科技有限公司 Method, system, equipment and medium for designing coal mine gas extraction hole position
CN117408085B (en) * 2023-12-12 2024-03-01 天津矿智科技有限公司 Method, system, equipment and medium for designing coal mine gas extraction hole position
CN117473792A (en) * 2023-12-22 2024-01-30 天津矿智科技有限公司 Method, system, equipment and medium for optimally designing gas extraction hole
CN117473792B (en) * 2023-12-22 2024-03-12 天津矿智科技有限公司 Method, system, equipment and medium for optimally designing gas extraction hole

Similar Documents

Publication Publication Date Title
CN116541904A (en) Gas extraction drilling design method
CN111305286B (en) Monitoring and verifying system and method for soil-rock double-element slope overall failure mode
US7096122B2 (en) Method for producing full field radial grid for hydrocarbon reservoir simulation
Shaw et al. Modification to the Suhara‐Fukuda method of network generation
CN108012608B (en) A kind of prepartion of land method based on GNSS
CN112214867A (en) Strip mine mining boundary and mining program collaborative optimization method under complex coal seam condition
CN113703068A (en) Method for determining uranium ore distribution
CN102455435A (en) Full-coverage automatic arrangement method for block-like template observation system
CN110159346B (en) Quasi-rectangular tunnel earth surface deformation prediction method based on non-uniform convergence mode
CN105069228B (en) A kind of method that spare via are added on spare cell
CN108804780A (en) A kind of method of three dimensional analysis pile-base supporting layer
CN110543660B (en) Low-impact development simulation method, system and related device
CN109558629A (en) The generation method and device of steel bar arrangement scheme
CN113010997A (en) Early warning treatment method for pile foundation impact hole-forming construction accidents in karst development area
CN102877471B (en) Method for carrying out early warning on instability collapse of vertical deep foundation pit
CN104866664A (en) Method for predicting construction difficulty of water seepage prevention curtain-grouting
CN107507179A (en) Rock And Soil quantitative analysis method based on GOCAD
CN109029383B (en) Ocean buoy station layout rationality evaluation method based on space analysis
CN114892688B (en) Three-dimensional design method and system for side slope anchor rod frame beam
CN116486028A (en) Self-adaptive inverse distance weighted interpolation method for three-dimensional geological structure modeling
CN115391871A (en) Commercial automatic layout method, electronic device, and storage medium
CN107292099A (en) A kind of dynamic based on grid rolls thickness analysis method
Kravchenko et al. Digital models of parametric structures
CN105785444B (en) Velocity field construction method
CN105335586B (en) Method for establishing automatic cloth tube model in design of heat exchanger

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