模型应用三

来自:    更新日期:早些时候
请教各位:如何运用三因素模型的组合收益率,来计算~

您好
运用三因素模型的组合收益率根本没有办法计算的,建议您现在做股票要快进快出,不要做中长线为好,因为现在的股票市场很不稳定,今天是超跌反弹,明天就是一次大跌,我国的股市已经进入下跌行情,加上诸多经济问题,股票市场今年几乎不会有很好的机会了,作为理财师我的经验告诉我,估计下星期还会下跌,所以您要注意资金安全,有什么问题可以继续问我,真诚回答,恳请您采纳!


以元宝山露天矿区地下水疏干井群优化设计为例。

元宝山露天煤矿于1954年发现,经过勘探,1987年由沈阳煤矿设计院完成露天矿初步设计,并于1990年10月15日正式开工建设。元宝山露天矿区,在历史上没有开采记载,其西南部有元宝山一、二、三、四井,东北部为老公营子和小风水沟井田,露天矿西南部现有两座地方小煤窑沿七煤露头向井田内开采。露天矿现已建成西、南两个排土场,一、二采区也正在建设中。采掘场南部边界形成于元宝山脚下,东部边界以F1断层为界;而西部边界则是沿着六煤层底板形成的。开采用一、二采区同时拉沟的方法,由南向北推进,北帮为工作帮。露天矿设计最终开采面积12.32 km2。露天开采储量54289万t,A+B级储量53265万t,设计规模为年产原煤500万t。

从1954年以来,地质、煤炭、水电等部门先后在该地区进行了大量的卓有成效的工作,取得了丰富的资料。1954~1955年,平庄矿务局地质队及东煤地质局107队在该地区进行了煤田地质普查工作;东煤地质局104队于1973年提交了“元宝山露天精查地质报告”,于1982年提交了“元宝山露天水文地质、工程地质(剥离物强度)勘察报告”;水电部东北电力设计院于1975年提交了“元宝山电厂新建工程供水水文地质勘察报告”,水电部南京水利科学研究院于1987年提交了“元宝山露天煤矿受英金河渗漏影响计算”;煤炭科学研究总院西安分院于1993年提交了“元宝山露天煤矿帷幕工程初步设计说明书”,“元宝山露天煤矿帷幕截流工程地质勘察报告”等成果。前人的工作积累了本区丰富的地质及水文地质资料,为今后的工作奠定了良好的基础。

但就在露天坑剥离建设过程中,丰富的第四系松散沉积含水层中的地下水给露天剥离带来了巨大的困难。目前正在剥离区外围大面积区高强度疏干地下水。自从1990年以来,先后施工疏干孔近120个,每天总排水量达40~50万m3左右。露天坑及其附近第四系含水层地下水位下降约8~20 m,采区内最大水位降深为27 m左右。区域内第四系含水层地下水流场已基本趋于稳定。但目前的地下水流场还远不能满足已经开挖的剥离区(一采区和二采区)的剥离要求,特别是二采区安全剥离的水位降深值应为20~50 m左右。因此,如要二采区安全出煤,还需对本区第四系地下水位疏降10~30 m左右。因此,目前的疏干方案、疏干工程和疏干水量都远不能满足矿山建设和生产的需要。选择最优秀的疏干工程以最小的总疏干水量达到疏干要求成为目前矿山建设急需解决的问题。

7.3.1 地下水系统及水文地质模型

7.3.1.1 研究区概况

(1)地理与交通

元宝山露天煤矿位于内蒙古自治区赤峰市东35 km,属赤峰市元宝山区建昌营子乡所辖。其地理坐标为东经119°17′55″~119°19′55″;北纬42°19′13″~42°22′21″。

矿区南部有专用铁路,在元宝山车站与叶赤线(叶柏寿—赤峰)接轨,矿区与赤峰及邻近旗、县均有三级公路联结,交通十分便利。

(2)地形与地貌

元宝山露天煤矿位于英金河河谷平原,英金河从矿区中部穿过,把矿区分为南、北两大部分。露天矿南部地处英金河右岸一级阶地,阶地宽500~4000 m,地面坡度为1‰~1.5‰,地面高程为472~482 m。露天矿北部位于英金河左岸一、二级阶地,地面坡度1‰~2‰,阶地宽500~5000 m,地面高程482~490 m,阶地表面有现代风积沙,沙丘呈波状分布。

总体来看,矿区为四周被低山丘陵所环绕的河流冲、洪积平原。周围丘陵高程一般在500~600 m。冲、洪积层形成本区极其富水的含水层,而丰富的巨厚层煤炭资源就隐伏在第四系含水层之下。

(3)气象与水文

该区为半干旱大陆性气候。其特征是:冬季长而干冷,夏季雨量集中,春秋季少雨雪多大风。

据赤峰市气象台资料,元宝山地区多年最高气温为42.5℃,最低气温为-31.4℃。冻结期一般为11月中旬至翌年3月末(平均气温-10.8℃,最低气温-27℃,最高气温7.6℃),最大冻结深度2.01 m。

本区多年平均相对湿度为49%,平均蒸发量为1867.1 mm,年平均降水量为372.34 mm(据1950~1994年资料),全年降水多集中在夏季6~8月,约占全年的68.55%。

英金河自西北向东南流经矿区,于东八家汇入老哈河,为老哈河左岸最大的支流。英金河发源于河北围场北部山区(七老头山)。流长194.6 km,流域面积10598 km2。历年最大洪峰流量2650 m3/s,最小流量为0.5 m3/s,多年平均流量为12.8 m3/s。河床宽度变化较大,在200~900 m之间。主流摆动对两岸侧蚀较强,洪水期常造成河岸坍塌。近年来由于上游水库蓄水和农业灌溉的发展,位于下游的露天矿区一带春秋冬季常常断流。该河枯水期、平水期、丰水期的流量分布与降水量的分布规律相同。在元宝山露天矿区以垂直渗透形式补给地下水。

老哈河由西南向东北流经河谷平原南部,距露天矿区3 km。该河发源于河北省平泉县七老头山脉的光头山,于昭乌达盟大兴乡海里吐附近与西拉木伦河汇合成西辽河。全长421.8 km。流域面积33076 km2。历年最大洪峰流量为9840 m3/s,最小流量为0 m3/s,多年平均流量为13.6 m3/s。研究区内河床宽度在500~1000 m之间。河床及河漫滩主要由沙、沙土及砾石组成,该河枯水期、平水期、丰水期的流量分布亦与降水量的分布规律相同。在元宝山露天矿区以垂直渗透形式补给地下水。

7.3.1.2 地质与水文地质条件

(1)矿区主要地层

A.上侏罗统杏园组:以灰白色中细砂岩为主,夹紫红色砂岩和泥岩,厚度大于100 m。中部为灰、灰绿色砂岩和砂砾岩夹黑色泥岩,厚度60~230 m。上部为灰绿、灰褐色厚层泥岩夹灰白色砂岩,厚度在200 m左右。分布于F1断层以东的断块中。

B.上侏罗统元宝山组:以灰白色中细砂岩为主,夹粗砂岩、泥岩和煤层,一般厚约340 m,含12个可采煤层,累计可采煤层平均厚度84.29 m,以五、六煤组为主要可采煤层。

C.第三系上新统(N2):底部为紫红色砂岩、泥砾岩和泥岩,不整合于元宝山含煤组地层之上,厚0~115 m,仅在露天矿南部的穹窿背斜处分布。上部为玄武岩、红土和砂砾石层,仅出现在矿区北部或覆盖于露天矿的南部紫红色砂砾岩之上。

D.第四系(Q):以现代冲积、洪积和冰水堆积物为主,由安山岩、玄武岩等砾石成分构成的圆砾、泥砾、卵石和砂组成。厚度14~85 m,一般为55 m,在整个元宝山盆地内均有分布。在矿区附近,厚度一般为14~60 m,自西向东由薄变厚。

(2)矿区主要地质构造特征

元宝山煤田为一断陷含煤盆地,受燕山构造变动的控制呈NNE—NE向狭长带展布,含煤盆地为一宽缓的复式向斜构造,由三个向斜和两个背斜组成,自东南向西北依次是:风水沟短轴向斜、五家背斜、南荒向斜、龙头山背斜和老窑短轴向斜。在向斜构造内赋存有可采煤层,向斜轴向NNE,地层倾角3°~5°。

总体而言,本区第四系下伏煤系地层比较平缓,主采煤层厚60 m左右,煤层产状及其赋存条件非常有利于煤炭资源的露天开采。

(3)矿区水文地质条件

第四系孔隙潜水含水层由冲积、洪积和冰水堆积的圆砾、砂砾、卵砾石、泥砾等组成,粒径为5~60 mm的占50%以上,大于60 mm的占20%,个别地段见漂石。砾石成分以安山岩、花岗岩为主,磨圆度较好,球度差。从西南往东北厚度增大,但变化趋势比较平缓,仅在基岩面的两阶地之间的阶坎处厚度变化比较大。

第四系地层由于成因不同,上部和下部地层渗透性有一定的差异,据以往勘探试验资料,上部渗透系数较大,约为256~710 m/d,下部较小约为16~146 m/d,但二者有密切的水力联系。水质属重碳酸钙镁型水,pH值大部分在7~8之间,固形物含量240~400 mg/L,钙镁总含量为260~365 mg/L,水温为8~11℃。

侏罗系孔隙裂隙弱含水层,由砂岩、砂砾岩、粉砂岩及煤层组成,在煤层中有少量的裂隙。一般厚度为50~150 m,平均厚度为113.9 m,根据抽水试验资料渗透系数为0.001~0.38 m/d。

由此可见,第四纪松散沉积潜水含水层是本区惟一的主要含水层。而其他基岩裂隙水可以忽略不计。

在天然条件下,本区地下水的流向和地表水流向一致,即由盆地的北西,南西向东南径流。水力梯度平缓(见图7.23)。近年来,随建昌营电厂水源地抽水及露天矿剥离疏干排水,使得第四系地下水位形成了以露天矿剥离区为中心的降落漏斗。地下水形成了从四周向漏斗中心汇流的新的径流条件。

第四系地下水的补给主要来自盆地内部季节性大气降水、流经盆地内部的英金河和老哈河的渗漏和来自北西、南西的上游侧向径流。特别应该注意的是目前地下水位降落漏斗已经越过英金河和老哈河向外扩展,所以河流对地下水的补给以渗入式为主,而不是注入式补给。由上述分析可见,研究区第四系地下水目前的主要补、径、排关系如图7.24所示。

(4)第四系地下水疏干现状

从1990年8月15日疏干工程开始至今,共投入运行13排,计121个疏干钻孔,有部分孔报废,部分孔停运。平均排水量达40万~50万m3/d(见表7.23)。采区内的漏斗中心在观5孔位置附近,水位标高在(438 m±1 m)变动。主要疏干钻孔见图7.25。

图7.23 矿区第四系地下水天然流场图

图7.24 元宝山矿区第四系地下水补、排关系示意图

与露天矿相邻的建昌营电厂采用群孔集中抽取第四系地下水作为电厂供水水源,其日抽水量在10万m3左右。由于两个水源地的长期抽水,区域地下水流场近几年基本趋于稳定,疏干排水结果使矿区水文地质条件发生了较大的变化。主要表现在:①沿英金河方向地下径流水基本被疏干截夺;②地下水水位低于英金河和老哈河的河水位,使二河成为“悬河”;③地下水从英金河和老哈河获得渗透补给。

(5)元宝山露天矿区第四系地下水流数学模型

根据前述元宝山露天煤矿水文地质条件,所选择模拟计算的主要含水层为位于英金河和老哈河冲积平原范围内的第四系潜水流含水层,面积约为210 km2。西北边界、南西边界和东南边界,一般离疏干区较远,可作为定水头边界处理,其余边界作为隔水边界。

图7.25 矿区目前主要疏干井群分布图

表7.23 1990~1995年疏干排水量统计表(万m3/mon)

将第四系含水层视为一个非均质,各向异性含水层,尽管第四系底板有一定的起伏,但地下水流仍可视为潜水二维非稳定流动。英金河和老哈河对潜水含水层以渗入方式补给。

根据元宝山矿区第四系潜水地下水特征及边界条件,文中建立了元宝山露天煤矿区第四系地下水流二维非均质,非稳定各向异性地下水运动数学模型。

地下水系统随机模拟与管理

式中:h——潜水地下水水位[L];

μ——给水度(量纲一);

Kxx,Kyy——第四系含水层x,y方向主渗透系数;

t——时间;

(x,y)——笛卡儿坐标;

Ω——地下水渗流区域;

Γ1——第一类边界条件;

Γ2——第二类边界条件;

h0(x,y,t)——初始水头分布;

ε———单位面积上的入渗补给强度[L3/(T·L2)],主要包括大气降水补给、河流入渗补给等;

W——源、汇项,本模型中主要反映了露天疏干和电厂取水等抽水量;

Z——含水层底板标高(L)。

对模型(7.3)采用迦辽金有限元方法进行求解,采用三角形单元剖分。平面上共剖分为1788个三角形单元,952个计算节点。见图7.26和图7.27。由于计算区域较大,所以图7.27是对露天剥离及井群疏干区域的放大。

图7.26 矿区地下水有限元计算剖分图

7.3.2 水文地质参数随机性及参数识别

研究区的水文地质条件和其他地下水系统类似,控制地下水流的主要参数亦具有极强的随机性和不确定性。就含水层介质结构而言,由于受其成因条件的限制,不管在垂向上或平面上,其沉积物性质变化较大,往往是砂、砾石及粉砂质黄土互为透镜体状产出。这种介质性质在空间上分布的随机性决定了含水层主要水文地质参数(Kx,Ky,μ)的随机性。其次,含水层的补给条件(河流入渗、大气降水及侧向径流)强烈地受到本区降水规律的影响。由于大气降水因素的随机性也决定了本区地下水的补给条件具有随机性。最后,地下水的主要排泄条件(电厂供水、露天坑疏水、农业用水等)都受到设备、用水量等多种人为因素的干扰和影响,所以排泄条件亦可视为随机因素。由此可见,影响本区地下水补给,径流和排泄的多种因素都具有随机性和不确定性。用任何一组确定的参数去刻画本区地下水系统的行为都欠准确。所以,引入随机理论和统计概念来研究该含水层系统的性质和对其进行规划管理则更具有实际意义。

图7.27 矿区地下水有限元计算剖分图

根据上述分析,采用研究区1991~1995年的水文地质资料(抽水资料、降水资料、地下水动态资料)分别求得了5组水文地质参数。调参过程采用了试估-校正法。即首先根据已有资料给出参数初值,运算地下水数学模型,求解地下水水位值,将所求结果与实测结果不断进行比较,修正参数,直到达到要求的拟合精度。并在假设检验的基础上,分别统计计算了每个参数的均值、方差及其概率分布形式。表1.1为15个参数区水文地质参数Kx,Ky及μ的反演值、均值和方差。通过假设检验,Kx,Ky及μ均服从[a,b]均匀分布。

在参数反演过程中,对下列问题进行了专门处理:

(1)初始流场。以每个模拟年份1月10日区域观测水位为基础,采用Kriging插值方法,插出每个节点的水位值,作为该年份模拟计算的初始流场。利用每年1~6月份水位观测资料作为调参拟合水位。

(2)降水量。模拟过程中的降水量采用赤峰气象站的实际观测资料。在疏干预测过程中采用1950~1995年月平均降水量资料。降水入渗系数取0.3。

(3)河流入渗补给。根据实际观测资料,英金河和老哈河的最大渗透量分别为4.12×104 m3/s 和1.6×104 m3/s。计算时用其总渗透量除以两河在计算区的面积得出单位面积的渗透量,加入到相应的面积单元中。农业取水和灌溉回渗因缺乏资料,在本模型中未进行专门考虑。

(4)由于潜水流为非线性偏微分方程,文中在利用响应系数法建立管理模型时,利用Boussinesq方程进行了近似线性处理,将其近似为线性问题考虑。

正如前面有关章节所述,元宝山露天矿目前生产和建设所遇到的最大问题就是第四系地下水的控制与管理问题。就这一问题前人已经做过大量的研究工作,并提出了建造防渗帷幕墙、疏干与回灌相结合等多种技术方案,其核心目标就是希望在保证矿山安全生产的同时,尽可能保护地下水资源,减少地下水位的大面积下降。但因经济、社会等多方面的原因,这些方案都未能付诸实施。目前仍以大面积井群疏干为主要防治水技术措施,且现有的疏干工程和设计能力无法满足生产要求。因此,基于现有工程状况提出新的既能满足矿井建设与生产的要求,又能确保矿井总疏水量最小的矿井疏干井群设计方案及孔位布置原则显得尤为必要和迫切。正是本着这一目的,以矿井剥离区四周水位降至设计标准为约束条件,以稳定总疏水量最小为目标函数,以疏干水量为决策变量建立和求解了随机地下水控制与管理模型。提出了在不同的约束条件置信度水平下总疏干水量及其疏干孔位优化设计的原则。

7.3.3 机会约束地下水管理模型建立

根据元宝山露天矿区第四系地下水疏干管理的约束条件,目标函数及决策变量,建立矿区第四系地下水疏水量优化设计的机会约束随机管理模型为:

地下水系统随机模拟与管理

式中符号意义同前,其中:n=1,即考虑了一个疏干阶段。m=96,即选择沿剥离区外围两排节点(总计96个)为疏干井的候选位置(决策变量),见图7.28。j=35,即沿剥离区边沿一周的节点(总计35个)为水位约束控制点。

将有关参数代入模型(7.4)式,并进行适当转换后,得下列管理模型:

地下水系统随机模拟与管理

式中:S(j)=H0(j)-ZL(j)-5;

i=1,2,…,96;

j=1,2,…,35。

7.3.4 模型求解

将本区第四系含水层随机水文地质参数的反演结果、随机分布形式及其他参数代入随机有限元模型,并采用Monte-carlo随机有限元求解技术解得随机管理模型的响应系数均值E[β(i,j)]和方差 r2(i,j)。采用Taylor展开随机地下水管理模型求解技术解得在不同随机约束置信度水平下的地下水总疏干水量及其疏干位置分布。计算结果见表 7.24。总疏干水量与约束置信度水平之间的关系见图7.29。

图7.28 可供选择的疏干井位分布图

表7.24 考虑目前开挖区范围条件下计算结果表

图7.29 总疏水量与置信度水平关系图

7.3.5 计算结果讨论

(1)从计算结果可知,由于水文地质参数的随机性,要达到疏干要求,其总疏水量与对约束条件满足的置信度水平有密切关系。随着约束条件置信度水平的降低,其总疏水量明显下降。而随约束条件置信度水平的提高,其总疏水量迅速增大。说明若要保证在水文地质参数出现不利于疏干进行的小概率事件时,仍能满足疏干要求,则总疏水量必然增加,这与理论分析及实际情况相一致。

(2)目前所具备的40万~45万m3/d的疏干能力显然太小,即使80%的约束置信度也不能满足。所以,总疏水量的增加是不可避免的。

(3)从疏干井的分布来说,主要集中于当前剥离区的东北和东南(第四系厚度较大区)。这也与水文地质条件分析结果一致。因为疏干区含水层基底总趋势是西高东低。若要保证整个疏干区地下水位的疏干,只要东部区能达到疏干要求,西部区水位自然可降至疏干要求(因含水层渗透性很好,地下水位将非常平缓)。因此,若将疏干孔布置在西部区必然是浪费和不必要的。

(4)为了检验管理结果的优越性和正确性,以目前的疏干水位降为约束条件代入管理模型进行了求解。如果考虑参数的方差为零(即为确定性模型),其计算结果如表7.25。如果考虑参数的方差及约束的置信度水平,则计算结果如表7.26所示。表中结果说明,对于确定性参数,如果按优化的疏干井群进行疏干,则总疏水量可比目前实际疏水量减少10万m3/d左右。如果按随机参数模型考虑,则目前的疏水量恰相当于约束条件置信度为(95%~100%)之间的计算水量。由此可见,管理模型的计算结果是符合实际情况的。

表7.25 元宝山露天矿疏干条件确定性模型优化水量及分配计算结果

表7.26 当前疏干条件优化计算结果表

续表

(5)为了进一步检验管理结果的正确性,我们将管理结果代入地下水模拟预测模型,进行了地下水疏干流场的模拟预测,预测的流场形态较好地反映了在不同约束条件置信度水平下对疏干要求的满足性。




模型应用三视频

相关评论:
  • 18942425180时间序列模型(三):MA模型
    倪习强2. 时间序列基础知识在深入MA模型前,务必掌握时间序列的基本概念,如单变量与多变量的区别,以及时序模型与机器学习的区别。《时间序列模型(一)》和《时间序列模型(二):AR模型》将为你提供扎实的背景知识。3. MA模型的实质与应用MA模型的精髓在于,它假设数据围绕均值波动,其中白噪声的线性组合决定了...

  • 18942425180地球物理探测三维模型的建立与应用
    倪习强图5-1-5 根据物探结果形成的煤火物探测量可视化模型 基于高密度电法电阻率180 Ω·m形成的着火区的实体模型。其中,红色部分是根据高密度电法测量结果生成的实体,蓝色部分是其外推的实体 (三)三维模型在煤火勘探及灭火中的应用 对乌达煤田着火区物探测量三维数字可视化模型的研究表明:(1)三维数字...

  • 18942425180毕业季:什么是计量模型?如何设计?如何应用?计量模型成功三要素!
    倪习强如模型识别、估计方法和软件应用;以及高质量的数据,作为模型推断的支撑。没有坚实的理论,模型就失去根基;没有适当的方法,模型无法转化为实用工具;而数据的质量直接影响模型的可信度和应用效果。因此,理论、方法和数据是计量模型成功不可或缺的三要素,它们共同塑造了模型的价值和影响。

  • 189424251803D模型软件在建筑设计中的应用有哪些?
    倪习强1. 概念设计:使用3D模型软件可以更加直观地展示建筑设计的概念或想法,减少建筑师和客户之间的沟通误差,提高效率。2. 空间布局和规划:通过3D模型软件,可以模拟建筑内外部空间的布局和规划,包括墙体、门窗、楼梯、电梯等部件的位置和大小等,进而预测建筑物的功能和实际效果。3. 可视化演示: 使用3D模型...

  • 18942425180蝴蝶模型三个定理
    倪习强总之,蝴蝶模型三个定理的研究,不仅深刻影响了气象学这一学科,而且也具有更广泛的意义。它探讨了科学建模中的局限性,提供了哲学层面上的思考。而这些思考,也可以应用到更多实际应用中。比如在数据采集、模拟计算等领域,更好地处理数据,更准确地模拟,无疑会为人们提供更好地决策基础。

  • 18942425180构建物理实验模型的应用有哪些
    倪习强3、生物物理学研究:生物物理学是研究生物系统的结构和功能的科学。构建物理实验模型在这个领域中有着广泛的应用,例如,科学家可以使用计算机模拟来模拟生物分子的行为和相互作用,从而更好地理解生命的机制。物理实验的作用 1、验证和发展理论:物理实验是验证物理理论的重要手段。通过实验,我们可以证实或...

  • 18942425180mvc是什么区块链,什么是mvc和mvvm
    倪习强MVC是一种设计模式,它强制性地把应用程序的输入、处理和输出分开。MVC把应用程序分成三个核心模块:模型、视图和控制器,它们分别担负不同的任务。 MVC处理过程: 首先控制器接受用户的请求,并决定应该调用哪个模型来进行处理;然后模型根据用户请求进行相应的业务逻辑处理,并返回数据;最后控制器调用相应的视图来格式化模型...

  • 18942425180因果分析(因果关系)的原理、方法论(模型)和应用
    倪习强因果关系,这个科学洞察的基石,在数据科学的殿堂中熠熠生辉。它不仅驱动着A\/B实验和异常分析的探索,还深刻影响着我们理解世界的方式。让我们一起走进这个神秘领域,探索因果分析的奥秘、原理与实践应用。首先,我们从基础概念出发。因果分析,就像一双解构现象的慧眼,揭示出事件背后的因果关系。它包含了...

  • 18942425180预测模型的建立及应用
    倪习强(二)证据层选择 根据对可可托海地区矿床控矿因素的分析,认为研究区的主要控矿因素取决于以下几个方面:有利地层岩性组合发育,岩浆岩脉群发育,线性构造、环形构造发育,有利地球化学异常,铁染和羟基蚀变带。(三)模型应用 本次工作以地质异常理论为指导,以中国地质调查局开发的矿产资源GIS评价系统...

  • 18942425180计量经济模型的应用主要包括哪些方面
    倪习强1、结构分析是对经济现象中变量之间相互关系的研究。它研究的是当一个变量或几个变量发生变化时会对其它变量以至经济系统产生什么样的影响。2、经济预测:计量经济学模型,是从经济预测,特别是短期预测而发展起来的。50年代与60年代的成功应用。70年代以来人们对计量经济学模型预测功能的置疑。3、政策评价...

  • 相关主题精彩

    版权声明:本网站为非赢利性站点,内容来自于网络投稿和网络,若有相关事宜,请联系管理员

    Copyright © 喜物网