Application of Moving Bed Biofilm Reactor and Fixed Bed Hybrid Biological Reactor for Oilfield Produced Water Treatment: Influence of Total Dissolved Solids Concentration

: This experimental paper deals with the development of a hybrid biological reactor for the treatment of a synthetic oilﬁeld produced water under an increase in total dissolved solids (TDS) concentration. To comply with strengthening regulations concerning produced water discharge and peculiar produced water compositions, a moving bed bioﬁlm reactor (MBBR) consisting in a combination of free activated sludge and moving bioﬁlm supports was compared to a ﬁxed bed hybrid biological reactor (FBHBR) consisting in a combination of free activated sludge and a ﬁxed bioﬁlm support. After a 216 days experimental period, the MBBR and the FBHBR were efﬁcient to treat a synthetic produced water with chemical oxygen demand (COD) removal rate above 90% under an increase in TDS concentrations from 1.5 to 20 g · L − 1 . Ecotoxicity measurements on freshwater and marine microorganisms revealed an absence of toxicity on treated waters. A decrease in bacterial diversity indices with respect to the inoculum was observed in both bioreactors. This suggests that the increase in TDS concentrations caused the predominance of a low number of bacterial species.


Introduction
Oil and gas activities generate a large amount of wastewater during crude oil and gas extraction.This wastewater, also called produced water (PW) is composed of formation water, initially present in the oil reservoir.Other sources of water may also be added in the reservoir.The resulting outcoming water is then produced water.The ageing of oil reservoirs is currently leading to an increase in PW flowrates [1].The water to oil ratio (WOR) quantifies this increase.Currently, the WOR is around 3:1, but its value is expected to grow [2].
Energies 2021, 14, 7297 2 of 20 Currently, the difficult access to water resources in water-stressed area has made oil producers to rethink their management of PW.Whether PW is discharged in the environment or reinjected in oil reservoirs, its quality should meet specific standards [8,9].Concerning PW discharge, regulations are becoming more stringent.It is also worth noting that these regulations are mainly site-specific.In the North Sea, the Convention for the Protection of Marine Environments in the North-East Atlantic (OSPAR) has set the average maximum discharge value for dispersed oil at 30 mg•L −1 [10].In the United States, the Environmental Protection Agency has set an oil and grease discharge value at 29 mg•L −1 and a daily maximum discharge limit value at 42 mg•L −1 .In a recent review, Zheng et al. (2016) summarized existing regulations of PW discharge all around the globe [11].Recently, regulations started to focus on dissolved hydrocarbons, which cannot be remove from PW with actual treatment technologies.For instance, OSPAR recommendation 2014/5 provides a list of Predicted No effect Concentrations (PNECs) for some compounds found in PW [12].This should force oil producers to adapt and optimize treatment facilities for a better quality of discharged PW.
PW are currently treated by mean of physico-chemical processes for the removal of dispersed hydrocarbons [13].However, these processes were not designed to remove dissolved hydrocarbons.Recently, research have been focused on biological treatments [14][15][16][17][18][19][20].Those treatments are considered as a good alternative to physico-chemical processes due to low operating costs.Furthermore, they are considered eco-friendly [21,22].In these processes, the action of a bacterial consortium in an aerated tank is responsible of the wastewater purification, by taking required nutrients and food for the biomass growth from the wastewater.However, few research reported the application of biological treatment for PW treatment.So far, conventional activated sludge (CAS) reactor, biological aerated filters (BAF), and membrane bioreactors (MBR) were proved to be efficient to treat PW [21][22][23][24][25].However, these processes were applied under operational conditions not representative of real field conditions.For instance, hydraulic retention times (HRT), i.e., the time that the water is in contact with the purifying biomass, are usually high in previous studies.This is not compatible with field applications where (i) PW flowrates are big and (ii) space allocated to treatments are limited especially on offshore fields.
Another challenge in PW treatment is the high salinity levels in these wastewaters.Measured with total dissolved solids concentration (TDS), this parameter can reach 300 g•L −1 in PW [26].The adverse effect of salinity on wastewater biological treatment is well documented in the literature [27,28].High salt content is responsible for cell membranes disruptions, enzymes denaturation and strong osmotic pressure that are lethal to conventional purifying microorganisms.This sensitivity to ionic strength and salt shock loads is responsible for decreases in removal efficiencies and volatile suspended solids.Furthermore, it also causes an increase in effluent suspended solids.One of the major issue in saline conditions is settling issue [29].Lefebvre and Moletta (2006) reported several effects of salinity on the property of activated sludge.The density of saline water being higher than fresh water, greater resistance to decantation occurs due to high buoyant forces [30].The lack of versatility of acclimated microorganisms is also pointed out.Once acclimated to saline conditions, it has been noted that bacteria easily lose their efficiency when salinity decreases.Working under hypersaline conditions either implies that the biomass used in the aerated tank has to undergo a well conducted acclimation time with a gradual increase in salinity or the use of halophilic consortia [19][20][21]26,30].
In the literature authors note obviously that the increase in salinity is associated with the decrease in COD removal.Indeed, the increase in salinity is responsible for biological treatment disturbance, as described before.In more details Sharghi et al. (2013) studied the effect of salinity on a sMBR under a HRT of 48 h [19].Their relatively low OLR (0.3-0.9 kg COD •m −3 •d −1 ) conjugated with the high HRT and the use of halophilic microorganisms led to high COD removal efficiencies (83%).It points out the feasibility of PW treatment at high salinity with biological processes.However, all studies reporting high COD removal efficiencies at high salinity were performed with synthetic produced water, suggesting that the experimental conditions were still soft enough to allow a good bacterial activity.In this context the study conducted by Pendashteh et al. (2012) revealed that the COD removal efficiency of a eMBR could decrease down to 17 % when the experiments were carried out with real produced water at high salinity (TDS = 250 g•L −1 ) [21].This result obtained under a HRT of 48 h, and with the use of halophilic microorganisms showed that there are still investigations needed to understand to what extent salinity can affect the biological treatment of PW.In their work, Pendashteh et al. (2012) only suggested that salinity could affect the biodegradation rate of the dissolved organic compounds.
To answer the challenge brought by the high TDS concentrations, attention should be paid to hybrid processes.These processes contain two or more treatments processes in one reactor.Synergistic effects of the combination of the two treatments may enhance the removal of recalcitrant compounds [31].In this context, the moving bed biofilm reactor (MBBR) and the fixed bed hybrid biological reactor (FBHBR) are promising to treat PW.Both technologies combine free biomass (activated sludge) and fixed biomass (biofilm) in the same reactor.In the MBBR, the packing material where the biofilm grows is allowed to move freely inside the bioreactor whereas in the FBHBR, the packing material is prevented from moving.These bioreactors have interesting advantages over conventional technologies.With the simultaneous presence of fixed and free biomass, higher biomass concentrations can be reached with respect to conventional technologies.Hybrid biological reactors can also be implemented with small footprints [32,33].Thus, those processes are of particular interest for offshore implementations.The presence of the biofilm allows the development of slow growing bacteria, which is beneficial for the removal of poorly biodegradable compounds [34].This type of reactor was proved to be mass transfer efficient and suitable for the treatment of domestic wastewater [35][36][37].Dong et al. (2011) were the first to study the treatment of PW in a hybrid biological reactor [38].Their experiments carried out in a MBBR filled with modified ceramic biocarriers, in order to increase the specific surface area, answered some questions about the development of hybrid biological reactors.Their comparison between a CAS reactor, a MBBR with ceramic biocarriers and a MBBR with modified ceramic biocarriers proved that the hybrid biological processes were better in terms of COD removal.When comparing the two MBBR, the modified ceramic biocarriers allowed a greater extent of biofilm colonization, which maintained high COD removal efficiencies when the HRT was reduced.This suggests a beneficial effect of the specific surface area in hybrid biological reactor.More recently, Hasanzadeh et al. (2020) used walnut shells as biocarriers in a MBBR treating PW [39].Interestingly, they used halophilic microorganisms to conduct the treatment under high salinity.In their experiments, COD removals were above 70% even under a salinity of 250 g•L −1 .This was attributed to a beneficial effect of the biocarriers.Authors pointed out that the accumulation of pollutants on the biocarriers allowed a gradual digestion by the biofilm under higher retention time, which maintained a good bacterial activity.
The aim of this experimental study is to demonstrate the efficiency of two hybrid biological reactors (namely a MBBR and a FBHBR) for the treatment of a synthetic PW.Attention was paid to the salinity by gradually increasing TDS concentrations from 1.5 to 20 g•L −1 .

Experimental Setup
The laboratory scale experimental device used for this study was described in a previous study [17].The experimental pilot is made of two identical cylindrical clear PVC columns as shown in Figure 1a.Both columns provide 30 L of working volume with an inner diameter of 190 mm and a total height of 1070 mm.As presented in Figure 1b, the left column was configured as a MBBR.The right column was configured as a FBHBR.The packing material in both the MBBR and the FBHBR was made of cylindrical solid packing rings (AnoxKaldnes ® , Lund, Sweden) as shown in Figure 2. The filling ratio was 33% v/v based on the results of previous studies [28,29].The packing rings were inserted in a closed 7-mm stainless steel grid basket.The activated sludge was separated from the treated water in a 5 L gravity settler.A peristaltic pump (Masterflex, Gelsenkirchen, Germany) was used to recirculate the settled sludge into the bioreactors.The recirculation rate was set at 200% of the feed rate to prevent sludge from accumulating in the settlers.The synthetic PW was also fed into the bioreactor by a peristaltic pump (Masterflex).Aeration was provided with compressed air and injected to the system through a porous membrane at the bottom of the column (Jäger).Flowrates were measured using a rotameter (Analyt, Müllheim, Germany).Air input ranged between 120 and 960 L•h −1 .
columns as shown in Figure 1a.Both columns provide 30 L of working volume with an inner diameter of 190 mm and a total height of 1070 mm.As presented in Figure 1b, the left column was configured as a MBBR.The right column was configured as a FBHBR.The packing material in both the MBBR and the FBHBR was made of cylindrical solid packing rings (AnoxKaldnes ® , Lund, Sweden) as shown in Figure 2. The filling ratio was 33% v/v based on the results of previous studies [28,29].The packing rings were inserted in a closed 7-mm stainless steel grid basket.The activated sludge was separated from the treated water in a 5 L gravity settler.A peristaltic pump (Masterflex, Gelsenkirchen, Germany) was used to recirculate the settled sludge into the bioreactors.The recirculation rate was set at 200% of the feed rate to prevent sludge from accumulating in the settlers.The synthetic PW was also fed into the bioreactor by a peristaltic pump (Masterflex).Aeration was provided with compressed air and injected to the system through a porous membrane at the bottom of the column (Jäger).Flowrates were measured using a rotameter (Analyt, Müllheim, Germany).Air input ranged between 120 and 960 L•h −1 .

Microbial Inoculum
The inoculum was withdrawn from a municipal wastewater treatment plant, in the activated sludge recycling line (Aix en Provence, France).Initially, the amount of collected activated sludge was set into two identical plastic batch reactors under gentle stirring.20 L of Clean AnoxKaldnes rings were inserted in the reactors to develop a biofilm for 30 columns as shown in Figure 1a.Both columns provide 30 L of working volume with an inner diameter of 190 mm and a total height of 1070 mm.As presented in Figure 1b, the left column was configured as a MBBR.The right column was configured as a FBHBR.The packing material in both the MBBR and the FBHBR was made of cylindrical solid packing rings (AnoxKaldnes ® , Lund, Sweden) as shown in Figure 2. The filling ratio was 33% v/v based on the results of previous studies [28,29].The packing rings were inserted in a closed 7-mm stainless steel grid basket.The activated sludge was separated from the treated water in a 5 L gravity settler.A peristaltic pump (Masterflex, Gelsenkirchen, Germany) was used to recirculate the settled sludge into the bioreactors.The recirculation rate was set at 200% of the feed rate to prevent sludge from accumulating in the settlers.The synthetic PW was also fed into the bioreactor by a peristaltic pump (Masterflex).Aeration was provided with compressed air and injected to the system through a porous membrane at the bottom of the column (Jäger).Flowrates were measured using a rotameter (Analyt, Müllheim, Germany).Air input ranged between 120 and 960 L•h −1 .

Microbial Inoculum
The inoculum was withdrawn from a municipal wastewater treatment plant, in the activated sludge recycling line (Aix en Provence, France).Initially, the amount of collected activated sludge was set into two identical plastic batch reactors under gentle stirring.20 L of Clean AnoxKaldnes rings were inserted in the reactors to develop a biofilm for 30

Microbial Inoculum
The inoculum was withdrawn from a municipal wastewater treatment plant, in the activated sludge recycling line (Aix en Provence, France).Initially, the amount of collected activated sludge was set into two identical plastic batch reactors under gentle stirring.20 L of Clean AnoxKaldnes rings were inserted in the reactors to develop a biofilm for 30 days.During this biofilm growth period, the activated sludge was fed with sodium acetate, ammonium chloride and potassium dihydrogen phosphate while respecting a C/N/P ratio of 100/5/1.The remaining sludge and colonized rings were equally put in the continuous bioreactor to configure a MBBR and a FBHBR.During the whole experiment, the sludge retention time was 20 days.

Composition of the Synthetic Produced Water
A synthetic PW was used to perform the experiments.However, the concentration of targeted pollutants may be low to insure the viability of the biomass.Therefore, dissolved nutrients were added to the PW to set values of Chemical Oxygen Demand (COD), Total Nitrogen (TN) and Total Phosphorous (TP) at 1600, 20.6 and 4 mg•L −1 , respectively.A detailed composition of the synthetic PW is summarized in Table 1 The salinity (expressed as TDS concentration) was gradually increased throughout the experiment.TDS concentrations were set at 1.5 g•L −1 from day 1 to day 100, 8 g•L −1 from day 101 to day 180 and 20 g•L −1 from day 181 to day 216.

Experimental Schedule
The experiment was divided into four steps.Each step consisted in an increase in the TDS concentration.Three TDS concentrations were successively studied: 1.5, 8 and 20 g•L −1 .The concentrations of pollutants were constant during the whole experimental period.The HRT was set at the constant value of 12 h.The OLR based on the working volume (30 L in both MBBR and FBHBR) was 3.2 kg COD •m −3 •d −1 .All these key parameters are summarized in Table 2.

Solids Concentration Measurements
Standards methods were applied to measure both free TSS and VSS concentrations [40].In both the MBBR and the FBHBR, the quantity of biofilm was assessed at each phase of the experimental period.The biofilm weight was determined using an experimental procedure described by Abtahi et al. [41].

Physical-Chemical Parameters
Dissolved oxygen (DO) concentrations, pH and temperature using oxygen electrodes (LDO10101, Hach, Loveland, CO, USA) and pH electrodes (PHC101, Hach, Loveland, CO, USA).It is also worth noting that all water samples were filtered using 0.45 µm polyether sulfone filters before analysis.COD and phenol were analyzed using Hach DR6000 analytical kits and spectrophotometric measurements.BTEX (toluene, m-xylene and oxylene) were measured by mean of a standardized chromatographic method (Headspace GC/MS) [42].PAHs (naphthalene, phenanthrene, benzo[a]pyrene) were measured by mean of gas chromatography coupled with mass spectrometry [43].

Evaluation of Water Ecotoxicity
Ecotoxicity tests were performed to evaluate the potential toxicity of both MBBR and FBHBR outlet waters.Samples were taken at the end of each step of the whole experimental period.The feed inlet and both MBBR and FBHBR outlet were tested.Standardized test used were the Microtox (Vibrio fischeri) acute toxicity test [44], the freshwater algae (Pseudokirchneriella subcapitata) acute toxicity test [45], the Daphnia magna acute toxicity test [46], and the Brachionus calyciflorus chronic toxicity test [47].When the increase in salinity was not compatible with the freshwater microorganisms, marine ecological tests were chosen to measure the ecotoxicity of the treated water.These were the growth inhibition of a marine algae (Phaeodactylum tricornutum) test [48] and the acute lethal toxicity of marine crustacea (Artemia salina) test [49].

Characterization of the Microbial Population
The evolution of the microbial population in both the CAS reactor and FBHBR was investigated between the different key stages of the experiment: at the beginning of the acclimation and after the end of each step.The extraction and characterization procedure was described in a previous study [17].

Evolution of Solids Concentrations
The evolution of free TSS in both reactors are represented in Figure 3b.During the first step (day 1-100) of the bioreactors operation (TDS = 1.5 g•L −1 ), a large decrease is observed from days 1 to 30 (from 10.23 g•L −1 to 1.18 g•L −1 in the FBHBR and from 7.5 g•L −1 to 0.8 g•L −1 in the MBBR.This decrease, due to the acclimation of the biomass, is the result of the toxicity of the PW compounds which is lethal for some bacterial species composing the activated sludge that we used to seed our bioreactors at t = 0.For instance, the analysis of the microbial population revealed differences between the inoculum and the population after the acclimation (disappearance of Planctomyces sp., Nitrospira sp., Dokdonella sp., Diaphorobacter sp. and a large part of unclassified species whereas Flavobacterium sp., Hydrogenophaga sp., and Dechloromonas sp.grew in all biomasses).Those observations will be further discussed in Section 3.4.
At the end of this first acclimation step, free TSS concentrations were 3.75 g•L −1 in the FBHBR and 2.94 g•L −1 in the MBBR.The second step started from day 101 to day 180, with a gradual increase of TDS concentration from 1.5 g•L −1 to 8 g•L −1 (day 101-145).At the end of the second step (TDS = 8 g•L −1 ), the free TSS concentrations were 2.39 g•L −1 and 2.19 g•L −1 in the FBHBR and the MBBR, respectively.This represents a decrease by 56% in the FBHBR and by 34% in the MBBR.An increase in TSS concentrations by 123% in the FBHBR and by 140% in the MBBR was observed during the third step, from day 181 to day 216 (TDS = 20 g•L −1 ).At the end of the process operation, free TSS concentrations were 5.33 g•L −1 and 5.26 g•L −1 in the FBHBR and the MBBR, respectively.The evolution of the free TSS concentrations revealed that there seems to be no significant differences between the FBHBR and the MBBR based on this parameter.
The evolution of free VSS concentrations in both reactors are represented in Figure 3c.The free VSS concentrations followed the same trend as the free TSS concentrations in both reactors.At the end of step I, VSS concentrations were 3.23 g•L −1 in the FBHBR and 2.56 g•L −1 in the MBBR.At the end of step II, free VSS concentrations decreased to 0.69 g•L −1 in the FBHBR and 1.16 g•L −1 in the MBBR.This represents a decrease by 79% in the FBHBR and by 55% in the MBBR.At the end of the bioreactors operation, free VSS concentrations increased by 333% to 2.99 g•L −1 in the FBHBR and by 237% to 3.91 g•L −1 in the MBBR.From the beginning of the process operation to day 127, there were no significant differences between the FBHBR and the MBBR concerning the free VSS concentrations.
When the TDS was increased to 5 g•L −1 , free VSS concentrations appeared to be bigger in the MBBR than in the FBHBR.The presence of the packing material in both MBBR and FBHBR allowed the to develop.In the MBBR, the biofilm grew on the internal surface of the Kaldnes r the FBHBR, the biofilm grew in the internal and external surface of the Kaldne There is no shear stress imposed to the external ring surface in the FBHBR, which the larger development of the biofilm.To better understand the evolution of solids concentrations, values of free VSS/TSS ratios are plotted against time in Figure 3d.During the first step (TDS = 1.5 g•L −1 ), free VSS/TSS ratios were high, reaching an average value of 0.8 in the FBHBR and 0.8 in Energies 2021, 14, 7297 8 of 20 the MBBR.These values are typical of biomasses grown in synthetic PWs and reflect the absence of accumulation of inert material in both reactors [24].As a comparison, Sambusiti et al. (2020) reported VSS/TSS ratios around 92% in a submerged MBR treating a synthetic PW.From the step II, free VSS/TSS ratios began to decrease.At the end of step II, VSS/TSS ratios were 0.29 in the FBHBR and 0.53 in the MBBR.At the end of process operation, an increase of free VSS/TSS ratios was observed, up to 0.6 in the FBHBR and 0.73 in the MBBR.Nonvisible in the evolution of free TSS and free VSS concentrations, the evolution of free VSS/TSS ratios seems to reflect the perturbations caused by the different increases in TDS concentrations.The decrease in VSS/TSS ratio in both bioreactors reflect the loss of bacteria with the increase of TDS concentrations.In the literature, salinity has been reported as a parameter that strongly affects biological treatments.Previous studies reported that salinity starts to affect biological treatment operation at chloride concentration above 5-8 g•L −1 [30,50].The increase in salinity is also responsible for a decrease in VSS concentrations and an increase in TSS concentrations.This is observed in both bioreactors after day 142.The phenomenon that could explain this increase is an accumulation of salts inside the suspended sludge.Nevertheless, the VSS/TSS ratio increased again at day 176 in the MBBR and day 188 in the FBHBR.The increase of the ratio at the end of bioreactors operations seems to show an acclimation of the biomass to the saline wastewater.The combined effects of the increase in biomass and the accumulation of inorganic solids are so responsible for the increase in free TSS concentrations.At the end of the process operation, the Figure 3d clearly shows that the VSS/TSS ratios are higher in the MBBR than in the FBHBR.Indeed, the average ratio is 0.67 in the MBBR and 0.47 in the FBHBR.This could be explained by the fact that under the same loading rate, the suspended sludge production is higher in the MBBR.In the FBHBR, the packing material offers a bigger specific area than in the MBBR.Therefore, the organic loading required for the biofilm growth is bigger in the FBHBR than in the MBBR.This should lead to a bigger biofilm extent in the FBHBR.It is also worth noting that we measured only the free VSS.The biofilm in both bioreactors also contains a part of volatile solids.Hence, as it will be discussed below, the biofilm has a bigger extent in the FBHBR than in the MBBR.This also explains why free VSS concentration are lower in the FBHBR at the end of the process operation.Indeed, the large part of the bacteria are present in the biofilm rather than in the free sludge.
The presence of the packing material in both MBBR and FBHBR allowed the biofilm to develop.In the MBBR, the biofilm grew on the internal surface of the Kaldnes rings.In the FBHBR, the biofilm grew in the internal and external surface of the Kaldnes rings.There is no shear stress imposed to the external ring surface in the FBHBR, which explain the larger development of the biofilm.
Results of biofilm quantification during the process operation are shown in Figure 4.When the salinity was 1.5 g•L −1 (step I), biofilm weights were 4.52 mg per carrier and 6.59 mg per carrier in the MBBR and FBHBR, respectively.A decrease in biofilm weights was observed when the salinity increased to 8 g•L −1 , especially in the MBBR (1.37 mg per carrier in the MBBR and 6.1 mg per carrier in the FBHBR).Finally, biofilm weights increased again at a salinity of 20 g•L −1 , reaching 3.64 mg per carrier in the MBBR and 12.69 mg per carrier in the FBHBR.Both bioreactors were filled with 10 L of carriers which corresponds to approximately 10,800 rings.Based on this number, the total biofilm weight should be 39.3 g for the MBBR and 137.1 g for the FBHBR.This suggests a significant biofilm development in the FBHBR compared to the MBBR and that the operational conditions are still satisfactory for biofilm growth.Finally, the biofilm weight is 3.5 times higher in the FBHBR than in the MBBR, which is consistent with the development of the biofilm on the external surface of the packing material.The presence of the biofilm brings interesting advantages such as a better resistance to loading shocks and toxic compounds which should be beneficial at higher salinity.
39.3 g for the MBBR and 137.1 g for the FBHBR.This suggests a significant biofilm de opment in the FBHBR compared to the MBBR and that the operational conditions are satisfactory for biofilm growth.Finally, the biofilm weight is 3.5 times higher in FBHBR than in the MBBR, which is consistent with the development of the biofilm on external surface of the packing material.The presence of the biofilm brings interes advantages such as a better resistance to loading shocks and toxic compounds w should be beneficial at higher salinity.

Chemical Oxygen Demand Removal
The time evolution of COD removal efficiencies in both reactors is represente Figure 5b.For better understanding, the increase in the TDS concentration presente Figure 3a is plotted again in Figure 5a.Large fluctuations of COD removal efficien values were observed during the step I of bioreactors operation.These fluctuations attributed to the acclimation of the activated sludge biomass to the new influent.At end of the first step, COD removal efficiencies reached 80% for the FBHBR and 74% the MBBR.During the step II, COD removal efficiencies increased again by 13% in FBHBR and by 22% in the MBBR, reaching values above 90%.The values remained h during the step III.At the end of the bioreactors operation, COD removal efficiencies w 93% for the FBHBR and 85% for the MBBR.These values reflect the high biodegradab abilities of the synthetic PW in both reactors.Furthermore, the evolution of the COD moval efficiencies shows that, once the biomass has been properly acclimated to the p ence of salt (end of step I), the TDS concentration does not affect the treatment efficie in both reactors.In the literature, there are few studies relating the treatment of PW w salinity variations.Our findings here are consistent with some studies reporting high moval efficiencies with synthetic PW in the presence of high TDS concentrations.A example, Sharghi et al. (2016) reported COD removal efficiencies from 94.6 to 81.6% w salinity increasing from 100 to 250 g•L −1 in a submerged membrane bioreactors [18].T OLR was 2.6 kgCOD•m −3 •d −1 with an HRT of 24 h.Their results were also obtained wi consortium of isolated microorganisms.Compared to these results, COD removal effic cies obtained in this study in the MBBR and the FBHBR seem to indicate that high rem

Chemical Oxygen Demand Removal
The time evolution of COD removal efficiencies in both reactors is represented in Figure 5b.For better understanding, the increase in the TDS concentration presented in Figure 3a is plotted again in Figure 5a.Large fluctuations of COD removal efficiencies values were observed during the step I of bioreactors operation.These fluctuations are attributed to the acclimation of the activated sludge biomass to the new influent.At the end of the first step, COD removal efficiencies reached 80% for the FBHBR and 74% for the MBBR.During the step II, COD removal efficiencies increased again by 13% in the FBHBR and by 22% in the MBBR, reaching values above 90%.The values remained high during the step III.At the end of the bioreactors operation, COD removal efficiencies were 93% for the FBHBR and 85% for the MBBR.These values reflect the high biodegradability abilities of the synthetic PW in both reactors.Furthermore, the evolution of the COD removal efficiencies shows that, once the biomass has been properly acclimated to the presence of salt (end of step I), the TDS concentration does not affect the treatment efficiency in both reactors.In the literature, there are few studies relating the treatment of PW with salinity variations.Our findings here are consistent with some studies reporting high removal efficiencies with synthetic PW in the presence of high TDS concentrations.As an example, Sharghi et al. ( 2016) reported COD removal efficiencies from 94.6 to 81.6% with salinity increasing from 100 to 250 g•L −1 in a submerged membrane bioreactors [18].Their OLR was 2.6 kg COD •m −3 •d −1 with an HRT of 24 h.Their results were also obtained with a consortium of isolated microorganisms.Compared to these results, COD removal efficiencies obtained in this study in the MBBR and the FBHBR seem to indicate that high removal efficiencies can be achieved at shorter HRT and with acclimated microorganisms but smaller TDS concentrations.

Ecotoxicity Assessments
Water toxicity was assessed at the end of step I and step II.As TDS concentration may interfere in the toxicity of freshwater microorganisms, saline microorganisms were preferred to freshwater microorganisms for the toxicity assessment at the end of step II.Results of toxicity assessments for the first step (TDS = 1.5 g•L −1 ) are summarized in Table 3. Results are here expressed as EC50 (%), meaning a minimum of % EC50 indicates a high toxicity of the sample to the tested microorganism, whereas a maximum of % EC 50 indicates an absence of toxicity of the sample to the tested microorganism.Here, values of % EC50 shows that the inlet water was highly toxic to Daphnia magna and moderately toxic to Brachionus calyciflorus, Pseudokirchneriella subcapitata and Vibrio fischeri.After treatment, an absence of toxicity is observed for all microorganisms.This absence of toxicity was due to the elimination of pollutants in both bioreactors.

Pollutant Removal Performance
Phenol removal efficiencies quickly reached satisfying levels (more than 96% in both bioreactors at day 42.During step II and step III, phenol removal efficiencies remained high, with values above 95% in both bioreactors.No significant differences were observed between both bioreactors.This suggests that a moderate salinity (TDS concentration up to 20 g•L −1 ) does not affect phenol removal efficiencies.Furthermore, both configurations proved their ability to remove this compound.
Removal efficiencies of toluene, o-xylene and m-xylene in both MBBR and FBHBR over 99% were observed at each step of the operation of the bioreactor.This suggests that the increase of TDS concentration does not affect the treatment efficiency of volatile organic compounds.The high volatility of these compounds combined with the aeration are suggesting that the main removal mechanism of these compounds from the PW is achieved by volatilization and/or adsorption [51].
Removal efficiencies of naphthalene, phenanthrene and benzo[a]pyrene were also assessed in both bioreactors.Here again, high removal efficiencies are observed for naphthalene and phenanthrene (over 99%).Concerning benzo[a]pyrene, a removal efficiency of 93% was observed in both bioreactors at the end of the first step.The removal efficiency then increased to 97% and 98% at the end of the second and third steps, respectively.This was due to different inlet concentrations of benzo[a]pyrene measured during the sample analyses after each step (outlet concentrations were below the quantification limit for both outlets at each step).Nevertheless, removal efficiencies for benzo[a]pyrene (above 90%), as well as for naphthalene and phenanthrene suggest that TDS concentrations up to 20 g•L −1 do not significantly affect the efficiency of both MBBR and FBHBR for removing PAHs.

Ecotoxicity Assessments
Water toxicity was assessed at the end of step I and step II.As TDS concentration may interfere in the toxicity of freshwater microorganisms, saline microorganisms were preferred to freshwater microorganisms for the toxicity assessment at the end of step II.Results of toxicity assessments for the first step (TDS = 1.5 g•L −1 ) are summarized in Table 3. Results are here expressed as EC50 (%), meaning a minimum of % EC50 indicates a high toxicity of the sample to the tested microorganism, whereas a maximum of % EC 50 indicates an absence of toxicity of the sample to the tested microorganism.Here, values of % EC50 shows that the inlet water was highly toxic to Daphnia magna and moderately toxic to Brachionus calyciflorus, Pseudokirchneriella subcapitata and Vibrio fischeri.After treatment, an absence of toxicity is observed for all microorganisms.This absence of toxicity was due to the elimination of pollutants in both bioreactors.
Table 3. Ecotoxicity measurements of inlet water (TDS = 1.5 g•L −1 ) and outlet water (MBBR and FBHBR) expressed as EC50 (%) (NT refers as Nont-Toxic).For the step II (TDS = 8 g•L −1 ), the toxicity of water samples was assessed on Vibrio fischeri, Phaeodactylum tricornutum (marine algae) and Artemia salina (crustaceans).Results are presented in Table 4. Results show here that the inlet water was moderately toxic to all three microorganisms.An absence of toxicity was observed for Vibrio fischeri and Artemia salina in the outlet water of both MBBR and FBHBR.However, a moderate toxicity is observed for Phaeodactylum tricornutum (% EC50 of 38.4 and 60.6 in the MBBR and FBHBR, respectively).These values coupled to the fact that all pollutants were efficiently removed from the PW suggests that the observed toxicity was not due to the pollutant but to another parameter.We suggest that in this case, at a TDS concentration of 8 g•L −1 , the salt concentration may be too low to use a marine algae for toxicity assessment purposes.The bacterial population was assessed along the MBBR and FBHBR operation.13 samples were analyzed: the inoculum, the MBBR biofilm and MBBR free sludge at the end of each steps under TDS concentrations of 1.5 g•L −1 , 8 g•L −1 and 20 g•L −1 , and the FBHBR biofilm and FBHBR free sludge at the end of each steps under TDS concentrations of 1.5 g•L −1 , 8 g•L −1 and 20 g•L −1 .The inoculum was sampled at the beginning of the bioreactors operation (day 0).Samples at TDS concentration of 1.5 g•L −1 were sampled on day 101, samples at TDS concentration of 8 g•L −1 were sampled on day 183 and samples at TDS concentration of 20 g•L −1 were sampled on day 216.

Microorganism
The number of OTU of each sample are represented in Figure 6b.In the inoculum, the number of OTU was 915.In the MBBR, the number of OTU was 437 for the biofilm and 301 for the free sludge at a TDS concentration of 1.5 g•L −1 .When the TDS concentration was set to 8 g•L −1 , the number of OTU was 473 for the biofilm and 225 for the free sludge.In the FBHBR, the number of OTU was 509 for the biofilm and 387 for the free sludge.When the concentration was set to 8 g•L −1 , the number of OTU decreased to 283 for the biofilm and 300 for the free sludge.At the end of the bioreactors operation, i.e., at a TDS concentration of 20 g•L −1 , the number of OTU were 230 and 163 in the MBBR biofilm and MBBR free sludge, respectively.In the FBHBR biofilm and FBHBR free sludge, the number of OTU was 239 and 178, respectively.These values suggest that the diversity of species decreased during the bioreactors operation.However, the decrease in the number of OTU along the time of operation of bioreactors is consistent the further increase of TDS concentration.There is a tendency of higher number of OTU in the biofilm in both bioreactors compared to the free sludge.This suggests that the diversity of species was moderately higher in the biofilm than in the free sludge.
Figure 6a represents the evolution of the Shannon-Weaver index in all biomasses during the bioreactors operation.The index of the inoculum was 5.56, suggesting a high diversity of species.A significant decrease of the indices was observed during the bioreactors operation.In the MBBR, the index was 3.83 in the biofilm and 3.5 for the free sludge at a TDS concentration of 1.5 g•L −1 .When the TDS concentration was set up to 8 g•L −1 , the Shannon-Weaver index decreased to 3.61 for the biofilm and 2.90 for the free sludge.In the FBHBR, the index was 3.55 for the biofilm and 3.2 for the free sludge at a TDS concentration of 1.5 g•L −1 .When the TDS concentration was set up to 8 g•L −1 , the index decreased to 2.94 for the biofilm and 3.1 for the free sludge.Finally, at the end of the bioreactors operation, Shannon-Weaver indices were 2.86 and 2.99 in MBBR biofilm and MBBR free sludge.In the FBHBR biofilm and FBHBR free sludge, Shannon-Weaver indices were 3.29 and 3.03.These values lead to the same conclusion than the evolution of the number of OTU, i.e., a significant speciation of species in both reactors with time as TDS concentrations increased.The increase in salinity is probably lethal to some bacteria communities and explain the decrease in the Shannon-Weaver index.
The evolution of the Simpson reciprocal index of all biomasses is presented in Figure 6c.The Simpson reciprocal index of the inoculum was 109.6.Here again, a significant decrease of the Simpson reciprocal index is observed for all biomasses during the bioreactors operation.In the MBBR, the Simpson reciprocal indices were 16.4 for the biofilm and 16.1 for the free sludge at a TDS concentration of 1.5 g•L −1 .When the TDS concentration was set up to 8 g•L −1 , Simpson reciprocal indices were 9.18 for the biofilm and 9.80 for the free sludge.In the FBHBR, the Simpson reciprocal indices were 8.7 in the biofilm and 10.5 in the free sludge at a TDS concentration of 1.5 g•L −1 ; When the TDS concentration was set up to 8 g•L −1 , the indices were 7.21 in the biofilm and 10 in the free sludge.Finally, at the end of bioreactors operation, Simpson reciprocal indices were 5.69 and 10.7 in the MBBR biofilm and MBBR free sludge, respectively.In the FBHBR biofilm and FBHBR free sludge, Simpson reciprocal indices were 13.87 and 11.1, respectively.
The number of OTU, the Shannon Weaver index and the Simpson reciprocal index all indicate that the microbial species diversity is significantly reduced over the time and in both MBBR and FBHBR, as the TDS concentration is slowly increased from 1.5 to 20 g•L −1 .This result is consistent with the fact that the increase in the TDS concentration is probably lethal to some bacterial species.Thus, a specialization appeared in both reactors, i.e., some bacterial species became more predominant in both MBBR and FBHBR.Even if the sludge retention time is not a limiting factor for biomass development in the biofilm, it appears that the PW composition has more influence on the bacterial development.Concerning the difference between the biofilm diversity and the free sludge diversity, the same conclusions as in the previous study are made here [17].Indeed, all the diversity indices seem to indicate that the diversity is bigger in the biofilm than in the free sludge.Lastly, the diversity of species is slightly higher in the MBBR than in the FBHBR.Figure 7 shows the results of 16S rRNA analysis of the bacterial communities at the genus level in the MBBR and the FBHBR when the TDS concentration was set to 1.5 g•L −1 .The represented sequences only represent a minimum of 1% of their total sequences.The percentage of unclassified genus in the inoculum was surprisingly high (53.6%).In the MBBR, a similarity is observed between the biofilm and the free sludge.
In the MBBR biofilm, four species were dominants (corresponding to 41% of the total sequences): Dechloromonas sp., Flavobacterium sp., Mycoplana sp. and Zooglea sp.In the MBBR free sludge, four species were also dominant: Dechloromonas sp., Hydrogenophaga sp., Flavobacterium sp. and Mycoplana sp.(corresponding to 47% of the total sequences).In   The represented sequences only represent a minimum of 1% of their total sequences.The percentage of unclassified genus in the inoculum was surprisingly high (53.6%).In the MBBR, a similarity is observed between the biofilm and the free sludge.
It is interesting to see that the most abundant species are Dechloromonas sp., Hydrogenophaga sp. and flavobacterium sp.In the literature, Dechloromonas sp were found to be able to degrade benzene, toluene, ethylbenzene and xylene [52].This seems to indicate an acclimation of microorganisms to the PW influent.The assessment of bacterial population at the genus level in the MBBR biofilm and free sludge as well as in the FBHBR biofilm and free sludge under a TDS concentration of 8 g•L −1 are shown in Figure 8. From this second assessment, there were 57%, 60%, 48% and 59 % of unclassified genus in the MBBR biofilm, MBBR free sludge, FBHBR biofilm and FBHBR free sludge, respectively.Some genus already present under a TDS concentration of 1.5 g•L −1 were found again in this second sequencing: Hydrogenophaga sp. in all biomasses, Thiothrix sp. in MBBR and FBHBR biofilm, Flavobacterium sp. in MBBR and FBHBR free sludge, Caulobacter sp. in FBHBR biofilm, Planctomyces sp. in FBHBR free sludge and Treponema sp. in MBBR biofilm.
With respect to the first sequencing, some gena developed further when the TDS concentration increased to 8 g•L −1 : Lewinella sp. and Pedobacter sp. in the MBBR and FBHBR free sludge, Nannocystis sp. in MBBR and FBHBR biofilm, Novosphingobium sp. and Sphingomonadales sp. in FBHBR free sludge, Haliscomenobacter sp. in FBHBR biofilm, A4b sp. and Mb2424 sp. in MBBR biofilm.The most preponderant new genus were Lewinella sp. and Pedobacter sp.
The proportion of Lewinella sp. were 15% in MBBR free sludge and 8% in FBHBR free sludge.Lewin (1970) originally described three novel marine species as Herpetosiphon sp.[53].Those strains were later sequenced and related to the genus Lewinella [54].Lewinella sp. are marine bacteria and originate from seawater.In the literature, there are only two articles relating the genus Lewinella to the degradation of hydrocarbon compounds [55,56]  The last DNA sequencing were performed at the end of the experiment at a TDS concentration of 20 g•L −1 .Results of the sequencing are represented in Figure 9.In this sequencing, there were 80%, 77%, 74%, 72% of unclassified genus in the MBBR biofilm, MBBR free sludge, FBHBR biofilm and FBHBR free sludge, respectively.There were two gena common to all tested samples, Hydrogenophaga sp. and Cytophagales sp. with abundances ranging from 4 to 12%.Lewinella sp. were found abundant at 4% only in MBBR free sludge.In the FBHBR biofilm, the genus Fusibacter is found abundant at 5.1 %.

Conclusions
This experimental study dealt with the influence of the increase in TDS concentration on treatment performances of a synthetic PW using two types of hybrid bioreactors, namely a MBBR and a FBHBR.The main conclusions obtained are:

•
Both bioreactors exhibited an increase in TSS concentrations, once the acclimation phase was over, with the increase in TDS concentrations.At the same time, the free VSS/TSS ratio tended to decrease in both bioreactors, which suggests an accumulation of inorganic solids in the free suspended sludge.

•
Both bioreactors were proved to be efficient to remove the COD from the influent as well as VOCs and PAHs.An absence of toxicity was noticed in the outlet water performing tests with different microorganisms.

•
A decrease in the bacterial diversity indices was observed with respect to the inoculum, leading to the predominance of a lower number of bacterial species.Despite a large part of unclassified gena, some gena, such as Lewinella sp.seem to indicate a logical shift of the bacterial community from freshwater bacteria towards saline bacteria.
In future work, increasing the TDS concentration to higher level will be of major importance to reach more realistic of real produced water treatment.

Energies 2021 ,Figure 3 .
Figure 3. Evolution of solids concentrations in both MBBR and FBHBR during the experim concentration (a), free TSS concentration (b), free VSS concentration (c) and free VSS/TSS r

Figure 3 .
Figure 3. Evolution of solids concentrations in both MBBR and FBHBR during the experiment: TDS concentration (a), free TSS concentration (b), free VSS concentration (c) and free VSS/TSS ratio (d).

Figure 4 .
Figure 4. Evolution of biofilm weight in both MBBR and FBHBR during the process operation.

Figure 4 .
Figure 4. Evolution of biofilm weight in both MBBR and FBHBR during the process operation.

Figure 5 .
Figure 5. Evolution of COD removal efficiency and phenol removal efficiency in both MBBR and FBHBR during the experiment: TDS concentration (a), COD removal (b) and phenol removal (c).

Figure 5 .
Figure 5. Evolution of COD removal efficiency and phenol removal efficiency in both MBBR and FBHBR during the experiment: TDS concentration (a), COD removal (b) and phenol removal (c).

Figure 6 .
Figure 6.Evolution of bacterial diversity indices in both MBBR and FBHBR: Shannon Weaver index (a), number of OTU (b) and Simspon reciprocal index (c).

Figure 7
Figure7shows the results of 16S rRNA analysis of the bacterial communities at the genus level in the MBBR and the FBHBR when the TDS concentration was set to 1.5 g•L −1 .The represented sequences only represent a minimum of 1% of their total sequences.The percentage of unclassified genus in the inoculum was surprisingly high (53.6%).In the MBBR, a similarity is observed between the biofilm and the free sludge.

Figure 7 .
Figure 7. Bacterial assessment at the genus level in the MBBR and the FBHBR at TDS = 1.5 g•L −1 .

Figure 7 .
Figure 7. Bacterial assessment at the genus level in the MBBR and the FBHBR at TDS = 1.5 g•L −1 .

Figure 8 .
Figure 8. Bacterial assessment at the genus level in the MBBR and FBHBR at TDS = 8 g•L −1 .

Figure 8 .
Figure 8. Bacterial assessment at the genus level in the MBBR and FBHBR at TDS = 8 g•L −1 .

Figure 9 .
Figure 9. Bacterial assessment at the genus level in the MBBR and the FBHBR at TDS = 20 g•L −1 (Un stands for unclassified).Figure 9. Bacterial assessment at the genus level in the MBBR and the FBHBR at TDS = 20 g•L −1 (Un stands for unclassified).

Figure 9 .
Figure 9. Bacterial assessment at the genus level in the MBBR and the FBHBR at TDS = 20 g•L −1 (Un stands for unclassified).Figure 9. Bacterial assessment at the genus level in the MBBR and the FBHBR at TDS = 20 g•L −1 (Un stands for unclassified).

Table 1 .
Composition of the synthetic produced water.

Table 2 .
Bioreactors operation during the experiment.