• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Recent progress on deep learning-based disruption prediction algorithm in HL-2A tokamak

    2023-09-05 08:48:08ZongyuYang楊宗諭YuhangLiu劉宇航XiaoboZhu朱曉博ZhengweiChen陳正威FanXia夏凡WulyuZhong鐘武律ZheGao高喆YipoZhang張軼潑andYiLiu劉儀
    Chinese Physics B 2023年7期
    關(guān)鍵詞:宇航

    Zongyu Yang(楊宗諭), Yuhang Liu(劉宇航), Xiaobo Zhu(朱曉博), Zhengwei Chen(陳正威), Fan Xia(夏凡),Wulyu Zhong(鐘武律),?, Zhe Gao(高喆), Yipo Zhang(張軼潑), and Yi Liu(劉儀)

    1Southwestern Institute of Physics,Chengdu 610043,China

    2Tsinghua University,Beijing 100084,China

    Keywords: macroinstabilities,tokamaks,neural networks,magnetic confinement and equilibrium

    1.Introduction

    Disruption places the tokamak-based fusion power plant a dilemma.The ignition condition and economical efficiency call for large-scale devices and high parameter plasma,[1]while in these conditions the risk of thermal load, electromagnetic force and runaway electron beam during disruptions could be unacceptable.[2]Disruption prediction and mitigation techniques take the responsibility to resolve this conflict.By accurately predicting the disruptions and mitigate them in advance,their harmful effects can be distinctly alleviated.[3,4]

    In experiment, disruptions could be induced by various of causes and event chains, e.g., mode locking, radiation collapse and MHD instabilities.[5]The complexity of their precursors makes the prediction a tough task.There are mainly 2 pathways to solve the problem.One is to design a dedicated framework to monitor disruption-related events and analysis the evolution of event chains,which leads to some representative programs like DECAF (disruption event characterization and forecasting).[6]The other way is to use data-driven algorithms, such as deep learning, to build an end-to-end model by learning from labeled disruptive and non-disruptive data.In recent years, many data-driven algorithms have been developed on various tokamaks, such as DIII-D, JET, C-Mod,J-TEXT and EAST.[7–14]

    The data-driven method has been tried in HL-2A at around 2010.A multi-layer perceptron model is built to predict the disruptions according to the evolution of signals from bolometer array.[15]However,the ability of machine learning to comprehensive analysis multi-diagnostic signals is not very mature then.The performance of algorithm is limited since only a part of the disruptions will be accompanied with precursors on bolometer array.

    Since 2018, deep learning method is tried in HL-2A.A disruption prediction algorithm based on 1.5 dimensional convolutional neural network and long short-term memory recurrent neural network (1.5D CNN+LSTM) is proposed and shows great performance by combining information from 25 disruption-related signals.[16]Therefore, a series of subsequent researches are implemented around this algorithm,such as the interpretation of algorithm,[17]real-time implementation in plasma control system (PCS),[18]cross-tokamak prediction and data-physics dual-driver algorithm.[19]These researches aim to make the disruption prediction algorithm a reliable, regular operated and operator friendly module in the PCS and provide a possible solution for future large-scale tokamaks.In this paper,we give an overview on the progresses about deep learning-based algorithm in HL-2A and discusses the possible pathways to further promote the algorithm.

    The rest part consists of 5 sections.Section 2 introduces the backbone of HL-2A’s disruption prediction algorithm.Then the backbone model is modified or combined with supplement modules to solve three key problems for disruption prediction in future fusion reactors, i.e., interpretability,real-time capacity and transferability.In Section 3, an integrated disruption prediction and mitigation system is developed in HL-2A.It realizes the real-time prediction and mitigation during online experiment, which could be a solution for the problem of real-time capacity.Section 4 gives the method and main results of model interpretation, which answers the question of interpretability.Section 5 focuses on the issue of transferability, i.e., how to develop an algorithm with limited data on future tokamaks.And three possible solutions will be introduced with preliminary implementation examples.Finally,a summary and outlook is given in Section 6.

    2.Disruption prediction backbone: 1.5D CNN+LSTM

    Disruption prediction task can also be seen as a pattern recognition of disruption precursors, which is a classical problem for machine learning techniques.Previous machine learning-based disruption prediction algorithms often use vectors consist of scalar features as their input.[7,8,10]Therefore they are not so good at dealing with high dimensional,temporal and multi-diagnostic signals, which should contain much more information.This problem is solved by deep learning recently.Fusion recurrent neural network(FRNN)introduces the profile data to promote its performance.[13]And the hybrid neural network in DIII-D proves that temporal information is pretty helpful for disruption prediction task.[14]In HL-2A,a disruption predictor based on 1.5D CNN+LSTM is developed,which takes signals from lots of different diagnostics with possibly different sample rates.Dedicated structure is designed to solve the numerical instability when training with data from different sources.Finally the algorithm gets an accuracy of 96.1%on testing set,which mainly benefits from the wide range of input signals.[16]

    In this part, Subsection 2.1 introduces the HL-2A disruption prediction dataset.Subsection 2.2 presents the data pipeline and neural network structure of HL-2A’s algorithm.Subsection 2.3 shows the accuracy of algorithm and discusses the design philosophy of 1.5D CNN+LSTM.

    2.1.Dataset

    Three aspects are selected here to give a global presentation of HL-2A’s disruption prediction dataset, i.e., shot numbers,input signal lists and typical disruption shots.

    Table 1 shows the shot numbers used for the training,validation and testing of neural network.Almost all the shots after 2013 are used without any intended selection,unless some of the input signals are lacked or obviously abnormal.Therefore the dataset consists of many types of disruptions that ever appeared in HL-2A,such as density limit,mode locking,and low safety factor on boundary.

    Table 1.Shot numbers used for training,validation,and testing of the HL-2A’s disruption predictor.

    Table 2 gives the main information of the algorithm’s 25 input signals.These signals are selected to cover the basic plasma parameters and disruption precursors as completely as possible.An important fact is that the signals have different statistical distributions and sample rates.This leads to a series of pre-processing methods and special neural network structures,which will be introduced in the subsequent part.

    Table 2.Signals in HL-2A used as the input of disruption predictor.Their sample rate,statistical features and physical meanings are also listed here.[16]

    Table 2.(Continued).

    Fig.1.The evolution of main plasma parameters, (a) plasma current, (b) line averaged electron density, (c) plasma stored energy measured by diamagnetic method, (d) magnetic perturbation measured by Mirnov probe, (e) spectrogram of the magnetic perturbation signal, during a m/n=2/1 mode locking-induced disruptive shot in HL-2A.The left subfigures plot the evolution during the whole shot.And the right subfigures plot the evolution close to disruption.A mode-locking process can be found on the Mirnov probe signal in the right bottom subfigure,which takes place at about 15 ms before the disruption.

    Figure 1 presents a mode locking-induced disruptive shot in HL-2A.The plasma current in HL-2A is normally 150 kA and the line averaged electron density is normally 0.5×1019–3×1019m?3.The time scale of mode locking process in HL-2A is around 15 ms.During disruptions, the slope of plasma current is around 10 MA/s, which leads to a current quench time of 5–20 ms.Many types of disruption precursors can be found in HL-2A’s experimental data, including but not limited to mode-locking, density limit, plasma displacement, radiation collapse and low safety factor.The variety and short timescale of disruption precursors make the disruption prediction a tough task.

    2.2.Data pipeline and neural network structure

    The working cycle of HL-2A’s disruption predictor is determined as 1 ms by the PCS.During each cycle, the algorithm takes the signals in past 20 ms from the channels given in Table 2.After that, pre-processing is taken to make these signals into structured and normalized matrix.Then the neural network takes the matrix as input and returns a possibility of disruption as output.Finally the PCS decides if disruption mitigation techniques should be triggered according to the output.

    The pre-processing consists of 3 steps.

    (a)Resampling all the signals into 100 kHz.

    (b) Reducing data from each channel in Table 2 by the channel’s mean value, and dividing them by standard deviation.

    (c)Replacing the values larger than 10 with 10 and values smaller than?10 with?10.

    After pre-processing,the input signal for each cycle will be uniformed into a 25 channels,20 ms long and 100 kHz matrix,whose size should be 2000×25.

    The structure of neural network is presented in Fig.2,which consists of convolutional layers,LSTM layers and fully connected layers.The cell number of the last fully connected layer is 2, therefore the neural network will output 2 values for each cycle, representing the possibility of disruption and non-disruption.

    Fig.2.Structure of 1.5D CNN+LSTM.This figure is adapted from Fig.3 in Ref.[16],where only the 1.5D CNN is presented.LSTM is merged into the figure after adaptation.

    2.3.The performance of 1.5D CNN+LSTM

    The performance of disruption prediction algorithm is often evaluated by two indicators, i.e., true positive rate (TPR)and true negative rate (FPR).The former refers to the possibility that the algorithm accurately gives an alarm during a disruptive shot.And the latter refers to the possibility that the algorithm gives a false alarm during a non-disruptive shot.

    In this research,the deep learning disruption predictor has a TPR of 92.2%and a FPR of 2.5%,which is a pretty good result.A key promotion is brought by its design to separately deal with signals from different diagnostics at first and merge them in the middle layer of neural network.The structure is more suitable for fusion data, which is an array consists of many heterogeneous parameters with various statistical distributions and waveforms, rather than a structured matrix like images or sentences.Table 3 gives the result of ablation experiment for this design,where 2D CNN means to merge signals from different diagnostics at the first layer of neural network,and 1D CNN means to deal with signals separately all the time.It is obvious that the 1.5D CNN gets the best performance.

    Table 3.Shot numbers used for training,validation,and testing of the HL-2A’s disruption predictor.

    3.Real-time disruption prediction and mitigation in HL-2A

    To implement the disruption prediction algorithm in the PCS, two main adjustments should be made on the version described in Section 2.Firstly, the principle to select input signals in Section 2 is to cover the disruption precursors as completely as possible.But not all these signals can be acquired in HL-2A’s PCS due to some engineering limitation of data acquisition system(DAS).Therefore the input signal list will be kind of different from the one in Table 2,which might be harmful to the algorithm’s performance.Secondly, some adjustments on neural network structure are required to accelerate its calculation speed into 1 ms per cycle.The detailed description about these updates can be found in Ref.[17].

    After that, the integrated disruption prediction and mitigation system can be built based on the real-time algorithm.Figure 3 shows the framework of this integrated system.The DAS gathers all the needed diagnostic signals and sends them to the PCS.These data are firstly used to do the position and shape control, in which process some secondary signals are produced.Then these secondary signals, together with the raw diagnostic signals, are sent to the disruption prediction module.In HL-2A, the PCS is mainly developed in C language while the disruption prediction algorithm is developed in Python, so a cross language interaction is required.Thus,there is a C-based disruption prediction module, which organizes the input data and calls the Python-based disruption prediction module.Finally, the prediction result is sent back to the C-based module to decide whether or not a trigger signal should be sent to the disruption mitigation system.

    To validate the reliability of the integrated system,open-loop online testing is tried in 382 shots during shot Nos.38650–39347.A TPR of 95.8% and a FPR of 22.5%are got during the online testing, which is acceptable but not as good as the result in Subsection 2.3 due to the limited input signal list and the unpredictable interference in DAS.

    Fig.3.Design of the real-time disruption prediction and mitigation system.[18]

    Fig.4.A demo shot of real-time disruption prediction and mitigation.A vertical displacement takes place after 800 ms due to the improper control parameters of the upper divertor coil,as indicated in the second subfigure.The output of algorithm increases rapidly when the vertical displacement amplitude becomes too large.Then the SMBI system is triggered to implement the disruption mitigation,as shown in the third subfigure.Finally a mitigated disruption occurs.Rapid decline can be seen on signals of plasma current, plasma density and bolometer radiation.[18]

    Fig.5.Comparison on the current quench and thermal quench processes between shot 39303(mitigated)and shot 39301(unmitigated).The thermal quench process is observed by upper bolometer array with 16 channels.[18]

    After the system validation, closed-loop testing is implemented.Figure 4 presents a demo shot, i.e., shot 39303,where a vertical displacement-induced disruption is predicted and then mitigated by SMBI, showing the system’s ability to real-timely predict and mitigate the coming disruptions.Figure 5 compares the thermal quench and current quench process between shot 39303 and an unmitigated disruption, shot 39301.During shot 39303, the thermal energy of plasma is released by radiation within a longer time range and the peak value is clearly lower than in the unmitigated one.[18]

    4.Prediction with interpretability and cause recognition

    In most situations,deep learning is considered as an endto-end black-box algorithm,which takes the disruption related signals as input and returns a probability of disruption without any explanation.Such an opaque system is difficult to accept for a fusion reactor with high safety requirements.In this section,an interpretation method is proposed to crack the black box.Subsection 4.1 gives the detail of the interpretation method.Subsection 4.2 shows the result of interpretation on single shot or from statistical point of view.Subsection 4.3 describes a preliminary disruption reason recognizer developed on top of the interpretation method.These results show that deep learning has the potential to realize high accuracy and interpretability at the same time.[17]

    4.1.Interpretation method

    The interpretation method is based on the 1.5-D CNN structure.In this structure a special node exists, as marked with red box in Fig.2,where the input signals have been dealt with by several neural network layers and will be concatenated in the next layer.On this node, a Gaussian noise could be added on data from each input channel to eliminate their contribution to the output.Then a distribution of correlation between each input channel and the algorithm’s output could be generated by perturbing the channels in turn and recording their contributions.

    4.2.Interpretation results

    Figure 6 gives the result of interpretation on 2 example shots.The 2 heatmaps are temporal evolution of correlation between algorithm’s output and each input channel.The brighter color means the corresponding input channel onyaxis is more related to disruption.In shot 35104, an obvious mode-locking phenomenon can be observed in the Mirnov probe signals before the disruption, as shown in the left bottom subfigures.And the rows in the left heatmap corresponding to the Mirnov probes are brighter as expected.In shot 35240,the plasma density gradually increases and approaches the Greenwald density limit before disruption.And the row corresponding to the electronic density in heat map turns to be the brightest one.The result of single shot interpretation seems to be consistent with the physical expectation.

    Fig.7.Averaged correlation distribution of input signals on each disruption type in HL-2A.Each subplot corresponds to a type of disruption.The x-axis is the input signal list and the y-axis is the averaged correlation between the input signal and a certain type of disruptions.[17]

    Statistical analysis isalso implemented to further validate the method.As a support, 613 shots in HL-2A are manually labeled into 6 disruption types, i.e., horizontal displacement,vertical displacement, radiation collapse, density limit, mode locking and low safety factor on plasma boundary.The averaged correlation distribution on each disruption type is calculated and shown in Fig.7.As expected,the most important signals for vertical displacement,lock mode,radiation,low safety factor on boundary and density limit induced disruptions areZEFIT,Mirnov probes,Pradiation/SoftX,qa,EFITandne,fir,respectively.The result of horizontal displacement seems to be kind of complex.It is suspected that other causes might also result in horizontal displacement,which calls for a further investigation.

    4.3.Disruption cause recognizer

    An important application of interpretable disruption prediction is the active disruption avoidance.As a sanity check of this technique routine, a disruption cause recognizer is developed based on the previous interpretation method.When a trigger signal is sent by disruption prediction algorithm,the recognizer will take the correlation distribution within the last 30 ms as input and try to classify the disruption types listed in Subsection 4.2.Note that the data from type of low safety factor on plasma boundary are too few to support the training of model,so this type is abandoned here.The rest dataset contains 605 shots and 5 disruption types.A naive Bayes classifier is trained on this dataset.It can classify the 5 types of disruptions with an accuracy of 71.2%in 5-fold cross validation.Figure 8 gives the algorithm’s confusion matrix during 10-fold cross validation.According to the result,it seems possible to develop a disruption cause recognizer and take measures to eliminate the disruption causes and avoid the disruptions in the future.

    Fig.8.Confusion matrix of the disruption cause recognizer.The x-axis represents the manually labeled disruption cause.And the y-axis refers to the predicted cause.So the numbers in the blocks on the diagonal line from top left to bottom right are the shot numbers that are correctly predicted,and vice versa.[17]

    5.Few-shot learning on future devices

    Since disruption prediction techniques are investigated to solve the problem of future fusion reactors, the challenge of few-shot learning becomes extremely important.Considering the unacceptable result of disruptions, future devices can only provide very limited amount of data to support the development of disruption prediction algorithm.[20]There are mainly two pathways to solve the problem within the deep learning paradigm.One is to develop a transferable algorithm, which learns data from existing tokamaks and can be adapted to future devices.Examples for this routine could be the adaptive learning algorithms,[21]FRNN with a glimpse of future device,[13]scenario adaptive learning[22]and open world learning.[23]In Subsection 5.1, a new method called device adversarial neural network (DANN) is tried and performs well when it is trained on HL-2A,J-TEXT and adapted to HL-2M.The other routine is to train a reliable algorithm with as fewer training data as possible,with the help of physical prior knowledge,some examples has been implemented in DIII-D[24]and J-TEXT.[25]Subsections 5.2 and 5.3 propose two possible physical inductive biases, i.e., plasma equilibrium equation and disruption related instabilities,that could be integrated into the disruption prediction algorithm.The future works in HL-2A will focus on these two inductive biases and preliminary works are implemented to validate the feasibility of these methods.

    5.1.Device adversarial neural network

    The researches in Refs.[13,14] prove that training with mixed dataset, which consists of a large amount of data from‘existing’ devices and a few data from ‘future’ devices is a possible way to develop a cross-tokamak algorithm.DANN makes a further step on top of this method,which uses adversarial training strategy to extract transferable features on multiple devices on purpose.Figure 9 demonstrates the structure of DANN.The main difference between DANN and normal neural networks is that it has two output branches.One branch will predict the disruptions as normal, and the other branch will recognize the source device of the input data.Note that there is a gradient reversal layer in the second branch, which will opposite the gradients back propagated from the second branch during the training of neural network.Therefore the neural network will try to predict the disruptions as accurately as possible and try to recognize the source device of data as inaccurately as possible.In other words it will recognize the disruptions by common features among the multiple devices.

    Fig.9.Structure of the device adversarial neural network.The array size of input, output and middle nodes are marked with red texts.The neural network layer parameters are marked with black texts.The position of gradient reversal layer, which is the key of DANN, is marked with blue texts.The idea is based on the research in Ref.[26].

    DANN is trained on HL-2A/J-TEXT/HL-2M mixed datasets.The main information of the mixed dataset is shown in Tables 4 and 5.The algorithm gets an accuracy of 88.9%with only 44 training shots from HL-2M.[19]However, the testing set provided by HL-2M is pretty limited now,27 shots.More data are required to implement a reliable enough test on DANN.

    Table 4.Shot numbers for training,validation,and testing of HL-2M’s DANN disruption predictor.D means disruptive and ND means non-disruptive.[19]

    Table 5.The input signal list of HL-2M’s DANN disruption predictor.[19]

    5.2.Parameter evolution equations

    Another possible way to relieve the lacking of training data on future devices is to introduce more physical information.In deep learning paradigm,the physical information can be seen as the prior knowledge in the neural networks.And the research in this chapter aims to make plasma equilibrium equations into the neural network’s prior knowledge.This purpose can be realized by two steps.Firstly,train a plasma parameter prediction algorithm, based on the parameter evolution equations.Secondly, embed the information hidden in parameter prediction algorithm into the disruption prediction algorithm.

    The plasma parameter prediction algorithm predicts the evolution of electron temperature (Te), electron density (Ne)and horizontal displacement (Dh) of HL-2A according to the control target,control actuators and the plasma state.Table 6 shows part of the input and output list of the algorithm, three empirical equilibrium equations are hidden in the design of the signal list,such as the equation for Te given here.The input and output lists of Ne and Dh are designed with similar principles.By training the neural network to do the parameter prediction task, the equation can be embedded into the neural network,and could possibly be brought into the disruption prediction algorithm in the future.

    Table 6.The hidden equation for electron temperature of the parameter prediction algorithm.

    The algorithm is based on a encoder–decoder model.The encoder extracts the characteristics of the inputs and converts the inputs into the intermediate state matrix.The decoder reconstructs the intermediate matrix and restores the encoder outputs to the signals to be predicted.The algorithm is trained by 1500 shots and tested with 200 shots.The mean absolute errors (MAE) of Ne, Te and Dh on the testing set are 0.0393 (baseline: 0.1036), 0.0359 (baseline: 0.0561) and 0.0506 (baseline: 0.0738) respectively, while the baseline is evaluated by the difference between actual parameters and the delayed parameters in the inputs.The typical result is shown in Fig.10, the output of algorithm is very close to the actual signals.It seems that the neural network model has correctly learned the plasma equilibrium equations.

    Fig.10.Panels(a),(b)and(c)are electron density,electron temperature and horizontal displacement,respectively.On the left is the comparison of algorithm inputs and target outputs,where an obvious delay of ?T can be observed.Here ?T is 30 ms for Ne and Te,while it is 10 ms for Dh.On the right is the comparison of predicted parameters and target outputs,they are much closer than the situation in the left subfigures.

    The second step, namely, to embed the prior knowledge in parameter prediction algorithm into the disruption prediction algorithm is still ongoing.A simple implementation is to concatenate the intermediate matrix given by encoder into the inputs of disruption prediction algorithm.With this auxiliary input, the AUC of disruption prediction algorithm could be promoted by 4%.

    5.3.Disruption related instabilities

    Plasma instabilities are also helpful prior knowledge for disruption prediction.[2]Therefore,researches on recognizing instabilities with neural network are also implemented in HL-2A.For example, an CNN+LSTM neural network is developed to recognize fishbone mode.The algorithm takes the signal of magnetic fluctuations measured by Mirnov coil during a 1 ms window with sample rate of 1 MHz,i.e.,a 1000 length vector.The vector is firstly dealt with by mean-std normalization.Then a Fourier transformation with frame length of 1000 and frame step of 100 is implemented to get a 501×91 spectrum.Since the fishbone in HL-2A usually has a frequency between 5–35 kHz,only the spectrum in this frequency range is reserved.Finally the 30×91 spectrum is fed into a neural network as shown in Fig.11.After training on 770 shots,the neural network can real-timely recognize if there is a fishbone instability in the plasma with an accuracy of 95.1%.Figure 12 presents the algorithm’s output during an example shot.Similar works has been implemented on tearing modes,sawtooth,ELM and other instabilities.Future works will focus on combining these modules into the disruption prediction algorithm and check if this could alleviate the algorithm’s demand for training data.

    Fig.12.An example shot of fishbone instability recognition algorithm.

    6.Summary and future works

    The deep learning-based disruption prediction algorithm has proved its powerful ability in HL-2A tokamak.With dedicated neural network design, optimized training strategy and novel interpretation method, the algorithm can realize high accuracy,good interpretability and priliminary cross-tokamak adaptation at the same time.Real-time disruption prediction and mitigation experiment is taken in HL-2A to validate the algorithm and the result proves its reliability.In general,deep learning can be considered as a sophisticated solution for disruption prediction in‘existing’tokamaks,both from a research and engineering perspective.Therefore the future works will go beyond the ‘prediction on existing tokamaks’.There are mainly 2 pathways, i.e., ‘prediction on future tokamaks’ and‘prevention on existing tokamaks’.

    To predict the disruption in‘future’tokamaks,the preliminary research in Section 5 will be carried on.Domain adaptation techniques and physical inductive bias will be introduced to alleviate the algorithm’s demand for training data.HL-2M,a newly built tokamak which has realized the first plasma in 2020, has promoted its plasma current to 1 MA recently and will provide a great testing platform for these researches.

    The other pathway is to prevent the disruptions based on interpretable disruption prediction algorithm.A correlation analysis has been realized on HL-2A’s algorithm in Section 4.However, the method needs to go beyond the correlation and realize causality analysis to find the root cause from a bunch of correlated phenomena.A cause-and-effect diagram needs to be designed with consideration of the callable control actuators.Therefore,the root cause can be matched to the control measures to support the disruption prevention.

    In general,deep learning provides a potential solution for the disruption prediction,prevention and mitigation on future tokamaks.The basic feasibility has been proved in HL-2A and many other tokamaks in recent years.Some bottleneck problems,such as cross tokamak adaptation,disruption cause recognition still exist, but are being investigated and preliminary solved.It can be expected that there will be more breakthroughs in this direction in the next few years.

    Acknowledgements

    Project supported by the National MCF R&D Program of China (Grant Nos.2018YFE0302100 and 2019YFE03010003).The authors wish to thank all the members at South Western Institute of Physics for providing data,technique assistance and co-operating during the experiment.

    猜你喜歡
    宇航
    《宇航計(jì)測(cè)技術(shù)》征稿簡(jiǎn)則
    狼蛛
    Steered coherence and entanglement in the Heisenberg XX chain under twisted boundary conditions*
    《宇航計(jì)測(cè)技術(shù)》2020年第6期目次
    我的收獲
    《宇航計(jì)測(cè)技術(shù)》征稿啟事
    宇航人 夢(mèng)想飛揚(yáng)
    2016年中國(guó)宇航發(fā)射將首超20次
    我的宇航夢(mèng)
    我的宇航夢(mèng)
    国产久久久一区二区三区| 身体一侧抽搐| 岛国毛片在线播放| 在线精品无人区一区二区三 | 亚洲av成人精品一二三区| 日韩制服骚丝袜av| 亚洲国产精品专区欧美| 国产真实伦视频高清在线观看| 中文欧美无线码| 又大又黄又爽视频免费| 一区二区三区四区激情视频| 一级av片app| h视频一区二区三区| 五月开心婷婷网| 久久99精品国语久久久| 深爱激情五月婷婷| 男人和女人高潮做爰伦理| 中文字幕亚洲精品专区| 女人久久www免费人成看片| 中文字幕av成人在线电影| 日韩伦理黄色片| av一本久久久久| 成年人午夜在线观看视频| 精品亚洲成国产av| 肉色欧美久久久久久久蜜桃| 国产精品秋霞免费鲁丝片| 校园人妻丝袜中文字幕| 特大巨黑吊av在线直播| 国产伦理片在线播放av一区| 久久国产乱子免费精品| 亚洲欧美日韩卡通动漫| 干丝袜人妻中文字幕| 亚洲av中文字字幕乱码综合| 欧美xxⅹ黑人| 99久久中文字幕三级久久日本| 黄色怎么调成土黄色| 亚洲精品成人av观看孕妇| 成人一区二区视频在线观看| 国产成人精品久久久久久| 亚洲精品乱久久久久久| 女性被躁到高潮视频| tube8黄色片| 久久国产精品大桥未久av | 亚洲内射少妇av| 日韩免费高清中文字幕av| 伊人久久国产一区二区| 日本黄大片高清| 日本wwww免费看| 色视频在线一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美日韩无卡精品| 久久97久久精品| 99热6这里只有精品| 国产色爽女视频免费观看| 伦理电影大哥的女人| 成人黄色视频免费在线看| 黄色怎么调成土黄色| 99热这里只有是精品在线观看| 高清在线视频一区二区三区| 一级毛片黄色毛片免费观看视频| 99热国产这里只有精品6| 午夜免费鲁丝| 热re99久久精品国产66热6| 国产精品国产av在线观看| 午夜福利高清视频| 免费久久久久久久精品成人欧美视频 | 亚洲色图综合在线观看| 纯流量卡能插随身wifi吗| 内射极品少妇av片p| 91久久精品电影网| 十分钟在线观看高清视频www | 亚洲高清免费不卡视频| 精品人妻一区二区三区麻豆| 亚洲欧美一区二区三区黑人 | 亚洲aⅴ乱码一区二区在线播放| 亚洲精品成人av观看孕妇| 午夜免费观看性视频| 91久久精品电影网| 国产欧美另类精品又又久久亚洲欧美| 大香蕉97超碰在线| 亚洲丝袜综合中文字幕| 亚洲精品,欧美精品| 久久影院123| 亚洲精品第二区| a级一级毛片免费在线观看| 插逼视频在线观看| 丰满人妻一区二区三区视频av| 久久97久久精品| 午夜精品国产一区二区电影| 色综合色国产| 亚洲怡红院男人天堂| 亚洲成人手机| 亚洲av成人精品一区久久| 黄色日韩在线| 亚洲精品乱久久久久久| 亚洲欧美一区二区三区国产| 不卡视频在线观看欧美| 色婷婷av一区二区三区视频| 国产精品嫩草影院av在线观看| 国产乱人偷精品视频| 人体艺术视频欧美日本| 美女国产视频在线观看| 国产 精品1| 国产精品国产av在线观看| 亚洲真实伦在线观看| 男女边吃奶边做爰视频| 啦啦啦啦在线视频资源| 国产成人freesex在线| 26uuu在线亚洲综合色| 91久久精品国产一区二区成人| 狂野欧美激情性bbbbbb| 国产毛片在线视频| 一区二区av电影网| 日本欧美国产在线视频| 国产色爽女视频免费观看| 80岁老熟妇乱子伦牲交| 欧美人与善性xxx| 国产欧美亚洲国产| 永久免费av网站大全| 亚洲av成人精品一二三区| 久久影院123| 毛片女人毛片| 最近2019中文字幕mv第一页| 亚洲成人一二三区av| 99久久人妻综合| 男人添女人高潮全过程视频| 亚洲真实伦在线观看| 国国产精品蜜臀av免费| 秋霞在线观看毛片| 九草在线视频观看| 精品午夜福利在线看| 国产精品一区www在线观看| 80岁老熟妇乱子伦牲交| 少妇人妻一区二区三区视频| 欧美zozozo另类| 欧美日韩亚洲高清精品| 精品熟女少妇av免费看| 久久久久人妻精品一区果冻| 蜜桃久久精品国产亚洲av| 不卡视频在线观看欧美| 国产欧美另类精品又又久久亚洲欧美| 狂野欧美激情性bbbbbb| 国产毛片在线视频| 国产精品一及| 国产片特级美女逼逼视频| 久热这里只有精品99| 男女啪啪激烈高潮av片| 亚洲,一卡二卡三卡| 97在线人人人人妻| 久久久亚洲精品成人影院| 噜噜噜噜噜久久久久久91| 在线免费观看不下载黄p国产| 青春草视频在线免费观看| 交换朋友夫妻互换小说| 高清不卡的av网站| 大又大粗又爽又黄少妇毛片口| 国产女主播在线喷水免费视频网站| 一本—道久久a久久精品蜜桃钙片| 好男人视频免费观看在线| 最近最新中文字幕大全电影3| 国产熟女欧美一区二区| 激情 狠狠 欧美| 中文字幕亚洲精品专区| 欧美精品人与动牲交sv欧美| 久久青草综合色| 久久久久久久大尺度免费视频| 久久精品国产亚洲av涩爱| 日韩精品有码人妻一区| 中文乱码字字幕精品一区二区三区| 日本欧美视频一区| 国产精品成人在线| 天天躁日日操中文字幕| 亚洲人成网站高清观看| 91午夜精品亚洲一区二区三区| 多毛熟女@视频| 啦啦啦啦在线视频资源| 国产男女超爽视频在线观看| 夫妻午夜视频| 久久久久久久久久久丰满| 交换朋友夫妻互换小说| 午夜福利视频精品| 国产午夜精品一二区理论片| 日韩不卡一区二区三区视频在线| 99久国产av精品国产电影| 小蜜桃在线观看免费完整版高清| 蜜臀久久99精品久久宅男| 亚洲欧美日韩卡通动漫| 一个人免费看片子| 亚洲激情五月婷婷啪啪| 99久久综合免费| 黄色欧美视频在线观看| 五月伊人婷婷丁香| 最后的刺客免费高清国语| tube8黄色片| 日产精品乱码卡一卡2卡三| 久久国产亚洲av麻豆专区| 精品久久久噜噜| 亚洲av国产av综合av卡| 插阴视频在线观看视频| 伊人久久精品亚洲午夜| 一本色道久久久久久精品综合| 久久久久网色| 欧美激情国产日韩精品一区| 自拍欧美九色日韩亚洲蝌蚪91 | 观看免费一级毛片| 日韩av不卡免费在线播放| 国产黄片视频在线免费观看| 日韩电影二区| 国产伦精品一区二区三区视频9| 欧美成人a在线观看| 国产在线一区二区三区精| 中文天堂在线官网| 欧美xxxx黑人xx丫x性爽| 99热这里只有是精品50| 欧美日韩视频精品一区| 一级毛片我不卡| 久久毛片免费看一区二区三区| 少妇人妻一区二区三区视频| 18禁在线无遮挡免费观看视频| 中文字幕人妻熟人妻熟丝袜美| 永久免费av网站大全| 嫩草影院新地址| 中文在线观看免费www的网站| 制服丝袜香蕉在线| 人人妻人人添人人爽欧美一区卜 | 男女边摸边吃奶| 国产高清有码在线观看视频| 亚洲人成网站高清观看| 亚洲欧美精品专区久久| 成人无遮挡网站| 久久久久久久亚洲中文字幕| 日本色播在线视频| 51国产日韩欧美| 亚洲精华国产精华液的使用体验| 中文欧美无线码| 免费观看的影片在线观看| 大又大粗又爽又黄少妇毛片口| 汤姆久久久久久久影院中文字幕| 欧美成人a在线观看| 夜夜爽夜夜爽视频| 久久精品国产自在天天线| videossex国产| 午夜视频国产福利| 久久99蜜桃精品久久| 一区二区三区免费毛片| 在线观看免费视频网站a站| 又大又黄又爽视频免费| 国产精品99久久99久久久不卡 | 欧美97在线视频| 欧美日韩精品成人综合77777| 国内精品宾馆在线| 人体艺术视频欧美日本| 国产成人a∨麻豆精品| 亚洲av中文字字幕乱码综合| 欧美日韩国产mv在线观看视频 | 国产精品一区二区性色av| 国产亚洲最大av| 欧美成人午夜免费资源| 熟女电影av网| 国产精品一及| 国内精品宾馆在线| av视频免费观看在线观看| 亚洲精品久久久久久婷婷小说| 成人漫画全彩无遮挡| 久热久热在线精品观看| 欧美激情国产日韩精品一区| 自拍偷自拍亚洲精品老妇| 男女国产视频网站| 在线观看三级黄色| 有码 亚洲区| 黄色怎么调成土黄色| 老熟女久久久| 欧美激情极品国产一区二区三区 | 男女边吃奶边做爰视频| 亚洲av不卡在线观看| 亚洲第一av免费看| 这个男人来自地球电影免费观看 | av在线观看视频网站免费| 中文欧美无线码| 18禁在线无遮挡免费观看视频| 久久毛片免费看一区二区三区| 99久久精品热视频| 日韩欧美精品免费久久| 99re6热这里在线精品视频| 最新中文字幕久久久久| 美女高潮的动态| 国产探花极品一区二区| 菩萨蛮人人尽说江南好唐韦庄| 最后的刺客免费高清国语| 中文资源天堂在线| 少妇人妻 视频| 最近手机中文字幕大全| 在线播放无遮挡| 亚洲熟女精品中文字幕| 久久国内精品自在自线图片| 亚洲色图av天堂| 一本—道久久a久久精品蜜桃钙片| 麻豆成人午夜福利视频| 久久毛片免费看一区二区三区| 交换朋友夫妻互换小说| 偷拍熟女少妇极品色| 久久精品久久久久久久性| h日本视频在线播放| 国产成人免费观看mmmm| 国产爽快片一区二区三区| 国产精品一区二区性色av| 国产精品爽爽va在线观看网站| 免费看av在线观看网站| 大陆偷拍与自拍| 有码 亚洲区| 免费看av在线观看网站| 日日啪夜夜撸| 免费观看无遮挡的男女| 国产淫语在线视频| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产最新在线播放| 99国产精品免费福利视频| 国产伦精品一区二区三区视频9| 日韩一区二区三区影片| 国产av精品麻豆| 亚洲精品国产色婷婷电影| 蜜桃亚洲精品一区二区三区| 黑丝袜美女国产一区| 多毛熟女@视频| 欧美极品一区二区三区四区| 日韩国内少妇激情av| 伦理电影免费视频| 国产在线一区二区三区精| freevideosex欧美| 久久国产亚洲av麻豆专区| 成年女人在线观看亚洲视频| 国产欧美另类精品又又久久亚洲欧美| 99热这里只有是精品在线观看| 精品国产一区二区三区久久久樱花 | 在线观看一区二区三区| 51国产日韩欧美| 99视频精品全部免费 在线| 午夜激情久久久久久久| 国产v大片淫在线免费观看| 纵有疾风起免费观看全集完整版| 国产黄片视频在线免费观看| 国产av一区二区精品久久 | 一区二区三区乱码不卡18| 午夜免费男女啪啪视频观看| 亚洲一级一片aⅴ在线观看| 高清在线视频一区二区三区| 精华霜和精华液先用哪个| 国产免费一区二区三区四区乱码| 国产精品99久久久久久久久| 色婷婷av一区二区三区视频| 免费看日本二区| 一个人看的www免费观看视频| 免费大片18禁| 国产一区有黄有色的免费视频| freevideosex欧美| 中文字幕精品免费在线观看视频 | 最黄视频免费看| 18禁裸乳无遮挡免费网站照片| 成人一区二区视频在线观看| 亚洲精品aⅴ在线观看| 欧美日韩国产mv在线观看视频 | 自拍偷自拍亚洲精品老妇| 日韩免费高清中文字幕av| 成人午夜精彩视频在线观看| 26uuu在线亚洲综合色| 狂野欧美激情性bbbbbb| 国产男女超爽视频在线观看| 高清不卡的av网站| 在线精品无人区一区二区三 | 亚洲精品一区蜜桃| 国产黄频视频在线观看| 国产人妻一区二区三区在| 大话2 男鬼变身卡| 最近手机中文字幕大全| 久久久久久久久久成人| 欧美日韩视频高清一区二区三区二| 男人狂女人下面高潮的视频| 中文字幕av成人在线电影| 在现免费观看毛片| 成人免费观看视频高清| 18禁在线播放成人免费| 中文在线观看免费www的网站| 免费高清在线观看视频在线观看| 夜夜骑夜夜射夜夜干| 日韩视频在线欧美| 国产成人91sexporn| 一区二区三区乱码不卡18| 久久久久久久久久成人| 亚洲精品色激情综合| 久久午夜福利片| 国产精品麻豆人妻色哟哟久久| 大陆偷拍与自拍| www.色视频.com| 久久精品夜色国产| 久久这里有精品视频免费| 我的老师免费观看完整版| 国产乱人视频| 日韩强制内射视频| 国产精品伦人一区二区| 精品少妇黑人巨大在线播放| 99热国产这里只有精品6| 水蜜桃什么品种好| 午夜免费男女啪啪视频观看| 黄片无遮挡物在线观看| 美女cb高潮喷水在线观看| 亚洲av不卡在线观看| 国产高清不卡午夜福利| 国产精品爽爽va在线观看网站| 在线免费十八禁| 另类亚洲欧美激情| 亚洲精品乱久久久久久| 水蜜桃什么品种好| 精品国产一区二区三区久久久樱花 | 日韩欧美一区视频在线观看 | 亚洲最大成人中文| 亚洲精品乱久久久久久| 国产大屁股一区二区在线视频| 两个人的视频大全免费| 亚洲欧美中文字幕日韩二区| 一级毛片黄色毛片免费观看视频| 最后的刺客免费高清国语| 亚洲国产日韩一区二区| 一二三四中文在线观看免费高清| 在线观看免费视频网站a站| 国产黄色免费在线视频| 啦啦啦啦在线视频资源| 国内揄拍国产精品人妻在线| 国产久久久一区二区三区| 少妇人妻 视频| 国产日韩欧美亚洲二区| 日日撸夜夜添| 美女中出高潮动态图| 亚洲精品国产色婷婷电影| 国产乱人偷精品视频| 亚洲自偷自拍三级| 亚洲av中文字字幕乱码综合| 777米奇影视久久| 小蜜桃在线观看免费完整版高清| 国产男人的电影天堂91| 小蜜桃在线观看免费完整版高清| 一本色道久久久久久精品综合| 中文字幕人妻熟人妻熟丝袜美| 免费观看a级毛片全部| av线在线观看网站| 欧美精品国产亚洲| 国产美女午夜福利| h日本视频在线播放| 少妇裸体淫交视频免费看高清| 亚洲第一av免费看| 久久久成人免费电影| 老女人水多毛片| 99热国产这里只有精品6| 高清欧美精品videossex| 国产爱豆传媒在线观看| 日本一二三区视频观看| 永久网站在线| 午夜福利影视在线免费观看| 在线天堂最新版资源| 欧美日韩国产mv在线观看视频 | 国产欧美另类精品又又久久亚洲欧美| 久久99热6这里只有精品| 国产一区二区在线观看日韩| 在线亚洲精品国产二区图片欧美 | 精品久久久噜噜| 国产成人免费观看mmmm| 伊人久久国产一区二区| 九色成人免费人妻av| 老女人水多毛片| 国产精品av视频在线免费观看| av黄色大香蕉| 成人国产麻豆网| 国内精品宾馆在线| 日日啪夜夜撸| 99热这里只有是精品在线观看| 国产真实伦视频高清在线观看| 亚洲人与动物交配视频| 亚洲欧美一区二区三区黑人 | 亚洲aⅴ乱码一区二区在线播放| 中文在线观看免费www的网站| 国产免费福利视频在线观看| 深爱激情五月婷婷| 熟妇人妻不卡中文字幕| 亚洲精品乱码久久久久久按摩| 99热网站在线观看| 国产人妻一区二区三区在| 高清黄色对白视频在线免费看 | 九九在线视频观看精品| 一本—道久久a久久精品蜜桃钙片| 亚洲国产精品成人久久小说| 在线精品无人区一区二区三 | 舔av片在线| 国产成人a区在线观看| 亚洲精品亚洲一区二区| 综合色丁香网| 国产在线一区二区三区精| 午夜日本视频在线| 毛片一级片免费看久久久久| 蜜桃在线观看..| 2018国产大陆天天弄谢| 超碰97精品在线观看| 天堂中文最新版在线下载| 在线观看美女被高潮喷水网站| 亚洲精品国产av蜜桃| 国产男女超爽视频在线观看| 啦啦啦视频在线资源免费观看| 久久国产精品男人的天堂亚洲 | 精品久久久久久久末码| 伊人久久精品亚洲午夜| 欧美一区二区亚洲| 久久女婷五月综合色啪小说| 高清欧美精品videossex| 全区人妻精品视频| 亚洲精品日本国产第一区| 亚洲欧美一区二区三区黑人 | 少妇人妻 视频| 亚洲精品国产av蜜桃| 少妇丰满av| 国产精品伦人一区二区| 3wmmmm亚洲av在线观看| 老司机影院成人| 亚洲av欧美aⅴ国产| 在线观看免费高清a一片| 两个人的视频大全免费| 大香蕉97超碰在线| 黄色日韩在线| 美女cb高潮喷水在线观看| 国产av码专区亚洲av| 亚洲精品第二区| 一区二区三区精品91| 亚洲美女搞黄在线观看| 亚洲伊人久久精品综合| 在线亚洲精品国产二区图片欧美 | 欧美老熟妇乱子伦牲交| 超碰av人人做人人爽久久| 亚洲在久久综合| 一区二区三区四区激情视频| av视频免费观看在线观看| 菩萨蛮人人尽说江南好唐韦庄| 欧美成人精品欧美一级黄| 91午夜精品亚洲一区二区三区| 国产精品久久久久久精品电影小说 | 亚洲av.av天堂| 久久久久久久亚洲中文字幕| 免费大片18禁| 亚洲精华国产精华液的使用体验| 国产国拍精品亚洲av在线观看| 免费少妇av软件| 噜噜噜噜噜久久久久久91| 中国美白少妇内射xxxbb| 亚洲精品久久午夜乱码| 午夜福利网站1000一区二区三区| 男人爽女人下面视频在线观看| 亚洲精品视频女| 久久久久性生活片| 日韩不卡一区二区三区视频在线| 国产成人a区在线观看| 欧美精品国产亚洲| 国产一区二区三区综合在线观看 | 日本黄色片子视频| 成人国产麻豆网| 老女人水多毛片| 黄色怎么调成土黄色| 国产在线免费精品| 黄色视频在线播放观看不卡| av在线观看视频网站免费| 久久久久久人妻| 黄片wwwwww| 国产大屁股一区二区在线视频| 亚洲精品国产色婷婷电影| 日韩国内少妇激情av| 亚洲欧美中文字幕日韩二区| 少妇人妻一区二区三区视频| 嘟嘟电影网在线观看| 亚洲精品aⅴ在线观看| 国产欧美日韩一区二区三区在线 | 人妻 亚洲 视频| 精品少妇久久久久久888优播| 亚洲精华国产精华液的使用体验| 国产黄片视频在线免费观看| 免费av不卡在线播放| 在线观看美女被高潮喷水网站| 九草在线视频观看| 男男h啪啪无遮挡| 亚洲激情五月婷婷啪啪| 99热这里只有是精品在线观看| 欧美日韩一区二区视频在线观看视频在线| 国产精品久久久久久精品电影小说 | 纵有疾风起免费观看全集完整版| 国产69精品久久久久777片| 国产一区二区三区av在线| 久久精品国产a三级三级三级| 国产亚洲最大av| 久久鲁丝午夜福利片| 国产人妻一区二区三区在| 国产午夜精品久久久久久一区二区三区| 女的被弄到高潮叫床怎么办| 国产成人免费观看mmmm| 联通29元200g的流量卡| av卡一久久| 国产熟女欧美一区二区| 五月天丁香电影| 国产伦精品一区二区三区视频9| 51国产日韩欧美| 亚洲国产av新网站| 亚洲第一av免费看| 毛片一级片免费看久久久久| 亚洲欧美一区二区三区黑人 | 日韩一区二区三区影片| 国产真实伦视频高清在线观看| 欧美一级a爱片免费观看看| 国国产精品蜜臀av免费| 国产探花极品一区二区| 久久久色成人| 亚洲精品自拍成人| 国产成人a区在线观看|