CN108446450B - Method for analyzing and calculating damage degree of building influenced by mining - Google Patents

Method for analyzing and calculating damage degree of building influenced by mining Download PDF

Info

Publication number
CN108446450B
CN108446450B CN201810154722.0A CN201810154722A CN108446450B CN 108446450 B CN108446450 B CN 108446450B CN 201810154722 A CN201810154722 A CN 201810154722A CN 108446450 B CN108446450 B CN 108446450B
Authority
CN
China
Prior art keywords
deformation
building
mining
membership
curvature
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.)
Active
Application number
CN201810154722.0A
Other languages
Chinese (zh)
Other versions
CN108446450A (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.)
General Technology Gropu Engineering Design Co ltd
Original Assignee
General Technology Gropu Engineering Design 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 General Technology Gropu Engineering Design Co ltd filed Critical General Technology Gropu Engineering Design Co ltd
Priority to CN201810154722.0A priority Critical patent/CN108446450B/en
Publication of CN108446450A publication Critical patent/CN108446450A/en
Application granted granted Critical
Publication of CN108446450B publication Critical patent/CN108446450B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Abstract

The invention belongs to the technical field of mineral engineering, and relates to an analysis and calculation method for the damage degree of a building influenced by mining. And selecting six influence factors of the horizontal deformation value of the building after coal seam mining, the plane area generating the horizontal deformation, the plane area generating the inclined deformation and the inclined deformation, the plane area generating the curvature deformation and the curvature deformation as the evaluation indexes of the damage degree of the building influenced by mining. The method for analyzing and calculating the damage degree of the building influenced by mining evaluates the damage degree of the building on the whole, can accurately analyze the influence damage condition of the coal seam on the building after mining, provides scientific basis for safe use of the building after mining, and is suitable for wide popularization and application.

Description

Method for analyzing and calculating damage degree of building influenced by mining
Technical Field
The invention belongs to the technical field of mineral engineering, and relates to an analysis and calculation method for the damage degree of a building influenced by mining.
Background
China is rich in coal resources, and not only a large amount of coal resources exist underground in plains, mountainous regions and other areas, but also a large amount of coal resources exist under a plurality of buildings. For a long time, China always adopts a strip mining or filling mining method to mine a large amount of coal resources under buildings. After the coal resources under the building are mined, the mining influence can be caused to the building, the building can be damaged in serious conditions, and even the building collapses, so that the life safety is threatened. At present, after coal and carbon resources under buildings in China are mined, influence damage to the buildings is mainly calculated through a probability integral method to obtain deformation values (including horizontal deformation epsilon, inclined deformation i and curvature deformation k) of the buildings, and according to a building damage grade table in building, water, railway and main roadway coal pillar reservation and coal briquetting mining procedures, the damage degree of the buildings is divided into I, II, III and IV grades, but only the deformation value of the buildings is considered, the damage degree of the buildings is not analyzed and evaluated from the damage range, and the damage degree of the buildings after coal seam mining under the buildings cannot be accurately reflected.
Disclosure of Invention
In order to make up for the defects of the prior art, the invention provides the method for analyzing and calculating the damage degree of the buildings influenced by mining, which is reasonable in analysis and convenient to calculate.
In order to achieve the purpose, the invention is realized by adopting the following technical scheme:
a method for analyzing and calculating the damage degree of a building influenced by mining is characterized by comprising the following steps: the analysis steps are as follows:
1) selecting six influence factors of horizontal deformation of a building deformation value and a plane area generating the horizontal deformation after coal seam mining, inclined deformation and a plane area generating the inclined deformation, curvature deformation and a plane area generating the curvature deformation as evaluation indexes of the damage degree of the building influenced by mining;
2) calculating the horizontal deformation epsilon, the inclined deformation i and the curvature deformation k of the building after coal seam mining by adopting a probability integration method, and respectively superposing a horizontal deformation epsilon contour map, an inclined deformation i contour map and a curvature deformation k contour map with a plan view of the building to obtain a plan area m for generating the horizontal deformation epsilonεPlane area m where inclined deformation i occursiAnd a plane area m where curvature deformation k occursk
3)Respectively constructing a membership function a of the horizontal deformation epsilon and a plane area m generating the horizontal deformation epsilonεMembership function b of (a), membership function c of oblique deformation i, and area m of plane generating oblique deformation iiMembership function d of (d), membership function e of curvature deformation k, plane area m of generating curvature deformation kkA membership function f;
membership function of horizontal deformation ε:
Figure 533927DEST_PATH_IMAGE001
in the formula: a is the horizontal deformation membership; epsilon is horizontal deformation generated by a building after coal seam mining, and the unit is mm/m (epsilon is used as the horizontal deformation value of the building, has positive or negative, and the absolute value is required to be taken when the epsilon range is determined); e is horizontal deformation which can be borne by the building and has a unit of mm/m;
plane area m generating horizontal deformation epsilonεMembership function of (a):
Figure 336798DEST_PATH_IMAGE002
in the formula: b is the area m of the plane where the horizontal deformation epsilon occursεDegree of membership of; m isεProducing a planar area for the building with a horizontal deformation of epsilon in m2(ii) a M is the plane area of the building, unit M2
Membership function of tilt deformation i:
Figure 475655DEST_PATH_IMAGE003
in the formula: c is the inclination deformation membership; i is the inclined deformation of the building after coal seam mining, and the unit is mm/m (i is used as the inclined deformation value of the building, has positive or negative values, and needs to take the absolute value when the range of i is determined); i is the inclination deformation which can be borne by the building, and the unit is mm/m;
plane area m of generating inclined deformation iiMembership function of (a):
Figure 59084DEST_PATH_IMAGE004
in the formula: d is the area m of the plane where the inclined deformation i occursiDegree of membership of; m isiProducing a planar area for the building with an oblique deformation of i, in m2(ii) a M is the plane area of the building, unit M2
Membership function of curvature deformation k:
Figure 246482DEST_PATH_IMAGE005
in the formula: e is the membership of curvature deformation k; k is the curvature deformation of the building after coal seam mining, unit 10-3M (k is used as a curvature deformation value of the building, has positive or negative values, and needs to take an absolute value when a k range is determined); k is the curvature deformation that the building can withstand, unit 10-3/m;
Plane area m generating curvature deformation kkMembership function of (a):
Figure 231756DEST_PATH_IMAGE006
in the formula: f is the area m of the plane where curvature deformation k occurskDegree of membership of; m iskProducing a planar area for a building with a curvature deformation of k, in m2(ii) a M is the plane area of the building, unit M2
In order to comprehensively evaluate the damage degree of the buildings affected by mining, the a, b, c, d, e and f are subjected to superposition calculation:
λ=βa+γb+δc+ηd+θe+μf
in the formula: lambda is the total membership of the damage degree of the building influenced by mining, beta is the weight of a, gamma is the weight of b, delta is the weight of c, eta is the weight of d, theta is the weight of e, mu is the weight of f, beta, gamma, delta, eta, theta and mu are directly selected according to the actual condition of the building, and are subjected to multiple superposition trial calculation and determined after comprehensive analysis;
4) after all the weights are selected through comprehensive analysis, calculating to obtain the total membership lambda of the damage degree of the building influenced by mining;
5) evaluating the damage condition of the building according to the total membership lambda of the damage degree of the building influenced by mining, wherein the larger the lambda value is, the larger the damage degree of the building is, and the smaller the damage degree of the building is otherwise;
6) and establishing an analysis and calculation model of the damage degree of the building influenced by mining.
Compared with the prior art, the invention has the advantages and positive effects that:
the method for analyzing and calculating the damage degree of the building influenced by mining evaluates the damage degree of the building as a whole, and can accurately analyze the damage condition of the coal seam after mining on the building. Compared with the previous research results, the method considers more influencing factors and enables the evaluation result to be more scientific. According to the specific evaluation result of the patent, the safety of the building can be accurately evaluated, a scientific basis is provided for whether the building can be safely used after coal seam mining, and the method is suitable for wide popularization and application.
Detailed Description
In order that the above objects, features and advantages of the present invention may be more clearly understood, the present invention will be further described with reference to specific embodiments. It should be noted that the embodiments and features of the embodiments of the present application may be combined with each other without conflict.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, however, the present invention may be practiced in other ways than those specifically described herein, and thus the present invention is not limited to the specific embodiments of the present disclosure.
Embodiment 1, this embodiment provides specific steps of a method for analyzing and calculating the damage level of a building affected by mining, as follows:
1) selecting six influence factors of horizontal deformation of a building deformation value and a plane area generating the horizontal deformation after coal seam mining, inclined deformation and a plane area generating the inclined deformation, curvature deformation and a plane area generating the curvature deformation as evaluation indexes of the damage degree of the building influenced by mining;
2) calculating the horizontal deformation epsilon, the inclined deformation i and the curvature deformation k of the building after coal seam mining by adopting a probability integration method, and respectively superposing a horizontal deformation epsilon contour map, an inclined deformation i contour map and a curvature deformation k contour map with a plan view of the building to obtain a plan area m for generating the horizontal deformation epsilonεPlane area m where inclined deformation i occursiAnd a plane area m where curvature deformation k occursk
3) Respectively constructing a membership function a of the horizontal deformation epsilon and a plane area m generating the horizontal deformation epsilonεMembership function b of (a), membership function c of oblique deformation i, and area m of plane generating oblique deformation iiMembership function d of (d), membership function e of curvature deformation k, plane area m of generating curvature deformation kkA membership function f;
membership function of horizontal deformation ε:
Figure 977733DEST_PATH_IMAGE001
in the formula: a is the horizontal deformation membership; epsilon is horizontal deformation of a building generated after coal seam mining, and the unit is mm/m; e is horizontal deformation which can be borne by the building and has a unit of mm/m;
plane area m generating horizontal deformation epsilonεMembership function of (a):
Figure 48457DEST_PATH_IMAGE002
in the formula: b is the area m of the plane where the horizontal deformation epsilon occursεDegree of membership of; m isεProducing a planar area for the building with a horizontal deformation of epsilon in m2(ii) a M is the plane area of the building, unit M2
Membership function of tilt deformation i:
Figure 39547DEST_PATH_IMAGE003
in the formula: c is the inclination deformation membership; i is the inclined deformation of the building after coal seam mining, and the unit is mm/m; i is the inclination deformation which can be borne by the building, and the unit is mm/m;
plane area m of generating inclined deformation iiMembership function of (a):
Figure 879327DEST_PATH_IMAGE004
in the formula: d is the area m of the plane where the inclined deformation i occursiDegree of membership of; m isiProducing a planar area for the building with an oblique deformation of i, in m2(ii) a M is the plane area of the building, unit M2
Membership function of curvature deformation k:
Figure 359987DEST_PATH_IMAGE005
in the formula: e is the membership of curvature deformation k; k is the curvature deformation of the building after coal seam mining, unit 10-3(ii)/m; k is the curvature deformation that the building can withstand, unit 10-3/m;
Plane area m generating curvature deformation kkMembership function of (a):
Figure 590111DEST_PATH_IMAGE006
in the formula: f is the area m of the plane where curvature deformation k occurskDegree of membership of; m iskProducing a planar area for a building with a curvature deformation of k, in m2(ii) a M is the plane area of the building, unit M2
In order to comprehensively evaluate the damage degree of the buildings affected by mining, the a, b, c, d, e and f are subjected to superposition calculation:
λ=βa+γb+δc+ηd+θe+μf
in the formula: λ is the total membership of the degree of destruction of the building affected by mining, β is the weight of a, γ is the weight of b, δ is the weight of c, η is the weight of d, θ is the weight of e, μ is the weight of f. Directly selecting beta, gamma, delta, eta, theta and mu according to the actual condition of the building, performing multiple superposition trial calculation, and determining after comprehensive analysis;
4) after all the weights are selected through comprehensive analysis, calculating to obtain the total membership lambda of the damage degree of the building influenced by mining;
5) evaluating the damage condition of the building according to the total membership lambda of the damage degree of the building influenced by mining, wherein the larger the lambda value is, the larger the damage degree of the building is, and the smaller the damage degree of the building is otherwise;
6) and establishing an analysis and calculation model of the damage degree of the building influenced by mining.
The above description is only a preferred embodiment of the present invention, and not intended to limit the present invention in other forms, and any person skilled in the art may apply the above modifications or changes to the equivalent embodiments with equivalent changes, without departing from the technical spirit of the present invention, and any simple modification, equivalent change and change made to the above embodiments according to the technical spirit of the present invention still belong to the protection scope of the technical spirit of the present invention.

Claims (1)

1. A method for analyzing and calculating the damage degree of a building influenced by mining is characterized by comprising the following steps: the analysis steps are as follows:
1) selecting six influence factors of horizontal deformation of a building deformation value and a plane area generating the horizontal deformation after coal seam mining, inclined deformation and a plane area generating the inclined deformation, curvature deformation and a plane area generating the curvature deformation as evaluation indexes of the damage degree of the building influenced by mining;
2) calculating the horizontal deformation epsilon, the inclined deformation i and the curvature deformation k of the building after coal seam mining by adopting a probability integration method, and summing the horizontal deformation epsilon contour map, the inclined deformation i contour mapThe curvature deformation k contour map is respectively superposed with the plane map of the building to obtain the plane area m generating the horizontal deformation epsilonεPlane area m where inclined deformation i occursiAnd a plane area m where curvature deformation k occursk
3) Respectively constructing a membership function a of the horizontal deformation epsilon and a plane area m generating the horizontal deformation epsilonεMembership function b of (a), membership function c of oblique deformation i, and area m of plane generating oblique deformation iiMembership function d of (d), membership function e of curvature deformation k, plane area m of generating curvature deformation kkA membership function f;
membership function of horizontal deformation ε:
Figure 952858DEST_PATH_IMAGE001
in the formula: a is the horizontal deformation membership; epsilon is horizontal deformation of a building generated after coal seam mining, and the unit is mm/m; e is horizontal deformation which can be borne by the building and has a unit of mm/m;
plane area m generating horizontal deformation epsilonεMembership function of (a):
Figure 848133DEST_PATH_IMAGE002
in the formula: b is the area m of the plane where the horizontal deformation epsilon occursεDegree of membership of; m isεProducing a planar area for the building with a horizontal deformation of epsilon in m2(ii) a M is the plane area of the building, unit M2
Membership function of tilt deformation i:
Figure 585145DEST_PATH_IMAGE003
in the formula: c is the inclination deformation membership; i is the inclined deformation of the building after coal seam mining, and the unit is mm/m; i is the inclination deformation which can be borne by the building, and the unit is mm/m;
plane area m of generating inclined deformation iiMembership function of (a):
Figure 108530DEST_PATH_IMAGE004
in the formula: d is the area m of the plane where the inclined deformation i occursiDegree of membership of; m isiProducing a planar area for the building with an oblique deformation of i, in m2(ii) a M is the plane area of the building, unit M2
Membership function of curvature deformation k:
Figure 272795DEST_PATH_IMAGE005
in the formula: e is the membership of curvature deformation k; k is the curvature deformation of the building after coal seam mining, unit 10-3(ii)/m; k is the curvature deformation that the building can withstand, unit 10-3/m;
Plane area m generating curvature deformation kkMembership function of (a):
Figure 514421DEST_PATH_IMAGE006
in the formula: f is the area m of the plane where curvature deformation k occurskDegree of membership of; m iskProducing a planar area for a building with a curvature deformation of k, in m2(ii) a M is the plane area of the building, unit M2
In order to comprehensively evaluate the damage degree of the buildings affected by mining, the a, b, c, d, e and f are subjected to superposition calculation:
λ=βa+γb+δc+ηd+θe+μf
in the formula: λ is total membership of the destruction degree of the building affected by mining, β is weight of a, γ is weight of b, δ is weight of c, η is weight of d, θ is weight of e, and μ is weight of f;
4) after all the weights are selected through comprehensive analysis, calculating to obtain the total membership lambda of the damage degree of the building influenced by mining;
5) evaluating the damage condition of the building according to the total membership lambda of the damage degree of the building influenced by mining, wherein the larger the lambda value is, the larger the damage degree of the building is, and the smaller the damage degree of the building is otherwise;
6) and establishing an analysis and calculation model of the damage degree of the building influenced by mining.
CN201810154722.0A 2018-02-23 2018-02-23 Method for analyzing and calculating damage degree of building influenced by mining Active CN108446450B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810154722.0A CN108446450B (en) 2018-02-23 2018-02-23 Method for analyzing and calculating damage degree of building influenced by mining

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810154722.0A CN108446450B (en) 2018-02-23 2018-02-23 Method for analyzing and calculating damage degree of building influenced by mining

Publications (2)

Publication Number Publication Date
CN108446450A CN108446450A (en) 2018-08-24
CN108446450B true CN108446450B (en) 2021-07-02

Family

ID=63192820

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810154722.0A Active CN108446450B (en) 2018-02-23 2018-02-23 Method for analyzing and calculating damage degree of building influenced by mining

Country Status (1)

Country Link
CN (1) CN108446450B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102279421A (en) * 2011-05-11 2011-12-14 北方工业大学 Slope rock mass stability evaluation method
CN102609782A (en) * 2012-01-10 2012-07-25 中国矿业大学(北京) Optimization method of simultaneous harvest and reclamation opportunity based on mining subsidence simulation and second ploughing rate
CN102918950A (en) * 2012-10-26 2013-02-13 中国矿业大学(北京) Method for determining reclamation construction sensitive area in coal mining subsidence land
CN103364533A (en) * 2013-06-27 2013-10-23 山东科技大学 Analytical method of exploiting baseplate disturbance destruction degree on coal mine confined aquifer
CN103670511A (en) * 2013-10-19 2014-03-26 煤炭工业济南设计研究院有限公司 Roof strata fracture degree analysis and calculation method
CN103726481A (en) * 2014-01-09 2014-04-16 辽宁工程技术大学 Repeatedly mined surface deformation predicting method
CN103902780A (en) * 2014-04-08 2014-07-02 中国矿业大学 Method for predicting deformation of solid-filled coal mining surfaces
JP2016212909A (en) * 2016-08-04 2016-12-15 前田建設工業株式会社 Building soundness evaluation device, building soundness evaluation method, and building soundness evaluation program
CN107092740A (en) * 2017-04-12 2017-08-25 乌鲁木齐城市轨道集团有限公司 A kind of old goaf residual deformation method for predicting of severe inclined thick coal seam

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102279421A (en) * 2011-05-11 2011-12-14 北方工业大学 Slope rock mass stability evaluation method
CN102609782A (en) * 2012-01-10 2012-07-25 中国矿业大学(北京) Optimization method of simultaneous harvest and reclamation opportunity based on mining subsidence simulation and second ploughing rate
CN102918950A (en) * 2012-10-26 2013-02-13 中国矿业大学(北京) Method for determining reclamation construction sensitive area in coal mining subsidence land
CN103364533A (en) * 2013-06-27 2013-10-23 山东科技大学 Analytical method of exploiting baseplate disturbance destruction degree on coal mine confined aquifer
CN103670511A (en) * 2013-10-19 2014-03-26 煤炭工业济南设计研究院有限公司 Roof strata fracture degree analysis and calculation method
CN103726481A (en) * 2014-01-09 2014-04-16 辽宁工程技术大学 Repeatedly mined surface deformation predicting method
CN103902780A (en) * 2014-04-08 2014-07-02 中国矿业大学 Method for predicting deformation of solid-filled coal mining surfaces
JP2016212909A (en) * 2016-08-04 2016-12-15 前田建設工業株式会社 Building soundness evaluation device, building soundness evaluation method, and building soundness evaluation program
CN107092740A (en) * 2017-04-12 2017-08-25 乌鲁木齐城市轨道集团有限公司 A kind of old goaf residual deformation method for predicting of severe inclined thick coal seam

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Collapses of underground cavities and soil-structure interactions: influences of the position of the structure relative to the cavity;Matthieu Caudron;Marwan Al Heib;Fabrice Emeriault;《The 12th International Conference of International Association for Computer Methods and Advances in Geomechanics》;20081231;第1-8页 *
采空塌陷区地表建筑物变形特征分析及破坏机理研究;翟淑花,张长敏,齐干,冉淑红;《城市地质》;20150331;第10卷;第59-63页 *

Also Published As

Publication number Publication date
CN108446450A (en) 2018-08-24

Similar Documents

Publication Publication Date Title
WO2019233114A1 (en) Method for evaluating risk of water damage from strata water inrush of overlying strata in mining area
CN106407545A (en) Soft soil foundation ditch safety evaluation method based on fuzzy comprehensive judgment method
CN107305700A (en) Judge coal and gas whether the method that can be protruded
CN106339798A (en) Soft rock inclined shaft supporting and protecting system stability evaluation method based on fuzzy comprehensive evaluation method
CN106651060B (en) A kind of method of the disconnected control natural gas pool enrichment degree of quantitative assessment
CN103258298A (en) Method for assessing safety of open pit coal mine
CN116597616B (en) Intelligent monitoring and early warning system for geological disasters in mining area
CN114841532A (en) Safety evaluation method and system for surface subsidence in shield excavation process
CN106251053A (en) The suitability assessment method of high voltage power transmission transmission tower is built on earth's surface, a kind of sinking land in coalmining areas
CN112907130A (en) Construction risk grade determination method, device and medium based on shield construction parameters
CN104992051A (en) Method and system for risk level evaluation of fuel gas polyethylene pipeline
CN108446450B (en) Method for analyzing and calculating damage degree of building influenced by mining
AU2021103225A4 (en) Embankment dam analysis to reduce flood impact using sensor modelling and gis techniques
CN106121721B (en) A kind of tight roof Face Ground Pressure Behavior stage division
CN107180306A (en) A kind of technical method of coal resources payable area block weighted superposition identification
CN106815445A (en) A kind of shallow buried coal seam section coal pillar LOAD FOR model and its computational methods
CN104008290B (en) Consider the fall of ground hidden danger stage division of two-way anisobaric stress field
CN113420457A (en) End slope stability analysis method for open pit coal mine steep wall mining
CN114329680A (en) Stability evaluation method for ore pillar dam body of underground reservoir in mining area and application of stability evaluation method
CN114135288B (en) Method for optimizing high-pressure water jet slotting pressure relief parameters of rock burst coal seam roadway
Zhao et al. Comprehensive optimization and engineering applications of thick residual coal re-mining methodology
CN112446592A (en) Rock burst risk grade evaluation method based on equivalent mining depth estimation
Gao et al. Stability analysis of Cavern No. 5 with the largest span in Heidong large ancient underground quarry
CN109723414A (en) A kind of critical water injection pressure calculation method of oil-gas anticlinal deposits crack elimination
Cao et al. Construction and Application of “Active Prediction-Passive Warning” Joint Impact Ground Pressure Resilience Prevention System: Take the Kuan Gou Coal Mine as an Example

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: No. 141, Dikou Road, Tianqiao District, Jinan City, Shandong Province

Applicant after: GENERAL TECHNOLOGY GROPU ENGINEERING DESIGN Co.,Ltd.

Address before: No. 141, Dikou Road, Tianqiao District, Jinan City, Shandong Province

Applicant before: COAL INDUSTRY JINAN DESIGN & RESEARCH Co.,Ltd.

GR01 Patent grant
GR01 Patent grant