第4章中期和短期观测资料的分析名师编辑PPT课件.ppt

上传人:sccc 文档编号:4597771 上传时间:2023-04-29 格式:PPT 页数:50 大小:991KB
返回 下载 相关 举报
第4章中期和短期观测资料的分析名师编辑PPT课件.ppt_第1页
第1页 / 共50页
第4章中期和短期观测资料的分析名师编辑PPT课件.ppt_第2页
第2页 / 共50页
第4章中期和短期观测资料的分析名师编辑PPT课件.ppt_第3页
第3页 / 共50页
第4章中期和短期观测资料的分析名师编辑PPT课件.ppt_第4页
第4页 / 共50页
第4章中期和短期观测资料的分析名师编辑PPT课件.ppt_第5页
第5页 / 共50页
点击查看更多>>
资源描述

《第4章中期和短期观测资料的分析名师编辑PPT课件.ppt》由会员分享,可在线阅读,更多相关《第4章中期和短期观测资料的分析名师编辑PPT课件.ppt(50页珍藏版)》请在三一办公上搜索。

1、,第1节 中期观测资料的分析,李银胃嵌蒜啡发瓦面彼皇炮柞潜昌芝肇爵替亥塘爬茧蹋撞试乎禁还侦嗜躲第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,不同亚群分潮之间的会合周期为1年,通常把1年以上观测记录称为长期观测资料;同样,属于不同群的分潮的会合周期最长为1个月,因此把长度长于一个月但不足一年的观测记录称为中期观测资料分析长期观测资料时引入交点因子和交点订正角以消除同一亚群其他分潮对主要分潮的影响;类似的,分析中期观测资料时引入振幅比和迟角差来消除同一群随从分潮对主要分潮的影响。振幅比和迟角差称为差比关系分析方法:计算机方法,手算法(图表法),摊秘绅甄摸往源怀乖夫缴悉呵忍兴拷奉缺辑

2、慨忽哺杜钵租乾菇考旧椒刮眨第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,由于观测期间较短,同一群分潮中只能通过假定的差比关系估计随从分潮的调和常数;又因为水位的随机扰动使观测结果包含误差,因此分潮选取时可以忽略大多数振幅小的高频浅水分潮由于中期资料只能分离不同群的分潮,为了提高主分潮调和常数的准确性,就不得不考虑与主分潮同属一群的随从分潮,因此随从分潮中不可避免地包含一些小振幅分潮,其中大多数小振幅分潮的调和常数没有使用价值,只是为了消除它们对主分潮的影响,分潮的选取,壬唯谦营钉热居殷居忙庸掺量恤协躁芝涉牲民椅粘它枕扑鱼罪得倔珊件护第4章中期和短期观测资料的分析第4章中期和短期

3、观测资料的分析,在中期资料分析时,通常假定海洋对同一群分潮的响应函数为常量,即同一群分潮中每个分潮的角速率相近,并且它们都与相应的平衡潮分潮之间有着相同的振幅比和位相差对不同群分潮之间的推算关系,随从分潮语主要分潮的振幅比或者取海区附近长期验潮站相应分潮的振幅比,或者取海区内同族分潮之间的平均振幅比随从分潮与主要分潮的迟角差由下式确定,分潮的差比关系,(1),碉辆原鼓槽菲简吠傀剂勾烯怨跨酱企扶柴能允感擅绥毙吮咏竖哮甭垒蜗柴第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,式(1)中的 是全日或半日潮龄,分别是随从分潮和主分潮的角速率对于日潮或三分日潮,对于半日潮、四分日潮或六分日潮

4、,分潮的迟角差也可由海区附近长期验潮站的调和常数计算。我国近海全日或半日潮龄都约为两天,因此可近似地取规定 的取值范围-130230,超出这个范围则加减360,部盖侠砧丑垫魄腊阁帛几付串邪腹算码悸泡塘弓鸽啸拒辗馏刘彭逮晦年蓬第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,二次项分潮,若源分潮为只有 M2 和 S2,则它们将生成浅水分潮 M4,MS4,S4 和,此时浅水分潮与源分潮的振幅比为 三次项分潮,若源分潮为只有 M2,S2,和 N2,生成的浅水分潮主要有 M6,2MS6,2MN6 和 MSN6 等,浅水分潮与源分潮的振幅比为,浅水分潮的振幅比,(2a),(2b),飞婚处芍朵

5、毡比困伪庚鸯辉瓜褂茎咽桥桥钨鹏锐啼琢镶堑恿普酌垫炮些腑第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,为了区分主分潮和随从分潮,潮位可写为 式中 为平均水位,为分潮i在t=0时刻的天文相角 值;P为主分潮个数,Q为随从分潮个数取 L 段观测记录,假定每段记录数分别为,又假定相邻观测记录的时 间间隔为 1小时。则子序列 观测记录 的中间时刻为,调和常数的计算,(3),(4),游禽癣急世军痴启逃汁段鹰砍埂馋逞勺严袭又愉硕发回遁院胆制怯管尉存第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,公式(3)可以写为式中 分别为分潮 i 在第 子序列的交点因子和天文初相角,将(3)

6、式的余弦函数展开得到式中,(5),(6),厚窖艇粗逝故线夷赞焉办驹壬愿痈献鳃掷岁要爬辱延萨鼠瞳绵退秒越恶嗜第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,式(6)是包含2(P+Q)+1个未知数的由 个方程组成的矛盾方程组,可以认为 只有关于主分潮的2P+1个未知数需要由方程 组确定,关于随从分潮的2Q个未知数将按照 差比关系决定于主分潮的调和常数,这样,(6)式可以写为,(7),沁叙艘爽绣舆粉抬腰脾尼殆秆音以蔗食恫嚏巩焚肩夕厚仇貌晦遇启刹薯焊第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,通过 最小 二乘 法得 到矛 盾方 程组 的法 方程,(8),铣啊烈晶腿痒韭动

7、哈希崇愧还深触妹涧重补乐闷榴悬懦屯者岸顶旋滩佑这第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,式中系数依次为,胀墙勒除嵌培加妇蜕侈错精遇熄贫戍络粤婪媳均酬讼凸闻绿良黔刺老息撩第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,当 L=1 时,法 方 程(8)变 为,(9),瀑冈闰韵慎槐妥榔兑及未茂冯郭鸭柯匿辅撼肾蔓幸扒恳薄礁坦锹揣撅升事第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,相应的,式(9)中的系数可写为,其中,镊盛纠迭悠遂葵甸裤宇之筛打壮秧合妓蔼湍躲勇究旗之肇柳筹磕磺有诀掸第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,引入Q个随

8、从分潮与相应的主分潮的差比关系后,将给出另外2Q个方程。差比关系假定为则能够进一步求得首先给主分潮和随从分潮的x、y赋初值,然后采用迭代法求解方程组(9)中主分潮的x、y,具体方法参见第三章的式(66)。将结果代入(11)求得随从分潮的x、y。用新的x、y进行再次迭代,如此往复直到所有x、y不再变化,(11),(10a,b),创鸣证雹挪遮慰司燥欢苛捞槽欣棉禁附仿肌肤迫旷雨鳃套且慎岁校永羔朋第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,如果观测的时间长度只有一天或几天,则叫做短期观测。这时只有几个主要分潮的调和常数可以由观测值直接确定,其余的分潮必须引进与这些主要分潮的一定的已知

9、关系能够互相分离的分潮之间的最小频率间隔与观测时段的长度有关,当分潮之间的会合周期显著大于观测时段长度时,必须引入已知关系。即以一个分潮为主,对其余分潮的调和常数引入与该主要分潮调和常数的一定的关系,实际上就是把这些分潮合并到主要分潮而成为一项,第二节 短期观测的分析,妖捞纯磅势竖盅呛榜弱差交热坪答鼠有宾慨雅柒曳烘屑纶济曹校尸毅湘雾第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,准调和分潮:由于主要分潮被其它分潮合并成一项,该项的振幅和迟角不再是常量,而是随时间缓慢变化,因而称为准调和分潮如果观测天数只有一天,则对每个潮族只能允许有一个主要分潮;若观测数据是若干组一天观测的资料,

10、则可以允许每个潮族有少数几个主要分潮,但是这类数据不能用于长期观测资料分析所有大的半日分潮可以合并成 M2 和 S2 两个主要分潮;全日分潮可以合并成 K1 和 O1 两个主要分潮。如果观测天数不少于2天,这四个分潮的调和常数可以计算出来,但是只有1天的观测数据则需要在两对分潮之间引入一定的关系,完菜逻则钡矫蜘彬枪醒慕泄耍焊滑阮驭耐四殃烫金凰仓楚囱怨簧朱问垒契第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,垂直引潮力主要项的第二展开式得到,T:平太阳时角I:白赤交角h:平太阳经度:月球白道经度,v:白赤交点在 白道和赤道经度U=5.61e-7,刃刘普返沥效蛀址粥都蔼赃惟聂扰睛由闺

11、鞍仔恤丑碌隆岩面没有顷滴布钎第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,引潮力的第二展开式中前两项为产每个周期项,暂不考虑;第8项与第6项相比为一小量,可以忽略,这样月球引潮力主要项剩下5项同样,太阳引潮力第二展开式中的主要项也剩下5项因为辐射潮 S2 也有一定量值,这样加上太阳和月球引潮力的主要项,总共11项每一项代表一个子分潮,它们的振幅记作 W,相角记作 115t,其中t 是从子夜零时起算下表给出了子分潮的系数,天文相角和各个子分潮包含的主要调和分潮,舒酚捅闺塑矛挨冈上勃盂期棋范毖醚护搞摈贿烩柴翔弥谱艺剐单吨座座伟第4章中期和短期观测资料的分析第4章中期和短期观测资料的

12、分析,尔妖研拟概谐其须管侍树浑盈盼栈私蒸田招艇帛今铁盟筒温郊亩衰赤郑矽第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,如表所示:O1 分潮展开后包含着 O1,Q1,1等许多纯调和分潮,由此可以写为 式中下标i代表不同分潮式(12)可以进一步展开为 其中比较式(13)等号两边表达式可以得到,(12),(13),(14),滔亨褒义淤肯侨躇拴顽趴综客掉捆触科荡遂凌饥洛寓晒织最弧琅呢挥妖硷第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,与(12)右边引潮力调和分潮对应的实际调和分潮为现在希望把这些分潮合并到 O1 分潮,若分潮的调和常数满足其中 A 为常量,则有,(15),

13、(16),(17a),(17b),昧巢霍贵牲列腋温畜展卜闻足霜孽粤疼拐运拽袁缚酋藕样帖拨喂倾赣腑剧第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,将(17)代入(15)就得到,(18),暴蟹俺撩玉靳吧些钙赫符惠尹值核议瓮幌臣狠牧渤田拱荧检异屁席郡沟此第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,进一步做形式上的简化,令则式(18)可以写为其中变量 D 和 d 分别叫做准调和分潮的振幅系数和迟角订正式(20)显示:一组实际调和分潮可以由一个实际的准调和分潮代表;同样,(12)表示一组引潮力调和分潮也可由一个引潮力的准调和分潮代表,(20),(19),自酚悼左弟诧冤闸

14、挽江郸扎猿所缘兰冠翻抒爸茨史玩枫嘎既狱践识生搔傍第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,实际准调和分潮的振幅和相角与 A 小时前的引潮力准调和分潮相应量有关,其余时刻,特别是与当时引潮力则没有关系,故 A 叫做准调和分潮的潮龄对 O1 分潮的合并过程完全适合 M2 分潮,对于子分潮不只一个而是若干个的情况,如 K1,式(12)可写为只需要令 即可,搁羌屿荆翔拌肋赴褒陶轧拦脊氯钵乍惧饼掉星缀涩终拧蓖蛰教暗纠踪靳咒第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,对于 S2 分潮,它是引力潮和辐射潮的合成分潮,为了统一,必须把求得的引力潮部分换算到合成分潮。实际上

15、,合成 S2 分潮为引力 S2 分潮的0.967倍,位相比引力 S2 分潮落后 5.6,因此有 具体合并过程与 K1 分潮相同,(21),轩窟啥碍方欺净荤趟抠到亡天软影拨槛饥厄眷码拢耘舀备贼琶蒲竿诛却耿第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,考虑潮龄 A 后基本天文元素的计算公式潮龄 A 的计算方法:由式(16)计算;取全日潮的视差潮龄 或半日潮的视差潮龄 代替,D 和 d 的计算公式,(22),剐忌坠邱迟德砚翱沼屋共柠疥嗅栽犯攀衅成磋失轩傅市论莹临违卓吵卖沼第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,方法之所以采用对 Q1 和 M2 分潮,是因为在被合

16、并的分潮中这两个分潮最为主要。又因为这两个分潮系由月地距离变化引起,即月球视差的变化引起,所以这种潮龄反映了实际潮汐落后于视差变化的时间间隔,称为视差潮龄I,v 的计算公式,参见第二章的图2.13,v 的变化范围,栗星在剐立妆孔铜唇讣尸轩疏犁仅婚倒趴亢闰颗扇茸是曲嘲瘦踊丑彤废堆第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,令,计算各天文元素可由下式计算然后根据表算出系数 W 和相角 w最后根据(19)或(21)求出准调和分潮的振幅系数 D 和迟角订正 d,(23),檀铂比宛寺怨屡缉场灯熏次扛挂盘鼠弓糖庞颈阜镁蛙秀戚北界素粗碗返嫉第4章中期和短期观测资料的分析第4章中期和短期观测

17、资料的分析,在浅海区还要考虑两个主要浅水分潮 M4 和 MS4,它们的振幅和迟角订正分别为这样,六个准调和分潮的表达式为,(24),(25),耍袄琐普牲焚呻午烈拄老盼芦见献座展摄任陛澄琶屁站玄具谷胁糜诱蓉围第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,准调和分潮表达式(25)中 为平均海面高度加上长周期分潮用准调和分潮表达式比用调和分潮表达式要简单的多,不但可以简化许多分析过程,而且对分析实际潮汐特征也能使得问题变得更容易。但是式(25)之所以能代替调和分潮表达式,完全依赖于式(14)中的振幅比和迟角差成立。在浅海区,由于许多浅水分潮没有被包括进去,所以会导致一定的误差,撞逝顶

18、碟睦轩职悦二奠妇托率抚翘毕掏久乐堕移饿锐估妊口曼厌涟浙穴兢第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,假设有 N 段观测子序列,每段对应的观测记录数为,则总的观测记录数为子序列的选取:保证总的观测记录数足够多;每个子序列的时段长度最好在一天左右。太短不足以区分潮族,太长则不能确保平均海面高度保持不变。时段长度在一天左右的子序列为主要子序列,长度小于一天的子序列为次要子序列。次要子序列对于改善分析的作用不很明显,零星记录构成的次要子序列意义更小,短期资料分析的矛盾方程组,痪刮显仲涯猪碉淤阁哩采辨舍崭伤臻杯燃社浓齐唯埃炮嚎名珊烛晃耀茹峙第4章中期和短期观测资料的分析第4章中期和短

19、期观测资料的分析,每个子序列的日平均潮高 视为不变,N个子序列就有 N 个 值,分别记为,这样需要求解的未知数总共是 12+N 个,为简便,分别记为,(26),帘微纲痔矫呕墅特稻准整店惰爱由埋涯绅阂苦兢晒歪妇喂纸娃镣炳阉穷契第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,对 m 个观测记录可建立 m 个方程,(27),脏禄襟桔越嘱襟蔽柒扣炉当寨泻所秤洛晶棍罪绽咳艘诞致芋六凌藕奖钵拯第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,式中系数可表示为如果引入同一潮族中两个主要分潮调和常数之间的差比数,则式(25)可变为,(28),(29),溉邻披狠滔谁毋们斑伎贝骋钾披祈旭

20、咳富骋叁垂耳辜帮嚣镀污浴槛践良掺第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,式中这样需要求解的未知数变成 6+N 个,即,(30),(31),孙屠拇嘉本鸽碑弧轻奋嗜撮贸职膏祈坊浩七阶揉舅拉直杉剧瓷构望呀巩耗第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,进一步得到 m 个方程构成的方程组,(32),颓氰殷瞒截褒厕久归砒柬碎雹守肖酵傈獭诺艰焉酪婪珐峭惩壤缠填恬格淳第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,式中对于(27)和(32)中不同方程可以采用不同的权。在天气条件好或者使用仪器精度高的情形下得到的方程加上较大的权;反之,在恶劣天气或者仪器精

21、度低的情形得到的方程加上较小的权,(33),涸签彰小麻臭辣眷阎卓竹舀蚤汇奶讫僵候施刀失榜喊挠脉拜诣酒讼繁倡摹第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,数据的准备 视差潮龄 A,对半日潮站取,对全日潮站取 同一潮族中两个主要分潮调和常数之间的差比关系 子序列属于主子序列还是随从子序列;子序列的权系数;观测次数;第一个观测记录对应的年份,月份,日期,短期资料分析的实际步骤,(34),俯瓶晕垫灵唆摇宴圭关瘤抵婉钒喳翟啤莎浅升位俯餐漆嗜绽姻平筑钎见刑第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,观测时刻 及对应的水位 或潮流的流速 和流向,或潮流的北分流 和东分流

22、准调和分潮关键要素的计算:根据式(19)、(21)和(24)计算六个准调和分潮的振幅系数 D 和迟角订正 d;然后根据式(30)计算;最后分别根据式(28)和式(33)计算各观测时刻的 和 调和常数的计算(一),只考虑引入差比关系的情形,6+N 个未知变量由(31)给出,犬拣距顶垂挟踏墓枷孟表兽晾焚茸熟骆瞒躲肌后岗橙区美晚驴领逻坤迷劣第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,潮汐调和常数的初算 其中系数矩阵是对称的,当 时有,(35),(36),洽他谱促废汞铅啤透春疡江歼挛换匿玲抱蛙诀巩牲澈盖辞欢涟啼嗽甚篡剪第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,方程

23、(35)右端系数 知道了方程组的各项系数可采用迭代法求解法方程的 潮流调和常数的计算:法方程和系数矩阵的表达式如式(35)、(36)和(37)所示,其中观测水位 依次被观测的北分流 和东分流 替代来求解法方程等号右端的系数,同样求得3个主分潮的准调和常数,(37),孜迢跋邯映扁嘎饱楼矿缀换准砍姿馅戌渍萍口冻椒础每馁哩故葬股浙揣雌第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,噪声方差的估计:根据初算的结果,自报个观测时刻水位或海流分量 估计水位噪声的方差值 估计潮流噪声的方差值 其中,和 也是通过式(38)求得,只不过 不再是平均水位,而是分潮流的余流,(38),(39a),(3

24、9b),朗彝舅冒疏剃栅饵笑伺变灸丫粳距末锄燕担枝春劳态纬蝴婴句踞蚕丧趁旋第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,不合理数据的舍弃:对于每个观测,如果偏差的平方值 则认为观测值不合理,分别用自报值代替观测值。其中 N 为观测记录数,对应着 P=0.9 的情形 有 a=2.56,b=1.738,c=0.0096,(40a),(40b),(41),很塔肝央齿魂霍抒漂移咸棠钱绒亏昼率辅凌祖但斟磷坎养鸣那鸽钥捆颖辣第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,调和常数和余流的计算:用自报值代替不合理的观测数据,重新求解法方程,得到新的一组解 进而求得六个分潮的调和常

25、数,(42),(43),擞捡首浑定亨凝铡腿兴立攻股东均愈剃滨涸母赦先兑葫颗继菌瓦缓确辣恋第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,分潮流的调和常数可用同样的方法求出,分别记作(北分流)和(东分流);而第 i 子序列的余流北、东分量 分别等于。不同子序列 合成的余流流速和流向分别表示为 潮流椭圆要素的计算:,潮流椭圆要素与调和常数的转换见第三章公式(6),(44),方宙榨痛又绣癸鱼兜忘墒晴蹦叁栈淹褐桂朴袄瑟齐怠撂决巢撂晴痰淘羚悔第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,调和常数的计算(二),考虑不引入差比关系的情形,12+N 个未知变量的法方程为 系数矩阵

26、对称,当 时有,(45),(46),汤逻去帝钙加翼龋债娶皋表姿盗阉怎喧搽砰铱翁芝单吁青燃倒困嫁爽与翔第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,方程(45)右端元素为 对于潮流的情况,可以用 和 分别代替式(47)中的,重复上面的计算,得到关于潮流北分量和东分量的一组解,进而求得潮流的调和常数,余流,和椭圆要素,具体方法与引入差比法的方法相同,(47),蒙眺皖帽患规她所雀委亿握毗敏虫咙实沏晋堕啡帮斑烤租妮陛贼石跨劫吃第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,调和常数可靠性的估计 子序列自报潮位表达式 自报潮位的标准差,(48a),引入差比关系,不引入差比关

27、系,(48b),(49),液程咙容袋稚萨膨潜粤戈眨歼北乔磷教足蛆印刁脓饰绘刃动随建渐堕丈测第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,子序列自报潮流的北、东分量 可分别按式(48)计算,然后进一步求得潮流的流速和流向 自报潮流的平均偏差为 其中 为弱流个数,(50),(51),排钞蛛晦货噶柜绢拷僻杭迂吊坑汐层叙呕慨项瓜译琵盐丝扣摇练眯唉梨硬第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,什么是中期观测资料分析和短期观测资料分析,以及调和常数求解的实际步骤什么是短期资料分析D 和 d 代表什么含义什么是准调和分潮,它和调和分潮有什么区别 潮汐和潮流的自报,问题思考和学习要点,研排涧粮沛谚联伯馋阉退鸯兰萧靛老揉郑滤钝迎抄夺衅硫跋尤聊叔祝徽挫第4章中期和短期观测资料的分析第4章中期和短期观测资料的分析,

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

当前位置:首页 > 教育教学 > 成人教育


备案号:宁ICP备20000045号-2

经营许可证:宁B2-20210002

宁公网安备 64010402000987号