CN109724570A - 地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法 - Google Patents
地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法 Download PDFInfo
- Publication number
- CN109724570A CN109724570A CN201811450752.2A CN201811450752A CN109724570A CN 109724570 A CN109724570 A CN 109724570A CN 201811450752 A CN201811450752 A CN 201811450752A CN 109724570 A CN109724570 A CN 109724570A
- Authority
- CN
- China
- Prior art keywords
- drop
- underground
- bank
- layer thickness
- water
- 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
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于陡坎跌水领域,具体涉及地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法。步骤1:在上游区打孔,记陡坎处为x轴的坐标原点,获得上游区打孔处距坐标原点的距离L;步骤2:测量上游区打孔处距坐标原点L处的含水层厚度hl、水平渗流的水力传导系数k、陡坎处水平渗流的水力传导系数kr和垂直渗流的水力传导系数kd;步骤3:根据步骤1和步骤2得到的L、hl、k、kr和kd的值通过式1计算地下跌水的跌水量q、跌水宽度R和坎上水层厚度h0。本发明客观模拟地下水渗流场的跌水现象,为地下水模拟与地下水评价提供数学模型,提高地下水计算与评价精度,可广泛应用于地下水资源评价与保护、地下水环境影响评价与污染防控、地下工程的水害评价与防治等领域。
Description
技术领域
本发明属于水文地质技术领域,特别是地下水流计算领域,具体涉及地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法。
背景技术
在盆地的山前多分布有平行于山体的大断层,第四系基底常因断层作用而产生大幅度的垂直位移,导致第四系潜水流经基底陡坎时由水平径流折转为垂直渗漏而形成地下跌水,使得跌水陡坎两侧的地下水位不连续,地下水位差可高达几十米甚至100m以上,跌水陡坎两侧上下游含水层的地下水已无直接的水力联系,仅发生着水量联系,即陡坎下游含水层将承接来自陡坎上游含水层的地下水径流,这意味着下游含水层的承接水量应等于上游含水层在陡坎处的排出水量,即跌水量。跌水量仅由上游含水层水平径流区与垂直渗漏区的渗透性和水力条件控制。
目前,针对陡坎跌水的解析模型研究较少,代表性的研究成果主要有“跌水量估算的半解析模型”(Wenping Li,Zhenying Liu等人,2011年)和“渗流场分析的解析模型”(A.R.Kacimov,2012),前者是基于均质各向异性含水层的剖面二维稳定流模型导出的,开发的单宽跌水量半解析解的模型包GWF,已作为地下水数值模拟大型软件MODFLOW的一个结构模块,并成功地应用于上层滞水、断层跌水的两个算例;后者是基于无界均质各向同性含水层的剖面二维稳定流模型,采用复变函数论的映射法导出的,该解析模型用于分析地下跌水渗流场有很好的效果。然而,二者均没有对非均质含水层的地下跌水解析模型进行研究,地下跌水陡坎两侧含水层的非均质是客观存在且非常显著的,不考虑地下跌水陡坎两侧含水层的非均质特性,将会导致计算分析得到的跌水量和渗流场不够准确,对地下跌水机理的解释不够清楚。
发明内容
针对现有技术中存在的现有陡坎跌水的解析模型未考虑陡坎两侧含水层的非均质特性而导致计算分析得到的跌水量和渗流场不够准确的问题,本发明提出了地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法,包括如下步骤:
步骤1:在上游区打孔,记陡坎处为x轴的坐标原点,获得上游区打孔处距坐标原点的距离L;
步骤2:测量上游区打孔处距坐标原点L处的含水层厚度hl、水平渗流的水力传导系数k、陡坎处水平渗流的水力传导系数kr和垂直渗流的水力传导系数kd;
步骤3:根据步骤1和步骤2得到的L、hl、k、kr和kd的值通过式Ⅰ计算地下跌水的跌水量q、跌水宽度R和坎上水层厚度h0,
其中:α和β为无量纲参数,且
地下跌水区渗流场模型的构造方法,包括如下步骤:
步骤1:采用如权利要求1所述的地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法,获得地下跌水的跌水量q、跌水宽度R和坎上水层厚度h0;
步骤2:利用二维流有限差分法及MODFLOW软件处理步骤1得到的跌水量q、跌水宽度R和坎上水层厚度h0,获得地下跌水区渗流场模型。
本发明具有以下有益效果:
本发明客观模拟地下水渗流场的跌水现象,为地下水模拟与地下水评价提供数学模型,提高地下水计算与评价精度,可广泛应用于地下水资源评价与保护、地下水环境影响评价与污染防控、地下工程的水害评价与防治等领域。
附图说明
图1为基底陡降引起地下跌水的纵向剖面结构模型图;
图2为地下跌水不同模拟情况的h0、R、q曲线图;
图3(a)为上游单元和断层垂直x方向图;
图3(b)为上游单元和断层垂直y方向图;
图4为地下跌水数值模拟区的位置与基本特征图;
图5为地下跌水数值模拟区的含水层结构剖面图;
图6为数值模拟区地下跌水模拟剖面图。
具体实施方式
1、地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法,其特征在于,包括如下步骤:
步骤1:在上游区打孔,记陡坎处为x轴的坐标原点,获得上游区打孔处距坐标原点的距离L;
步骤2:测量上游区打孔处距坐标原点L处的含水层厚度hl、水平渗流的水力传导系数k、陡坎处水平渗流的水力传导系数kr和垂直渗流的水力传导系数kd;
步骤3:根据步骤1和步骤2得到的L、hl、k、kr和kd的值通过式Ⅰ计算地下跌水的跌水量q、跌水宽度R和坎上水层厚度h0,
其中:α和β为无量纲参数,且
具体的,式1还可写为:
具体的,L的数值由钻孔的位置到陡坎的位置确定,当L确定后hl的数值能够直接测得,k、kr和kd为常数能够测量。
基底陡降引起地下跌水的纵向剖面如图1所示,选择陡坎处为坐标x轴的原点,选择陡坎下游含水层的基底为潜水位H的起点,潜水位H与坐标轴x的交汇点用符号O表示。陡坎上游含水层基底高度为O,上游距坐标原点L处的潜水位为Hl、含水层厚度为hl,陡坎上的潜水位为HO、含水层厚度为hO,陡坎下的潜水位为Hr,地下跌水宽度为R、跌水高度为D。显然有:h=H-B,D=B-Hr。
地下跌水区渗流场模型的构造方法,包括如下步骤:
步骤1:获得地下跌水的跌水量q、跌水宽度R和坎上水层厚度h0;
步骤2:利用二维流有限差分法及MODFLOW软件处理步骤1得到的跌水量q、跌水宽度R和坎上水层厚度h0,获得地下跌水区渗流场模型。
Wenping Li,Zhenying Liu等人(2011年)基于断层陡坎单宽跌水量q的半解析模型,研究了断层陡坎跌水量Q的有限差分的数值解,开发了MODFLOW的GWF软件包。该GWF包包括三个子程序:GWFAL,GWFRP和GWFFM,前两个子程序分别用于分配存储器和数据准备;子程序GWFFM用来计算地下水穿越断层的跌水量Q,并将该跌水量由断层上游节点分配给断层下游节点。
其中,基于断层陡坎单宽跌水量q的解析模型,得到了断层陡坎跌水量Q的有限差分的数值解:假设断层上游单元的剖分间隔为Δx或Δy,如图3所示,则上游单元节点到断层的距离为Δx/2(或Δy/2),故基于式Ⅰ的断层陡坎上游单元的跌水量Qx和Qy数值解可写为式2和式3:
利用式Ⅱ或式Ⅲ得到断层陡坎跌水量Q,采用平面二维流有限差分法及MODFLOW软件获得地下跌水区渗流场模型。
本发明客观模拟地下水渗流场的跌水现象,为地下水模拟与地下水评价提供数学模型,提高地下水计算与评价精度,可广泛应用于地下水资源评价与保护、地下水环境影响评价与污染防控、地下工程的水害评价与防治等领域。
下面结合地下跌水渗流精进行计算分析,给定不同的L、hl、k、kr和kd的值,利用式Ⅰ计算出地下跌水不同模拟情况的无量纲综合参数α和β以及地下跌水的跌水量q、跌水宽度R和坎上水层厚度h0,根据地下跌水不同模拟情况绘制表1。
通过表1绘制的不同距离L、不同含水层厚度hl、不同参数(k、kr和kd)的q、R和h0的曲线,如图2所示。
从图2中曲线①③⑤⑦可以看出,坎上水层厚度h0特征为:L越长、hl越低时,h0值越小;kr或kd值越小时,h0值越大。
从图2中曲线①④⑤⑧可以看出,跌水宽度R特征为:L越长、hl越低时,R值越小;kr值越小时,R值越小,kd值越小时,R值越大。
从图2中曲线②④⑥⑧可以看出,单宽跌水量q特征为:L越长、hl越低时,q值越小;kr或kd值越小时,q值越大。
表1地下跌水不同模拟情况的q、R、h0计算成果表
以下给出本发明的具体实施方式,需要说明的是本发明并不局限于以下具体实施例,凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。
实施例1:
地下跌水数值模拟区为东西长70000m、南北宽30000m的矩形区域如图4所示。如图5所示,模拟区的南部山区与中部隆起区主要为第三系基岩,为不透水岩层;南盆地和倾斜平原区及穿越隆起区的四个过水通道(即为图4的①、②、③、④)均为第四系松散堆积层,其中赋存孔隙潜水;中部隆起区北侧断层将潜水含水层的基底大幅垂直错落,并导致断层两侧潜水位不连续,四个过水通道在断层处的潜水位下跌高达200~20m。
本渗流场数值模拟采用规则网格剖分,剖分间距为200m×200m,南北剖分为150行,东西剖分为350列,共52500个单元,其中有效计算单元39983个。模拟区为单一潜水含水层,断层两侧含水层的水力传导系数kx=ky90m/day,含水层给水度μ=0.08~0.20,跌水区松散层的水力传导系数kr=kd=5m/day。断层陡坎跌水量Q依据式Ⅱ或式Ⅲ给出,采用平面二维流有限差分法及MODFLOW软件模拟了地下跌水区渗流场,如图6所示,模拟的渗流场准确地反映了断层跌水现象。
Claims (2)
1.地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法,其特征在于,包括如下步骤:
步骤1:在上游区打孔,记陡坎处为x轴的坐标原点,获得上游区打孔处距坐标原点的距离L;
步骤2:测量上游区打孔处距坐标原点L处的含水层厚度hl、水平渗流的水力传导系数k、陡坎处水平渗流的水力传导系数kr和垂直渗流的水力传导系数kd;
步骤3:根据步骤1和步骤2得到的L、hl、k、kr和kd的值通过式Ⅰ计算地下跌水的跌水量q、跌水宽度R和坎上水层厚度h0,
其中:α和β为无量纲参数,且
2.地下跌水区渗流场模型的构造方法,其特征在于,包括如下步骤:
步骤1:采用如权利要求1所述的地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法,获得地下跌水的跌水量q、跌水宽度R和坎上水层厚度h0;
步骤2:利用二维流有限差分法及MODFLOW软件处理步骤1得到的跌水量q、跌水宽度R和坎上水层厚度h0,获得地下跌水区渗流场模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811450752.2A CN109724570B (zh) | 2018-11-30 | 2018-11-30 | 地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811450752.2A CN109724570B (zh) | 2018-11-30 | 2018-11-30 | 地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109724570A true CN109724570A (zh) | 2019-05-07 |
CN109724570B CN109724570B (zh) | 2020-05-19 |
Family
ID=66295883
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811450752.2A Active CN109724570B (zh) | 2018-11-30 | 2018-11-30 | 地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109724570B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111965720A (zh) * | 2020-08-19 | 2020-11-20 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | 一种基于地-井联合获取水力传导系数的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL7807471A (nl) * | 1978-07-11 | 1980-01-15 | Skadoc 77 I O B V | Werkwijze en inrichting voor het bepalen van het pro- fiel van de bodem van een water. |
JP2009300346A (ja) * | 2008-06-17 | 2009-12-24 | Aichi Tokei Denki Co Ltd | 水流計測システム及び水流計測方法 |
CN109213966A (zh) * | 2018-09-25 | 2019-01-15 | 深圳市工勘岩土集团有限公司 | 悬挂式截水帷幕绕渗涌水量计算方法 |
-
2018
- 2018-11-30 CN CN201811450752.2A patent/CN109724570B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL7807471A (nl) * | 1978-07-11 | 1980-01-15 | Skadoc 77 I O B V | Werkwijze en inrichting voor het bepalen van het pro- fiel van de bodem van een water. |
JP2009300346A (ja) * | 2008-06-17 | 2009-12-24 | Aichi Tokei Denki Co Ltd | 水流計測システム及び水流計測方法 |
CN109213966A (zh) * | 2018-09-25 | 2019-01-15 | 深圳市工勘岩土集团有限公司 | 悬挂式截水帷幕绕渗涌水量计算方法 |
Non-Patent Citations (3)
Title |
---|
DEY S.: "Free overall in open channels: state-of-the-art review", 《FLOW MEASUREMENT AND INSTRUMENTATION》 * |
陈小攀等: "基于堰流的矩形明渠跌水量水研究", 《灌溉排水学报》 * |
陈小攀等: "矩形明渠跌水量水试验研究", 《四川大学学报( 工程科学版)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111965720A (zh) * | 2020-08-19 | 2020-11-20 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | 一种基于地-井联合获取水力传导系数的方法 |
CN111965720B (zh) * | 2020-08-19 | 2023-05-23 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | 一种基于地-井联合获取水力传导系数的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109724570B (zh) | 2020-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106884405B (zh) | 一种无资料地区溃堤型山洪灾害分析评价方法 | |
Huang et al. | A physical similarity model of an impulsive wave generated by Gongjiafang landslide in Three Gorges Reservoir, China | |
Cheng et al. | Flood risk zoning by using 2D hydrodynamic modeling: A case study in Jinan City | |
Mansuri et al. | Effect of location and angle of cutoff wall on uplift pressure in diversion dam | |
Wu et al. | Nonmonotonic and spatial‐temporal dynamic slope effects on soil erosion during rainfall‐runoff processes | |
Wang et al. | Hydraulic barrier function of the underground continuous concrete wall in the pit of subway station and its optimization | |
Zhang et al. | Effects of initial step height on the headcut erosion of bank gullies: A case study using a 3D photo-reconstruction method in the dry-hot valley region of southwest China | |
Shrestha et al. | Model-based estimation of land subsidence in Kathmandu Valley, Nepal | |
Li et al. | Seepage analysis in a fractured rock mass: The upper reservoir of Pushihe pumped-storage power station in China | |
Lee et al. | Study on inlet discharge coefficient through the different shapes of storm drains for urban inundation analysis | |
Armanuos et al. | Impact of inclined double-cutoff walls under hydraulic structures on uplift forces, seepage discharge and exit hydraulic gradient | |
CN112685874A (zh) | 一种降雨入渗下非饱和黄土边坡渗流控制方程计算方法 | |
Rodríguez et al. | Three‐dimensional flow in centered pool‐riffle sequences | |
CN106680454A (zh) | 一种具拦沙坝已治理崩岗土壤侵蚀模数测算方法 | |
Gartner | Stream power: Origins, geomorphic applications, and GIS procedures | |
CN101845815A (zh) | 一种软土地层含承压水基坑突涌塑性破坏计算处理方法 | |
CN108643246B (zh) | 基坑施工中承压水降水对周边环境影响的预测方法 | |
CN109724570A (zh) | 地下跌水的跌水量、跌水宽度、坎上水层厚度的计算方法 | |
Zhu et al. | Experimental and numerical investigations of influence on overland flow and water infiltration by fracture networks in soil | |
Alahiane et al. | Practical method proposed to estimate silting’s rate in small and hillside dams | |
Adib et al. | Flood inundation modeling for Kota Tinggi catchment by combination of 2D hydrodynamic model and flood mapping approach | |
Thomas | Simulating the hydrologic impact of distributed flood mitigation practices, tile drainage, and terraces in an agricultural catchment | |
Aberle et al. | Directional dependency of flow resistance in an unlined rock blasted hydropower tunnel | |
Wang et al. | A new method for determining weathering rates in weathering pits | |
CN105912861B (zh) | 基于gis的非参数概念性水文模型的计算方法 |
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 |