Invention content
In order to solve the deficiencies in the prior art, the present invention proposes a kind of sinking land in coalmining areas basic data informationization analysis exhibition
Show method.
The present invention adopts the following technical scheme that:
A kind of sinking land in coalmining areas basic data informationization analysis methods of exhibiting, described method includes following steps:
(1) acquisition and standardization of sinking land in coalmining areas difference basic data:
Different basic datas include InSAR data extraction and dem data update, satellite image extracting data pools zone model
Enclose, the acquisition of manual measurement data, the extraction of subsided land from coal mine administration project data, the acquisition of mining production digging related data
And analysis, land use data obtain and analysis;
(2) Overlaying analysis between different data sources and secondary development:
Overlaying analysis including sinking land in coalmining areas deformation data Yu mining production mining operation data;Sinking land in coalmining areas deformation
The Overlaying analysis of data, depression ponding data and administration project data;Sinking land in coalmining areas deformation data, depression ponding data and soil
Ground utilizes the Overlaying analysis of present situation;The Overlaying analysis of sinking land in coalmining areas deformation data, depression ponding data;
(3) extraction and analysis of spatial data:
All kinds of initial data, superposition of data and secondary development data make it have unification by processing in step (1), (2)
Attribute structure, database can allowed arbitrarily to call, and simple calculations can be carried out, to realize intuitive human-computer interaction;
(4) scientific information of sinking land in coalmining areas basic data:
By the way that step (1), (2), (3) data are concentrated on a platform realizing the unified display of data base.
In the step (1), different basic datas include InSAR data extraction and dem data update:InSAR data source
After improvement, registration and coherent analysis based on two rail difference, simultaneously using multi-baseline interference accumulation of phase stacking method
Landform phase is eliminated after being overlapped with known DEM terrain datas, can be sought the accumulation deformation quantity in corresponding observation period, be led to
Being superimposed for cumulative deformation and original DEM is crossed, and can get newest DEM present situations, dem data has coordinate and elevation information etc.
Attribute structure can be directly used for database;
Satellite image extracting data pools zone range:Select sinking land in coalmining areas wet season sequence satellite remote sensing for many years shadow
Picture, image coordinate system is corrected to it is consistent with present situation dem data, pass through GIS software extract different times ponding range line, will
Ponding range vector quantization, imparting attribute structure information can be used for spatial database calling and analyzes pools zone range reason with this
And rule;
The acquisition of manual measurement data:Mistake null point is made up, skin is used using GPS combination lining ropes for small range ponding region
Raft ship carries out manual measurement, and a wide range of waters is measured using Water depth measuring instrument, and measurement data carries attribute information, can be direct
For spatial database, for the amendment of pools zone cumulative deformation and the update of dem data;
The extraction of subsided land from coal mine administration project data:All subsided land from coal mine administration project related datas in area are carried out
Analysis is aggregated, data content includes:Administer area's range, as-built drawing, control investment volume, ginseng after topographic map, improvement before improvement
It builds unit, restore land type and area etc., by the data vector after summarizing and assigning attribute structure makes database facilitate tune
With;
The acquisition and analysis of mining production digging related data:Mining production digging related data includes:Mining area range, coal
Charcoal resource actual output, Subsidence Area actual observation surface strata movement parameter, mining engineering plan etc., by mine mining area range,
Goaf range is extracted one by one, is overlapped with present situation accumulated deformation data, this two kinds of data are mutually corrected, scarce mending-leakage is looked into,
Ensure that Subsidence Area data cover is comprehensive, and eliminates the influence of other ground deformation factors.Such data generally pass through vector quantization and
With attribute information, database can be directly used for;
Land use data obtains and analysis:Regional land use type shown in present landuse map reflection, will
Above-mentioned different range boundary line carries out overlapping extraction with present landuse map, you can obtains present status of land utilization within the scope of this
Statistical data, present landuse map attribute structure information is complete, can directly invoke its various information.
In the step (1), the InSAR data impact analysis carries out D-InSAR differential SAR Interferometry analyses first,
Selected Space Baseline and interference pair are interfered, landform phase difference, filtering, solution twine processing;
Then centainly range line is collapsed using interpolation accumulation stack result extraction;
Later, while obtaining sinking land in coalmining areas totality deformation data using multi-baseline interference accumulation of phase stacking method,
Obtain changing character;
Finally, multidate SAR Strength Changes detection methods are utilized after the range of satellite image extracting data pools zone, are carried out
Pools zone identifies.
In the step (2), the Overlaying analysis and two based on MSPS Prediction of Collapse with Magnetic Resonance Imaging analysis softwares between different data sources
Secondary exploitation:
The Overlaying analysis of sinking land in coalmining areas deformation data and mining production mining operation data:Ensure coal-mining subsidence data
The influence that all standing, accuracy, the exclusive PCR factor for verifying deformation data are brought, and adopted for analysis and summary different type
Dig engineering, the generation development law of different coal mining conditions coal-mining subsidence;Realize that different administrative regions, different barghs are adopted
The analysis summary of coal depression present situation;Divide left over by history subsided land range;Improve the accurate of bargh's ten thousand tons of subsidence rates of coal mining
Property;The anti-reasonability for pushing away mine coal-mining subsidence parameter;
The Overlaying analysis of sinking land in coalmining areas deformation data, depression ponding data and administration project data:Show administration project
With the relative position relation of Subsidence Area;Summarize the governance model rule of different depression types and ponding situation Subsidence Area;Realization is collapsed
The management and regulation effect for administering progress sunkenly summarize demonstration;Also can be the deployment guideline direction of administration project;
The Overlaying analysis of sinking land in coalmining areas deformation data, depression ponding data and present status of land utilization:Divide administrative region, divide
Sink deformation area damage to land type, area ownership feature situation are reflected in bargh, plough especially in the sinking land in coalmining areas
Deformation failure and ponding situation, and with this evaluate coal-mining subsidence to agricultural production, urban traffic, water conservancy projects destroying infection,
Economic loss situation caused by refining coal-mining subsidence;
The Overlaying analysis of sinking land in coalmining areas deformation data, depression ponding data:Carry out rainfall statistics, floods analysis, converges
Water analysis can system the analysis and summary sinking land in coalmining areas ponding rule, calculate it is estimated floods possibility, and be the sinking land in coalmining areas
Administration project direction.
In the step (4), the scientific information of sinking land in coalmining areas basic data includes:
The information-based displaying of sinking land in coalmining areas distribution:Showed in the form of closed line on striograph or other geographical map sources
Mining area range and sinking land in coalmining areas distribution, and depression area can be counted by administrative area or bargh, form intuitive chart;
The information-based displaying of coal-mining subsidence degree and influence:On striograph or other geographical map sources in the form of closed line
Show different stage depression depth area's distribution and depression pools zone range, and can count different by administrative area or bargh
Depression rank, ponding degree distribution area and account for damage land area situation, form intuitive chart;
The information-based displaying of subsided land from coal mine administration project:On striograph or other geographical map sources in the form of closed line
Show the distribution situation of sinking land in coalmining areas administration project, and shows single engineering unit of taking part in building, investment, work in the form of list
The basic informations such as Cheng Jindu, regulation effect are realized and restore land type, economy by administrative area or bargh's statistics administration project
Environmental benefit.
Use the advantageous effects that technical solution as above obtains for:
1. data base is unified, the present situation realization of former each department, every profession and trade data format disunity, one-sidedness is changed
The all standing of basic data.Each department's data are concentrated on a platform, inherently a kind of performance of uniform format.
2. the promotion of data utilization ratio.Multiple portions such as the governments such as mineral products, soil, house property, planning, urban transportation, water conservancy
Door is required for using for reference the data of other relevant departments when carrying out the work or specifying decision, and the realization of this patent is by various data
Carry out fusion and secondary development, multilayer, Layering manifestation work(can solve the demand that data are used for reference mutually, greatly improve
The utilization ratios of data.
3. the figure management combined using geography information, management information, GIS auxiliary phases, coal mine subsidence land is believed substantially
Breath carries out fusion and background process greatly, so that dull number is converted into figure, table, and by being combined the face of realizing with network service
Scientific information-based displaying to user.
Specific implementation mode
It is described further in conjunction with 1 to 7 pair of specific implementation mode of the invention of attached drawing:
A kind of sinking land in coalmining areas basic data informationization analysis methods of exhibiting, described method includes following steps:
(1) acquisition and standardization of mining collapse area difference basic data:
Different basic datas include InSAR data extraction and dem data update:InSAR data source (uses interferometry skill
The synthetic aperture radar of art or data) after improvement, registration and coherent analysis based on two rail difference, using more bases
Line interferometric phase accumulates stacking method and (digital elevation model passes through limited landform altitude number with known DEM terrain datas
The factually existing digitized simulation to ground surface or terrain.) be overlapped after eliminate landform phase, can seek tired in corresponding observation period
Product deformation quantity, by being superimposed for cumulative deformation and original DEM, and available newest DEM as-is datas, dem data have
Coordinate and elevation information attribute structure, are directly used in database.
Satellite image extracting data pools zone range:Select sinking land in coalmining areas wet season sequence satellite remote sensing for many years shadow
Picture, image coordinate system is corrected to it is consistent with present situation dem data, pass through GIS software extract different times ponding range line, will
Ponding range vector quantization, assign attribute structure information for spatial database call with this analyze pools zone range reason and
Rule.
The acquisition of manual measurement data:Mistake null point is made up, skin is used using GPS combination lining ropes for small range ponding region
Raft ship carries out manual measurement, and a wide range of waters is measured using Water depth measuring instrument, and measurement data carries attribute information, directly uses
In spatial database, for the amendment of pools zone cumulative deformation and the update of dem data.
The extraction of subsided land from coal mine administration project data:All subsided land from coal mine administration project related datas in area are carried out
Analysis is aggregated, data content includes:Administer area's range, as-built drawing, control investment volume, ginseng after topographic map, improvement before improvement
It builds unit, restore land type and area etc., by the data vector after summarizing and assigning attribute structure makes database identify simultaneously
It calls.
The acquisition and analysis of mining production digging related data:Mining production digging related data includes:Mining area range, coal
Charcoal resource actual output, goaf data, Subsidence Area actual observation surface strata movement parameter, mining engineering plan etc., by mine
Mountain mining area range, goaf range are extracted one by one, are overlapped with present situation accumulated deformation data, by mining production digging dependency number
It is mutually corrected according to, two kinds of data of present situation accumulated deformation data, looks into scarce mending-leakage, ensure that Subsidence Area data cover is comprehensive, and eliminate it
He influences ground deformation factor.Such data generally pass through vector quantization and have attribute information, can be directly used for database.
Land use data obtains and analysis:Regional land use type shown in present landuse map reflection, will
Above-mentioned different range boundary line carries out overlapping extraction with present landuse map, obtains the statistics of present status of land utilization within the scope of this
Data, present landuse map attribute structure information is complete, can directly invoke its various information.
(2) Overlaying analysis between different data sources and secondary development:
The Overlaying analysis of sinking land in coalmining areas deformation data and mining production mining operation data:Ensure coal-mining subsidence data
The influence that all standing, accuracy, the exclusive PCR factor for verifying deformation data are brought, and adopted for analysis and summary different type
Dig engineering, the generation development law of different coal mining conditions coal-mining subsidence;Realize that different administrative regions, different barghs are adopted
The analysis summary of coal depression present situation;Divide left over by history subsided land range;Improve the accurate of bargh's ten thousand tons of subsidence rates of coal mining
Property;The anti-reasonability for pushing away mine coal-mining subsidence parameter.
The Overlaying analysis of sinking land in coalmining areas deformation data, depression ponding data and administration project data:Show administration project
With the relative position relation of Subsidence Area;Summarize the governance model rule of different depression types and ponding situation Subsidence Area;Realization is collapsed
The management and regulation effect for administering progress sunkenly summarize demonstration;For the deployment guideline direction of harnessing project.
The Overlaying analysis of sinking land in coalmining areas deformation data, depression ponding data and present status of land utilization:Divide administrative region, divide
Sink deformation area damage to land type, area ownership feature situation are reflected in bargh, plough especially in the sinking land in coalmining areas
Deformation failure and ponding situation, and with this evaluate coal-mining subsidence to agricultural production, urban traffic, water conservancy projects destroying infection,
Economic loss situation caused by refining coal-mining subsidence;
The Overlaying analysis of sinking land in coalmining areas deformation data, depression ponding data:Carry out rainfall statistics, floods analysis, converges
Water analysis, can system the analysis and summary sinking land in coalmining areas ponding rule, calculate it is estimated floods possibility, and be the sinking land in coalmining areas
Harnessing project direction.
(3) extraction and analysis of spatial data:
All kinds of initial data, superposition of data and secondary development data make it have unification by processing in step (1), (2)
Attribute structure, database can allowed arbitrarily to call, and simple calculations can be carried out, to realize intuitive human-computer interaction.
Such as:Somewhere Subsidence Area area is 10000m on coal-mining subsidence present situation figure2, which includes affiliated mine
Deng;Area of Subsidence Area depression 1m or more is 3000m at this on coal-mining subsidence degree figure2, face of the Subsidence Area depression less than 1m
Product is 7000m2, attribute includes depression grade, ground deformation degree etc.;Subsidence Area includes ploughing at this on present landuse map
2000m at ground 42, 8000m at forest land 22, each plot attribute information includes land type, area, soil ownership etc.;Pass through
Extraction, analysis and the plus and minus calculation of spatial database, you can know the cultivated area of region yields 1m or more at this, forest land face
The information such as long-pending, affiliated administrative region can also equally obtain the damage degree ploughed in the area, the damage degree in forest land etc..
(4) scientific information of sinking land in coalmining areas basic data:
By the way that step (1), (2), (3) data are concentrated on a platform realizing the unified display of data base.
The information-based displaying of sinking land in coalmining areas distribution:Showed in the form of closed line on striograph or other geographical map sources
Mining area range and sinking land in coalmining areas distribution, and depression area can be counted by administrative area or bargh, form intuitive chart.
The information-based displaying of coal-mining subsidence degree and influence:On striograph or other geographical map sources in the form of closed line
Show different stage depression depth area's distribution and depression pools zone range, and can count different by administrative area or bargh
Depression rank, ponding degree distribution area and account for damage land area situation, form intuitive chart.
The information-based displaying of subsided land from coal mine administration project:On striograph or other geographical map sources in the form of closed line
Show the distribution situation of sinking land in coalmining areas administration project, and shows single engineering unit of taking part in building, investment, work in the form of list
The basic information of Cheng Jindu, regulation effect are realized and restore land type, economy by administrative area or bargh's statistics administration project
Environmental benefit.
InSAR data impact analysis carries out D-InSAR differential SAR Interferometry analyses first, selectes Space Baseline and interference
It is right, interfered, landform phase difference, filtering, solution twine processing;Then centainly model is collapsed using interpolation accumulation stack result extraction
Contour;Later, it while obtaining sinking land in coalmining areas totality deformation data using multi-baseline interference accumulation of phase stacking method, obtains
Changing character;Finally, multidate SAR Strength Changes detection side is utilized after the range of satellite image extracting data pools zone
Method carries out pools zone identification.
The selected a certain area in southwest in Shandong Province, carries out D-InSAR differential SAR Interferometry analyses, is less than to meeting time reference line first
200 days, less than 34 of 300 meters interference pair of Space Baseline are selected, are interfered, landform phase difference, filtering, solution twine processing,
Specifically use two rail difference processing methods.
It is extracted using interpolation accumulation stack result and centainly collapses range line, selection collapses 5cm range line drawings, sees and extracted
Range wire shaped.
The apparent region of depression is taken from 5cm range line selections are collapsed, is obtained using multi-baseline interference accumulation of phase stacking method
While the totality deformation data of the sinking land in coalmining areas, changing character is obtained.
Follow-up 5cm ranges line drawing is used for as this mining area InSAR monitoring results using more baselines accumulation stack result.
Stacking method correlation analysis is accumulated according to more baselines, is supervised more than InSAR deformation using non-interpolation accumulation stack result for extracting
" the null value area " of survey ability extracts, and null value area is included in 5cm range lines;It is collapsed using interpolation accumulation stack result for extracting
5cm range lines.
Deformation results are superimposed to accumulation, the 5cm automatically extracted by deformation magnitude collapses line of demarcation, it can be seen that in mining area model
Monitor that the activity that collapses in multiple goafs, gob collapse spatially show as funnel shaped, draw with mining area exploitation in enclosing
The surface subsidence space characteristics risen meet preferable.
While more baseline accumulations are superimposed upon to obtain mining area totality deformation data, it can also reflect that mining area deformation occurs well
Process.From the apparent region of multiple sedimentations that accumulation deformation results are selected, and multiple characteristic points are chosen, carries out multiple regions respectively
The changing character of sedimentation is analyzed, and the period is by taking in April, 2013 in October, 2015 as an example.
By taking the A1 of certain point as an example, its sequential change procedure is analyzed, point A1 settling characteristics are apparent, and subsiding center more connects
It is continuous.Its deformation generating process is further analyzed, mining area deformation is relatively stable between June in April, 2013-, and in June, 2013-, September occurred
The trend of deformation aggravation;Deformation is more steady between in September, 2013 in June ,-2014, between in June, 2014-September deformation further add
It is acute;Deformation is steady between in September, 2014 in May ,-2015, and October in May, 2015-, deformation tendency is further exacerbated by;Therefore, the area
In domain between in October, -2015 in April, 2013, deformation was happened between annual 5, June -9, October, be rendered as it is non-linear by
Increased process is walked, as shown in Figure 6.
Multidate SAR Strength Changes detection methods are utilized after the range of satellite image extracting data pools zone, carry out ponding
Area identifies.On SAR images, black or dead color is generally presented in water body, pools zone can be identified according to this feature,
More scape images are regarded more using multidate SAR Strength Changes detection method, accuracy registration, multidate filter, geocoding
With radiation correction etc. processing.According to pools zone sequential intensity map, observation judges pools zone range in a certain period of time.
Data concentrate on a platform realizing the unified display of data base, as shown in Figure 2-5, form mining and collapse
Fall into ground dynamic monitoring management system centralized displaying.
The accuracy of InSAR monitoring results passes through 60 monitoring points randomly selected in thousands of working faces in a certain mining area
For, eyeball and InSAR monitoring results are compared respectively, monitoring time 20141005,20141122,20141207,
20141221,20150110,20150131,20150214,20150322,20150523 and 20150627.Corresponding SAR
Data obtaining time 20140921,20141202,20150212,20150401,20150519,20150730.
When carrying out precision test, it should ensure that the date of the two is consistent as possible, or the difference time is little.Therefore it selects
20141005-20141207 and 20140921-20141202,20141005-20150214 and 20140921-20150212 and
The data of tri- periods of 20141005-20150523 and 20140921-20150519 carry out precision test respectively.
As shown in fig. 7, as can be seen from the figure observation and InSAR value variation tendencies are more consistent, actual observed value with
The average error of InSAR values is 0.013m, coincidence loss requirement.
Certainly, described above is only that presently preferred embodiments of the present invention is answered the present invention is not limited to enumerate above-described embodiment
When explanation, any technical person familiar with the field is all equivalent substitutes for being made, bright under the guidance of this specification
Aobvious variant, all falls within the essential scope of this specification, ought to be protected by the present invention.