更全的杂志信息网

Time Sequence Change-Point Model of Electrostatic State Parameters of Aircraft Engine

更新时间:2016-07-05

0 Introduction

Condition monitoring and fault diagnosis of aircraft engines are effective means to guaranteeing air safety and reduce maintenance costs.Engine monitoring system (EMS),prognostics and health management[1](PHM)have been recently further developed toward condition monitoring technology to realize condition-based maintenance(CBM)[2].Although traditional monitoring and diagnosis technologies hold advantages in aviation safety and economic efficiency,they are limited by the principle of monitoring means with many defects.

In bore-detection based on internal damage detection,the light source is limited by the endoscopic diameter.The internal structure of the engine is complex,and it presents an irregular geometric shape.Under the illumination of point light sources,many shadow areas will be easily generated in the image.These drawbacks hinder the extraction of image features and the subsequent diagnosis[34].

In vibration monitoring and diagnostic technology,abnormal vibration(i.e.,vibration generated under normal engine operation)is accompanied by strong background signals,making abnormal signals difficult to capture accurately[56].In addition,measured vibration signals can reflect only the overall vibration of engine,and the abnormal vibration is concealed in the macroscopic system-level vibration.Therefore,this technology is not ideal for fault localization.

In gas path performance parameter monitoring technology[78],the deficiency of this method lies in its poor reaction capability to short-term performance change of the engine[9].Moreover,the system structure of the aircraft engine is extremely complex and the operating environment is harsh,leading to poor fault diagnosis capacity and localizing depth.

生产计划的功能是完成公司月计划的接收,并以此进行全厂生产计划的分解,生成制丝批次计划及卷包开台计划的排产与调度,指导和监督各生产车间展开生产。

Gas-path electrostatic monitoring is an approach for intelligent health management of online monitoring early-stage faults.An increasing amount of literature on its methodologies,applications,and experiments has been presented in recent years.Vatazhin et al.[10]presented the theoretical method,laboratory modeling and simulation of electrostatic monitoring for engine.The electrostatic monitoring technology was also used in engine diagnostics through monitoring the overall electrostatic charge level.Considering the great potential of electrostatic monitoring,Fisher et al.[11]conducted an on-line monitoring experiment on gas path debris using electrostatic sensors.They reported that the electrostatic monitoring technology would be an important PHM tool.Powrie et al.[12] applied this technology to the PHM system of aero-engines.The inlet debris monitoring system and exhaust debris monitoring system (EDMS),using the electrostatic monitoring,were also developed and applied to a certain type of joint-fighter.Wilcox et al.[13]investigated the application in industrial gas turbines for gas path condition monitoring and discussed the sensor installation issues.Addabbo et al.[1415]presented a theoretical modeling of an electrostatic gas path debris detection system and conducted an experimental validation.A similar research was conducted by the RMS Center of Nanjing University of Aeronautics and Astronautics in China.Wen et al.[1617]optimized the electrostatic sensor design and conducted a simulated experimental study with a simulation test bench.Fu et al.[1819]conducted a verified electrostatic monitoring experiment on a certain type of turbojet engine.Yin et al.[20]conducted a verified electrostatic monitoring experiment on a certain type of civil turbofan engine accordingly,and found a blade case rubbing fault and combustion fault.

Combining gas-path electrostatic with the change-point analysis method[21]of time sequence AR(P)model,we diagnosed the occurrence of a corresponding fault and further localized the damaged component and the damage degree.The proposed diagnosis method has the following advantages:Guaranteed air safety,reduced cost,and improved service efficiency of aircraft,and so on.

1 Formulation of Change-Point Model

The change-point model is defined as″one or a number of points that suddenly change in the model″.This problem involves a sequence of samples in a chronological order.At an unknown time t,the mathematical features or statistical distribution of these samples suddenly change.Thus,t is called the change point at this time.For convenience,in the unary linear regression changepoint model,we assume that:X is an independent variable,Y is a dependent variable,and n times of observed values are taken;i.e.,(Xi,Yi),i= 1,2,…,n.The model is set when i=1,2,…,m-1,(Xi,Yi)follows the linear regression model,that is

When imm+1,…;(Xi,Yi)follows linear regression

2)大力推行“项目驱动、案例教学、1+2+3课程改革”等教学模式,实施“1+2+3”校企深度融合课程改革,即在专业课程中遴选8-10门课程,每一门课程至少两位教师(在校教师和企业资深工程师)协作完成教学,至少有三种以上的考核方式(项目答辩、小组讨论、期中测试、期末测试、实验考核等)完成课程考核。

where eiis a random error in the model,a1and a2 are the constant terms,b1and b2the independent variable coefficients of the equation,and Yiis the i-th valuation of the dependent variable.At least one of a1=a2and b1=b2is false.At point m,if the coefficients of this regression equation are unequal,m is called the regressive change point.

2 Time Sequence Change-Point Model

2.1 Difference calculation

Most time series are non-stationary,and they could not accord with the precise demand to analyze aero-engine condition by change-point model of stationary time series.To realize a refined analysis of the aircraft engine state,stationary processing of the sample data must be conducted.At present, non-stationary time sequences are analyzed by using certain factor de-composition method,and difference method is a convenient and effective technique for information extraction.After the original data are stabilized with difference method,the stable time sequence change-point model(we take AR (P)as an example)can be used to analyze the identification of change points.

For any p-order auto-regression AR (P)process

(2)Establish time sequence AR (k)model

2.1.5 精密度试验 精密吸取混合对照品溶液10 μL,按照“2.1.1”项下条件重复进样6次,测定10种黄酮的峰面积,结果金丝桃苷、朝藿定B、朝藿定A、朝藿定C、淫羊藿苷、木犀草素、槲皮素、川陈皮素、山柰酚、宝藿苷I峰面积的RSD值分别为1.31%、0.48%、0.25%、0.63%、0.44%、1.81%、1.43%、1.11%、0.41%、0.57%,表明该条件下仪器精密度良好。

The discrete-time sequence of p order difference is equal to p order derivation of the continuous time series.According to Cramer decomposition theorem,certain information in the sample sequence{xt}can be sufficiently extracted from the p order difference.Taking first-order difference as an example

This equation indicates that the essence of the first-order difference is an auto-regression process,that is,historical data{xt}after onephase delay are used as independent variables to explain the change in date value{xt1}in the current phase.Thus

②The residual sum of squares of the AR(k)model recurred:k=1,2,…,m

where Bis the backward shift operator,BmXt=Xtm,m the time span,▽the backward difference operator,Xt=Xt-Xt1=(1-B)Xt,d the times of difference.The AR (P)model can be obtained through the deformation of Eq.(4).Under suitable difference order operation,certain formation contained in the data can be sufficiently extracted by applying Cramer decomposition theorem.

绿色矿业属于国民经济行业序列之一,绿色矿山则是实现绿色矿业的中心环节和主要载体。绿色矿山的实施主体是矿业企业,绿色矿业实施主体一般包括政府主管部门、矿业企业、矿业协会、社区组织等,比之前者涵盖范围更广。如果说绿色矿业是矿业发展的“线”的维度,那么绿色矿山则是这一线条上的“点”,也即绿色矿山是绿色矿业的基础经营单位,绿色矿业是绿色矿山的成果集成和延伸。

2.2 Model analysis

Comparisons of Figs.2(b,c)and Figs.7(b,c)show that the activity level and event rate after the engine starts and stabilizes tend to be 0.However,when the change point occurs,a corresponding activity level AL is maintained at a stable level 0.5—0.8(nc),and the proportion of the corresponding event rate reaches 50%,which is far greater than the stable status under normal operation of the engine.Actual operation of the engine shows that normal background noise signals are at milli-volt levels,but obvious abnormal pulse amplitude appears within the period of the 216th—217th test run and even volt-level signals appear.

Ifф= (φ1,…,φp)τ≠ ф′= Eq.(6)is called the discrete auto-regression time sequence change-point model,where m is the change point in the model,coefficientsф=(φ1gression parameters,and residual error atis the white noise series.

2.3 Estimation and inspection of change points

(1)Roll investigated sample sequence

First,the sample length n and the investigation interval length n0are determined.That is,the sample data sequence is set as{Xt},t=1,2,… ,N,the total data sequence length n is the length of the phase I sample,the interval n0(n0≤n≤N)is rolled to divide the total data sequence into several subsequences:The first phase {x1,x2,…,xn},the second phase {xn+1,xn+2,…,xn+n},…,the m-th phase {x(m1)n+1,x(m1)n+2,…,x(m1)n+n};the″rolling″analysis is conducted on the subsequences.

where xtis the observed value at time t;pthe order of auto-regression;φ1,…,φpare the auto-regression coefficients;andεtis the white noise series.

① When k=1,2,…,m-1,kn,parametersφk1,φk2,…,φkkof AR(k)model are calculated as follows

In Eq.(7),pjis the autocorrelation function and is estimated through the samples as follows

上述利用技术中,钢铁生产再利用和建材化利用可以消纳大量钢渣,是未来钢渣处置利用的重要途径,但各自存在如下问题:钢渣中含有P元素,会导致钢铁冶炼过程P富集,加量过多会对产品造成负面影响,因此再利用程度有限;钢渣中含有大量游离氧化物,在建材应用过程中游离氧化物会遇水体积膨胀1到2倍,内部膨胀容易导致路面开裂,制品破碎与粉化。从国外发达国家钢渣资源化处置成熟经验看,建材化利用可以有效解决钢渣利用率低的难题,但关键要解决钢渣体积稳定性差的问题。因此,探讨钢渣稳定化技术,将有助于从源头提高钢渣利用难和利用率低的难题。

第二,司法实践模糊。以涉外贴牌加工中的商标侵权问题为例,由于我国自贸区拥有自身独有的特点,这为其发展带来一系列问题。在有其他相关商标注册人时,自贸区内加工的委托方尚未取得商标注册权的产品是否能在中国大陆销售。由于市场监管的放松,在自贸区内进行加工极易造成知识产权的侵权问题。上述问题若无法厘清,实践中的执法问题就很难解决。我国《商标法》规定,若委托方不在中国取得注册权,那么无论是否在中国销售,自贸区的商品都不能进行委托加工,否则即属于商标侵权。加之由于自贸区自身特点所带来的问题,若不结合实际情况其硬性规定就会出现偏颇,所以在司法实践中认定涉外贴牌加工是否构成侵权也存有异议。

③ The amount of information BIC(k),k=1,2,…,m of the AR (k)model is calculated as follows

The solved model is named AR (P),and its parameters areφp1,φp2,…,φpp.

(3)Test judgment

Difference test is conducted on neighboring models.If the orders are the same and the coefficients are approximate,then the difference is insignificant;otherwise,Fstatistical judgment will be used

真的是“道高一尺、魔高一丈”?继“宗派色彩”“山头主义”以及“码头文化”之后,新形势下,有的地方,又出来了个“酒局圈”。“酒局圈”是个什么怪胎?近日,河南省纪委和组织部联袂出手,以治理“酒局圈”为切入口,在全省集中开展整治“帮圈文化”专项排查工作,我们才知道这个“酒局圈”究竟是何等货色!

After the first-order difference of the original time sequence data (Fig.3){Qt}of sensor in Fig.1recorded as {Qt},Qt= Qt- Qt1,{Qt}is as shown in Fig.4.

1.3 诱导培养和增殖培养条件 温度(25±2)℃,光照时间12~14 h/d,光照强度1 500~3 000 lx,光质为T4灯管发出的光。

The confidence levelαis provided in advance.The Fαvalue,which satisfies Fα(P(F≥Fα)=α),is searched through the F distribution table.When F>Fα,the subsequences of the two phases will be significantly different.

We consider one single change point as an example.If the subsequences of two neighboring phases are significantly different,t0is the terminal point of the early-phase subsequence,and the moments in section[t0-n0,t0+n0]are taken as alternate change-point moments.When an alternative change-point moment is taken as the boundary,the total data sequence is divided into two segments for modeling.The statistical quantity F will then be calculated,and then the place with the most significant before-after difference value will be taken as the change-point estimation.

3 Model Application and Discussions

3.1 Data source

The data were provided by aparticular turbojet engine test.This experiment was used to detect electrostatic particles in engine exhaust.To avoid the destruction of engine structure,the electrostatic sensor was installed at the outer bracket of the engine tailpipe.This test bench was applicable to all life tests for turbojets,and mainly used for overall performance,applicability performance,long performance and lifespan,etc.The test bench consisted of a thrust test bench system,a fuel oil supply system,and an electric control system.This experimental engine started a 200hlifespan performance test in July 2011,and a 40htest was added later,resulting in a total 240htest).The interval of the test was 1h.Electrostatic monitoring started from the 100th phase of life test and ended in the 240th phase.The first hour was removed because collection line was connected at the ground and it effectively monitored 139phases.These data were collected on November 12,2011between 08:00a.m.—10:00a.m.,and the searching interval was 1 min.The sampling site was as shown in Fig.1.

Fig.1 Schematic of actual sampling site

(1)Background signals

According to Refs.[11-12],there are two important feature parameters in engine electro-static signal analysis,namely,activity level and event rate,which are calculated as follows:

①Activity level(AL):The activity level reflects the quantity of small particles(such as soot particles and fraction-lets)that continuously appear within a period of time T,and it is a measurement of ratio-frequency components of signals.Activity level is defined as follows

where N is the number of samples within time T,and its activity level is the root-mean-square value.

② Event rate(ER):The event rate is used to measure the number of abnormal particles in airflow within unit time.Typical events include abnormal large particles generated by component fault and large soot particles generated by incomplete combustion.Event rate is expressed as follows

where M is the number of events within the time T,N the total number of samples within time T.The physical significance represented by the event rate is within a certain time interval(1s)of the electrostatic sensor.The percentage of the number of points exceeding K times of charge of current interval AL value in the total number of sampled points.

In ignition phase after the experiment starts,a large quantity of positive and negative ions were generated inside the pipeline,resulting in a drastic change in the induced voltage on the sensor,as shown in Fig.2(a),where Arepresents the amplitude of original data.After approximately 5s,the combustion tends to stabilize,and so does the induced voltage on the sensor.

Fig.2 Signal and its feature parameters in startup phase

In Fig.2(b),when the engine starts in the ignition phase,large soot particles generated by incomplete combustion are contained in tail gas,resulting in a major change in the activity level in the first 5s.When the AL charge reaches 30pc,the signal event rate ER in Fig.2(c)ascends accordingly,and the largest proportion reaches 10%.When the engine is under normal operating status(after 5s),combustion generates soot particles with minimal size(nanometer level).These soot particles are the main components of garpath charged particles.Moreover,the combustion process in the combustion chamber of the engine under normal operation is very stable.Thus,the generated soot particles are relatively stable,and the corresponding activity level and event rate tend to be 0.In the analysis chart in this paper,the activity level parameters are uniformly expressed by red solid dots″·″,and event rates are expressed by blue solid square dots“■”.

(3)Abnormal signal

This paper uses typical data monitored from 8:00to 10:00AM on November 12,2011.The value on each sampled point represents the electrostatic signal of the tail gas that passes through the sensor in last 1min.The data after difference stationary processes are taken as an example to establish the time sequence AR (P)change-point model

where Qis the tail-gas electrostatic induction signal detected by the sensor,where mis the change point in the model,coefficientsф=represent auto-regression parame-ters,and the residual error atis the white noise series.The model significance is shown in Eq.(6).Fig.3shows the original data curves of the electrostatic signals.

Fig.3 Original data curves acquired by sensor

3.2 Stabilizations method

where A0and A1are the residual sums of squares corresponding to subsequences in two different phases(calculation is shown in Eq.(9)),s is the difference of numbers of subsequence parameters in two different phases,and r is the number of high-order model parameters.

As shown in Figs.3,4,the original sequence Qthas linear ascending and descending tendencies,and sequence {Qt}after difference has been stablizied nearby a fixed value.

Fig.4 Original data first-order difference

3.3 Change-point searching algorithm flow

The change-point searching algorithm in Section 3can rapidly and accurately detect the time position at which quantitative change of electrostatic data occurs.The specific searching steps are as shown in Fig.5.

3.4 Test results

The test level is set asα=0.05.The parameters are calculated according to Eqs.(7—12)and then used to test the existence of the change points.The calculation results are shown in Table 1.

Fig.5 Change-point searching flowchart

Table 1 Change-point searching results

(1)Within 5sin the engine ignition phase,the induced voltage amplitude experiences obvious change;the corresponding activity level is 30pc and event rate is approximately 10%.After it stabilizes,both the activity level and event rate of the induced signals return to a value near 0.

Fig.6 Carbon deposit in fuel spray nozzle

Fig.7 Electrostatic signals and features at typical change-point time

Random dynamic data in time sequence are arranged according to time order.According to before-after correlation of dynamic data,a convenient and feasible AR (P)(p-order auto regression model)model is used to establish the time sequence change-point model.Its descriptions are as follows:A sequence of observed values{xt},t=1,2,…,N meet the following auto regression model

ATB混合料碾压应遵循“紧跟慢压、高频低幅”原则,压路机应紧跟摊铺机后进行碾压,避免混合料热量散失过多,初压使用13T双钢轮振动压路机静压2遍;复压采用13T振动压路机振动碾压3~4遍,碾压温度控制为100~110℃;终压使用11T钢轮压路机静压2~3遍以消除轮印,碾压温度不低于80℃。

4 Conclusions

The time sequence change-point model established in this paper rapidly and effectively detects the time when quantitative change occurs.

Change-point searching results show that change points appear five times at 8:11a.m.,8:26a.m.,8:48a.m.,9:17a.m.,and 9:41a.m.Meanwhile,the engine produces a loud,abnormal sound.Detection by test personnel after disassembly revealed that large-particle carbon deposits appear in the combustion chamber(Fig.6).The real-world scenario verifies that this paper can perform real-time monitoring of change status of gas-path charging level of the engine to provide early warning of the initial fault status.Fig.7 shows the data and signal features acquired at a typical change-point time 08:11a.m.

(2)The activity level at a typical changepoint time is maintained at 0.5—0.8 (nc)and the event rate reaches 50%,both of which are far greater than the level under stable working conditions.

生态经济首先是一个系统工程。有人把生态经济亦称为循环经济,这就更加形象地说明了他是一个系统工程,其重要目的是达到资源的循环利用,这样一方面可以达到低消耗高产出的要求,另一方面又可以达到低排放(甚至零排放)、低污染的要求。循环经济还需要运用一系列的清洁技术。经济发展严格遵循自然客观规律。一方面要严格按自然生态系统的规律办事,不损害自然生态系统。另一方面又要求给予自然生态系统积极的回报,以不断优化自然生态系统。还要求尽量利用可再生能源,以更好地节约自然资源,减少污染。

(3)Based on the above searching algorithm,results show that five change-point times,namely,8:11a.m.,8:26a.m.,8:48a.m.,9:17a.m.,and 9:41a.m.,are detected,and the corresponding jumping degrees are 0.020 189,0.014 452,0.083 381,0.147 074,and 0.099 978.Normal background noise signals of the engine are at millivolt level,but volt-level signals appear at changepoint time.This finding indicates that excess soot particles are present at change-point time,and these abnormalities change the charging level of the charged particles in the gas path.The engine disassembly report shows that large-particle carbon deposits can be found in the combustion chamber,and change-point results of electrostatic induction signals determine the corresponding fault reflection.

(4)Difference method can concisely and effectively extract certain information.After the original data are processed by the first-order differential formula,linear trend terms of this data sequence are effectively eliminated,thereby facilitating the extraction of information while the sample sequence is smoothed.

㉑㊸S.M.Cretney and J.M.Masson,Principles of Family Law,Sweet and Maxwell Publisher,1997,p.642,p.639.

(5)As an advanced early fault detection method,the change-point model can reflect the working conditions of the engine.This model is beneficial to improving monitoring techniques and realizing fault prediction and health management.

Acknowledgements

The work is supported by the Initial Scientific Research Fund (No.2015QD02S),the Foundation Research Funds for the Central Universities (No.3122016A004,3122017027),which are highly appreciated by the authors.

References:

[1] GUO Y M,CAI X B,ZHANG B Z.Review of prognostics and health management technology[J].Computer Measurement & Control,2008,16(9):1213-1216.

[2] ROEMER M J,KACPRZYNSKI G J,SCHOELLER M H.Improved diagnostic and prognostic assessments using health management information fusion[C]//2001IEEE Systems Readiness Technology Conference.Valley Forge,PA:IEEE,2001:365-377.

[3] YU H,ZUO H F,HUANG C Q.Advanced endoscopy and fault testing of aeronautic engine[J].Aviation Engineerging & Mainienance,2002(2):20-22.

[4] LI C Y,SHI H,YAO H Y.Research on borescope image properties[J].Testing Equipment & Trchnology,2006,239(5):38-40.

[5] CHEN G,LI C G,WANG D Y.Nonlinear dynamic analysis and experiment verification of rubbing faults of rotor-ball bearing-support-stator coupling system for aero-engine [J].Journal of Aerospace Power,2008,23(7):1304-1311.

[6] CHEN G.Vibration modeling and analysis for dualrotor aero-engine[J].Journal of Vibration Engineering,2011,24(6):619-632.

[7] VODOPIANOV V,NIKITIN V.A new approach to GPA-system for gas turbine engine[C]//The 39th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit.Huntsville,Alabama:AIAA,2003:4986.

[8] ROMESSIS C,MATHIOUDAKIS K.Bayesian network approach for gas path fault diagnosis[J].Journal of Engineering for Gas Turbines and Power,2006,128(1):64-72.

[9] LI Y G.A gas turbine diagnostic approach with transient measurements[J].Proceedings of the Institution of Mechanical Engineers,Part A:Journal of Power and Energy,2003,217(2):169-177.

[10]VATAZHIN A,STARIK A,KHOLSHCHEVNIKOVA E.Electric charging of soot particles in aircraft engine exhaust plumes[J].Fluid Dynamics,2004,39(3):384-392.

[11]FISHER C.Data and information fusion for gas path debris monitoring [C]// Aerospace Conference,2001,IEEE Proceedings.[S.l.]:IEEE,2001:3017-3022.

[12]POWRIE H,NOVIS A.Gas path debris monitoring for F-35Joint Strike Fighter propulsion system PHM[C]// Aerospace Conference.[S.l.]:IEEE,2006:8.

[13]WILCOX M,RANSOM D,HENRY M,et al.Engine distress detection in gas turbines with electrostatic sensors[C]//ASME Turbo Expo 2010:Power for Land,Sea,and Air.Marine Glasgow,UK:[s.n.],2010:39-51.

[14]ADDABBO T,FORT A,MUGNAINI M,et al.Theoretical modeling of an electrostatic gas-path debris detection system with experimental validation[C]∥Proceedings of 2015IEEE Sensors Applications Symposium.Zadar,Croatia:[s.n.],2015:13-15.

[15]ADDABBO T,FORT A,GARGIN R,et al.Theoretical characterization of a gas path debris detection monitoring system based on electrostatic sensors and charge amplifiers[J].Measurement,2015,64:138-146.

[16]WEN Z H,ZUO H F,LI Y H.Gas path debris electrostatic monitoring technology and experiment[J].Journal of Aerospace Power,2008,23(12):2321-2326.

[17]WEN Z,MA X,ZUO H F.Characteristics analysis and experiment verification of electrostatic sensor for aero-engine exhaust gas monitoring [J].Measurement,2014,47(1):633-644.

[18]FU Y,ZUO H F,LIU P P,et al.Gas path electrostatic sensors monitoring and comparison experiment on turbojet engine[J].Transactions of Nanjing University of Aeronautics & Astronautics,2013,30(4):361-365.

[19]FU Y,YIN Y B,WANG H,et al.Study on detection for aero-engine abnormal condition based on permutation entropy of electrostatic signal[J].Automation&Instrumentation,2017(11):1-5,19.

[20]YIN Y,CAI J,ZUO H,et al.Experimental investigation on electrostatic monitoring technology for civil turbofan engine [J].Journal of Vibroengineering,2017,19(2):967-987.

[21]FU Y.Recognition for change-point of traffic flow binary linear regressions based on projective transform method [J].Journal of Computer Applications,2010,30(1):263-262.

Fu Yu,Wei Dongdong,Zuo Hongfu,Feng Zhengxing
《Transactions of Nanjing University of Aeronautics and Astronautics》2018年第1期文献

服务严谨可靠 7×14小时在线支持 支持宝特邀商家 不满意退款

本站非杂志社官网,上千家国家级期刊、省级期刊、北大核心、南大核心、专业的职称论文发表网站。
职称论文发表、杂志论文发表、期刊征稿、期刊投稿,论文发表指导正规机构。是您首选最可靠,最快速的期刊论文发表网站。
免责声明:本网站部分资源、信息来源于网络,完全免费共享,仅供学习和研究使用,版权和著作权归原作者所有
如有不愿意被转载的情况,请通知我们删除已转载的信息 粤ICP备2023046998号