土力学1第二章土的渗透性和渗流问题.ppt

上传人:sccc 文档编号:5298850 上传时间:2023-06-23 格式:PPT 页数:91 大小:6.16MB
返回 下载 相关 举报
土力学1第二章土的渗透性和渗流问题.ppt_第1页
第1页 / 共91页
土力学1第二章土的渗透性和渗流问题.ppt_第2页
第2页 / 共91页
土力学1第二章土的渗透性和渗流问题.ppt_第3页
第3页 / 共91页
土力学1第二章土的渗透性和渗流问题.ppt_第4页
第4页 / 共91页
土力学1第二章土的渗透性和渗流问题.ppt_第5页
第5页 / 共91页
点击查看更多>>
资源描述

《土力学1第二章土的渗透性和渗流问题.ppt》由会员分享,可在线阅读,更多相关《土力学1第二章土的渗透性和渗流问题.ppt(91页珍藏版)》请在三一办公上搜索。

1、清华大学土木水利学院 岩 土 工 程 研 究 所,张 丙 印,土力学1之第二章,土的渗透性和渗流问题,本章提要 本章特点 学习要点,第二章:土的渗透性和渗流问题,土的渗透性和渗透规律平面渗流及流网渗透力与渗透变形 有严格的理论(水流的一般规律)有经验性规律(散粒多孔介质特性)注意对物理概念和意义的把握注意把握土是散粒多孔介质这一特点,作业,本章作业:2-1 2-32-42-52-7,2.1 概述2.2 土的渗透性与渗透规律 2.3 平面渗流与流网2.4 渗透力与渗透变形,第二章:土的渗透性和渗流问题,2.1 概述,土体中的渗流,土是一种碎散的多孔介质,其孔隙在空间互相连通。当饱和土中的两点存在

2、能量差时,水就在土的孔隙中从能量高的点向能量低的点流动,水在土体孔隙中流动的现象称为渗流 土具有被水等液体透过的性质称为土的渗透性,透水层,不透水层,土石坝坝基坝身渗流,2.1 概述-渗流问题,防渗体,坝体,渗流问题:1.渗流量?2.渗透破坏?3.渗透力?,工程实例,浸润线:水从土坝(或土堤)迎水面,经过坝体向下游渗透所形成的自由水面和坝体横剖面的相交线。,板桩围护下的基坑渗流,2.1 概述-渗流问题,渗流问题:1.渗流量?2.渗透破坏?3.渗水压力?,透水层,不透水层,基坑,板桩墙,工程实例,渗流问题:1.渗流量Q?2.降水深度?,透水层,不透水层,天然水面,水井渗流,漏斗状潜水面,Q,2.

3、1 概述-渗流问题,渗流问题:1.渗流量?2.地下水影响 范围?,渠道、河流渗流,2.1 概述-渗流问题,原地下水位,渗流时地下水位,2.1 概述-渗流问题,降雨入渗引起的滑坡,渗流问题:1.渗透力?2.入渗过程?,事故实例,渗流量扬压力渗水压力渗透破坏渗流速度渗水面位置,挡水建筑物 集水建筑物 引水结构物 基础工程地下工程边坡工程,渗透特性变形特性强度特性,土的渗透特性,2.1 概述-土渗流特性,2.1 概述 2.2 土的渗透性与渗透规律 2.3 平面渗流与流网2.4 渗透力与渗透变形,第二章:土的渗透性和渗流问题,2.2 土的渗流性与渗透规律,水头与水力坡降土的渗透试验与达西定律渗透系数的

4、测定及影响因素层状地基的等效渗透系数,土的渗透性与渗透规律,渗流的驱动能量,反映渗流特点的定律,土的渗透性,地基的渗透系数,位置水头:到基准面的竖直距离,代表单位重量的液体从基准面算起所具有的位置势能压力水头:水压力所能引起的自由水面的升高,表示单位重量液体所具有的压力势能测管水头:测管水面到基准面的垂直距离,等于位置水头和压力水头之和,表示单位重量液体的总势能在静止液体中各点的测管水头相等,2.2 土的渗流性与渗透规律,位置、压力和测管水头,2.2 土的渗流性与渗透规律,水流动的驱动力,水往低处流,水往高处“跑”,位置:使水流从位置势能高处流向位置势能低处,流速:水具有的动能压力:水所具有的

5、压力势能,也可使水流发生流动,压力势能:,动能:,总能量:,称为总水头,是水流动的驱动力,单位重量水流的能量:,2.2 土的渗流性与渗透规律,水流动的驱动力-水头,渗流中的水头与水力坡降,2.2 土的渗流性与渗透规律,A,B,L,渗流为水体的流动,应满足液体流动的三大基本方程:连续性方程、能量方程、动量方程,2.2 土的渗流性与渗透规律,总水头:单位重量水体所具有的能量,位置水头Z:水体的位置势能(任选基准面)压力水头u/w:水体的压力势能(u孔隙水压力)流速水头V2/(2g):水体的动能(对渗流多处于层流0),渗流的总水头:,渗流问题的水头,也称测管水头,是渗流的总驱动能,渗流总是从水头高处

6、流向水头低处,A点总水头:,2.2 土的渗流性与渗透规律,水力坡降,B点总水头:,二点总水头差:反映了两点间水流由于摩阻力造成的能量损失,水力坡降 i:单位渗流长度上的水头损失,达西渗透试验,2.2 土的渗流性与渗透规律,1856 年达西(Darcy)在研究城市供水问题时进行的渗流试验,或:,其中,A是试样的断面积,达西定律,2.2 土的渗流性与渗透规律,达西定律:在层流状态的渗流中,渗透速度v与水力坡降i的一次方成正比,并与土的性质有关 层流是流体的一种流动状态。流体在管内流动时,其质点沿着与管轴平行的方向作平滑直线运动。渗透系数k:反映土的透水性能的比例系数,其物理意义为水力坡降i1时的渗

7、流速度,单位:cm/s,m/s,m/day渗透速度v:土体试样全断面的平均渗流速度,也称假想渗流速度,其中,Vs为实际平均流速,孔隙断面的平均流速,达西定律的适用范围,2.2 土的渗流性与渗透规律,适用条件:层流(线性流动),岩土工程中的绝大多数渗流问题,包括砂土或一般粘土,均属层流范围在粗粒土孔隙中,水流形态可能会随流速增大呈紊流状态,渗流不再服从达西定律。可用雷诺数进行判断:,Re5时层流Re 200时紊流200 Re 5时为过渡区,达西定律的适用范围,2.2 土的渗流性与渗透规律,在纯砾以上的很粗的粗粒土如堆石体中,在水力坡降较大时,达西定律不再适用,此时:,两种特例,对致密的粘性土,存

8、在起始水力坡降i0?,ii0,v=k(i-i0),渗透系数的测定方法,常水头试验法 变水头试验法,井孔抽水试验 井孔注水试验,2.2 土的渗流性与渗透规律,室内试验方法 野外试验方法,室内试验方法-常水头试验法,试验条件:h,A,L=const 量测变量:体积V,t,适用土类:透水性较大的砂性土,2.2 土的渗流性与渗透规律,室内试验方法-变水头试验法,试验条件:h变化 A,a,L=const 量测变量:h,t,适用土类:透水性较小 的粘性土,2.2 土的渗流性与渗透规律,土样,A,L,Q,水头测管,开关,a,在tt+dt时段内:,2.2 土的渗流性与渗透规律,室内试验方法-变水头试验法,室内

9、试验方法小结,2.2 土的渗流性与渗透规律,常水头试验,变水头试验,条件,已知,测定,公式,取值,h=const,h变化,h,A,L,V,t,重复试验后,取均值,a,A,L,h,t,不同时段试验,取均值,适用,粗粒土,粘性土,2.2 土的渗流性与渗透规律,现场测定法抽水试验,试验条件:Q=const 量测变量:r=r1,h1=?r=r2,h2=?,优点:可获得现场较为可靠的平均渗透系数 缺点:费用较高,耗时较长,2.2 土的渗流性与渗透规律,现场测定法抽水试验,是土中孔隙直径大小的主要影响因素因由粗颗粒形成的大孔隙可被细颗粒充填,故土体孔隙的大小一般由细颗粒所控制。因此,土的渗透系数常用有效粒

10、径d10来表示,如哈臣公式:,2.2 土的渗流性与渗透规律,土的性质 水的性质,粒径大小及级配孔隙比矿物成分结构,渗透系数的影响因素,是单位土体中孔隙体积的直接度量对于砂性土,常建立孔隙比e与渗透系数k之间的关系,如:,2.2 土的渗流性与渗透规律,渗透系数的影响因素,土的性质 水的性质,粒径大小及级配孔隙比矿物成分结构,对粘性土,影响颗粒的表面力不同粘土矿物之间渗透系数相差极大,其渗透性大小的次序为高岭石伊里石蒙脱石;当粘土中含有可交换的钠离子越多时,其渗透性将越低塑性指数Ip综合反映土的颗粒大小和矿物成份,常是渗透系数的参数,2.2 土的渗流性与渗透规律,渗透系数的影响因素,土的性质 水的

11、性质,粒径大小及级配孔隙比矿物成分结构,影响孔隙系统的构成和方向性,对粘性土影响更大在宏观构造上,天然沉积层状粘性土层,扁平状粘土颗粒常呈水平排列,常使得k水平k垂直在微观结构上,当孔隙比相同时,凝聚结构将比分散结构具有更大的透水性,2.2 土的渗流性与渗透规律,渗透系数的影响因素,土的性质 水的性质,粒径大小及级配孔隙比矿物成分结构,2.2 土的渗流性与渗透规律,渗透系数的影响因素,土的性质 水的性质,粒径大小及级配孔隙比矿物成分结构,2.2 土的渗流性与渗透规律,渗透系数的影响因素,水的动力粘滞系数:温度,水粘滞性,k饱和度(含气量):封闭气泡对k影响很大,可减少有效渗透面积,还可以堵塞孔

12、隙的通道,土的性质 水的性质,粒径大小及级配孔隙比矿物成分结构,2.2 土的渗流性与渗透规律,层状地基的等效渗透系数,等效渗透系数,确立各层土的ki 根据渗流方向确定等效渗流系数,天然土层多呈层状,多个土层用假想单一土层置换,使得其总体的透水性不变,2.2 土的渗流性与渗透规律,等效渗透系数:,层状地基的水平等效渗透系数,层状地基的垂直等效渗透系数,2.2 土的渗流性与渗透规律,等效渗透系数:,算例说明,按层厚加权平均,由较大值控制类似串联电路电阻?,层厚倒数加权平均,由较小值控制类似并联电路电阻?,层状地基的等效渗透系数,2.2 土的渗流性与渗透规律,层状地基的等效渗透系数,2.2 土的渗流

13、性与渗透规律,水平渗流情形,垂直渗流情形,条件,已知,等效,公式,2.2 土的渗流性与渗透规律,小 结,水头与水力坡降渗透试验与达西定律渗透系数的测定及影响因素层状地基的等效渗透系数,总水头=位置水头+压力水头 水头是渗流的驱动力,达西定律渗透系数、渗透速度达西定律的适用条件,常水头试验 变水头试验 抽水试验 渗透系数影响因素,水平等效渗透系数 垂直等效渗透系数,2.1 概述 2.2 土的渗透性与渗透规律 2.3 平面渗流与流网2.4 渗透力与渗透变形,第二章:土的渗透性和渗流问题,平面渗流的基本方程及求解 流网的绘制及应用,平面问题:渗流剖面和产生渗流的条件沿某一个方向不发生变化,则在垂直该

14、方向的各个平面内,渗流状况完全一致。对平面问题,常取dy=1m单位宽度的一片来进行分析,2.3 平面渗流与流网,h=h(x,z),v=v(x,z),与时间无关,稳定渗流:流场不随时间发生变化的渗流,平面稳定渗流,2.3 平面渗流与流网,渗流的连续性方程,单位时间流入单元的水量:,渗流的连续性方程:,单位时间内流出单元的水量:,连续性条件:,渗流的运动方程,达西定律:,渗流的连续性方程:,渗流的运动方程:,2.3 平面渗流与流网,特例:各向同性均质土体 kx=kz,Laplace方程,描述渗流场内水头的分布,是平面稳定渗流的基本方程,2.3 平面渗流与流网,H1,H2,2,3,4,渗流的边界条件

15、,课堂讨论:土石坝渗流问题的边界条件,渗流域,不透水层,2.3 平面渗流与流网,H1,H2,2,3,4,水头边界条件 在边界1上给定水头,2,流速边界条件 在边界2上给定法向流速,4,3,渗出面 在边界3上H=z,vn0,自由水面*在边界4上H=z,vn=0,渗流的边界条件,平面稳定渗流问题描述,运动方程:,边界条件:,水头边界条件 在边界1上给定水头,2,流速边界条件 在边界2上给定法向流速,4,3,渗出面 在边界3上h=z,vn0,自由水面*在边界4上h=z,vn=0,2.3 平面渗流与流网,或:,2.3 平面渗流与流网,数学解析法或近似解析法:求取渗流运动方程在特定边界条件下的理论解,或

16、者在一些假定条件下,求其近似解数值解法:有限元、有限差分、边界元法等,近年来得到迅速地发展电比拟试验法:利用电场来模拟渗流场,简便、直观,可以用于二维问题和三维问题流网法:简便快捷,具有足够的精度,可分析较复杂断面的渗流问题,渗流分析的方法,流速势或势函数:,则有:,满足达西定律的渗流问题是一个势流问题,势 函 数,势函数的特性:等势面是等水头面 两条等势面的势值差同其水头差成正比,2.3 平面渗流与流网,存在函数,称为流函数,且有:,流线的方程:,流 函 数,2.3 平面渗流与流网,定义:流线是流场中的曲线,在这条曲线上所有各质点的流速矢量都和该曲线相切,即表示该点的渗流的方向,亦即一颗颗渗

17、流水流经的轨迹,设上式右边为某一函数全微分,即,性质一:流线互不相交,在同一条流线上,流函数的值为一常数,流函数的性质,2.3 平面渗流与流网,性质二:两条流线流函数的差值等于其间通过的流量,a和b为两流线间的过水断面,a(x,z),b(x-dx,z+dz),2.3 平面渗流与流网,1)势函数和流函数均满足拉普拉斯方程2)势函数和流函数正交,一点两线的斜率互成负倒数3)势函数和流函数是互为共轭的调和函数,两者均完备 地描述了同一个渗流场4)当对调边界条件时,势函数和流函数两组曲线可互换,势函数与流函数,流网及其特性,流线和等势线正交流网中应使相邻流线间的流函数差和相邻等势线间的势函数(水头)差

18、不变流网中每一网格的边长比为常数,通常取为1,在流场中,流线和等势线(等水头线)组成的网格称为流网,+d,+d,v,s,l,2.3 平面渗流与流网,流网的画法,1)确定边界条件:边界流线和首尾等势线2)研究水流的方向:流线的走向3)判断网格的疏密大致分布4)初步绘制流网的雏形:正交性、曲边正方形5)反复修改和检查,2.3 平面渗流与流网,要点:边界条件、正交性、曲边正方形、多练习,a,b,c,d,e,f,g,h,流网的应用,2.3 平面渗流与流网,测管水头 h,确定孔压,确定流速,确定流量,水力坡降,流道数,小 结,平面渗流的基本方程及求解流网的绘制及应用,连续性方程 运动方程 边界条件 求解

19、方法,势函数 流函数 流网及特性 流网的画法 流网的应用,2.3 平面渗流与流网,2.1 概述 2.2 土的渗透性与渗透规律 2.3 平面渗流与流网 2.4 渗透力与渗透变形,第二章:土的渗透性和渗流问题,渗透力临界水力坡降渗透变形(渗透破坏),2.4 渗透力与渗透变形,渗透力-试验观察,h=0 静水中,土骨架会受到浮力作用。h0 水在流动时,水流受到来自土骨架的阻力,同时流动的孔隙水对土骨架产生一个摩擦、拖曳力。,渗透力j:渗透作用中,孔隙水对土骨架的作用力,方向与渗流方向一致,2.4 渗透力与渗透变形,渗透力-试验观察,土粒,渗 流,渗透力 j:体积力,渗透力j:单位土体内土骨架所受到的渗

20、透水流的拖曳力,2.4 渗透力与渗透变形,渗透力-受力分析,W=L sat L(+w),P1=whw,P2=wh2,R=?,R+P2=W+P1,R+wh2=L(+w)+whw,R=L,土水整体受力分析-静水,2.4 渗透力与渗透变形,渗透力-受力分析,W=L sat L(+w),P1=whw,P2=wh1,R=?,R+P2=W+P1,R+wh1=L(+w)+whw,R=L-wh,土水整体受力分析-渗流,2.4 渗透力与渗透变形,渗透力-受力分析,R=L-wh,土水整体受力分析-对比,静水中的土体,渗流中的土体,向上渗流存在时,滤网支持力减少,R=L,减少的部分由谁承担?,总渗透力:J=wh,渗

21、透力j:单位土体内土骨架所受到的渗透水流的拖曳力,j=J/V=wh/L=wi,向上渗流存在时,滤网支持力减少。当滤网支持力为零时的水力坡降称为临界水力坡降icr,它是土体开始发生流土破坏时的水力坡降:,2.4 渗透力与渗透变形,渗透力-受力分析,渗透力-受力分析,R=L-wh=0,临界水力坡降,icr=h/L=/w,土水隔离受力分析,2.4 渗透力与渗透变形,渗透力-受力分析,渗透力-受力分析,R=L-wh,土骨架受力分析:有效重量:W=L总渗透力:J=Lj滤网的反力:R,孔隙水受力分析:水压力:P1=whw P2=wh1总渗透力:J=J水重+浮力反力:Ww=Vvw+Vsw=Lw,孔隙水受力平

22、衡,j=wi,土骨架受力平衡,渗透力的性质,物理意义:单位土体内土骨架所受到的渗透水流的拖曳力,它是一种体积力大小:j=wi方向:与水力坡降方向一致作用对象:土骨架,2.4 渗透力与渗透变形,2.4 渗透力与渗透变形,土工建筑物及地基由于渗流作用而出现的变形或破坏称为渗透变形或渗透破坏。渗透变形是土工建筑物发生破坏的常见类型基本类型:管涌 流土 接触流土 接触冲刷,渗透变形,单一土层渗透变形的两种基本型式,2.4 渗透力与渗透变形,渗透变形-流土,流土:在向上的渗透作用下,表层局部范围内的土体或颗粒群同时发生悬浮、移动的现象。任何类型的土,只要水力坡降达到一定的大小,都可发生流土破坏,渗流,2

23、.4 渗透力与渗透变形,渗透变形 管涌,原因,内因:有足够多的粗颗粒形成大于细粒直径的孔隙外因:渗透力足够大,在渗流作用下,一定级配的无粘性土中的细小颗粒,通过较大颗粒所形成的孔隙发生移动,最终在土中形成与地表贯通的管道,渗流,过程演示,1.在渗透水流作用下,细颗粒在粗颗粒形成的孔隙中移动流失,2.孔隙不断扩大,渗流速度不断增加,较粗颗粒也相继被水带走,3.形成贯穿的渗流通道,造成土体塌陷,2.4 渗透力与渗透变形,流土与管涌的比较,流土,土体局部范围的颗粒同时发生移动,管涌,只发生在水流渗出的表层,只要渗透力足够大,可发生在任何土中,破坏过程短,导致下游坡面产生局部滑动等,现象,位置,土类,

24、历时,后果,土体内细颗粒通过粗粒形成的孔隙通道移动,可发生于土体内部和渗流溢出处,一般发生在特定级配的无粘性土或分散性粘土,破坏过程相对较长,导致结构发生塌陷或溃口,2.4 渗透力与渗透变形,Fs:安全系数1.52.0,i:允许坡降,i icr:土体处于稳定状态,i=icr:土体处于临界状态,i icr:土体发生流土破坏,工程设计:,流土可能性的判别,在自下而上的渗流逸出处,任何土,包括粘性土和无粘性土,只要满足渗透坡降大于临界水力坡降这一水力条件,均要发生流土:,2.4 渗透力与渗透变形,土是否会发生管涌,取决于土的性质:粘性土(分散性土例外)属于非管涌土无粘性土中发生管涌必须具备相应的几何

25、条件和水力条件,管涌可能性的判别,2.4 渗透力与渗透变形,较均匀土(Cu10),几何条件 水力条件,无粘性土管涌的判别,级配,孔隙及细粒,判定,非管涌土,粗颗粒形成的孔隙小于细颗粒,不均匀土(Cu10),不连续,连续,d0=0.25d20,细粒含量35%,细粒含量25%,细粒含量=25-35%,d0 d3,d0 d5,d0=d3-d5,管涌土,过渡型土,非管涌土,非管涌土,管涌土,过渡型土,P(%),lgd,骨架,充填料,发生管涌的必要条件:粗颗粒所构成的孔隙直径大于细颗粒直径,2.4 渗透力与渗透变形,几何条件 水力条件,无粘性土管涌的判别,渗透力能够带动细颗粒在孔隙间滚动或移动。可用管涌

26、临界水力坡降表示,0 5 10 15 20 25 30 35,1.51.00.50,icr,Cu,流土,过渡,管涌,伊斯托敏娜(苏),中国学者,Cu 20时,icr=0.25-0.30,考虑安全系数后:i=0.10-0.15,渗透变形的防治措施,减小i:上游延长渗径 下游减小水压增大i:下游增加透水 盖重,2.4 渗透力与渗透变形,改善几何条件:设反滤层等 改善水力条件:减小渗透坡降,防治流土,防治管涌,2.1 概述 2.2 土的渗透性与渗透规律 2.3 平面渗流与流网 2.4 渗透力与渗透变形,第二章:土的渗透性和渗流问题,小 结,2.4 渗透力与渗透变形,工程实例渗流问题,土的渗透性及渗透

27、规律,二维渗流及流网,渗透力与渗透变形,渗流中的水头与水力坡降渗透试验与达西定律渗透系数的测定及影响因素层状地基的等效渗透系数,平面渗流的基本方程及求解流网的绘制及应用,渗透力:概念与计算渗透变形:类型、条件、防治,本章作业:2-1 2-32-42-52-7,请注意对第一章和第二章的复习!,作业,土坝,高90m,长1000m,1975年建成,次年6月失事渗透破坏:冲蚀 水力劈裂,Teton坝失事现场现状,原因,土石坝坝基坝身渗流破坏实例,失事原因研究结论,土石坝坝基坝身渗流破坏实例,九江大堤决口,1998年8月7日13:10发生管涌险情,很快形成宽62m的溃口,堤基管涌,焦点词汇:豆腐渣工程,

28、原因,土石坝坝基坝身渗流破坏实例,沟后面板砂砾石坝,位于青海省,高71米,长265米,建于1989年。1993年8月7日突然发生溃坝,是现代碾压堆石坝垮坝的先例。,溃坝原因:面板止水失效,下游坝体排水不畅,造成坝坡失稳,土石坝坝基坝身渗流破坏实例,广州京广广场基坑塌方,基坑渗流破坏,珠海祖国广场基坑失事,基坑渗流破坏,西藏易贡巨型滑坡,时间:2000年4月9日约20时规模:滑坡体自相对高差近3330m的雪峰阳坡滑下,历时约10分钟,滑程8km。堆积体长、宽各约2500m,平均厚60m,最厚100m,体积约2.8亿-3.0亿m3。地质:滑坡堆积体80%以上是砂性土险情:堵塞易贡藏布江成堰塞湖,湖

29、水面积22km2,湖长17km,水位以每天0.5-0.6m的速度上涨,湖水无下泄通道,预计6月底湖水将上涨至堆积体顶,拦存湖水将达40亿-60亿m3,降雨入渗引起的滑坡,西藏易贡巨型高速滑坡,降雨入渗引起的滑坡,西藏易贡巨型高速滑坡,降雨入渗引起的滑坡,湖水每天上涨50cm!,易贡巨型滑坡现场,降雨入渗引起的滑坡,预案一:加强监测,上、下游移民工作,库满自溢漫顶溃口预案二:在堆积体最低处开渠引流,水库溢流漫顶溃口。预案三:在右岸山体垭口开溢洪道,改造堆积体成坝,堰塞湖成库,易贡滑坡处理预案,最终采用预案二,降雨入渗引起的滑坡,武警部队等700多名抢险人员,奋战33天,累计开挖土石方135.5万m3,有效降低了堆积体过水高程24.1m,减少拦存湖水约20亿m3施工队伍6月4日撤离;6月8日6时40分,泄水渠过水,至11日2时50分堆积体溃决;11日21时,滑坡体拦存的湖水按预定方案完全下泄 泄洪时易贡湖库容量已达30多亿立方米,易贡滑坡处理结果,降雨入渗引起的滑坡,

展开阅读全文
相关资源
猜你喜欢
相关搜索
资源标签

当前位置:首页 > 建筑/施工/环境 > 农业报告


备案号:宁ICP备20000045号-2

经营许可证:宁B2-20210002

宁公网安备 64010402000987号