Disclosure of Invention
The invention aims to solve the defects in the prior art, and provides a photosynthetic organism carbon fixation evaluation system under the condition of coking flue gas.
In order to achieve the purpose, the invention adopts the following technical scheme that the photosynthetic organism carbon fixation evaluation system under the condition of coking flue gas comprises:
The carbon kinetic energy evolution monitoring module acquires injection data of the coking flue gas, extracts a time node with the largest mutation amplitude in each gas injection period according to the carbon dioxide migration rate in a unit volume of a microalgae culture system in the gas injection process, constructs a trend track and generates a carbon kinetic energy evolution track;
The air pressure permeation threshold extraction module invokes the carbon kinetic energy evolution track, selects an inlet static pressure change curve in an injection period, monitors the fluctuation range of the curve in unit time before and after an inflection point, and combines the occurrence frequency and the duration time of a plurality of period lower limit values to generate a carbon migration trigger boundary group;
the spectral response adsorption evaluation module calls the carbon migration triggering boundary group, extracts the adsorption change trend, judges the reaction continuity under the continuous wave band, marks the adsorption active section in the curve, and combines the duration and the response change performance for evaluation to generate the illumination response adsorption section;
and the flux adaptation change comparison module calls the illumination response adsorption section, and performs trend comparison analysis according to a carbon fixation flux data curve of unit volume in a coverage time period and in combination with synchronous response in a carbon kinetic energy evolution track in coking flue gas injection response, so as to generate carbon fixation adaptation change information.
As a further scheme of the invention, the carbon kinetic energy evolution track comprises carbon migration critical time point distribution, a gas injection period evolution trend curve andThe carbon migration triggering boundary group comprises an interface pressure limit value set, pressure fluctuation threshold value classification and gas injection period pressure activation mapping, the illumination response adsorption section comprises an active spectrum bandwidth section, an adsorption intensity distribution label and illumination section duration classification, and the carbon fixation adaptation change information comprises a flux trend response offset set, an adsorption behavior response comparison group and an adaptation difference classification result.
As a further aspect of the present invention, the carbon kinetic energy evolution monitoring module includes:
The injection parameter acquisition submodule acquires injection data of the coking flue gas, acquires carbon dioxide migration rate in a unit volume in a corresponding time period in the microalgae culture system, monitors partial pressure change rate of carbon dioxide in a gas injection period, establishes a ratio sequence according to a numerical ratio between the carbon migration rate in the unit volume and the partial pressure change rate, and generates a ratio sequence integration value;
the critical inflection point identification submodule extracts a time node with the largest mutation amplitude in each curve based on the ratio sequence integrated value, judges whether the direction is ascending and is in a local fluctuation high position, integrates the time sequence distribution of the marking nodes, and obtains a migration inflection point distribution section;
And the trend track construction submodule calls the migration inflection point distribution section, performs difference operation on migration inflection point position time and gas injection starting time under the same period, classifies and maps trend change directions to corresponding period groups according to load grades, establishes a change evolution path under injection conditions, and generates a carbon kinetic energy evolution track.
As a further aspect of the present invention, the pneumatic permeation threshold extraction module includes:
The pressure curve extraction submodule calls the carbon kinetic energy evolution track, acquires an inlet static pressure change curve in an injection period where each time point is located according to a marked migration inflection point time point, intercepts a symmetrical time section taking the inflection point as a center on each curve, establishes a periodic static pressure fluctuation data set and generates a static pressure change sequence group;
The fluctuation section identification submodule judges whether the difference sign of adjacent pressure values in unit time is positive continuously according to each sequence in the static pressure change sequence group, calculates the length of a time segment with positive continuity, extracts the pressure value corresponding to the starting moment of the section, and obtains a starting pressure set of a continuous rising section;
the boundary limit classifying sub-module calls the continuous ascending segment starting pressure set, counts the frequency of occurrence in a plurality of injection periods, records the duration range of each pressure value, constructs a combined characteristic sequence according to the frequency and the duration, classifies the combined characteristic sequence according to the stability sequence of the combined characteristic sequence, and then classifies the combined characteristic sequence into a boundary condition set to generate a carbon migration triggering boundary group.
As a further aspect of the present invention, the spectral response adsorption assessment module includes:
the flux density acquisition submodule acquires the carbon migration trigger boundary group, and acquires photon flux density and alga surface unit area under a plurality of wave bands in each periodThe adsorption rate increases the data, establishes a photon flux density and adsorption rate corresponding relation matrix, extracts a change curve of the corresponding relation in a spectrum band section, and generates a spectrum response curve group;
The adsorption trend identification submodule extracts a band section meeting continuous rising according to the spectrum response curve group, calculates an adsorption intensity response variation value in the continuous band section, screens and marks a band with the variation value being greater than a continuous response reference intensity difference, and generates an adsorption active section set;
and the response section evaluation submodule calls the adsorption active section set, extracts the photon flux density change quantity of the corresponding time segment of the section and the adsorption quantity change trend of unit time, constructs an active wave band adsorption curve set, and establishes an active response wave band time coverage map to obtain the illumination response adsorption section.
As a further aspect of the present invention, the formula for calculating the adsorption strength response variation value in the continuous band section specifically includes:
;
Wherein, theRepresents the variation value of the adsorption response,Represents the firstNormalized value of CO2 adsorption rate variation corresponding to each wave band,Represents the firstNormalized value of adsorption duration per unit area in each band,Represents the firstThe center wavelength of the individual wavelength bands,Represents the firstNormalized values of photon flux density corresponding to the individual bands,Is the total number of extracted bands.
As a further aspect of the present invention, the flux adaptation change control module includes:
Based on the illumination response adsorption section, the carbon fixation data extraction submodule collects continuous data curves of the carbon fixation flux of unit volume in a time period, constructs a flux change sequence, classifies and integrates each section of data by combining with timestamp information, establishes a flux performance set and generates a carbon fixation flux data set;
The trend response comparison submodule calls the carbon fixation flux data set, extracts synchronous flux change rate and kinetic energy fluctuation change rate sequences, calculates a fluctuation trend difference value after the two types of sequences are matched according to time indexes, and calculates and obtains a flux response offset value;
And the adaptation performance merging sub-module extracts corresponding fluctuation interval data according to the flux response offset value, screens continuous sections with the variation amplitude smaller than a response tolerance threshold, performs time statistics on the screened sections, and constructs a flux response behavior classification set to obtain carbon fixation adaptation variation information.
As a further aspect of the present invention, the formula for calculating the flux response offset value is specifically:
;
Wherein, theIndicating the value of the flux response offset,Represent the firstThe carbon fixation flux per unit volume at the moment,Represents the maximum value of the carbon fixation flux in the same period,Represent the firstThe kinetic energy value of carbon per unit volume at the moment,Represents the maximum value of the kinetic energy of carbon in the same period,Representing the total number of time points involved in the matching.
As a further aspect of the present invention, the system further includes:
The carbon fixation stable fluctuation indicating module invokes the carbon fixation adaptation change information, extracts amplitude values corresponding to five time points of each period according to the carbon fixation rate change amplitude of each stage, performs extremely poor extraction on the difference value, marks the most severe period of fluctuation, counts the fluctuation duration ratio in the period, and takes the classification summary of the ratio sections as fluctuation indicating sections to generate a carbon fixation stable fluctuation section;
the carbon fixation stable fluctuation section comprises a stability grade marking section, a fluctuation continuous proportion section and a periodic fluctuation structure identification tag.
As a further aspect of the present invention, the carbon sequestration stabilization fluctuation indicating module includes:
The amplitude difference extraction submodule acquires the carbon fixation adaptation change information, extracts speed amplitude values corresponding to five fixed time points in each period, calculates the difference between every two time points, extracts the maximum value and the minimum value, performs the range operation between the difference values, gathers the range output values of the periods, and generates a carbon fixation fluctuation range set;
the fluctuation severe identification submodule calls the carbon fixation fluctuation extremely difference set, screens the period which is larger than the fluctuation identification standard extremely difference, marks the period as a high variation period section, and extracts the start and stop time of the marked period and the carbon fixation response curve segment to obtain a high variation period index group;
And the period section gathering sub-module is used for counting the proportion of the fluctuation duration of the fixed carbon rate in each period to the total operation period according to the corresponding time range of each period based on the high variation period index group, classifying the intervals between the duty ratio data and the set fluctuation duration ratio threshold value, and establishing a fixed carbon stable fluctuation section.
Compared with the prior art, the invention has the advantages and positive effects that:
According to the invention, through real-time tracking of the abrupt change amplitude and change direction of the carbon dioxide migration rate in the coking flue gas, the carbon migration inflection point is accurately captured, the boundary pressure for starting the carbon migration reaction is determined and extracted based on the static pressure change fluctuation range at the inflection point, and the migration pressure threshold of the gas-liquid interface is clearly defined, so that the coking flue gas can be more accurately mastered in the actual operation processThe critical condition of migration to the microalgae system is achieved by fusing the photon flux density and the microalgae under the spectral bandThe response relation of the adsorption rate is clear, the corresponding characteristics of the light condition and the algae adsorption reaction in the continuous wave band section are clear, the specific spectral adsorption response section is determined, the effective distinction and evaluation of the carbon fixation efficiency under different light conditions are realized, the trend contrast analysis is carried out by combining the carbon fixation flux data and the carbon kinetic energy evolution track at different stages, the adaptability change of the carbon fixation process to the coking flue gas injection period can be accurately judged, the extremely poor analysis is carried out on the variation amplitude of the carbon fixation rate in different periods, the fluctuation section with the most severe carbon fixation reaction is clear, the biological carbon fixation process evaluation of the coking flue gas is more targeted, the dynamic handling capability of the carbon fixation reaction stability in the operation process is enhanced, the integral evaluation precision is improved, and an effective basis is provided for the technical optimization and policy formulation of the carbon emission reduction effect.
Detailed Description
The technical scheme of the invention is described below with reference to the accompanying drawings.
In embodiments of the invention, words such as "exemplary," "such as" and the like are used to mean serving as an example, instance, or illustration. Any embodiment or design described herein as "exemplary" is not necessarily to be construed as preferred or advantageous over other embodiments or designs. Rather, the term use of an example is intended to present concepts in a concrete fashion. Furthermore, in embodiments of the present invention, the meaning of "and/or" may be that of both, or may be that of either, optionally one of both.
In the embodiments of the present invention, "image" and "picture" may be sometimes used in combination, and it should be noted that the meaning of the expression is consistent when the distinction is not emphasized. "of", "corresponding (corresponding, relevant)" and "corresponding (corresponding)" are sometimes used in combination, and it should be noted that the meaning of the expression is consistent when the distinction is not emphasized.
In embodiments of the present invention, sometimes a subscript such as W1 may be written in a non-subscript form such as W1, and the meaning of the expression is consistent when de-emphasizing the distinction.
In order to make the technical problems, technical solutions and advantages to be solved more apparent, the following detailed description will be given with reference to the accompanying drawings and specific embodiments.
Referring to fig. 1, a photosynthetic organism carbon fixation evaluation system under the condition of coking flue gas, the system comprises:
the carbon kinetic energy evolution monitoring module acquires injection data of coked smoke, the injection data comprise volume flow rate in unit time, inlet pressure and duration time of injection period, partial pressure change rate is called according to carbon dioxide migration rate in unit volume in a microalgae culture system in the gas injection process, ratio judgment is carried out on the two, a time node with the largest mutation amplitude in each gas injection period is extracted, and the node is acquiredJudging whether an adsorption upper limit feature exists according to the migration rate change direction and the change direction, if so, marking the adsorption upper limit feature as an inflection point, constructing a trend track, acquiring change comparison of the track under a plurality of flue gas load periods, and generating a carbon kinetic energy evolution track;
The method comprises the steps that an air pressure permeation threshold extraction module calls a carbon kinetic energy evolution track, an inlet static pressure change curve in an injection period is selected according to a pressure data segment corresponding to each inflection point, the fluctuation range of the curve in unit time before and after the inflection point is monitored, the air-liquid interface pressure which is continuously increased for the first time and causes migration rate change is extracted and recorded as a limit value, and a carbon migration trigger boundary group is generated by combining the occurrence frequency and duration of a plurality of period lower limit values;
the spectral response adsorption evaluation module calls a carbon migration trigger boundary group, and photon flux density and alga surface unit area under a plurality of wave bands are obtained in the period according to the corresponding gas injection periodExtracting adsorption change trend and judging reaction continuity under continuous wave bands under the condition of increasing adsorption rate, marking active adsorption sections in a curve, and evaluating by combining duration and response change performance to generate illumination response adsorption sections;
The flux adaptation change comparison module invokes an illumination response adsorption section, and performs trend comparison analysis according to a carbon fixation flux data curve of unit volume in a coverage time period and in combination with synchronous response in a carbon kinetic energy evolution track in coking flue gas injection response to generate carbon fixation adaptation change information;
The carbon fixation stable fluctuation indicating module calls carbon fixation adaptation change information, extracts amplitude values corresponding to five time points in each period according to the carbon fixation rate change amplitude of each stage, performs extremely poor extraction on the difference value, marks the most severe period of fluctuation, counts the fluctuation duration ratio in the period, and gathers the classification of the ratio section as a fluctuation indicating section to generate a carbon fixation stable fluctuation section;
The carbon kinetic energy evolution track comprises carbon migration critical time point distribution, gas injection period evolution trend curve andThe carbon migration triggering boundary group comprises an interface pressure limit value set, pressure fluctuation threshold value classification and gas injection period pressure activation mapping, the illumination response adsorption section comprises an active spectrum bandwidth section, an adsorption intensity distribution label and illumination section duration classification, the carbon fixation adaptation change information comprises a flux trend response offset set, an adsorption behavior response comparison group and an adaptation difference classification result, and the carbon fixation stable fluctuation section comprises a stability grade marking section, a fluctuation duration section and a period fluctuation structure identification label.
Referring to fig. 2 and 3, the carbon kinetic energy evolution monitoring module includes an injection parameter acquisition sub-module, a critical inflection point identification sub-module, and a trend track construction sub-module;
The method comprises the steps that an injection parameter acquisition submodule acquires injection data of coked flue gas, the injection data comprise volume flow rate in unit time, inlet pressure and duration of an injection period, carbon dioxide migration rate in unit volume in a corresponding time period in a microalgae culture system is acquired, meanwhile, partial pressure change rate of carbon dioxide in the gas injection period is monitored, a ratio sequence is established according to numerical ratio between the carbon migration rate in unit volume and the partial pressure change rate, and ratio sequence groups of all gas injection periods are screened to generate a ratio sequence integration value;
The method comprises the steps of acquiring volume flow, inlet pressure and cycle duration of coking flue gas in a gas injection cycle, firstly arranging a standard flow sensor and a pressure acquisition device at a flue gas discharge port of a coking section, measuring volume flow data in different gas injection time periods through a standard orifice plate flowmeter or an ultrasonic gas flowmeter, and recording the volume flow data as sampling frequency per minuteSynchronously collecting the instantaneous pressure of the inlet end pipe section at the corresponding time point through a pressure transmitterDuration of gas injection cycleThen the automatic statistics is obtained by setting the start-stop signal flag bit, for example, in a typical coking exhaust period, the gas injection starting time is 08:00:00, and the gas injection ending time is 08:30:00, then the period is thatIn seconds, the carbon dioxide migration rate in unit volume in the microalgae culture system is collected in the gas injection process, and the carbon dioxide migration rate in the photosynthetic bioreactor is required to be setMeasuring point of concentration difference, recording inlet concentration by laser infrared gas sensorAnd outlet concentrationBy unit volumeAccording to the formula:
;
deriving migration rateWherein all concentration units are converted into mol/m3, the volume flow is calculated in m3/s, and the partial pressure change rate is calculated by the inlet pressure value difference in each second and is recorded asThen constructing a ratio sequence at each momentRepeating the operation for all the gas injection cycles to generateThe sequence set is classified and screened according to the period, and the period in the classification is corresponding to each periodThe array is stored in a structured database to form multi-period ratio integration, and finally, a ratio sequence integration value is established;
The critical inflection point identification submodule extracts a time node with the largest mutation amplitude in each curve according to the ratio change curve in each period based on the ratio sequence integrated value, acquires the carbon migration rate change direction corresponding to the node, judges whether the direction is ascending and is in a local fluctuation high position, marks the ascending and high fluctuation direction as migration upper limit trigger points if the judgment is met, and integrates the time sequence distribution of the marked nodes to obtain a migration inflection point distribution section;
based on the ratio change curve in each period in the ratio sequence integrated value, the method needs to be carried out for eachAll time points in the sequence are processed by the first-order difference to obtain instantaneous increment valueSearching forIndex position with maximum absolute valueThe time corresponding to the index is the node with the maximum mutation amplitude, at the moment, trend judgment is carried out on the two ratios before and after the time point, if the ratio meets the requirementAnd is also provided withAnd the current valueAbove the average of all points in the first 30 secondsThe point is determined to be at a local fluctuation high level and marked as a migration upper limit trigger point, for example, if,,The point is marked, the operation is repeated on each periodic sequence, all marked nodes are integrated, the corresponding time index sets are extracted, a cross-periodic time sequence mapping table is constructed, and finally a migration inflection point distribution section is formed;
The trend track construction submodule calls a migration inflection point distribution section, performs difference operation on migration inflection point position time and gas injection starting time in the same period according to corresponding gas injection period information, calculates distribution trend of difference among periods, maps trend change directions to corresponding period groups according to load grade classification, establishes a change evolution path under injection conditions, and generates a carbon kinetic energy evolution track;
Invoking gas injection period data to which each time node belongs in the migration inflection point distribution section, and setting the moment of each migration inflection pointFor the initial time of the gas injection periodPerforming a difference operation, and recording asCounting the distribution of the difference under different periods, such as 90 seconds, 120 seconds and 150 seconds after gas injection at the inflection point in a certain period to obtain a sequenceComparing in different periodsMean of (2)And standard deviation ofCalculating the direction of the variation trend, if a certain period groupThe trend is continuously increased and the standard deviation interval is larger than 5 seconds, the trend direction is mapped to be 'fluctuation increment', and the load grades of each period are classified, for example, the volume flow is classifiedAnd (3) classifying to be high-load, making mapping marks by combining trend directions, finally constructing trend identifiers and track path sets of each period group, completing the summarization of evolution paths of inflection point distribution trends among all periods, and generating a carbon kinetic energy evolution track.
Referring to fig. 2 and 4, the barometric permeation threshold extraction module includes a pressure curve extraction sub-module, a fluctuation section identification sub-module, and a boundary limit classification sub-module;
The pressure curve extraction submodule calls a carbon kinetic energy evolution track, acquires an inlet static pressure change curve in an injection period where each time point is located according to a marked migration inflection point time point, intercepts a symmetrical time section taking the inflection point as a center on each curve, extracts pressure data points in the section, forms a continuous pressure change sequence in a time sequence arrangement manner, establishes a periodic static pressure fluctuation data set and generates a static pressure change sequence group;
When a migration inflection point time point of a mark in a carbon kinetic energy evolution track is called, each inflection point index in an established carbon kinetic energy evolution track database is required to be obtainedAnd extracting and processing original inlet static pressure data of the gas injection period of each inflection point, wherein the static pressure data is from a high-precision pressure sensor arranged at the front end of the gas injection pipe section, and the sampling frequency is set to be 1Hz. At each inflection point timeFor the center, extend for 30 seconds before and after the center, form a symmetrical time windowAll data points within the time window are extracted from the original pressure sequence and recorded asA total of 61 data points. Each pressure curve sequence is required to mark a period number, an inflection point number and a time window index in a data table, and a structural mapping relation is established, if the inflection point number is #3 and the period is the 5 th round of gas injection period, the corresponding entry is 'C5_M3'. And (3) arranging window pressure sequences corresponding to all inflection points in time sequence and storing the window pressure sequences in batches into a time sequence set to form a static pressure change sequence group, wherein if the static pressure of the period of the 3 rd inflection point is 121.3kPa, 122.0kPa and 123.1kPa to 125.0kPa respectively under a certain coking working condition, the complete sequence is used as a single record input static pressure fluctuation data set for subsequent discrimination.
The fluctuation section identification submodule judges whether the difference sign of adjacent pressure values in unit time is positive continuously according to each sequence in the static pressure change sequence group, calculates the time segment length which is positive continuously, screens the section sections which meet continuous rising for the first time and the duration time is larger than the threshold value of the set fluctuation time, extracts the pressure value corresponding to the starting moment of the section, integrates the section starting points which meet the condition to form an index set, and obtains a continuous rising section starting pressure set;
According to each time sequence in the static pressure change sequence groupJudging difference sign of adjacent pressure points with unit time of 1 secondIf continuously presentRegarding as continuous rising trend of pressure, recording start point, end point and duration of each continuous positive difference section during traversal, setting fluctuation time threshold to 5 seconds, i.e. at least 5 continuous time points with difference greater than 0, for example, satisfying within 10-16 seconds in a sequenceThe segment is marked as a candidate segment, the initial time point is 10 seconds, and the initial pressure value is. Repeating the judgment on all inflection point window sequences, screening out fragments meeting the condition for the first time, and extracting the corresponding initial pressure valueIf a sequence continuously rises at t=14s and continues until t=21s and the pressure increases from 122.3kPa to 124.1kPa, the segment is marked and recorded. And uniformly merging all initial pressure values meeting the continuous positive increasing condition to form a continuous rising section initial pressure set, recording inflection points and period indexes of the initial pressure set, and establishing a one-to-one corresponding form structure for classifying operation in the next stage.
The boundary limit value classifying submodule calls a continuous ascending section starting pressure set, counts the frequency of occurrence in a plurality of injection periods, records the duration range of each pressure value, constructs a combined characteristic sequence according to the frequency and the duration, classifies the combined characteristic sequence according to the stability sequence of the combined characteristic sequence, and then classifies the combined characteristic sequence into a boundary condition set to generate a carbon migration triggering boundary group;
Calling each pressure value in the starting pressure set of the continuous rising sectionCounting the frequency of each value occurring in multiple periodsAnd calculate the corresponding paragraph durationConstructing combination feature pairs by termsThen, stability sorting treatment is carried out on all feature pairs, and stability weights are setToSorting and classifying the sizes in order, and thenThose values above the set stability reference value (e.g., set to 100) are classified as "stability boundary segments", those below which are classified as "secondary boundary segments", if a starting pressure value of 122.8kPa occurs in 6 cycles and the average duration is 20 secondsMeets the standard of the stabilized segment. And (3) sorting boundary condition sets according to the sorting result, merging to form a carbon migration triggering boundary group, wherein each record in the content structure is attached with four fields of 'initial pressure value', 'frequency', 'duration' and 'classification label', and the four fields are used as important basic parameters for judging the active window of the carbon response in the coking flue gas period of the system.
Referring to fig. 2 and 5, the spectral response adsorption evaluation module includes a flux density acquisition sub-module, an adsorption trend identification sub-module, and a response section evaluation sub-module;
the flux density acquisition submodule acquires a carbon migration trigger boundary group, and acquires photon flux densities and surface unit areas of algae surfaces in a plurality of wave bands in each period according to the covered gas injection periodThe adsorption rate increases the data, establishes a photon flux density and adsorption rate corresponding relation matrix, extracts a change curve of the corresponding relation in a spectrum band section, and generates a spectrum response curve group;
After the gas injection period covered in the carbon migration triggering boundary group is acquired, the illumination control parameters of the corresponding period and the sensor data in the bioreactor are required to be extracted according to the period number. In the illumination system, each gas injection period is provided with an adjustable multi-band light source, the band range is set to 400nm to 700nm, each 20nm is a segment, and 16 bands are divided. Measuring photon flux density per unit area using an radiometer in each bandThe unit is converted into a relative value through normalization processing so as to eliminate the influence of energy input of different wave bands. Synchronously, the gas exchange layer sensor array embedded on the surface of the reactor is utilized to acquire the unit timeAdsorption rateThe algae volume area is treated as a normalization factor to obtainAnd record the corresponding adsorption durationIn seconds, the time unit, normalized for normalization comparison. Will be in each bandThe three groups of data are correspondingly stored in a three-dimensional matrix, and the dimensions of the matrix respectively correspond to the band sequence number, the flux value and the adsorption response parameter. Under practical experimental conditions, e.g. at band 7 (520 nm)、、Corresponding points are recorded in a matrix, 16 change curves are generated according to the sequence numbers of the wave bands, and the transverse axes of the curves are the wavelengthsThe vertical axis is the adsorption rate of a unit area, and finally a spectrum response curve group is generated, so that a basic data set is provided for the identification of the subsequent adsorption reaction trend;
The adsorption trend identification submodule judges whether the directions of the adjacent wave bands corresponding to the adsorption rate change are continuously consistent or not according to the spectrum response curve group, extracts the wave band sections meeting continuous rising, calculates adsorption strength response variation values in the continuous wave band sections, screens and marks wave bands with variation values larger than continuous response reference strength difference, and generates an adsorption active section set;
the formula for calculating the adsorption strength response variation value in the continuous wave band section is specifically as follows:
;
Wherein, theRepresents the variation value of the adsorption response,Represents the firstNormalized value of CO2 adsorption rate variation corresponding to each wave band,Represents the firstNormalized value of adsorption duration per unit area in each band,Represents the firstThe center wavelength of the individual wavelength bands,Represents the firstNormalized values of photon flux density corresponding to the individual bands,Is the total number of the extracted wave bands;
According to each curve data in the spectrum response curve group, firstly judging the change direction of the adsorption rate between adjacent wave bands, namely calculating the normalized adsorption rate difference value under the adjacent wave bandsIf the difference is positive, the adsorption response of the current section is determined to be in a continuous increasing trend. For example, in three continuous wave bands of 520nm, 540nm and 560nm, the normalized adsorption rates are respectively、、Adsorption duration per unit area is respectively、、The central wavelength of the wave band is respectively、、(In nm) corresponding normalized photon flux densities of、、. Substituting the parameters into the adsorption response variation value formula:
;
first, the sum of the adsorption rate and duration product is calculated:
;
and calculating the average value of the products of the wavelength and the flux:
;
substituting the formula to obtain:
;
if the continuous response reference intensity difference is set to 10.0, the band segment calculation resultGreater than the baseline value, the adsorption response activity determination criterion is satisfied, and therefore the band segment will be screened into the adsorption activity segment set. After screening according to the same flow in all spectral response curves, all the output satisfying the formula is reservedContinuous wave band segment greater than reference value, recording its start-stop wave band number, its period and its correspondent positionAnd numerical values and the like, which are marked as adsorption active section sets and serve as basic data sets for subsequent section reaction intensity assessment and map generation. The structure of the set can be set as a band section ID, a start wavelength, a stop wavelength, an air injection period ID and an adsorption response intensity value, so that graphical display and data backtracking processing are facilitated.
The response section assessment submodule calls an adsorption active section set, extracts photon flux density variation quantity of a section corresponding to a time segment and adsorption quantity variation trend in unit time according to the marked wave section, constructs an active wave section adsorption curve set, judges adsorption reaction expression level according to fluctuation amplitude and adsorption increment of the curve, establishes an active response wave section time coverage map after grouping and classifying, and obtains an illumination response adsorption section;
And (3) calling each section index in the adsorption active section set, extracting a photon flux density actual measured value sequence and a unit time adsorption quantity sequence in a corresponding time window, constructing an adsorption trend curve set in the section, wherein the abscissa of the curve is time, the ordinate is adsorption rate, and then classifying and dividing according to the maximum fluctuation amplitude and the adsorption total increment of each curve. The fluctuation range grade limit is set to be 0.15, the adsorption increment grade limit is set to be 0.4, if the fluctuation of a certain curve is 0.18 and the increment is 0.45, the grade is classified as a high reaction grade, and if the fluctuation is 0.09 and the increment is 0.22, the grade is classified as a low reaction grade, and the corresponding wave band sections are marked respectively. Mapping all the classified active sections onto corresponding gas injection period time shafts according to time start-stop positions, integrating all the marks to generate an active response wave band time coverage map, recording wave band time span and response intensity by the map in a Gantt mode, and finally obtaining the illumination response adsorption section.
Referring to fig. 2 and 6, the flux adaptation change comparison module includes a carbon sequestration data extraction sub-module, a trend response comparison sub-module, and an adaptation performance merging sub-module;
The carbon fixation data extraction submodule is used for collecting continuous data curves of carbon fixation flux of unit volume in a time period according to the marked time period based on the illumination response adsorption section, extracting periodic carbon fixation data segments based on a time axis, constructing a flux change sequence, classifying and integrating each segment of data by combining with time stamp information, and establishing a flux behavior expression set to generate a carbon fixation flux data set;
after the marked time period in the illumination response adsorption section is obtained, the corresponding time interval in each gas injection period is subjected to partition identification, and the carbon fixation flux of unit volume recorded in the microalgae reactor at intervals of seconds is readThe flux is calculated by the CO2 concentration difference in combination with the reactor gas exchange capacity, typically in mol/m3.multidot.s. In practical implementation, the acquisition system willConcentration at inlet and outlet、Flux calculation uses the formulaWhereinFor the gas injection flow rate,For the effective volume of the culture system, if,,,Then. All of each gas injection periodTime series data are ordered according to acquisition time and are archived in cooperation with system time stamps to form time-flux pairs, e.g. Then slicing the data according to the coverage time of the illumination response adsorption section, classifying according to the cycle number, and extracting all the data under the section if the illumination response section in the cycle C1 is t=300 s to t=600 sAnd are classified into flux sequencesRepeating the operation to form all periodic data sets, establishing a structured carbon sequestration behavior database, and summarizing to obtain a carbon sequestration flux data set;
The trend response comparison submodule calls a carbon fixation flux data set, and according to flux curves of each time period, synchronous response information marked in a carbon kinetic energy evolution track is compared, synchronous flux change rate and kinetic energy fluctuation change rate sequences are extracted, change trend difference value calculation is carried out after the two sequences are matched according to time indexes, and flux response offset values are calculated and obtained;
the formula for calculating the acquired flux response offset value is specifically as follows:
;
Wherein, theIndicating the value of the flux response offset,Represent the firstThe carbon fixation flux per unit volume at the moment,Represents the maximum value of the carbon fixation flux in the same period,Represent the firstThe kinetic energy value of carbon per unit volume at the moment,Represents the maximum value of the kinetic energy of carbon in the same period,Representing the total number of time points involved in the matching;
invoking each cycle in the fixed carbon flux data setFlux sequence and invoking kinetic energy value sequence of unit time corresponding to synchronous gas injection period in carbon kinetic energy evolution trackThe two sequences need to be aligned according to the time stamp to ensure eachCorresponding to the moment of timeMatching. And then respectively normalizing the flux value and the kinetic energy value at each moment j, wherein the normalization rule is as follows:、 wherein、. Let m be the total number of time points of the sequence, in each gas injection period, the following formula is executed:
;
Setting 5 time points in a certain period, wherein the time points are respectively as follows:
maximum value of;
Maximum value of;
Then after normalization:
;
;
the difference sequence is: averaging absolute value:
;
recording the calculation result into the period index, and forming a flux response offset value set for all periods according to the calculation result, wherein the flux response offset value set is used as basic evaluation data of flux-kinetic energy matching relations in different load stages;
the adaptation performance merging sub-module extracts corresponding fluctuation interval data according to the flux response offset value, screens continuous sections with the variation amplitude smaller than the response tolerance threshold, counts the duration of the screened sections, extracts the response stability duration, establishes an adaptation grade label in combination with the response matching strength, and constructs a flux response behavior grading set to obtain carbon fixation adaptation variation information;
From recordings in flux response offset valuesExtracting the corresponding fluctuation interval of each period, namely, within a continuous time periodThe instantaneous difference of (2) does not exceed the given response tolerance threshold, the threshold is set to 0.05, in the previous exampleThe time period for continuously satisfying the condition is t2 to t5, the corresponding duration is 4s, and the effective response period is recorded. After screening and duration extraction are carried out on all sections meeting the condition in each period, a response stable duration sequence is constructed, a response grade standard is defined according to the response stable duration sequence, and a grade threshold is set to be low adaptation <3s, medium adaptation 3-6s and high adaptation >6s, and the medium adaptation grade is adopted in the previous example. And adding an adaptation grade label to each section, then summarizing the adaptation grade labels into a response expression list, forming a flux response behavior grading set according to grade ordering, and outputting carbon fixation adaptation change information.
Referring to fig. 2 and 7, the carbon fixation stable fluctuation indicating module includes an amplitude difference extracting sub-module, a fluctuation intensity identifying sub-module, and a period section collecting sub-module;
The amplitude difference extraction submodule acquires carbon sequestration adaptation change information, extracts speed amplitude values corresponding to five fixed time points in each period according to the carbon sequestration speed change amplitude data in each period, calculates a difference value between every two time points, extracts a maximum value and a minimum value, performs a difference operation between the difference values, summarizes the difference output values of the periods, and generates a carbon sequestration fluctuation difference set;
When acquiring the rate change amplitude data of each period in the fixed carbon adaptation change information, firstly, confirming the selection rules of five fixed time points in each period, generally taking equidistant time points in each period as references, for example, the total period length is 300 seconds, then the five time points are respectively set as 60s, 120s, 180s, 240s and 300s, and recorded asTo the point of. Reading a carbon sequestration rate value per unit volume at a corresponding time point from a carbon sequestration behavior databaseTo the point ofFor example, the data collected during a certain period isSequentially calculating a sequence of differences between adjacent time pointsExtracting the maximum difference from the aboveMinimum differencePerforming polar difference operation, namely, polar difference = maximum difference-minimum difference = 0.003- (-0.002) =0.005, recording the polar difference output value of the current period as 0.005, processing the 5-point flux data of all the periods according to the same method, summarizing the polar difference output values under each period to form a standard data set, and establishing a carbon fixation fluctuation polar difference set;
The fluctuation severe recognition submodule calls a carbon fixation fluctuation extremely difference set, screens periods larger than fluctuation recognition standard extremely difference according to each period extremely difference value, marks the periods as high variation period sections, extracts start and stop time and carbon fixation response curve segments of the marked periods, constructs a high fluctuation curve distribution index set and obtains a high variation period index set;
And calling each item of data in the fixed carbon fluctuation range, setting the fluctuation identification standard range to be 0.004 mol/m & lt 3 & gt.s, screening the range values of all the periods according to the standard, classifying the periods into high variation periods if the range value of a certain period is more than 0.004, and conforming to screening conditions if the range value of the period is 0.005 in the previous example. And extracting the cycle numbers, the starting and ending time points and the corresponding carbon fixation response curves of all the high variation cycles to form a high fluctuation identification index. Taking the period C3 as an example, the fluctuation time ranges from t=60 s to t=300 s, and the response curve segments are markedAnd the curve segment is encoded into C3_W1, all periodic curve index number sets meeting the conditions are summarized, and a high-variation periodic index group is established, and the index group is used for periodic classification and stability zone drawing in the follow-up.
The period section collecting sub-module is used for counting the proportion of the fluctuation duration of the carbon fixation rate in each period to the total operation period according to the corresponding time range of each period, classifying the intervals between the duty ratio data and a set fluctuation duration ratio threshold value, integrating the classified time periods to construct a multi-level fluctuation time band spectrum, and establishing a carbon fixation stable fluctuation section;
based on the time ranges recorded for each period in the high variation period index set, the duration of the fixed carbon rate fluctuation (i.e. in a higher rate segment than the period average fluctuation) per period is countedAnd with the total period of timePerforming ratio operation to obtain fluctuation continuous duty ratio. For example, if the total duration of the C3 period is 300s and the fluctuation period lasts 75s, then. Setting the fluctuation continuous ratio threshold value to be 0.2 and 0.4 as grade boundaries, ifThen labeled "low fluctuation zone", ifThen "medium fluctuation zone",Is a "high fluctuation section". And collecting the fluctuation level and the time period of each period, and establishing a multi-level fluctuation time distribution map covering all the gas injection periods to finally form a carbon fixation stable fluctuation section for representing the stability expression period of the system under different working conditions.
It should be understood that the term "and/or" is merely an association relationship describing the associated object, and means that three relationships may exist, for example, a and/or B, and may mean that a exists alone, while a and B exist alone, and B exists alone, wherein a and B may be singular or plural. In addition, the character "/" herein generally indicates that the associated object is an "or" relationship, but may also indicate an "and/or" relationship, and may be understood by referring to the context.
In the present invention, "at least one" means one or more, and "a plurality" means two or more. "at least one of" or the like means any combination of these items, including any combination of single item(s) or plural items(s). For example, at least one (a, b, or c) of a, b, c, a-b, a-c, b-c, or a-b-c may be represented, wherein a, b, c may be single or plural.
It should be understood that, in various embodiments of the present invention, the sequence numbers of the foregoing processes do not mean the order of execution, and the order of execution of the processes should be determined by the functions and internal logic thereof, and should not constitute any limitation on the implementation process of the embodiments of the present invention.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
It will be clearly understood by those skilled in the art that, for convenience and brevity of description, specific working procedures of the apparatus, device and unit described above may refer to corresponding procedures in the foregoing method embodiments, which are not repeated herein.
In the several embodiments provided by the present invention, it should be understood that the disclosed apparatus, device and method may be implemented in other manners. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of the units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another device, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention may be embodied essentially or in a part contributing to the prior art or in a part of the technical solution, in the form of a software product stored in a storage medium, comprising several instructions for causing a computer device (which may be a personal computer, a server, a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. The storage medium includes a U disk, a removable hard disk, a read-only memory (ROM), a random access memory (random access memory, RAM), a magnetic disk, an optical disk, or other various media capable of storing program codes.
The foregoing is merely illustrative of the present invention, and the present invention is not limited thereto, and any person skilled in the art will readily recognize that variations or substitutions are within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.