DNV认证严格吗?

DNV认证严格吗?,第1张

不知道你问的是工厂认证?还是产品认证?也或者是质量体系的认证?

DNV是运国际著名船级社,运作相当规范,如果只是从职业化的角度说,DNV的认证是很严格的,因为其专业,敬业,职业,所以相对严格。

但是如果你说的严格是指容不容易通过认证,那么我觉得DNV不严格,因为只要你满足认证的要求,他们不会难为你,绝不会和政府机关一样刁难人!比如CCS那种,即使你满足条件,也会找麻烦,就是为了捞油水!

DET NORSKE VERITAS 挪威船级社

―――成立于1864年的一个独立基金会,提供风险管理服务的国际性机构

DNV致力于保护生命,财产和环境!

DET NORSKE VERITAS缩写DNV,成立于1864年,总部位于挪威的奥斯陆,是一个独立的基金会组织,自其成立之日起,便以保护生命、财产和安全为最高宗旨。DNV的铁锚与天秤标志,象征着安全、可靠与公正,蓝绿双色象征着海洋、绿地与环保,在世界各地享有极高的信誉。DNV现有员工6200多人,在世界130个国家设有330多个分支机构或办事处。在其发展初期,从事于船舶的检验,迄今为止DNV已成为全球最大的船舶入级检验机构,DNV的服务范围也已经由传统的业务扩展到诸如海上石油平台认证,风险及安全的分析控制管理,制造及服务行业的管理认证,国际质量分级系统,及相关的软件开发等领域。虽然DNV的开展的业务范围已经今非昔比,但他还是被约定俗成地被称为挪威船级社。

中国国家***近年出访挪威,都分别参观了DNV总部。目前有二十多个工业国家的政府授权DNV从事管理体系认证工作。DNV颁发的证书在全球均获得广泛信任。DNV也是欧盟成员国认可的“权威机构”和指定提供CE产品标志认证服务的机构。DNV的服务网络遍及全球,在工业及服务领域,我们提供最广泛的管理认证服务。我们开展的服务从传统的ISO9000到QS9000, VDA61, OHSAS18001,SA8000等。在全球近上千家认证机构中,DNV被认为是最具公信力的认证机构之一。在全球的认证市场中,DNV占近9%的市场份额,共发出近5万张认证证书。在船舶方面,全世界25%的新造船入DNV船级;130个国家授权DNV作船舶的调查和认证; DNV同时制定许多船舶的建造规范。在近海工业方面,超过200个近海安装设备及100个移动平台入DNV级;DNV是中国政府授权颁发海上油气田安全生产证的机构。中国海域中有75%的油气生产设施是由DNV检验发证。

DNV的认证业务始终秉持着技术领先的优势,在世界上获得了超过80多个国家的授信认可,在全球已发出各类管理体系证书65000张。作为最早进入大中国地区的国际认证机构,DNV挪威船级社迄今为止已在大中国地区成立13个办事处,颁发了近5000多张证书。

uid就是用户编码是你注册的时候在数据库里面自动生成的。

每个uid对应一个用户名就这样一般用用户名登陆,将手机摄像头对着条码扫描,成功后摄像机的uid就会显示在数据栏里。可以点击搜索,手机将会自动显示局域网内的摄像机uid。

一、海上油田风、浪、流、潮年极值要素计算

以秦皇岛32-6油田为例,海域风、浪、流、潮等年极值要素的分析计算,采用数值计算后报方式,以197\01993年间影响该海域的台风及强寒潮天气过程下的历史海面气压场及天气图作为基础资料,采用中国海油与青岛海洋大学合作开发的数值计算模式,即海面风场计算模式(包括二层诊断模式、行星边界层的动力学模式、台风模式)、浅水浪流耦合海浪数值计算模式、水位和流场的三维模式,以1/8°网格距划分计算区域,由海面风场数值计算模式计算出每年影响该海域的强寒潮及台风天气过程下各格点的风速。由浪流耦合海浪数值模式、水位场和流场的三维模式计算出各格点的波高、周期、水位、风增、减水等要素,作为推算设计环境要素的样本序列。然后依据样本序列采用WEIBULL分布推算工程点各重现期设计波要素值。

(一)强天气过程个例选取

在1970-1993年间的24年里,有120多个强天气过程掠袭渤海,但能够影响秦皇岛(QHD)32-6油田海区风、浪、流及水位产生较大变化的,只有82个天气过程,平均每年约有3或4个天气过程。它们不可能都造成年极值,只有逐个对天气过程进行模拟计算后,才能确定哪些是可取的天气过程。

(二)模拟结果衍生的其他结果

110m高度处的各种平均风速

根据海洋工程要求,依据API标准,将10m高度处平均风速的年极值按下列诸式换成同一高度的各种平均风速年极值:

1h平均风速=30min平均风速/103

3h平均风速=1h平均风速/105

10min平均风速=1h平均风速×107

1min平均风速=1h平均风速×120

3sec平均风速=1h平均风速×160

2.波浪各要素间的关系

波浪模拟给出的是有效波高Hs及对应的周期Ts。

a.假定单个波高H与有效波高Hs满足瑞利分布:

中国海洋石油高新技术与实践

F(H)是波高≥H的累积概率,则最大单个波波高Hmax与Hs有如下关系:

中国海洋石油高新技术与实践

N构成Hs的波序列波个数。N是一定时间内大浪处于平稳状态的个数。我们都知道,大浪的平稳过程持续时间不是很长,在此时间内若波的平均周期 ,取N=800个波,则持续时间τ=10s×800=8000s=222h,如果这个时间内的海浪满足平稳过程的要求,则下式为:

中国海洋石油高新技术与实践

最大单个波高Hmax按式(17-13)计算。

b.特征周期:

平均周期 是谱峰周期S(ωp)对应的圆频率); 谱峰周期

最大波的周期Tmax=106Ts(据Goda提供的关系式)。

(三)多年一遇高、低水位

多年一遇极值水位的推算,对于具有长期水位观测的地点来说,可利用实测资料的年极值高、低水位和年极值风暴增、减水,依长期极值分布求得。对于不具备长期水位观测的地点,多年一遇风增水、风减水可利用对天气过程模拟得到的年极值风增水和风减水推算出来,而多年一遇的高水位和低水位只能通过多年一遇的风增水和风减水与特征水位(例如最高天文潮位和最低天文潮位)的线形组合求得。

在渤海海区,主要水位站由实测潮位资料获得的多年一遇高、低水位和多年一遇风增(减)水,以及最高天文潮位、最低天文潮位如表17-5所示。

选用以下形式的线形组合,利用多年一遇风增、减水和最高天文、最低天文潮位推算出多年一遇高、低水位:

多年一遇高水位=072m(多年一遇风增水+最高天文潮位)

多年一遇低水位=091m(多年一遇风减水+最低天文潮位)

表17-5 渤海潮位观测站多年一遇高、低水位及天文最高、最低潮位

其结果与实测水位推算的多年一遇高、低水位较接近,表17-6列出了100年一遇的结果作为比较。因此,对于不具备长期实测水位的地点,可利用这种形式的线性组合求得接近实际的多年一遇高、低水位。

为了能利用上述方法求得多年一遇高、低水位,首先需要获得可靠的多年一遇风增、减水,表17-7给出了据模拟的强天气过程所得的风暴增水、风暴减水推算的秦皇岛和塘沽多年一遇风暴增水、风暴减水与用实测资料推算的结果,发现两种结果相当接近,从而可以利用上述线性组合获得QHD32-6油田可信的多年一遇高、低水位。

表17-6 渤海主要水位观测站高、低水位及风减水推算值

表17-7 秦皇岛、塌沽实测与计算出的风暴增减水的对比结果

最高(低)天文潮位,利用在渤海潮波数值模拟得到的8个主要分潮调和常数(M2,S2, N2,K2,K1,D1,P1,Q1)以及邻近地点的长期分潮Sa的调和常数,推算19年潮位获得。潮位表达式为:

中国海洋石油高新技术与实践

根据主要分潮调和常数,利用费拉基尔斯基方法推算理论深度基准面和BPF方法推算近最低潮面。为了保持统一取理论深度基准面为海图深度基准面。

二、不同重现期环境要素设计参数推算

依据QHD32-6油田工程设计需要,运用上述1970-1993年的强天气过程数值计算出风、浪、水位、海流极值的样本序列,采用Weibull极值概率分布推算多年一遇的设计参数。

Weibul1概率分布:

中国海洋石油高新技术与实践

式中:a,b,c为待定参数,a>0,为位置参数;b>0为尺度参数;c>0为形状参数。a=0时,上式为二参数Weibul1分布。对上式取对数可得:

中国海洋石油高新技术与实践

将上式移项得:

中国海洋石油高新技术与实践

在Weibull机率格纸上,ln(x-a)为横坐标,ln[-1n(1-F(x))]为纵坐标。

改写上式得:

中国海洋石油高新技术与实践

由于1n(x-a)中必须x>a,取a=Xn,即以序列中最小要素为位置参数,然后用最小二乘法对E、D两参数进行参数估计。这样多年一遇环境要素设计参数由下式求得:

中国海洋石油高新技术与实践

从而就可以计算出QHD32-6油田海域风速、波浪、海流及潮位多年一遇工程设计参数。

三、风、浪、潮、流联合概率计算技术

风、浪、潮、流等近岸灾害性海洋动力环境直接引起近海、海岸工程、河口海岸城市的灾害性破坏,以及人员伤亡,从而造成巨大的经济损失,给近岸带灾害动力环境的优化评估方法,首先必须考虑它的“同时出现”,“同时作用”于近海和海岸工程。因为一场台风(或风暴)过程,其最大的危害就是各种灾害性动力环境同时出现所引起的后果。但是如何考虑“同时出现”特点,国内外传统的评估方法都是分别对各种灾害动力环境进行概率分析,选用不同重现期的动力因素,再把它们叠加起来,作为近海和海岸工程极端海况的评估方法,并用作“设计标准”。例如在海岸工程设计中分别采用五十年一遇的波高和五十年一遇的潮位作为设计标准,而在近海工程设计中分别采用百年一遇波高,百年一遇风速和百年一遇海流作为设计标准。很明显,近岸带灾害动力环境各因素的某种概率条件下的极值组合在一起同时出现的事件,是一个小概率事件,用以作为设计标准,是非常保守的,特别是我国一些油田属于“边际油田”,采用过高的灾害性动力环境作为设计标准,导致过高的经济投入,将会使很多油田失去开发价值,因此传统的评估方法是不可取的。

关于风、浪、潮、流联合出现的极端海况国际上广泛使用的规范、DNV规范等都作了相应规定。我国海洋石油开发工程将规范作为行业标准,此规范针对极端海况提出了三种不同的解决途径:

第一,分别采用百年一遇的风速、百年一遇的波高和百年一遇的海流速度作为设计标准;

第二,采用百年一遇的波高以及与此伴生的风速和流速;

第三,采用联合概率为百年一遇的同时再现的风、浪、流值。

API规范指出,第一种方法过于保守,建议采用第二种方法。但又在规范中指出,“与百年一遇波高伴生的风、流值”的“伴生”这一词是含混不清的。第三种方法,由于联合概率是非单一解,使用起来有不少困难。因此在海洋工程界至今仍大多采用第二种方法。

本项研究的出发点认为,研究同时出现于一次大风或台风过程中的风、浪、流联合概率是正确而实用的途径,因为第一种方法实质上是将风、浪、流视为互相独立的随机变量,各自采用百年一遇值的概率作为三者概率的乘积,这是不合理的。第二种方法的缺点,一为API规范指出的“associated”一词是“含混不清”;二为百年一遇的波高及伴生的风速和流速(二者皆有相应的概率水平)三者同时出现的联合概率必然为百年一遇以上的重现期。本研究的若干结论也证明了这一点。

(一)联合概率随机模拟方法

多维联合概率的推求,实际上就是求解公式(17-20):

中国海洋石油高新技术与实践

式(17-20)只有当随机变量均为正态分布时才可解析求解,对非高斯函数,具有不同相关性的多维随机变量的联合概率密度函数,用随机模拟法是行之有效的。

随机模拟法是基于现实资料和一些假定,通过计算机来重复某些过程的方法。常用的模拟求解法是蒙特卡洛法,但M-C方法求解较小失效概率时耗时巨大,且结果精度不高,必须寻找新的方法。重点抽样法是减少机时、降低方差的一种行之有效的方法,其基本原理是集中对分布的最重要的区域抽样,即对失效概率做出主要贡献的那部分抽样,而不是扩展到整个定义域均匀抽样。

图17-1 二维联合概率示意图

为设计点坐标;fx(x)为联合概率密度; 为权函数密度

为了说明重点抽样法的特点,图17-1以二维情况示意,图中贯穿x1、x2平面的曲线,代表模拟联合概率状态的联合概率曲线,曲线右边是超过联合概率的分布;左边是低于联合概率的部分,曲面上可找到一个距离原点最近的点,称为“设计点”, 为其坐标。重点抽样法的特点,在于围绕设计点附近抽样;其另一个特点是引入了权密度函数hy(x),通过它将模拟引向以设计点为中心的区域,这样才能达到降低模拟结果方差的目的。

公式(17—20)的计算公式为:

中国海洋石油高新技术与实践

式中:

概率密度函数fx(x)代表 权密度函数特殊情况下真正的联合概率密度。

由上式可见,权密度并不重要,因此,联合概率的期望值可写为:

中国海洋石油高新技术与实践

式中:N为模拟次数; 为第I次模拟向量。

优点在于无须考虑基本随机变量的分布类型而适用于原始空间。对具有相关关系的随机变量,基本随机变量转换为独立的标准正态变量向量是比较困难的。但是,实际上联合概率的计算使用原始分布,任何非高斯分布是不会影响到加权样本的。在使用任何非高斯分布和相关的随机变量时,需进行下列变换:

中国海洋石油高新技术与实践

式中:Fxi(ixi)是基本随机变量x的原始累积分布函数;Φ-1()是标准高斯累积分布函数的反函数。假定Z为标准正态,则联合概率密度函数可写为:

中国海洋石油高新技术与实践

式中:Zi为式(17-23)的计算结果;φ()为标准正态密度函数;φn(Z,R´)为均值0、标准差为1的多维标准高斯密度函数;R´为 构成的修正相关矩阵; 是由相关系数 系列定义的值。

中国海洋石油高新技术与实践

对每一对边缘分布,式(17-25)可迭代求解。

应该指出,边缘分布密度和协方差对非高斯随机变量的联合密度并非为一对一定义,然而对大多数情况可利用的信息仅仅局限于边缘分布和协方差。因此,任何适合的模式,只要不与资料矛盾,都是可采用的,其使用范围依据变量间的相关系数,而无须考虑对任何模式的数学依据。

(二)联合概率随机模拟软件开发

基于重点抽样法的联合概率随机模拟JOPAP软件包括两部分。

1.设计点的计算

这样随机点的抽样能在以设计点为中心的有效区域内进行,可以减少计算时间,提高模拟效率。

2.重点抽样并计算联合概率

在执行主程序前,需针对所要计算联合概率的实际问题的不同,设计出不同的联合概率模型的基本方程,称之为极限状态方程。模拟的极限状态方程为:

中国海洋石油高新技术与实践

调整不同的设置值(A1,A2,A3)即可得出不同的联合概率,如将Ai分别赋予各自Xi变量独立统计时所得的特征值(如百年一遇值、五十年一遇值),则可得各变量Xi超过各自百年一遇(五十年一遇)值联合出现的概率。本研究已经完成了上述软件的研制工作。

(三)联合概率随机模拟非单一解的求解方法及其工程应用

最大波高及相应同时出现的风速、流速联合概率计算,对于固定式海上结构由荷载分析可知,当采用最大波高与相应同时出现的风速、流速作为随机模拟样本时,其联合作用下产生的荷载效应对平台结构最为不利,如以下计算结果所示。

1.求各要素单独的百年一遇极值

波高、风速、流速的分布通常符合某种概率分布,如冈贝尔分布(Gumbell distribution)、威布尔分布(Weibul1 distribution)、对数正态分布(Log-Normal distribution)等。因此可根据已有多次大风过程中波高、风速、流速的数据资料,进行上述分布的拟合。可根据最大波高或选取相应的一组风浪流组合,从而形成以波为主及以流为主和以风为主的风、浪、流序列。每个序列都进行冈贝尔分布、威布尔分布及对数正态分布的适线计算。通过分布曲线拟合优度检验后,采取一定判别标准来进行判别,最终可以确定每一种组合中相应的风速、波高及流速的分布形式。这里使用柯尔莫哥洛夫检验法,即常用的K-S检验法进行检验。并以频率绝对离差和离差均方和最小为判别标准。所得各要素分布形式的计算结果见表17-8。

由表17-8中还可看出一种似乎“矛盾”的现象,即以流最大序列统计所得的百年一遇海流值,反而小于以风最大序列统计所得的百年一遇海流值,而且以浪最大序列统计所得的百年一遇波高值,小于以风最大序列统计所得的百年一遇波高值。这一现象可以解释为海流速度序列(波高序列)的方差较小,因此尽管以流(浪)为主海流(波高)原始序列值都大于以风为主的海流(波高)原始序列值,但其概率分布曲线较平缓,所以百年一遇值反而低于以风为主的海流(波高)百年一遇值。

表17-8 风、浪、流各自的分布形式及百年一遇极值

2.联合概率随机模拟分析方法

风、浪、流三种环境要素之间是相关而不是独立的,联合概率法以“风暴”或“台风”过程中同时出现的风、浪、流作为随机分析的基本系列,从而得到某组风、浪、流组合及相应的概率水平。即求解公式:

中国海洋石油高新技术与实践

式中:g(x)<0表示失效域。

对于非高斯且相关的多维随机变量,用随机模拟法求解该式可以收到很好的效果。随机模拟法是基于现实资料的某些特征和假定,通过计算机来实现某些过程的方法,蒙特卡洛法是其中的一种。它的基本思想是建立一个概率模型,使其数学期望等于结构的失效概率,然后对这个模型进行随机抽样,最后用子样均值估计失效概率。但求解较小联合概率时,直接用蒙特卡洛法需耗费巨大机时,且会产生较大误差。因此蒙特卡洛法只适用于失效概率不太小的情况。而“重点抽样法”可提高蒙特卡洛法抽样的效率,它的基本原理是集中对分布中的最重要区域进行抽样,即对联合概率做出主要贡献的部分抽样,而不是扩展到整个定义域均等的抽样。在用随机模拟法求解公式(17-26)的过程中,存在非单一解问题,即应用“重点抽样法”的随机模拟程序,以海洋环境要素作用下风浪流的分布形式、均值、方差及相关系数矩阵为输入条件,经过计算,可得到很多组相当于联合概率为百年一遇的海洋环境要素组合。解决非单一解问题,需有一定的控制条件,一般以响应最大为判断标准。若响应以 Y表示,则Y=f(H,Uc,Uw),其中H为波高,Uc为流速,Uw为风速。H、Uc、Uw的取样将影响 Y值大小,因此通常根据资料取以波高为主、以流为主、以风为主的相应的风浪序列进行模拟。应用“重点抽样法”的随机模拟程序,可得到的是分别以流为主、以浪为主、以风为主的风浪流组合。用平台响应判断标准,取最大响应对应的组合为最终解,由此可得到相当于百年一遇的海洋环境联合设计标准。计算结果见表17-9。

表17-9 联合概率方法下的风、浪、流值

导管架平台风、浪流单因素分布形式、百年一遇值及响应值见表17-10。

表17-10 导管架平台风、浪流单因素分布形式、百年一遇值及响应值

(四)最大流速及相应同时出现的波高、风速联合概率计算

对于自由浮动的半潜式平台,平台因波浪引起的稳定运动而承受的荷载最大,因此采用波高最大时及相应风速、流速的组合作为联合概率选择准则是合适的。

对于约束较强的半潜式平台,如计算结果所示,平台受到的环境荷载在这种情况下最大,见表17-11。

表17-11 浮式平台海洋环境联合设计标准

(五)最大波高和最大流速相应同时出现的联合概率

海洋立管(Riser)是海洋工程结构系统中的重要部件,也是薄弱易损的构件。计算结果表明必须采用以波浪或海流为控制条件的风、浪、流联合概率作为设计标准。见图17-2和图17-3以及表17-12各种方法的比较。

图17-2 以流速为主百年一遇曲线

图17-3 以波高为主百年一遇曲线

表17-12 各种方法的比较

(六)最大风速及相应的波高、流速联合概率

由于自升式平台重要组成部分迎风面积较大,因此,最大风速及相应同时出现的波高、流速可以导致最大的响应,现以实例计算见表17-13。

表17-13 联合概率百年一遇风浪流组合、相应倾覆力矩与极值响应法的比较

(七)结论

a.按照预计计划完成了“非高斯过程,具有不同相关性的多维随机变量联合概率”随机模拟的软件开发(JOPAP软件开发)。随机模拟技术具有收敛快(比M-C法快5~10倍)、误差小(最大相对误差不超过10%)的优点。

b.随机模拟软件可使用于各类灾害性动力环境经常采用的概率分布模式(如极值I、Ⅱ、Ⅲ型分布,WEIBULL分布,对数正态分布,复合极值分布等),可满足各类工程防灾的概率分析。

c.本研究针对不同工程结构,从最不利的灾害环境组合产生的响应出发,进行不同类别同时出现的灾害环境资料抽样,从而达到了将联合概率随机模拟的非单一解转化为固定解的目的。

d.针对美国API规范以及中国和其他规范,对同时出现灾害动力环境评估方法中存在的问题,本项目具体解决了:①对固定式平台,建议采用以波高为控制条件的风、浪、流联合概率作为设计标准;②对自升平台,建议采用以风速为控制条件的风、浪、流联合概率作为设计标准;③对平台立管在海流较大的海区,必须采用以波浪或海流为控制条件的风、浪、流联合概率作为设计标准。

e.本项目成果解决了美国API最新版本(中国海洋石油行业规范)中“含糊不清”的一些提法,在国际上首先以联合概率的方法,解决了不同工程结构设计标准的合理确定方法。

DNV GL的历史可上溯至1864年,当年DNV作为一家会员制组织在奥斯陆成立。挪威的海事保险俱乐部联合起来建立一套统一的规范和程序用于评估船舶保险的风险,目的是为挪威籍商船提供可靠并统一的评级和课税。

当时挪威航运业正经历着快速发展并摆脱传统束缚的羁绊,因此正需要建立一个新兴的全国统一的海事保险市场。三年以后在德国,600多个船东、船厂和保险商聚集在汉堡证券交易所大厅,召开GL作为一家新的设在汉堡的非盈利性组织的成立大会。 0571-85860756

GL的成立是基于对于透明的愿望。商人、船东和保险商通常对于船舶状况的信息了解甚少。作为一家独立的船级社,GL的成立是为了评估船舶的质量并告知利益相关者评估结论。GL的第一本国际船舶登记册从1868年开始记录了273艘船的报告。到1877年,登记的船舶数量增长了数十倍,验船师网络也因此扩展迅速。

DNV的船队也得以迅速增加。第一批代理即当时称作永久性验船师,在很多国家被任命,以服务于在国外的挪威籍船舶。1870年代出现了蒸汽船舶,极大地改变了船级业务和对验船师的工作和能力要求。

DNV和GL在一开始就开始了合作。DNV委员会的记录表明从1868年9月开始就计划建立一个两家船级社的共同入级登记。这一计划最终没能实现,同样在1891年双方有类似的讨论希望能互相认可证书并建立共同的船舶登记册。

社会驱动力

社会日益成为私营占主导地位行业的重要利益相关方。1891年起Samuel Plimsoll提出的装载线在英籍船舶上被强制执行,以保护英国沿海的海员生命。装载线于1907年在挪威被强制实施。

1912年的泰坦尼克号惨案将海上安全问题呈现在公众面前,成为公众关切。国际船级社在讨论船舶安全的议题中发挥了重要作用。然而,在海上安全和人命的第一次国际大会上只有两名船级社的代表,他们分别是GL的董事总经理Carl Pagel 和来自DNV的Johannes Bruun。

战争的代价                  杭州中国产品CCC认证网

对GL来说,第一次世界大战是一次严厉的打击。国际关系被斩断,外国旗的船舶纷纷更换船级社。在二次大战之前的战争间隙GL又得以一定程度的发展。二次大战以后德国经济的恢复促进了GL的发展。

第一次世界大战以后,由帆船向蒸汽船舶的过渡为船级社行业的技术和技能的发展带来了深刻的变化。过时的规范已经不能适应当时的造船方法。在1920年至1940年期间,DNV在技术上保持独立,并且建立了优先发展工程、建造和设计的企业文化。随之而来的第二次世界大战,DNV几乎面临分裂的局面。

欢迎分享,转载请注明来源:浪漫分享网

原文地址:https://hunlipic.com/meirong/10455327.html

(0)
打赏 微信扫一扫微信扫一扫 支付宝扫一扫支付宝扫一扫
上一篇 2023-11-06
下一篇2023-11-06

发表评论

登录后才能评论

评论列表(0条)

    保存