CN102918950A - 一种采煤沉陷地复垦施工敏感区的确定方法 - Google Patents

一种采煤沉陷地复垦施工敏感区的确定方法 Download PDF

Info

Publication number
CN102918950A
CN102918950A CN2012104178432A CN201210417843A CN102918950A CN 102918950 A CN102918950 A CN 102918950A CN 2012104178432 A CN2012104178432 A CN 2012104178432A CN 201210417843 A CN201210417843 A CN 201210417843A CN 102918950 A CN102918950 A CN 102918950A
Authority
CN
China
Prior art keywords
surface units
mining
construction
distortion
reclamation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN2012104178432A
Other languages
English (en)
Other versions
CN102918950B (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.)
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
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 China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN201210417843.2A priority Critical patent/CN102918950B/zh
Publication of CN102918950A publication Critical patent/CN102918950A/zh
Application granted granted Critical
Publication of CN102918950B publication Critical patent/CN102918950B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Or Calibration Of Command Recording Devices (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明涉及一种采煤沉陷地复垦施工敏感区的确定方法,属于采矿技术、土地利用和土地复垦领域,该办法包括:地面单元的划分、地面单元栅格化;根据开采阶段进行沉陷预计,获取各个开采阶段的地表变形等值线;通过GIS的插值与拟合功能,获取各个开采阶段的地表的水平变形栅格图与曲率变形栅格图;实现地表变形等值线到变形栅格图的转换;根据抗水平变形的阈值和抗曲率变形的阈值,确定各个开采阶段复垦施工的敏感地面单元。该发明方法简单易懂,易于指导实际复垦施工的实施,避免井工开采煤矿山边采边复过程中施工时机及区域选择不当导致的复垦失败。适用于单一水平(近水平)煤层。是井工开采煤矿区边采边复技术的重要组成。

Description

一种采煤沉陷地复垦施工敏感区的确定方法
技术领域
本发明属于采矿技术、土地利用和土地复垦领域,特别涉及高潜水位-平原矿区煤炭开采下正在沉陷土地的边采边复工作,主要通过地面单元的划分,定量化的确定正在沉陷土地进行复垦施工过程中的敏感地面单元,科学的指导复垦施工。
背景技术
人多地少,耕地资源十分稀缺,这是我国的基本国情。随着经济的高速发展,对能源的需求与日俱增,煤炭作为我国的最主要能源,占一次能源消费的70%左右,由于我国90%以上的煤炭产量来自于井工开采,且多采用走向长壁全部垮落法开采,土地不可避免地产生下沉,造成大量土地的沉陷损毁。高潜水位矿区主要分布于我国的华东及东北平原区,开采所损毁的土地大多为优质的高产农田。对高潜水位矿区损毁土地的复垦是保证耕地数量,保障粮食安全的重要举措。
以往的治理多侧重于沉陷稳定后再复垦,根据该类地区多年的复垦经验,在现有技术水平下复耕率只有30%-40%,这也就意味着大面积的优质农田将会变为水域,直接影响耕地的数量。造成这一现象的原因有很多,一方面,特殊的地质采矿条件使得这一地区沉陷具有典型性,高潜水位矿区大多为平原矿区,所在区域大多为地势平坦、土壤资源丰富、肥力较高的优质基本农田区,由于本身地表地势平坦且潜水位较高,而同时煤层赋存条件又特别好,很多矿区煤层数量多、煤层厚度大,这就造成沉陷后积水的现象严重,大量的表土资源沉入水底;另一方面,由于是平原地区,缺乏固体充填物,目前主要采用煤矸石、粉煤灰、建筑垃圾等作为塌陷区治理的充填材料,与地表巨大的沉陷体积相比,此类充填物是车水杯薪。且珍贵的土壤资源作为多年风化沉积形成的资源,从某种意义上来说是不可再生的资源,复垦时缺少充填材料与表土覆盖材料,造成恢复土地和耕地的比例低。
今年来,相关学者开展“动态复垦”与“超前复垦”的研究,研究人员希望通过提前采取复垦措施来达到保护耕地的目的。由于“动态复垦”与“超前复垦”不可避免的涉及到地面沉陷前施工或者对正在沉陷的土地进行复垦施工。因此,复垦施工的对象(如耕地区、水利设施、道路等)不可避免的要遭受动态沉陷,如拉伸、倾斜、曲率变形的影响。导致很多工程在施工过程中或者施工后遭受二次破坏,致使复垦失败。如何合理的确定复垦的时间与区域,既确保复垦耕地率最大,又避免在动态开采阶段的敏感区域(受水平变形、曲率变形区域)进行施工,成为井工开采煤矿区边采边复技术发展的瓶颈,直接关系到复垦工程的成功与否。
地理信息系统(Geographic Information System,GIS),它是一种特定的十分重要的空间信息系统。它是在计算机硬、软件系统支持下,对整个或部分地球表层空间中的有关地理分布数据进行采集、储存、管理、运算、分析、显示和描述的技术系统。GIS与测绘学和地理学有着密切的关系。大地测量、工程测量、矿山测量、地籍测量、航空摄影测量和遥感技术为GIS中的空间实体提供各种不同比例尺和精度的定位数。GIS的技术优势在于它的数据综合、模拟与分析评价能力,可以得到常规方法或普通信息系统难以得到的重要信息,实现地理空间过程演化的模拟和预测。
1965年,刘宝琛、廖国华出版了《煤矿地表移动的基本规律》一书,将概率积分法全面引入中国矿区地表移动沉陷预计,目前,概率积分法也成为中国最成熟应用范围最广泛的煤矿开采沉陷预计理论。应用该理论,可获得煤炭开采过程中各阶段的地表下沉值及相应的水平变形、曲率变形、倾斜变形等各种变形值。
发明内容
本发明的目的是为解决上述难题,提出一种采煤沉陷地复垦施工动态敏感区域的确定方法,本发明能确定正在沉陷土地在采取复垦措施过程中的施工敏感区域,避免井工开采煤矿山边采边复过程中施工时机及区域选择不当导致的复垦失败。
本发明方法的流程图如图1所示。
本发明提出的采煤沉陷地复垦施工动态敏感区域的确定方法,包括地面单元的划分、地面单元栅格化、各个阶段地表变形等值线的获取、施工敏感单元的确定等几个步骤,具体叙述如下:
1)划分地面单元,并栅格化地面单元。
11)将地面划分为方格地面单元,方格地面单元各边长度为H/20所得数的10位取整;H为平均采深,单位为米;
12)将方格地面单元栅格化,栅格化单位取值范围为10m-100m;
2)根据开采阶段进行沉陷预计,获取各个开采阶段的地表变形等值线(包括水平变形等值线与曲率变形等值线);开采阶段以月或年为单位进行划分;
3)通过GIS的插值与拟合功能,获取各个开采阶段的地表的水平变形栅格图与曲率变形栅格图;实现地表变形等值线到变形栅格图的转换;
4)将获取的水平变形栅格图与曲率变形栅格图分别通过叠加运算,赋值给地面单元,获取每个地面单元各个开采阶段的变形值:地面单元的变形值取所在地面单元内栅格单元的最大值:
Si=Max(E1,E2,E3...En)
Qi=Max(F1,F2,F3...Fn)
式中,Si为地面单元i的水平变形,En为地面单元i内第n个栅格单位的水平变形值;
Qi为地面单元i的曲率变形,Fn为地面单元i内第n个栅格单位的曲率变形值;
5)确定各个开采阶段复垦施工的敏感地面单元,若Si≧Bj,则该地面单元为施工敏感单元,否则为非敏感单元,Bj为对应的施工对象的抗水平变形的阈值(施工对象的抗水平变形的阈值由《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规程》确定);若Qi≧Gj,则该地面单元为施工敏感地面单元,否则为非敏感单元,Gj为对应的施工对象的抗曲率变形的阈值(施工对象的抗曲率变形的阈值由《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规程》确定);对于各开采阶段中的每个地面单元,只要Si或者Qi中任意一个超过阈值,则该地面单元为该开采阶段的复垦施工敏感地面单元。
本发明主要具有以下技术优点:
本发明能确定正在沉陷土地的施工敏感地面单元,避免了由于在敏感区域进行施工导致的复垦失败的情况。由于采用了GIS与开采沉陷预测结合的方式,能定量的确定各个开采阶段地面的敏感单元,方法简单易懂,易于操作,是井工开采煤矿区边采边复技术的重要组成。
附图说明
图1为本发明方法流程图。
图2第6开采阶段耕地复垦施工敏感地面单元。
图3第6开采阶段沟渠及堤坝复垦施工敏感地面单元
具体实施方式
本发明根据划分的地面单元,以及分阶段进行的沉陷预计,确定地面单元任意开采阶段进行复垦施工工作时的敏感单元,使得井工煤矿区边采边复工作进行时的施工敏感区域的获取成为可能,避免在敏感区域施工导致的复垦失败。结合实施例详细说明如下:
本发明提出的采煤沉陷地复垦施工动态敏感区域的确定方法,首先获取各个阶段的水平变形等值线与曲率变形等值线,经过栅格化处理后,获取每个地面单元各个开采阶段的变形值,通过与《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规程》中对应的施工对象的抗变形的阈值,确定复垦施工的敏感单元。该方法是利用GIS的栅格数据及空间分析功能,定量确定地面复垦施工敏感单元的方法,适用于井工开采煤矿区的边采边复工作。该方法的流程图如图1所示,包括以下步骤:
1)划分地面单元,并栅格化地面单元。
11)将地面划分为方格地面单元,由H/20确定,方格地面单元各边长度为H/20所得数的10位取整;H为平均采深,单位为米;
12)将方格地面单元栅格化,栅格化单位取值范围为10m-100m;
2)根据开采阶段进行沉陷预计,获取各个开采阶段的地表变形等值线(包括水平变形等值线与曲率变形等值线);开采阶段以月或年为单位进行划分;
3)通过GIS的插值与拟合功能,获取各个开采阶段的地表的水平变形栅格图与曲率变形栅格图;实现地表变形等值线到变形栅格图的转换;
4)将获取的水平变形栅格图与曲率变形栅格图分别通过叠加运算,赋值给地面单元,获取每个地面单元各个开采阶段的变形值:地面单元的变形值取所在地面单元内栅格单元的最大值:
Si=Max(E1,E2,E3...En)
Qi=Max(F1,F2,F3...Fn)
式中,Si为地面单元i的水平变形,En为地面单元i内第n个栅格单位的水平变形值;
Qi为地面单元i的曲率变形,Fn为地面单元i内第n个栅格单位的曲率变形值;
5)确定各个开采阶段复垦施工的敏感地面单元,若Si≧Bj,则该地面单元为施工敏感单元,否则为非敏感单元,Bj为对应的施工对象的抗水平变形的阈值(施工对象的抗水平变形的阈值由《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规程》确定);若Qi≧Gj,则该地面单元为施工敏感地面单元,否则为非敏感单元,Gj为对应的施工对象的抗曲率变形的阈值(施工对象的抗曲率变形的阈值由《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规程》确定);对于各开采阶段中的每个地面单元,只要Si或者Qi中任意一个超过阈值,则该地面单元为该开采阶段的复垦施工敏感地面单元。
实施例
本实施例为某高潜水位矿区,煤层平均厚度9.0m,埋藏深度800m;地势平坦,地表自然标高+43.0~+44.50m,地下水埋深约3.0m左右,地面坡度绝大部分为0~2°之间。
本实施例的一种采煤沉陷地复垦施工动态敏感区域的确定方法的步骤如下:
1)划分地面单元,并栅格化地面单元。
11)将地面划分为方格地面单元,方格地面单元各边长度为H/20所得数的10位取整;H为平均采深,单位为米;由此确定本实施例地面单元尺寸为40m×40m,共7875个地面单元。
12)栅格化地面单元,栅格化单位选取为10m×10m;则每个地面单元有16个栅格单元。
2)根据开采阶段进行沉陷预计,获取各个开采阶段的地表变形等值线(包括水平变形等值线与曲率变形等值线);开采阶段以月为单位进行划分;共获得12期开采阶段的水平变形等值线与曲率变形等值线。
3)通过GIS的插值与拟合功能,获取各个开采阶段的地表的水平变形栅格图与曲率变形栅格图;实现地表变形等值线到变形栅格图的转换;
4)将获取的水平变形栅格图与曲率变形栅格图分别通过叠加运算,赋值给地面单元,获取每个地面单元各个开采阶段的变形值:地面单元的变形值取所在地面单元内栅格单元的最大值:
Si=Max(E1,E2,E3...En)
Qi=Max(F1,F2,F3...Fn)
式中,Si为地面单元i的水平变形,En为地面单元i内第n个栅格单位的水平变形值;
Qi为地面单元i的曲率变形,Fn为地面单元i内第n个栅格单位的曲率变形值;
5)确定各个开采阶段复垦施工的敏感地面单元,若Si≧Bj,则该地面单元为施工敏感单元,否则为非敏感单元,Bj为对应的施工对象的抗水平变形的阈值(施工对象的抗水平变形的阈值由《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规程》确定);若Qi≧Gj,则该地面单元为施工敏感地面单元,否则为非敏感单元,Gj为对应的施工对象的抗曲率变形的阈值(施工对象的抗曲率变形的阈值由《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规程》确定);对于各开采阶段中的每个地面单元,只要Si或者Qi中任意一个超过阈值,则该地面单元为该开采阶段的复垦施工敏感地面单元。
本实施例中,由于耕地区施工过程中主要考虑拉伸变形导致的裂缝,因此对于复垦为耕地的只选取Bj=5mm/m;对于进行砖和混凝土的沟渠及鱼塘堤坝建设的取Bj=2.5mm/m,Gj=12km;由此,确定了各个开采阶段地面施工敏感单元,以第6开采阶段为例,耕地复垦施工敏感地面单元共51个,如图2所示;复垦为沟渠与鱼塘堤坝的复垦施工敏感地面单元共320个,如图3所示。在第6开采阶段进行复垦工作时,应当避免在上述地面单元进行对应的复垦施工。图2、图3中,虚线为开采工作面,阴影部分为施工敏感地面单元。
可见,本发明能定量的划分地面单元,清晰明了的确定任意开采阶段地面复垦施工的敏感地面单元,为井工开采煤矿区边采边复技术的实施提供了重要的保障。从而避免由于在施工敏感区域建设造成的复垦失败与重复建设。

Claims (1)

1.一种采煤沉陷地复垦施工敏感区的确定方法,通过开采沉陷与GIS的结合,合理划分地面单元,科学划定开采阶段。由此确定各开采阶段地面复垦施工对象的敏感施工区域。本发明所提供的地面单元划分方法、栅格化单元、开采阶段划分及各开采阶段复垦施工敏感区域的确定,简单易懂,易于指导实际复垦施工的实施,避免施工时间与区域选择不当导致的复垦失败。适用于单一水平(近水平)煤层;该方法包括以下步骤:
1)划分地面单元,并栅格化地面单元;
11)将地面划分为方格地面单元,方格地面单元各边长度为H/20所得数的10位取整;H为平均采深,单位为米;
12)将方格地面单元栅格化,栅格化单位取值范围为10m-100m;
2)根据开采阶段进行沉陷预计,获取各个开采阶段的地表变形等值线;地表变形等值线包括水平变形等值线与曲率变形等值线;开采阶段以月或年为单位进行划分;
3)通过GIS的插值与拟合功能,获取各个开采阶段的地表的水平变形栅格图与曲率变形栅格图;实现地表变形等值线到变形栅格图的转换;
4)将获取的水平变形栅格图与曲率变形栅格图分别通过叠加运算,赋值给地面单元,获取每个地面单元各个开采阶段的变形值:地面单元的变形值取所在地面单元内栅格单元的最大值:
Si=Max(E1,E2,E3...En)
Qi=Max(F1,F2,F3...Fn)
式中,Si为地面单元i的水平变形,En为地面单元i内第n个栅格单位的水平变形值;
Qi为地面单元i的曲率变形,Fn为地面单元i内第n个栅格单位的曲率变形值;
5)确定各个开采阶段复垦施工的敏感地面单元,若Si≧Bj,则该地面单元为施工敏感单元,否则为非敏感单元,Bj为对应的施工对象的抗水平变形的阈值;若Qi≧Gj,则该地面单元为施工敏感地面单元,否则为非敏感单元,Gj为对应的施工对象的抗曲率变形的阈值;对于各开采阶段中的每个地面单元,若Si或者Qi中任意一个超过阈值,则该地面单元为该开采阶段的复垦施工敏感地面单元。
CN201210417843.2A 2012-10-26 2012-10-26 一种采煤沉陷地复垦施工敏感区的确定方法 Active CN102918950B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210417843.2A CN102918950B (zh) 2012-10-26 2012-10-26 一种采煤沉陷地复垦施工敏感区的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210417843.2A CN102918950B (zh) 2012-10-26 2012-10-26 一种采煤沉陷地复垦施工敏感区的确定方法

Publications (2)

Publication Number Publication Date
CN102918950A true CN102918950A (zh) 2013-02-13
CN102918950B CN102918950B (zh) 2015-04-01

Family

ID=47634001

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210417843.2A Active CN102918950B (zh) 2012-10-26 2012-10-26 一种采煤沉陷地复垦施工敏感区的确定方法

Country Status (1)

Country Link
CN (1) CN102918950B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104564069A (zh) * 2015-01-09 2015-04-29 中国矿业大学(北京) 一种基于方格网法的地面动态沉陷预测与复垦方法
CN105248007A (zh) * 2015-11-06 2016-01-20 中国矿业大学(北京) 黄河泥沙复垦采煤沉陷地交替多层多次充填土壤重构方法
CN105509708A (zh) * 2015-11-21 2016-04-20 西安科技大学 一种黄土沟壑区下煤层开采地表下沉系数的确定方法
CN108446450A (zh) * 2018-02-23 2018-08-24 煤炭工业济南设计研究院有限公司 一种受采动影响的建筑物破坏程度的分析计算方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101176400A (zh) * 2006-11-08 2008-05-14 中国矿业大学(北京) 一种矿山塌陷区土地复垦新方法
CN101422094A (zh) * 2007-11-02 2009-05-06 煤炭科学研究总院唐山研究院 平原矿区厚煤层开采塌陷地动态预复垦方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101176400A (zh) * 2006-11-08 2008-05-14 中国矿业大学(北京) 一种矿山塌陷区土地复垦新方法
CN101422094A (zh) * 2007-11-02 2009-05-06 煤炭科学研究总院唐山研究院 平原矿区厚煤层开采塌陷地动态预复垦方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
曹化平等: "基于GIS的矿山开采沉陷信息可视化应用", 《测绘工程》, vol. 19, no. 3, 30 June 2010 (2010-06-30) *
赵晓东等: "开采沉陷预测的GIS面元栅格数字化模型及其应用", 《煤炭学报》, vol. 25, no. 5, 31 October 2000 (2000-10-31) *
郝长胜: "基于GIS的煤矿开采沉陷预测分析系统研究", 《煤炭技术》, vol. 29, no. 6, 30 June 2010 (2010-06-30) *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104564069A (zh) * 2015-01-09 2015-04-29 中国矿业大学(北京) 一种基于方格网法的地面动态沉陷预测与复垦方法
CN105248007A (zh) * 2015-11-06 2016-01-20 中国矿业大学(北京) 黄河泥沙复垦采煤沉陷地交替多层多次充填土壤重构方法
CN105509708A (zh) * 2015-11-21 2016-04-20 西安科技大学 一种黄土沟壑区下煤层开采地表下沉系数的确定方法
CN108446450A (zh) * 2018-02-23 2018-08-24 煤炭工业济南设计研究院有限公司 一种受采动影响的建筑物破坏程度的分析计算方法
CN108446450B (zh) * 2018-02-23 2021-07-02 通用技术集团工程设计有限公司 一种受采动影响的建筑物破坏程度的分析计算方法

Also Published As

Publication number Publication date
CN102918950B (zh) 2015-04-01

Similar Documents

Publication Publication Date Title
Li et al. A review on anthropogenic geomorphology
Xiao et al. Dynamic subsidence simulation and topsoil removal strategy in high groundwater table and underground coal mining area: a case study in Shandong Province
CN102609782B (zh) 基于采煤沉陷模拟和复耕率的边采边复时机的优选方法
Hu et al. Review on current status and challenging issues of land subsidence in China
Wu et al. Zonation of the landslide hazards in the forereservoir region of the Three Gorges Project on the Yangtze River
Li et al. Innovation for sustainable mining: Integrated planning of underground coal mining and mine reclamation
Jemec Auflič et al. The variety of landslide forms in Slovenia and its immediate NW surroundings
Marschalko et al. Subsidence map of underground mining influence for urban planning: an example from the Czech Republic
Li et al. Mining scheme design for super-high water backfill strip mining under buildings: a Chinese case study
CN103643956B (zh) 一种采煤塌陷区煤矸石充填复垦表层覆土厚度设计方法
CN111798052B (zh) 一种高潜水位采煤沉陷积水区三维空间信息动态预测方法
CN102918950B (zh) 一种采煤沉陷地复垦施工敏感区的确定方法
Pan et al. Mechanism of sinkhole formation during groundwater-level recovery in karst mining area, Dachengqiao, Hunan province, China
CN102889083B (zh) 基于gis栅格单元的采煤沉陷地表土剥离时空确定方法
CN103235853A (zh) 用于采煤沉陷地边采边复的基于土方平衡的基塘布局方法
Kılıç et al. Microzonation of Zeytinburnu region with respect to soil amplification: a case study
Wang et al. Three-dimensional analyses of the sliding surface distribution in the Huangtupo No. 1 riverside sliding mass in the Three Gorges Reservoir area of China
CN102877846B (zh) 基于单元法的边采边复表土剥离时机确定方法
Hsieh et al. Late Holocene episodic river aggradation along the Lao-nong River (southwestern Taiwan): An application to the Tseng-wen Reservoir Transbasin Diversion Project
Wang et al. Surface subsidence and its reclamation of a coal mine locating at the high groundwater table, China
Zhao et al. Discrete element analysis of deformation features of slope controlled by karst fissures under the mining effect: a case study of Pusa landslide, China
Atiq et al. Dam Site Identification Using Remote Sensing and GIS (A case study Diamer Basha Dam Site)
CN103742151B (zh) 应用于边采边复技术的土壤重构方法
Marschalko et al. Relevance of building site categories implementation into the land use plan in underground mining area in the region of Petrvald (Czech Republic)
CN103226738B (zh) 一种基于预期复垦耕地率和复垦成本的复垦时机优选方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant