CN109902416B - Ground subsidence multivariable mathematical relationship calculation method - Google Patents
Ground subsidence multivariable mathematical relationship calculation method Download PDFInfo
- Publication number
- CN109902416B CN109902416B CN201910174224.7A CN201910174224A CN109902416B CN 109902416 B CN109902416 B CN 109902416B CN 201910174224 A CN201910174224 A CN 201910174224A CN 109902416 B CN109902416 B CN 109902416B
- Authority
- CN
- China
- Prior art keywords
- relation
- relationship
- unit
- accumulated
- ground subsidence
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The application discloses a ground subsidence multivariable mathematical relationship, which comprises a relationship I, a relationship II and a relationship III; the expressions of the relation I, the relation II and the relation III are respectively as follows: l (L) cn ＝4×10 ‑5 P n 3 ‑0.003P n 2 +0.079P n +0.873；S cn ＝‑0.086P n 3 +3.562P n 2 ‑21.899P n +444.09；The application can identify the response relation between the underground water exploitation amount and the ground subsidence, and through verification, the correlation of the relation I reaches 0.9954, the correlation of the relation II reaches 0.9889, the correlation of the relation III reaches 0.8712, and the accuracy is better.
Description
Technical Field
The application belongs to the field of geological disaster prevention and control management, and particularly relates to a ground subsidence multivariable mathematical relationship calculation method.
Background
Groundwater resources are important water resources for human beings, and reasonable management and control are needed. In order to maintain the requirements of economy on water resources in a part of areas, underground water is excessively and intensively exploited, so that a large-range water level dropping funnel is generated in urban water sources, even geological disasters such as ground subsidence and karst collapse are caused, the quality of the underground water in the part of areas is continuously deteriorated, pollution is serious, and immeasurable loss is brought to social economy and people's life.
Super-mining of groundwater is the most significant cause of ground subsidence. The response relationship between the underground water exploitation amount and the ground settlement amount (the central accumulated settlement amount, the settlement area and the funnel area) is not clear at home and abroad. And the response relation between the underground water exploitation amount and the ground subsidence is clear, so that the method has important significance for identifying the ground subsidence mechanism and controlling the ground subsidence.
Disclosure of Invention
The technical problem to be solved by the application is to provide a ground subsidence multivariable mathematical relationship calculating method capable of clearly and accurately identifying the response relationship between the underground water exploitation amount and the ground subsidence.
In order to solve the technical problems, the application adopts the following technical scheme: a calculation method of a ground subsidence multivariable mathematical relationship is provided, and a relationship is established according to underground water exploitation quantity, central accumulated subsidence quantity, subsidence area and subsidence funnel area;
the expression of the relation I is as follows:
L cn ＝4×10 -5 P n 3 -0.003P n 2 +0.079P n +0.873 relationship I;
L cn represents the accumulated sedimentation quantity of the center of the nth year, and the unit is m; p (P) n Represents the accumulated exploitation amount of groundwater in the nth year, and the unit is hundred million m 3 。
Further, the method also comprises a relation pi, wherein the expression of the relation pi is as follows:
S cn ＝-0.086P n 3 +3.562P n 2 -21.899P n +444.09 relationship n;
S cn represents the sedimentation area of the nth year, and the unit is km 2 ；P n Represents the accumulated exploitation amount of groundwater in the nth year, and the unit is hundred million m 3 。
Further, the method also comprises a relation pi I, wherein the expression of the relation pi I is as follows:
P n represents the accumulated exploitation amount of groundwater in the nth year, and the unit is hundred million m 3 ；S cn Represents the sedimentation area of the nth year, and the unit is km 2 ；L cn Represents the accumulated sedimentation quantity of the center of the nth year, and the unit is m; s is S ln Represents the area of the dropping funnel of the nth year, and the unit is km 2 。
The beneficial effects of the application are as follows:
the application can identify the response relation between the underground water exploitation amount and the ground subsidence, and through verification, the correlation of the relation I reaches 0.9954, the correlation of the relation II reaches 0.9889, and the correlation of the relation III reaches 0.8712, so that the application can clearly show the mathematical relation between the underground water exploitation amount and the central accumulation subsidence amount, the mathematical relation between the underground water exploitation amount and the subsidence area, and the mathematical relation between the underground water exploitation amount and the central accumulation subsidence amount, the subsidence area and the subsidence funnel area, and has better accuracy.
Drawings
FIG. 1 is a diagram of Fuyang city for years L cn And P n Is a visual graph of (1);
FIG. 2 is a diagram of Fuyang city for years cn And P n Is a visual graph of (1);
FIG. 3 is a diagram of Fuyang city for yearsAnd->Is a visual graph of (c).
Detailed Description
In order to make the technical problems, technical solutions and advantages to be solved more apparent, the following detailed description will be given with reference to the accompanying drawings and specific embodiments. It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other.
Aiming at the defect that the response relationship between the underground water exploitation amount and the ground settlement amount is not clear at home and abroad at present, the application provides a ground settlement multivariable mathematical relationship calculation method with better accuracy. The mathematical relationship includes a relationship I, the expression of which is as follows:
L cn ＝4×10 -5 P n 3 -0.003P n 2 +0.079P n +0.873 relationship I;
L cn represents the accumulated sedimentation quantity of the center of the nth year, and the unit is m; p (P) n Represents the accumulated exploitation amount of groundwater in the nth year, and the unit is hundred million m 3 。
In one embodiment, the method further comprises a relationship pi, wherein the expression of the relationship pi is as follows:
S cn ＝-0.086P n 3 +3.562P n 2 -21.899P n +444.09 relationship n;
S cn represents the sedimentation area of the nth year, and the unit is km 2 ；P n Represents the accumulated exploitation amount of groundwater in the nth year, and the unit is hundred million m 3 。
In one embodiment, the method further comprises a relationship pi I, wherein the relationship pi I has the following expression:
P n represents the accumulated exploitation amount of groundwater in the nth year, and the unit is hundred million m 3 ；S cn Represents the sedimentation area of the nth year, and the unit is km 2 ；L cn Represents the accumulated sedimentation quantity of the center of the nth year, and the unit is m; s is S ln Represents the area of the dropping funnel of the nth year, and the unit is km 2 。
The application relates to a relation I, a relation II and a relation III, which are ground subsidence multivariable mathematical relation calculation methods with better accuracy finally confirmed by a great deal of creative labor through a great deal of verification, analysis and modification by designing hundreds of mathematical models through years of researches.
Table 1 is the measured data of groundwater production, central accumulated sedimentation, sedimentation area, and dropping funnel area from 1996 to 2016 in Fuyang;
TABLE 1
Verification of the application with years of measured data from Fuyang CityAccuracy of relation I, relation II and relation III, L of Fuyang city for years cn And P n 、S cn And P n 、And->As shown in figures 1 to 3, the correlation of the relation I reaches 0.9954, the correlation of the relation II reaches 0.9889, and the correlation of the relation III reaches 0.8712, which shows that the application can clearly show the mathematical relationship between the accumulated groundwater yield and the central accumulated settlement, the mathematical relationship between the accumulated groundwater yield and the settlement area, and the mathematical relationship between the accumulated groundwater yield and the central accumulated settlement, the settlement area and the settlement funnel area, and has better accuracy.
It should be understood that the examples and embodiments described herein are for illustrative purposes only and are not intended to limit the present application, and that various modifications or changes in light thereof will be suggested to persons skilled in the art and are to be included within the spirit and purview of this application.
Claims (3)
1. A method for calculating a multi-variable mathematical relationship of ground subsidence is characterized by comprising the following steps: establishing a relation according to the underground water exploitation quantity, the central accumulated settlement quantity, the settlement area and the dropping funnel area;
the expression of the relation I is as follows:
L cn ＝4×10 -5 P n 3 -0.003P n 2 +0.079P n +0.873 relationship I;
L cn represents the accumulated sedimentation quantity of the center of the nth year, and the unit is m; p (P) n Represents the accumulated exploitation amount of groundwater in the nth year, and the unit is hundred million m 3 。
2. A method of calculating a ground subsidence multivariate mathematical relationship according to claim 1, wherein: also included is a relationship pi, the expression of which is as follows:
S cn ＝-0.086P n 3 +3.562P n 2 -21.899P n +444.09 relationship n;
S cn represents the sedimentation area of the nth year, and the unit is km 2 ；P n Represents the accumulated exploitation amount of groundwater in the nth year, and the unit is hundred million m 3 。
3. A method of calculating a ground subsidence multivariate mathematical relationship according to claim 1 or 2, wherein: also included is a relationship pi I, the expression of which is as follows:
P n represents the accumulated exploitation amount of groundwater in the nth year, and the unit is hundred million m 3 ；S cn Represents the sedimentation area of the nth year, and the unit is km 2 ；L cn Represents the accumulated sedimentation quantity of the center of the nth year, and the unit is m; s is S ln Represents the area of the dropping funnel of the nth year, and the unit is km 2 。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910174224.7A CN109902416B (en) | 2019-03-08 | 2019-03-08 | Ground subsidence multivariable mathematical relationship calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910174224.7A CN109902416B (en) | 2019-03-08 | 2019-03-08 | Ground subsidence multivariable mathematical relationship calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109902416A CN109902416A (en) | 2019-06-18 |
| CN109902416B true CN109902416B (en) | 2023-09-12 |
Family
ID=66946861
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910174224.7A Active CN109902416B (en) | 2019-03-08 | 2019-03-08 | Ground subsidence multivariable mathematical relationship calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109902416B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1776458A (en) * | 2005-11-17 | 2006-05-24 | 上海交通大学 | Method for determining underground water extractable amount according to annual ground setting and pump-out |
| RU2376469C1 (en) * | 2008-07-21 | 2009-12-20 | Учреждение Российской академии наук Центр геофизических исследований Владикавказского научного центра РАН и Правительства Республики Северная Осетия-Алания (ЦГИ ВНЦ РАН и РСО-А) | Device for determination of ground surface subsidence |
| KR100982448B1 (en) * | 2010-03-03 | 2010-09-16 | 한국지질자원연구원 | Ground subsidence prediction system and predicting method using the same |
| CN102494667A (en) * | 2011-11-30 | 2012-06-13 | 同济大学 | Characterizing method of land subsidence |
| CN104091063A (en) * | 2014-07-01 | 2014-10-08 | 天津市水文水资源勘测管理中心 | Method for using regression analysis for determining groundwater control level |
-
2019
- 2019-03-08 CN CN201910174224.7A patent/CN109902416B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1776458A (en) * | 2005-11-17 | 2006-05-24 | 上海交通大学 | Method for determining underground water extractable amount according to annual ground setting and pump-out |
| RU2376469C1 (en) * | 2008-07-21 | 2009-12-20 | Учреждение Российской академии наук Центр геофизических исследований Владикавказского научного центра РАН и Правительства Республики Северная Осетия-Алания (ЦГИ ВНЦ РАН и РСО-А) | Device for determination of ground surface subsidence |
| KR100982448B1 (en) * | 2010-03-03 | 2010-09-16 | 한국지질자원연구원 | Ground subsidence prediction system and predicting method using the same |
| CN102494667A (en) * | 2011-11-30 | 2012-06-13 | 同济大学 | Characterizing method of land subsidence |
| CN104091063A (en) * | 2014-07-01 | 2014-10-08 | 天津市水文水资源勘测管理中心 | Method for using regression analysis for determining groundwater control level |
Non-Patent Citations (1)
| Title |
|---|
| 地面沉降与地下水开采相关关系及沉降预测;杨丽萍;周志华;;南水北调与水利科技(第03期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109902416A (en) | 2019-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Murray | State-scale perspective on water use and production associated with oil and gas operations, Oklahoma, US | |
| Fan et al. | Industrialization from scratch: The “Construction of Third Front” and local economic development in China's hinterland | |
| Coppola Jr et al. | A neural network model for predicting aquifer water level elevations | |
| CN107331132B (en) | A kind of method and system of Urban Fires hidden danger dynamic prediction monitoring | |
| CN109236273B (en) | Dynamic data processing method for oil field development and production | |
| CN112990739B (en) | Phosphorite area water environment risk assessment method based on multi-medium pollution characteristic identification | |
| He et al. | Assessment of ecological vulnerability of resource-based cities based on entropy-set pair analysis | |
| CN101266299A (en) | Method for forecasting oil gas utilizing earthquake data object constructional features | |
| Xu et al. | China Sponge City database development and urban runoff source control facility configuration comparison between China and the US | |
| CN110469299A (en) | A kind of exploitation of oil-extracting well water injection takes effect effect evaluation method | |
| CN104899602A (en) | User cluster analysis system based on K-means algorithm | |
| CN105545265A (en) | Large-area tight gas reservoir exploitation method | |
| CN111191817B (en) | Bus network topology division method based on transfer passenger flow | |
| CN109902416B (en) | Ground subsidence multivariable mathematical relationship calculation method | |
| CN106382102B (en) | A kind of overflow early warning method based on clustering algorithm | |
| CN107605474B (en) | Method and device for predicting gas formation yield while drilling | |
| Minh et al. | Assessment of groundwater level and quality: A case study in O Mon and Binh Thuy districts, Can Tho City, Vietnam | |
| CN108442921A (en) | It is a kind of to consider that the oil well output of time-varying and interlayer interference splits a point method | |
| Valdes‐Abellan et al. | Validating the KAGIS black‐box GIS‐based model in a Mediterranean karst aquifer: Case of study of Mela aquifer (SE Spain) | |
| CN108119121A (en) | A kind of horizontal well refracturing potentiality well rapid screening method | |
| Gosselain et al. | A large-scale stream benthic diatom database | |
| CN107818189B (en) | Method for calculating recoverable resource amount of shale gas | |
| CN106803139A (en) | Help stall oil well economic life time Benefit Evaluation Method long | |
| Belis | Palaeoenvironmental reconstruction of Lago di Albano (Central Italy) during the Late Pleistocene using fossil ostracod assemblages | |
| CN111485868A (en) | Coal bed gas field power utilization reserve estimation method based on development scheme |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |